Lines Matching refs:i
242 sal_uInt16 i; in TriDiagGS() local
253 for (i = 1; i < n; i++) in TriDiagGS()
254 { if ( fabs(diag[i-1]) < MACH_EPS ) /* Wenn ein diag[i] = 0 */ in TriDiagGS()
256 lower[i] /= diag[i-1]; /* gung. */ in TriDiagGS()
257 diag[i] -= lower[i] * upper[i-1]; in TriDiagGS()
263 for (i = 1; i < n; i++) /* Vorwaertselimination */ in TriDiagGS()
264 b[i] -= lower[i] * b[i-1]; in TriDiagGS()
268 i=j; in TriDiagGS()
269 b[i] = ( b[i] - upper[i] * b[i+1] ) / diag[i]; in TriDiagGS()
376 sal_uInt16 i; in ZyklTriDiagGS() local
391 for (i = 1; i < n-2; i++) in ZyklTriDiagGS()
392 { diag[i] -= lower[i] * upper[i-1]; in ZyklTriDiagGS()
393 if ( fabs(diag[i]) < MACH_EPS ) return(2); in ZyklTriDiagGS()
394 temp = 1.0 / diag[i]; in ZyklTriDiagGS()
395 upper[i] *= temp; in ZyklTriDiagGS()
396 ricol[i] = -lower[i] * ricol[i-1] * temp; in ZyklTriDiagGS()
402 for (i = 1; i < n-2; i++) in ZyklTriDiagGS()
403 lowrow[i] = -lowrow[i-1] * upper[i-1]; in ZyklTriDiagGS()
408 for (temp = 0.0, i = 0; i < n-2; i++) in ZyklTriDiagGS()
409 temp -= lowrow[i] * ricol[i]; in ZyklTriDiagGS()
416 for (i = 1; i < n-1; i++) in ZyklTriDiagGS()
417 b[i] = ( b[i] - b[i-1] * lower[i] ) / diag[i]; in ZyklTriDiagGS()
419 for (temp = 0.0, i = 0; i < n-2; i++) in ZyklTriDiagGS()
420 temp -= lowrow[i] * b[i]; in ZyklTriDiagGS()
426 i=j; in ZyklTriDiagGS()
427 b[i] -= upper[i] * b[i+1] + ricol[i] * b[n-1]; in ZyklTriDiagGS()
454 sal_uInt16 i; in NaturalSpline() local
463 for (i=0;i<n;i++) { in NaturalSpline()
464 h[i]=x[i+1]-x[i]; in NaturalSpline()
465 if (h[i]<=0.0) { delete[] a; delete[] h; return 1; } in NaturalSpline()
467 for (i=0;i<n-1;i++) { in NaturalSpline()
468 a[i]=3.0*((y[i+2]-y[i+1])/h[i+1]-(y[i+1]-y[i])/h[i]); in NaturalSpline()
469 b[i]=h[i]; in NaturalSpline()
470 c[i]=h[i+1]; in NaturalSpline()
471 d[i]=2.0*(h[i]+h[i+1]); in NaturalSpline()
509 for (i=0;i<n-1;i++) c[i+1]=a[i]; in NaturalSpline()
536 for (i=0;i<n;i++) { in NaturalSpline()
537 b[i]=(y[i+1]-y[i])/h[i]-h[i]*(c[i+1]+2.0*c[i])/3.0; in NaturalSpline()
538 d[i]=(c[i+1]-c[i])/(3.0*h[i]); in NaturalSpline()
562 sal_uInt16 i,im1,nm1; //integer in PeriodicSpline() local
570 for (i=0;i<=nm1;i++) if (x[i+1]<=x[i]) return 2; // muss streng nonoton fallend sein! in PeriodicSpline()
579 c[1]=c[1]-3.0*((y[i]-y[0])/(x[1]-x[0])); in PeriodicSpline()
583 for (i=1;i<=nm1;i++) { in PeriodicSpline()
584 im1=i-1; in PeriodicSpline()
585 hl=x[i]-x[im1]; in PeriodicSpline()
586 hr=x[i+1]-x[i]; in PeriodicSpline()
590 a[im1]=3.0*((y[i+1]-y[i])/hr-(y[i]-y[im1])/hl); in PeriodicSpline()
607 for (i=0;i<=nm1;i++) c[i+1]=a[i]; in PeriodicSpline()
610 for (i=0;i<=nm1;i++) { in PeriodicSpline()
611 hl=x[i+1]-x[i]; in PeriodicSpline()
612 b[i]=(y[i+1]-y[i])/hl; in PeriodicSpline()
613 b[i]=b[i]-hl*(c[i+1]+2.0*c[i])/3.0; in PeriodicSpline()
614 d[i]=(c[i+1]-c[i])/hl/3.0; in PeriodicSpline()
644 sal_uInt16 i; in ParaSpline() local
653 for (i=0;i<n;i++) { in ParaSpline()
654 deltX=x[i+1]-x[i]; deltY=y[i+1]-y[i]; in ParaSpline()
657 T[i+1]=T[i]+sqrt(delt); in ParaSpline()
731 sal_uInt16 i; in CalcSpline() local
740 for (i=0;i<rPoly.GetSize();i++) { in CalcSpline()
741 Pt=rPoly.GetPoint(i); in CalcSpline()
742 if (i==0 || Pt!=P0) { in CalcSpline()
832 sal_uInt16 i; // aktuelles Teilpolynom in Spline2Poly() local
843 i=0; in Spline2Poly()
844 while (i<n) { // n Teilpolynome malen in Spline2Poly()
845 t=tv[i]+Step; in Spline2Poly()
848 bEnde=t>=tv[i+1]; in Spline2Poly()
849 if (bEnde) t=tv[i+1]; in Spline2Poly()
850 dt1=t-tv[i]; dt2=dt1*dt1; dt3=dt2*dt1; in Spline2Poly()
851 x=long(ax[i]+bx[i]*dt1+cx[i]*dt2+dx[i]*dt3); in Spline2Poly()
852 y=long(ay[i]+by[i]*dt1+cy[i]*dt2+dy[i]*dt3); in Spline2Poly()
863 i++; // naechstes Teilpolynom in Spline2Poly()