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

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

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

revision 1063 by gross, Wed Mar 14 00:15:24 2007 UTC revision 1064 by gross, Tue Mar 27 06:21:02 2007 UTC
# Line 40  void Finley_Assemble_integrate(Finley_No Line 40  void Finley_Assemble_integrate(Finley_No
40      Finley_resetError();      Finley_resetError();
41      if (nodes==NULL || elements==NULL) return;      if (nodes==NULL || elements==NULL) return;
42      /* set some parameter */      /* set some parameter */
43                                                                                                                                                      jac=Finley_ElementFile_borrowJacobeans(elements,nodes,FALSE,Finley_Assemble_reducedIntegrationOrder(data));
     if (data_type==FINLEY_ELEMENTS) {  
         reducedIntegrationOrder=FALSE;  
     } else if (data_type==FINLEY_FACE_ELEMENTS)  {  
         reducedIntegrationOrder=FALSE;  
     } else if (data_type==FINLEY_CONTACT_ELEMENTS_1)  {  
         reducedIntegrationOrder=FALSE;  
     } else if (data_type==FINLEY_CONTACT_ELEMENTS_2)  {  
         reducedIntegrationOrder=FALSE;  
     } else {  
        Finley_setError(TYPE_ERROR,"Finley_Assemble_integrate: integration of data is not possible.");  
     }  
   
     /* get access to jacobean */  
     jac=Finley_ElementFile_borrowJacobeans(elements,nodes,FALSE,reducedIntegrationOrder);  
   
44      if (Finley_noError()) {      if (Finley_noError()) {
45    
46          /* check the shape of the data  */          /* check the shape of the data  */
47          if (! numSamplesEqual(data,jac->ReferenceElement->numQuadNodes,elements->numElements)) {          if (! numSamplesEqual(data,jac->ReferenceElement->numQuadNodes,elements->numElements)) {
48             Finley_setError(TYPE_ERROR,"Finley_Assemble_integrate: illegal number of samples of integrant kernel Data object");             Finley_setError(TYPE_ERROR,"Finley_Assemble_integrate: illegal number of samples of integrant kernel Data object");
49          }          }
   
   
50          /* now we can start */          /* now we can start */
51    
52          if (Finley_noError()) {          if (Finley_noError()) {
Legend:
Removed from v.1063  
changed lines
  Added in v.1064
  ViewVC Help
Powered by ViewVC 1.1.26