/[escript]/trunk-mpi-branch/finley/src/ShapeFunctions.c
ViewVC logotype

Diff of /trunk-mpi-branch/finley/src/ShapeFunctions.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1306 by ksteube, Tue Sep 18 05:51:09 2007 UTC revision 1307 by gross, Tue Sep 18 08:53:11 2007 UTC
# Line 519  void Finley_Shape_Tet10(int NumV,double* Line 519  void Finley_Shape_Tet10(int NumV,double*
519      S(8,i) = (1.-x-y-z)*4.*z;      S(8,i) = (1.-x-y-z)*4.*z;
520      S(9,i) = 4.*x*z;      S(9,i) = 4.*x*z;
521      S(10,i)= 4.*y*z;      S(10,i)= 4.*y*z;
522    
523      DSDV(1,1,i)= -3.+4.*x+4.*y+4.*z;      DSDV(1,1,i)= -3.+4.*x+4.*y+4.*z;
524      DSDV(1,2,i)= -3.+4.*x+4.*y+4.*z;      DSDV(1,2,i)= -3.+4.*x+4.*y+4.*z;
525      DSDV(1,3,i)= -3.+4.*x+4.*y+4.*z;      DSDV(1,3,i)= -3.+4.*x+4.*y+4.*z;
526    
527    
528      DSDV(2,1,i)= -1.+4.*x;      DSDV(2,1,i)= -1.+4.*x;
529      DSDV(2,2,i)=   0.;      DSDV(2,2,i)=   0.;
530      DSDV(2,3,i)=   0.;      DSDV(2,3,i)=   0.;
531    
532      DSDV(3,1,i)=   0.;      DSDV(3,1,i)=   0.;
533      DSDV(3,2,i)=  -1.           +4.*y;      DSDV(3,2,i)=  -1.           +4.*y;
534      DSDV(3,3,i)=   0.;      DSDV(3,3,i)=   0.;
535    
536      DSDV(4,1,i)=   0.;      DSDV(4,1,i)=   0.;
537      DSDV(4,2,i)=   0.;      DSDV(4,2,i)=   0.;
538      DSDV(4,3,i)=  -1.                        +4.*z;      DSDV(4,3,i)=  -1.                        +4.*z;
539    
540      DSDV(5,1,i)=   4. -8.*x -4.*y -4.*z;      DSDV(5,1,i)=   4. -8.*x -4.*y -4.*z;
541      DSDV(5,2,i)=      -4.*x;      DSDV(5,2,i)=      -4.*x;
542      DSDV(5,3,i)=      -4.*x;      DSDV(5,3,i)=      -4.*x;
543      DSDV(6,1,i)=                      4.*y;  
544      DSDV(6,2,i)=         4.*x;      DSDV(6,1,i)=             4.*y;
545        DSDV(6,2,i)=       4.*x;
546      DSDV(6,3,i)=   0.;      DSDV(6,3,i)=   0.;
547      DSDV(7,1,i)=                   -4.*y;  
548        DSDV(7,1,i)=            -4.*y;
549      DSDV(7,2,i)=   4. -4.*x -8.*y -4.*z;      DSDV(7,2,i)=   4. -4.*x -8.*y -4.*z;
550      DSDV(7,3,i)=                   -4.*y;      DSDV(7,3,i)=            -4.*y;
551    
552      DSDV(8,1,i)=                                -4.*z;      DSDV(8,1,i)=                                -4.*z;
553      DSDV(8,2,i)=                                -4.*z;      DSDV(8,2,i)=                                -4.*z;
554      DSDV(8,3,i)=   4. -4.*x -4.*y -8.*z;      DSDV(8,3,i)=   4. -4.*x -4.*y -8.*z;
555      DSDV(9,1,i)=                                   4.*z;  
556        DSDV(9,1,i)=                                 4.*z;
557      DSDV(9,2,i)=   0.;      DSDV(9,2,i)=   0.;
558      DSDV(9,3,i)=         4.*x;      DSDV(9,3,i)=         4.*x;
559    
560      DSDV(10,1,i)=  0.;      DSDV(10,1,i)=  0.;
561      DSDV(10,2,i)=                                  4.*z;      DSDV(10,2,i)=                                  4.*z;
562      DSDV(10,3,i)=                      4.*y;      DSDV(10,3,i)=                      4.*y;

Legend:
Removed from v.1306  
changed lines
  Added in v.1307

  ViewVC Help
Powered by ViewVC 1.1.26