Lines Matching refs:h
456 double* h; in NaturalSpline() local
462 h=new double[n+1]; 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()
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()
479 a[0] =a[0]*h[1]/(h[0]+h[1]); in NaturalSpline()
480 a[n-2]=a[n-2]*h[n-2]/(h[n-1]+h[n-2]); in NaturalSpline()
481 d[0] =d[0]-h[0]; in NaturalSpline()
482 d[n-2]=d[n-2]-h[n-1]; in NaturalSpline()
483 c[0] =c[0]-h[0]; in NaturalSpline()
484 b[n-2]=b[n-2]-h[n-1]; in NaturalSpline()
488 a[0] =a[0]-1.5*((y[1]-y[0])/h[0]-Marg0); in NaturalSpline()
489 a[n-2]=a[n-2]-1.5*(MargN-(y[n]-y[n-1])/h[n-1]); in NaturalSpline()
490 d[0] =d[0]-h[0]*0.5; in NaturalSpline()
491 d[n-2]=d[n-2]-h[n-1]*0.5; in NaturalSpline()
494 a[0] =a[0]-h[0]*Marg0*0.5; in NaturalSpline()
495 a[n-2]=a[n-2]-h[n-1]*MargN*0.5; in NaturalSpline()
498 a[0] =a[0]+Marg0*h[0]*h[0]*0.5; in NaturalSpline()
499 a[n-2]=a[n-2]-MargN*h[n-1]*h[n-1]*0.5; in NaturalSpline()
500 d[0] =d[0]+h[0]; in NaturalSpline()
501 d[n-2]=d[n-2]+h[n-1]; in NaturalSpline()
508 if (error!=0) { delete[] a; delete[] h; return error+2; } in NaturalSpline()
517 c[0]=c[1]+h[0]*(c[1]-c[2])/h[1]; in NaturalSpline()
518 c[n]=c[n-1]+h[n-1]*(c[n-1]-c[n-2])/h[n-2]; in NaturalSpline()
522 c[0]=1.5*((y[1]-y[0])/h[0]-Marg0); in NaturalSpline()
523 c[0]=(c[0]-c[1]*h[0]*0.5)/h[0]; in NaturalSpline()
524 c[n]=1.5*((y[n]-y[n-1])/h[n-1]-MargN); in NaturalSpline()
525 c[n]=(c[n]-c[n-1]*h[n-1]*0.5)/h[n-1]; in NaturalSpline()
532 c[0]=c[1]-Marg0*h[0]*0.5; in NaturalSpline()
533 c[n]=c[n-1]+MargN*h[n-1]*0.5; 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()
541 delete[] h; in NaturalSpline()