int spcoef int  evflag,
double  u,
double *  us
 

Definition at line 338 of file vecmatop.cpp.

00339 {
00340   double u1 = u;
00341   double u2 = u*u;
00342   double u3 = u*u*u;
00343   if(       evflag==EVFLAG_VL) {
00344          us[0] = (  1 - 3*u1 + 3*u2 -   u3)/6.0;
00345          us[1] = (  4        - 6*u2 + 3*u3)/6.0;
00346          us[2] = (  1 + 3*u1 + 3*u2 - 3*u3)/6.0;
00347          us[3] = (                  +   u3)/6.0;
00348   } else if(evflag==EVFLAG_FD) {
00349          us[0] = (- 3 + 6*u1 - 3*u2)/6.0;
00350          us[1] = (    -12*u1 + 9*u2)/6.0;
00351          us[2] = (  3 + 6*u1 - 9*u2)/6.0;
00352          us[3] = (             3*u2)/6.0;
00353   } else if(evflag==EVFLAG_SD) {
00354          us[0] = (  6 - 6*u1 ) / 6.0;
00355          us[1] = (-12 +18*u1 ) / 6.0;
00356          us[2] = (  6 -18*u1 ) / 6.0;
00357          us[3] = (      6*u1 ) / 6.0;    //assert(0); //TODO;
00358   }  
00359   return 0;
00360 }


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