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

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

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

revision 616 by elspeth, Wed Mar 22 02:46:56 2006 UTC revision 1028 by gross, Wed Mar 14 00:15:24 2007 UTC
# Line 32  Line 32 
32    
33  void Finley_Assemble_interpolate(Finley_NodeFile *nodes, Finley_ElementFile* elements,escriptDataC* data,escriptDataC* interpolated_data) {  void Finley_Assemble_interpolate(Finley_NodeFile *nodes, Finley_ElementFile* elements,escriptDataC* data,escriptDataC* interpolated_data) {
34    double* local_data=NULL,*S=NULL,*data_array;    double* local_data=NULL,*S=NULL,*data_array;
35    index_t dof_offset,*resort_nodes;    index_t dof_offset,*resort_nodesi, NN, NS;
36    dim_t q,i,NS_DOF,NN_DOF,numNodes,e;    dim_t q,i,NS_DOF,NN_DOF,numNodes,e, numQuad;
37      dim_t numComps=getDataPointSize(data);
38      index_t id[MAX_numNodes], *resort_nodes;
39      type_t data_type=getFunctionSpaceType(data);
40    type_t type;    type_t type;
41      Finley_resetError();
42    #define NODES 0    #define NODES 0
43    #define DOF 1    #define DOF 1
44    #define REDUCED_DOF 2    #define REDUCED_DOF 2
45    if (nodes==NULL || elements==NULL) return;    if (nodes==NULL || elements==NULL) return;
46    dim_t NN=elements->ReferenceElement->Type->numNodes;    NN=elements->ReferenceElement->Type->numNodes;
47    dim_t NS=elements->ReferenceElement->Type->numShapes;    NS=elements->ReferenceElement->Type->numShapes;
48    dim_t numComps=getDataPointSize(data);    numQuad=elements->ReferenceElement->numQuadNodes;
   type_t data_type=getFunctionSpaceType(data);  
   dim_t numQuad=elements->ReferenceElement->numQuadNodes;  
   index_t id[NN];  
49    for (i=0;i<NN;i++) id[i]=i;    for (i=0;i<NN;i++) id[i]=i;
   Finley_resetError();  
50    
51    /* set some parameter */    /* set some parameter */
52    

Legend:
Removed from v.616  
changed lines
  Added in v.1028

  ViewVC Help
Powered by ViewVC 1.1.26