Changeset 17056

Timestamp:

12/24/13 01:19:21 (11 years ago)

Author:

Mathieu Morlighem

Message:

BUG: fixed final nodal function derivatives P1xP3

File:

• issm/trunk-jpl/src/c/classes/Elements/PentaRef.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/classes/Elements/PentaRef.cpp

@@ -1986 +1986 @@
                         dbasis[NUMNODESP1xP3*0+10] = - (27./32.) *(zeta-1)*(zeta+1./3.)*(zeta+1);
                         dbasis[NUMNODESP1xP3*1+10] = ((9.*SQRT3)/32.)  *(zeta-1.)*(zeta+1./3.)*(zeta+1);
-                        dbasis[NUMNODESP1xP3*2+10] = -gauss->coord2 *(27./16.) * 2.* zeta *( zeta + (1./3.) + ( zeta - 1. )*( zeta + 1. ));
+                        dbasis[NUMNODESP1xP3*2+10] = -gauss->coord2 *(27./16.) *( 2.* zeta *( zeta + (1./3.)) + ( zeta - 1. )*( zeta + 1. ));
                         /*Nodal function 12*/   
                         dbasis[NUMNODESP1xP3*0+11] = 0.;
                         dbasis[NUMNODESP1xP3*1+11] = -((9.*SQRT3)/16.) *(zeta-1.)*(zeta+1./3.)*(zeta+1);
-                        dbasis[NUMNODESP1xP3*2+11] = -gauss->coord3 *(27./16.)* 2.* zeta *( zeta + (1./3.) + ( zeta - 1. )*( zeta + 1. ));
+                        dbasis[NUMNODESP1xP3*2+11] = -gauss->coord3 *(27./16.)*( 2.* zeta *( zeta + (1./3.)) + ( zeta - 1. )*( zeta + 1. ));
                         return;
                 default: