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

Diff of /trunk/finley/src/Assemble_setNormal.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 45  void Finley_Assemble_setNormal(Finley_No Line 45  void Finley_Assemble_setNormal(Finley_No
45        sign=1;        sign=1;
46    }    }
47    /* check the dimensions of normal */    /* check the dimensions of normal */
48    if (! (numDim==numDim_local || numDim-1==numDim_local)) {    if (numDim==numDim_local || numDim==numDim_local-1) {
49         Finley_ErrorCode=TYPE_ERROR;         Finley_ErrorCode=TYPE_ERROR;
50         sprintf(Finley_ErrorMsg,"Cannot calculate normal vector");         sprintf(Finley_ErrorMsg,"Cannot calculate normal vector");
51    } else if (! isDataPointShapeEqual(normal,1,&(numDim))) {    } else if (! isDataPointShapeEqual(normal,1,&(numDim))) {
# Line 68  void Finley_Assemble_setNormal(Finley_No Line 68  void Finley_Assemble_setNormal(Finley_No
68            {            {
69               local_X=dVdv=NULL;               local_X=dVdv=NULL;
70               /* allocation of work arrays */               /* allocation of work arrays */
71               local_X=THREAD_MEMALLOC(NS*numDim,double);               local_X=(double*) THREAD_MEMALLOC(NS*numDim*sizeof(double));
72               dVdv=THREAD_MEMALLOC(numQuad*numDim*numDim_local,double);               dVdv=(double*) THREAD_MEMALLOC(numQuad*numDim*numDim_local*sizeof(double));
73               if (!(Finley_checkPtr(local_X) || Finley_checkPtr(dVdv) ) ) {               if (!(Finley_checkPtr(local_X) || Finley_checkPtr(dVdv) ) ) {
74                         /* open the element loop */                         /* open the element loop */
75                         #pragma omp for private(e,q,normal_array) schedule(static)                         #pragma omp for private(e,q,normal_array) schedule(static)
# Line 91  void Finley_Assemble_setNormal(Finley_No Line 91  void Finley_Assemble_setNormal(Finley_No
91  }  }
92  /*  /*
93   * $Log$   * $Log$
94   * Revision 1.2  2004/12/14 05:39:30  jgs   * Revision 1.3  2004/12/15 03:48:45  jgs
95   * *** empty log message ***   * *** empty log message ***
96   *   *
  * Revision 1.1.1.1.2.2  2004/11/24 01:37:13  gross  
  * some changes dealing with the integer overflow in memory allocation. Finley solves 4M unknowns now  
  *  
  * Revision 1.1.1.1.2.1  2004/11/15 00:12:42  gross  
  * small bug fixed  
  *  
97   * Revision 1.1.1.1  2004/10/26 06:53:57  jgs   * Revision 1.1.1.1  2004/10/26 06:53:57  jgs
98   * initial import of project esys2   * initial import of project esys2
99   *   *

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

  ViewVC Help
Powered by ViewVC 1.1.26