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

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

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

revision 3220 by jfenwick, Thu Sep 23 06:08:56 2010 UTC revision 3221 by jfenwick, Wed Sep 29 01:00:21 2010 UTC
# Line 47  void Dudley_Assemble_getSize(Dudley_Node Line 47  void Dudley_Assemble_getSize(Dudley_Node
47    // now we look up what type of elements we need based on the function space of element_size    // now we look up what type of elements we need based on the function space of element_size
48    // if it is DUDLEY_REDUCED_ELEMENTS or DUDLEY_REDUCED_FACE_ELEMENTS then we have single quad point    // if it is DUDLEY_REDUCED_ELEMENTS or DUDLEY_REDUCED_FACE_ELEMENTS then we have single quad point
49    
 //  refElement=Dudley_ReferenceElementSet_borrowReferenceElement(elements->referenceElementSet,Dudley_Assemble_reducedIntegrationOrder(element_size));  
 //  shape= refElement->BasisFunctions;  
 //  numQuad=shape->numQuadNodes;  
   
50    if (Dudley_Assemble_reducedIntegrationOrder(element_size))    if (Dudley_Assemble_reducedIntegrationOrder(element_size))
51    {    {
52      numQuad=1;      numQuad=1;
# Line 61  void Dudley_Assemble_getSize(Dudley_Node Line 57  void Dudley_Assemble_getSize(Dudley_Node
57    }    }
58    
59    NN=elements->numNodes;    NN=elements->numNodes;
 //  NS=refElement->BasisFunctions->Type->numShapes;  
60    NS=elements->numDim+1;    NS=elements->numDim+1;
 //  NVertices=refElement->BasisFunctions->Type->numVertices;  
61    NVertices=elements->numDim+1;    NVertices=elements->numDim+1;
62        
     
63    /* check the dimensions of element_size */    /* check the dimensions of element_size */
64    
65    if (! numSamplesEqual(element_size,numQuad,elements->numElements)) {    if (! numSamplesEqual(element_size,numQuad,elements->numElements)) {

Legend:
Removed from v.3220  
changed lines
  Added in v.3221

  ViewVC Help
Powered by ViewVC 1.1.26