/[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 97 by jgs, Tue Dec 14 05:39:33 2004 UTC revision 100 by jgs, Wed Dec 15 03:48:48 2004 UTC
# Line 98  void Finley_Assemble_interpolate(Finley_ Line 98  void Finley_Assemble_interpolate(Finley_
98         {         {
99            local_data=NULL;            local_data=NULL;
100            /* allocation of work arrays */            /* allocation of work arrays */
101            local_data=THREAD_MEMALLOC(NS*numComps,double);            local_data=(double*) THREAD_MEMALLOC(NS*numComps*sizeof(double));
102            if (! Finley_checkPtr(local_data)) {            if (! Finley_checkPtr(local_data)) {
103    
104          /* open the element loop */          /* open the element loop */
# Line 112  void Finley_Assemble_interpolate(Finley_ Line 112  void Finley_Assemble_interpolate(Finley_
112                             i=elements->Nodes[INDEX2(resort_nodes[dof_offset+q],e,NN)];                             i=elements->Nodes[INDEX2(resort_nodes[dof_offset+q],e,NN)];
113                             data_array=getSampleData(data,i);                             data_array=getSampleData(data,i);
114                             Finley_copyDouble(numComps,data_array,local_data+q*numComps);                             Finley_copyDouble(numComps,data_array,local_data+q*numComps);
                         }  
115                          break;                          break;
116                   case DOF:                   case DOF:
117                          for (q=0;q<NS_DOF;q++) {                          for (q=0;q<NS_DOF;q++) {
# Line 138  void Finley_Assemble_interpolate(Finley_ Line 137  void Finley_Assemble_interpolate(Finley_
137    
138            }            }
139        THREAD_MEMFREE(local_data);        THREAD_MEMFREE(local_data);
140       } /* end of parallel region */          }
141          }
142    }    }
143    #undef NODES    #undef NODES
144    #undef DOF    #undef DOF
# Line 146  void Finley_Assemble_interpolate(Finley_ Line 146  void Finley_Assemble_interpolate(Finley_
146  }  }
147  /*  /*
148   * $Log$   * $Log$
149   * Revision 1.2  2004/12/14 05:39:30  jgs   * Revision 1.3  2004/12/15 03:48:45  jgs
150   * *** empty log message ***   * *** empty log message ***
151   *   *
  * Revision 1.1.1.1.2.2  2004/11/24 01:37:12  gross  
  * some changes dealing with the integer overflow in memory allocation. Finley solves 4M unknowns now  
  *  
  * Revision 1.1.1.1.2.1  2004/10/27 23:39:21  gross  
  * bug in finley interpolation fixed. GradTest is running now  
  *  
152   * Revision 1.1.1.1  2004/10/26 06:53:57  jgs   * Revision 1.1.1.1  2004/10/26 06:53:57  jgs
153   * initial import of project esys2   * initial import of project esys2
154   *   *

Legend:
Removed from v.97  
changed lines
  Added in v.100

  ViewVC Help
Powered by ViewVC 1.1.26