/[escript]/trunk/finley/src/ElementFile_jacobeans.c
ViewVC logotype

Diff of /trunk/finley/src/ElementFile_jacobeans.c

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

revision 3258 by jfenwick, Thu Jan 28 02:03:15 2010 UTC revision 3259 by jfenwick, Mon Oct 11 01:48:14 2010 UTC
# Line 120  Finley_ElementFile_Jacobeans* Finley_Ele Line 120  Finley_ElementFile_Jacobeans* Finley_Ele
120            Finley_setError(SYSTEM_ERROR,"Finley_ElementFile_borrowJacobeans: 1D does not support local dimension 0.");            Finley_setError(SYSTEM_ERROR,"Finley_ElementFile_borrowJacobeans: 1D does not support local dimension 0.");
121               } else if (refElement->numLocalDim==1) {               } else if (refElement->numLocalDim==1) {
122                    if (out->numSides==2) {                    if (out->numSides==2) {
123                         Assemble_jacobeans_1D(nodes->Coordinates,out->numQuadTotal,shape->QuadWeights,                         Finley_Assemble_jacobeans_1D(nodes->Coordinates,out->numQuadTotal,shape->QuadWeights,
124                                              shape->Type->numShapes,self->numElements,numNodes,self->Nodes,                                              shape->Type->numShapes,self->numElements,numNodes,self->Nodes,
125                                              shape->dSdv,basis->Type->numShapes,dBdv, out->DSDX,out->volume,self->Id);                                              shape->dSdv,basis->Type->numShapes,dBdv, out->DSDX,out->volume,self->Id);
126    
# Line 137  Finley_ElementFile_Jacobeans* Finley_Ele Line 137  Finley_ElementFile_Jacobeans* Finley_Ele
137               } else if (refElement->numLocalDim==1) {               } else if (refElement->numLocalDim==1) {
138                    if (out->BasisFunctions->Type->numDim==2) {                    if (out->BasisFunctions->Type->numDim==2) {
139                       if (out->numSides==1) {                       if (out->numSides==1) {
140                          Assemble_jacobeans_2D_M1D_E2D(nodes->Coordinates,out->numQuadTotal,shape->QuadWeights,                          Finley_Assemble_jacobeans_2D_M1D_E2D(nodes->Coordinates,out->numQuadTotal,shape->QuadWeights,
141                                                        shape->Type->numShapes,self->numElements,numNodes,self->Nodes,                                                        shape->Type->numShapes,self->numElements,numNodes,self->Nodes,
142                                                        shape->dSdv,basis->Type->numShapes,dBdv,                                                        shape->dSdv,basis->Type->numShapes,dBdv,
143                                                        out->DSDX,out->volume,self->Id);                                                        out->DSDX,out->volume,self->Id);
144                       } else if (out->numSides==2) {                       } else if (out->numSides==2) {
145                          Assemble_jacobeans_2D_M1D_E2D_C(nodes->Coordinates,out->numQuadTotal,shape->QuadWeights,                          Finley_Assemble_jacobeans_2D_M1D_E2D_C(nodes->Coordinates,out->numQuadTotal,shape->QuadWeights,
146                                                          shape->Type->numShapes,self->numElements,numNodes,self->Nodes,                                                          shape->Type->numShapes,self->numElements,numNodes,self->Nodes,
147                                                          shape->dSdv,basis->Type->numShapes,dBdv,                                                          shape->dSdv,basis->Type->numShapes,dBdv,
148                                                          out->DSDX,out->volume,self->Id);                                                          out->DSDX,out->volume,self->Id);
# Line 151  Finley_ElementFile_Jacobeans* Finley_Ele Line 151  Finley_ElementFile_Jacobeans* Finley_Ele
151                        }                        }
152                    }  else if (out->BasisFunctions->Type->numDim==1) {                    }  else if (out->BasisFunctions->Type->numDim==1) {
153                       if (out->numSides==1) {                       if (out->numSides==1) {
154                          Assemble_jacobeans_2D_M1D_E1D(nodes->Coordinates,out->numQuadTotal,shape->QuadWeights,                          Finley_Assemble_jacobeans_2D_M1D_E1D(nodes->Coordinates,out->numQuadTotal,shape->QuadWeights,
155                                                        shape->Type->numShapes,self->numElements,numNodes,self->Nodes,                                                        shape->Type->numShapes,self->numElements,numNodes,self->Nodes,
156                                                        shape->dSdv,basis->Type->numShapes,dBdv,                                                        shape->dSdv,basis->Type->numShapes,dBdv,
157                                                        out->DSDX,out->volume,self->Id);                                                        out->DSDX,out->volume,self->Id);
158                       } else if (out->numSides==2) {                       } else if (out->numSides==2) {
159                          Assemble_jacobeans_2D_M1D_E1D_C(nodes->Coordinates,out->numQuadTotal,shape->QuadWeights,                          Finley_Assemble_jacobeans_2D_M1D_E1D_C(nodes->Coordinates,out->numQuadTotal,shape->QuadWeights,
160                                                          shape->Type->numShapes,self->numElements,numNodes,self->Nodes,                                                          shape->Type->numShapes,self->numElements,numNodes,self->Nodes,
161                                                          shape->dSdv,basis->Type->numShapes,dBdv,                                                          shape->dSdv,basis->Type->numShapes,dBdv,
162                                                          out->DSDX,out->volume,self->Id);                                                          out->DSDX,out->volume,self->Id);
# Line 168  Finley_ElementFile_Jacobeans* Finley_Ele Line 168  Finley_ElementFile_Jacobeans* Finley_Ele
168                    }                    }
169               } else if (refElement->numLocalDim==2) {               } else if (refElement->numLocalDim==2) {
170                    if (out->numSides==1) {                    if (out->numSides==1) {
171                       Assemble_jacobeans_2D(nodes->Coordinates,out->numQuadTotal,shape->QuadWeights,                       Finley_Assemble_jacobeans_2D(nodes->Coordinates,out->numQuadTotal,shape->QuadWeights,
172                                             shape->Type->numShapes,self->numElements,numNodes,self->Nodes,                                             shape->Type->numShapes,self->numElements,numNodes,self->Nodes,
173                                             shape->dSdv,basis->Type->numShapes,dBdv,                                             shape->dSdv,basis->Type->numShapes,dBdv,
174                                             out->DSDX,out->volume,self->Id);                                             out->DSDX,out->volume,self->Id);
# Line 186  Finley_ElementFile_Jacobeans* Finley_Ele Line 186  Finley_ElementFile_Jacobeans* Finley_Ele
186               } else if (refElement->numLocalDim==2) {               } else if (refElement->numLocalDim==2) {
187                    if (out->BasisFunctions->Type->numDim==3) {                    if (out->BasisFunctions->Type->numDim==3) {
188                       if (out->numSides==1) {                       if (out->numSides==1) {
189                          Assemble_jacobeans_3D_M2D_E3D(nodes->Coordinates,out->numQuadTotal,shape->QuadWeights,                          Finley_Assemble_jacobeans_3D_M2D_E3D(nodes->Coordinates,out->numQuadTotal,shape->QuadWeights,
190                                                        shape->Type->numShapes,self->numElements,numNodes,self->Nodes,                                                        shape->Type->numShapes,self->numElements,numNodes,self->Nodes,
191                                                        shape->dSdv,basis->Type->numShapes,dBdv,                                                        shape->dSdv,basis->Type->numShapes,dBdv,
192                                                        out->DSDX,out->volume,self->Id);                                                        out->DSDX,out->volume,self->Id);
193                       } else if (out->numSides==2) {                       } else if (out->numSides==2) {
194                          Assemble_jacobeans_3D_M2D_E3D_C(nodes->Coordinates,out->numQuadTotal,shape->QuadWeights,                          Finley_Assemble_jacobeans_3D_M2D_E3D_C(nodes->Coordinates,out->numQuadTotal,shape->QuadWeights,
195                                                          shape->Type->numShapes,self->numElements,numNodes,self->Nodes,                                                          shape->Type->numShapes,self->numElements,numNodes,self->Nodes,
196                                                          shape->dSdv,basis->Type->numShapes,dBdv,                                                          shape->dSdv,basis->Type->numShapes,dBdv,
197                                                          out->DSDX,out->volume,self->Id);                                                          out->DSDX,out->volume,self->Id);
# Line 200  Finley_ElementFile_Jacobeans* Finley_Ele Line 200  Finley_ElementFile_Jacobeans* Finley_Ele
200                        }                        }
201                    }  else if (out->BasisFunctions->Type->numDim==2) {                    }  else if (out->BasisFunctions->Type->numDim==2) {
202                       if (out->numSides==1) {                       if (out->numSides==1) {
203                          Assemble_jacobeans_3D_M2D_E2D(nodes->Coordinates,out->numQuadTotal,shape->QuadWeights,                          Finley_Assemble_jacobeans_3D_M2D_E2D(nodes->Coordinates,out->numQuadTotal,shape->QuadWeights,
204                                                        shape->Type->numShapes,self->numElements,numNodes,self->Nodes,                                                        shape->Type->numShapes,self->numElements,numNodes,self->Nodes,
205                                                        shape->dSdv,basis->Type->numShapes,dBdv,                                                        shape->dSdv,basis->Type->numShapes,dBdv,
206                                                        out->DSDX,out->volume,self->Id);                                                        out->DSDX,out->volume,self->Id);
207                       } else if (out->numSides==2) {                       } else if (out->numSides==2) {
208                          Assemble_jacobeans_3D_M2D_E2D_C(nodes->Coordinates,out->numQuadTotal,shape->QuadWeights,                          Finley_Assemble_jacobeans_3D_M2D_E2D_C(nodes->Coordinates,out->numQuadTotal,shape->QuadWeights,
209                                                          shape->Type->numShapes,self->numElements,numNodes,self->Nodes,                                                          shape->Type->numShapes,self->numElements,numNodes,self->Nodes,
210                                                          shape->dSdv,basis->Type->numShapes,dBdv,                                                          shape->dSdv,basis->Type->numShapes,dBdv,
211                                                          out->DSDX,out->volume,self->Id);                                                          out->DSDX,out->volume,self->Id);
# Line 217  Finley_ElementFile_Jacobeans* Finley_Ele Line 217  Finley_ElementFile_Jacobeans* Finley_Ele
217                    }                    }
218               } else if (refElement->numLocalDim==3) {               } else if (refElement->numLocalDim==3) {
219                    if (out->numSides==1) {                    if (out->numSides==1) {
220                       Assemble_jacobeans_3D(nodes->Coordinates,out->numQuadTotal,shape->QuadWeights,                       Finley_Assemble_jacobeans_3D(nodes->Coordinates,out->numQuadTotal,shape->QuadWeights,
221                                             shape->Type->numShapes,self->numElements,numNodes,self->Nodes,                                             shape->Type->numShapes,self->numElements,numNodes,self->Nodes,
222                                             shape->dSdv,basis->Type->numShapes,dBdv,                                             shape->dSdv,basis->Type->numShapes,dBdv,
223                                             out->DSDX,out->volume,self->Id);                                             out->DSDX,out->volume,self->Id);

Legend:
Removed from v.3258  
changed lines
  Added in v.3259

  ViewVC Help
Powered by ViewVC 1.1.26