/[escript]/branches/domexper/dudley/src/Assemble_setNormal.c
ViewVC logotype

Diff of /branches/domexper/dudley/src/Assemble_setNormal.c

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

revision 3151 by jfenwick, Fri Sep 3 03:49:32 2010 UTC revision 3152 by jfenwick, Fri Sep 3 05:48:31 2010 UTC
# Line 38  void Dudley_Assemble_setNormal(Dudley_No Line 38  void Dudley_Assemble_setNormal(Dudley_No
38    NN=elements->numNodes;    NN=elements->numNodes;
39    numDim=nodes->numDim;    numDim=nodes->numDim;
40    
41    numQuad=reference_element->Parametrization->numQuadNodes;    numQuad=reference_element->BasisFunctions->numQuadNodes;
42    numDim_local=reference_element->Parametrization->Type->numDim;    numDim_local=reference_element->BasisFunctions->Type->numDim;
43    NS=reference_element->Parametrization->Type->numShapes;    NS=reference_element->BasisFunctions->Type->numShapes;
44    
45    
46        
# Line 76  void Dudley_Assemble_setNormal(Dudley_No Line 76  void Dudley_Assemble_setNormal(Dudley_No
76                            /* gather local coordinates of nodes into local_X: */                            /* gather local coordinates of nodes into local_X: */
77                            Dudley_Util_Gather_double(NS,&(elements->Nodes[INDEX2(0,e,NN)]),numDim,nodes->Coordinates,local_X);                            Dudley_Util_Gather_double(NS,&(elements->Nodes[INDEX2(0,e,NN)]),numDim,nodes->Coordinates,local_X);
78                            /*  calculate dVdv(i,j,q)=local_X(i,n)*DSDv(n,j,q) */                            /*  calculate dVdv(i,j,q)=local_X(i,n)*DSDv(n,j,q) */
79                            Dudley_Util_SmallMatMult(numDim,numDim_local*numQuad,dVdv,NS,local_X,reference_element->Parametrization->dSdv);                            Dudley_Util_SmallMatMult(numDim,numDim_local*numQuad,dVdv,NS,local_X,reference_element->BasisFunctions->dSdv);
80                            /* get normalized vector:  */                            /* get normalized vector:  */
81                            normal_array=getSampleDataRW(normal,e);                            normal_array=getSampleDataRW(normal,e);
82                            Dudley_NormalVector(numQuad,numDim,numDim_local,dVdv,normal_array);                            Dudley_NormalVector(numQuad,numDim,numDim_local,dVdv,normal_array);

Legend:
Removed from v.3151  
changed lines
  Added in v.3152

  ViewVC Help
Powered by ViewVC 1.1.26