/[escript]/trunk/dudley/src/Assemble_gradient.c
ViewVC logotype

Diff of /trunk/dudley/src/Assemble_gradient.c

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

revision 3489 by jfenwick, Mon Oct 11 01:48:14 2010 UTC revision 3490 by caltinay, Wed Mar 30 02:24:33 2011 UTC
# Line 33  void Dudley_Assemble_gradient(Dudley_Nod Line 33  void Dudley_Assemble_gradient(Dudley_Nod
33      register dim_t e, q, l, s, n;      register dim_t e, q, l, s, n;
34      register __const double *data_array;      register __const double *data_array;
35      register double *grad_data_e;      register double *grad_data_e;
36      dim_t numNodes = 0, numShapesTotal = 0, numComps, NN = 0, numDim = 0, numShapesTotal2 = 0, numQuad = 0;      dim_t numNodes = 0, numShapesTotal = 0, numComps, NN = 0, numDim = 0, numQuad = 0;
37      type_t data_type = getFunctionSpaceType(data);      type_t data_type = getFunctionSpaceType(data);
38      bool_t reducedIntegrationOrder = FALSE;      bool_t reducedIntegrationOrder = FALSE;
39      Dudley_ElementFile_Jacobeans *jac = NULL;      Dudley_ElementFile_Jacobeans *jac = NULL;
# Line 86  void Dudley_Assemble_gradient(Dudley_Nod Line 86  void Dudley_Assemble_gradient(Dudley_Nod
86      numShapesTotal = jac->numShapes;      numShapesTotal = jac->numShapes;
87      numQuad = jac->numQuad;      numQuad = jac->numQuad;
88      localGradSize = sizeof(double) * numDim * numQuad * numComps;      localGradSize = sizeof(double) * numDim * numQuad * numComps;
     numShapesTotal2 = elements->numShapes;  
89      /* check the dimensions of data */      /* check the dimensions of data */
90    
91      if (!numSamplesEqual(grad_data, numQuad, elements->numElements))      if (!numSamplesEqual(grad_data, numQuad, elements->numElements))

Legend:
Removed from v.3489  
changed lines
  Added in v.3490

  ViewVC Help
Powered by ViewVC 1.1.26