/[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 3195 by jfenwick, Fri Sep 3 05:48:31 2010 UTC revision 3196 by jfenwick, Wed Sep 22 01:18:52 2010 UTC
# Line 24  Line 24 
24  #include <omp.h>  #include <omp.h>
25  #endif  #endif
26    
27    #include "ShapeTable.h"
28    
29  /**************************************************************/  /**************************************************************/
30    
31    
# Line 32  void Dudley_Assemble_setNormal(Dudley_No Line 34  void Dudley_Assemble_setNormal(Dudley_No
34    index_t sign;    index_t sign;
35    Dudley_ReferenceElement* reference_element=NULL;    Dudley_ReferenceElement* reference_element=NULL;
36    dim_t e,q, NN, NS, numDim, numQuad, numDim_local;    dim_t e,q, NN, NS, numDim, numQuad, numDim_local;
37      bool_t reduced_integration;
38      const double* dSdv=0;
39    if (nodes==NULL || elements==NULL) return;    if (nodes==NULL || elements==NULL) return;
40      switch (elements->numDim)
41      {
42      case 2: dSdv=&(DTDV_2D[0][0]); break;
43      case 3: dSdv=&(DTDV_3D[0][0]); break;
44      default:
45        dSdv=&(DTDV_1D[0][0]); break;
46      }
47    Dudley_resetError();    Dudley_resetError();
48    reference_element=Dudley_ReferenceElementSet_borrowReferenceElement(elements->referenceElementSet, Dudley_Assemble_reducedIntegrationOrder(normal));    reference_element=Dudley_ReferenceElementSet_borrowReferenceElement(elements->referenceElementSet, Dudley_Assemble_reducedIntegrationOrder(normal));
49    NN=elements->numNodes;    NN=elements->numNodes;
50    numDim=nodes->numDim;    numDim=nodes->numDim;
51      reduced_integration = Dudley_Assemble_reducedIntegrationOrder(normal);
52    numQuad=reference_element->BasisFunctions->numQuadNodes;  //  numQuad=reference_element->BasisFunctions->numQuadNodes;
53    numDim_local=reference_element->BasisFunctions->Type->numDim;    numQuad=(!reduced_integration)?(elements->numDim+1):1;
54    NS=reference_element->BasisFunctions->Type->numShapes;  //  numDim_local=reference_element->BasisFunctions->Type->numDim;
55      numDim_local=elements->numDim;
56    //  NS=reference_element->BasisFunctions->Type->numShapes;
57      NS=elements->numDim+1;
58    //fprintf(stderr,"\nnumQuad=%d,%d numDim_local=%d,%d NS=%d,%d\n", numQuad, ((!reduced_integration)?(elements->numDim+1):1), numDim_local, elements->numDim, NS, numQuad);
59    
60        
61    /* set some parameter */    /* set some parameters */
62    
63    sign=1;    sign=1;
64    /* check the dimensions of normal */    /* check the dimensions of normal */

Legend:
Removed from v.3195  
changed lines
  Added in v.3196

  ViewVC Help
Powered by ViewVC 1.1.26