/[escript]/trunk/esys2/finley/src/finleyC/Assemble_gradient.c
ViewVC logotype

Diff of /trunk/esys2/finley/src/finleyC/Assemble_gradient.c

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

revision 96 by jgs, Tue Oct 26 06:53:54 2004 UTC revision 97 by jgs, Tue Dec 14 05:39:33 2004 UTC
# Line 104  void Finley_Assemble_gradient(Finley_Nod Line 104  void Finley_Assemble_gradient(Finley_Nod
104            {            {
105                 local_X=local_data=dVdv=dvdV=Vol=d_datadv=NULL;                 local_X=local_data=dVdv=dvdV=Vol=d_datadv=NULL;
106            /* allocation of work arrays */            /* allocation of work arrays */
107            local_X=(double*) THREAD_MEMALLOC(NS*numDim*sizeof(double));            local_X=THREAD_MEMALLOC(NS*numDim,double);
108            local_data=(double*) THREAD_MEMALLOC(NS*numComps*sizeof(double));            local_data=THREAD_MEMALLOC(NS*numComps,double);
109            dVdv=(double*) THREAD_MEMALLOC(numQuad*numDim*numDim*sizeof(double));            dVdv=THREAD_MEMALLOC(numQuad*numDim*numDim,double);
110            dvdV=(double*) THREAD_MEMALLOC(numQuad*numDim*numDim*sizeof(double));            dvdV=THREAD_MEMALLOC(numQuad*numDim*numDim,double);
111            Vol=(double*) THREAD_MEMALLOC(numQuad*sizeof(double));            Vol=THREAD_MEMALLOC(numQuad,double);
112            d_datadv=(double*) THREAD_MEMALLOC(numQuad*numComps*numDim*sizeof(double));            d_datadv=THREAD_MEMALLOC(numQuad*numComps*numDim,double);
113            if (!(Finley_checkPtr(local_X) || Finley_checkPtr(dVdv) || Finley_checkPtr(dvdV) || Finley_checkPtr(Vol) || Finley_checkPtr(d_datadv) || Finley_checkPtr(local_data) ))  {            if (!(Finley_checkPtr(local_X) || Finley_checkPtr(dVdv) || Finley_checkPtr(dvdV) || Finley_checkPtr(Vol) || Finley_checkPtr(d_datadv) || Finley_checkPtr(local_data) ))  {
114              /* open the element loop */              /* open the element loop */
115                  #pragma omp for private(e,i,q,data_array) schedule(static)                  #pragma omp for private(e,i,q,data_array) schedule(static)
# Line 165  void Finley_Assemble_gradient(Finley_Nod Line 165  void Finley_Assemble_gradient(Finley_Nod
165  }  }
166  /*  /*
167   * $Log$   * $Log$
168   * Revision 1.1  2004/10/26 06:53:57  jgs   * Revision 1.2  2004/12/14 05:39:30  jgs
169   * Initial revision   * *** empty log message ***
170     *
171     * Revision 1.1.1.1.2.1  2004/11/24 01:37:12  gross
172     * some changes dealing with the integer overflow in memory allocation. Finley solves 4M unknowns now
173     *
174     * Revision 1.1.1.1  2004/10/26 06:53:57  jgs
175     * initial import of project esys2
176   *   *
177   * Revision 1.5  2004/08/26 12:03:52  gross   * Revision 1.5  2004/08/26 12:03:52  gross
178   * Some other bug in Finley_Assemble_gradient fixed.   * Some other bug in Finley_Assemble_gradient fixed.

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

  ViewVC Help
Powered by ViewVC 1.1.26