int lagcoef int  evflag,
double  u,
double *  us
 

Definition at line 564 of file vecmatop.cpp.

00565 {
00566   //double u1 = u;
00567   double u2 = u*u;
00568   double u3 = u*u*u;
00569   if(       evflag==EVFLAG_VL) {
00570          us[0] = (u3-3*u2+2*u)/(-6.0);
00571          us[1] = (u3-2*u2-u+2)/(2.0);
00572          us[2] = (u3-u2-2*u)  /(-2.0);
00573          us[3] = (u3-u)       /(6.0);
00574   } else if(evflag==EVFLAG_FD) {
00575          us[0] = (3*u2-6*u+2)/(-6.0);
00576          us[1] = (3*u2-4*u-1)/(2.0);
00577          us[2] = (3*u2-2*u-2)/(-2.0);
00578          us[3] = (3*u2-1)    /(6.0);
00579   } else if(evflag==EVFLAG_SD) {
00580          assert(0); //TODO
00581   }
00582   return 0;
00583 }


Generated on Sun Dec 4 18:13:13 2005 for common by  doxygen 1.4.5