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

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

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

revision 102 by jgs, Wed Dec 15 07:08:39 2004 UTC revision 123 by jgs, Fri Jul 8 04:08:13 2005 UTC
# Line 43  void Finley_Assemble_RobinCondition(Finl Line 43  void Finley_Assemble_RobinCondition(Finl
43    double *EM_S=NULL,*EM_F=NULL,*V=NULL,*dVdv=NULL,*dSdV=NULL,*Area=NULL;    double *EM_S=NULL,*EM_F=NULL,*V=NULL,*dVdv=NULL,*dSdV=NULL,*Area=NULL;
44    double time0;    double time0;
45    Assemble_Parameters p;    Assemble_Parameters p;
46    maybelong *index_row=NULL,*index_col=NULL;    index_t *index_row=NULL,*index_col=NULL,color;
47    int dimensions[ESCRIPT_MAX_DATA_RANK],e,q,color;    dim_t dimensions[ESCRIPT_MAX_DATA_RANK],e,q;
48    
49    if (nodes==NULL || elements==NULL) return;    if (nodes==NULL || elements==NULL) return;
50    if (S==NULL && isEmpty(F)) return;    if (S==NULL && isEmpty(F)) return;
# Line 54  void Finley_Assemble_RobinCondition(Finl Line 54  void Finley_Assemble_RobinCondition(Finl
54    if (Finley_ErrorCode!=NO_ERROR) return;    if (Finley_ErrorCode!=NO_ERROR) return;
55    
56    /*  get a functionspace */    /*  get a functionspace */
57    int funcspace=UNKNOWN;    type_t funcspace=UNKNOWN;
58    updateFunctionSpaceType(funcspace,d);    updateFunctionSpaceType(funcspace,d);
59    updateFunctionSpaceType(funcspace,y);    updateFunctionSpaceType(funcspace,y);
60    if (funcspace==UNKNOWN) return; /* all  data are empty */    if (funcspace==UNKNOWN) return; /* all  data are empty */
# Line 128  void Finley_Assemble_RobinCondition(Finl Line 128  void Finley_Assemble_RobinCondition(Finl
128          V=THREAD_MEMALLOC(p.NN*p.numDim,double);          V=THREAD_MEMALLOC(p.NN*p.numDim,double);
129          dVdv=THREAD_MEMALLOC(p.numDim*p.numElementDim*p.numQuad,double);          dVdv=THREAD_MEMALLOC(p.numDim*p.numElementDim*p.numQuad,double);
130          Area=THREAD_MEMALLOC(p.numQuad,double);          Area=THREAD_MEMALLOC(p.numQuad,double);
131          index_row=THREAD_MEMALLOC(p.NN_row,maybelong);          index_row=THREAD_MEMALLOC(p.NN_row,index_t);
132          index_col=THREAD_MEMALLOC(p.NN_col,maybelong);          index_col=THREAD_MEMALLOC(p.NN_col,index_t);
133          if (! ( Finley_checkPtr(EM_S) || Finley_checkPtr(EM_F) || Finley_checkPtr(index_row) || Finley_checkPtr(index_col) ||          if (! ( Finley_checkPtr(EM_S) || Finley_checkPtr(EM_F) || Finley_checkPtr(index_row) || Finley_checkPtr(index_col) ||
134                  Finley_checkPtr(V) || Finley_checkPtr(dVdv) || Finley_checkPtr(Area))) {                  Finley_checkPtr(V) || Finley_checkPtr(dVdv) || Finley_checkPtr(Area))) {
135             /*  open loop over all colors: */             /*  open loop over all colors: */
136             for (color=0;color<elements->numColors;color++) {             for (color=elements->minColor;color<=elements->maxColor;color++) {
137                /*  open loop over all elements: */                /*  open loop over all elements: */
138                #pragma omp for private(e,q) schedule(static)                #pragma omp for private(e,q) schedule(static)
139                for(e=0;e<elements->numElements;e++){                for(e=0;e<elements->numElements;e++){
# Line 223  void Finley_Assemble_RobinCondition(Finl Line 223  void Finley_Assemble_RobinCondition(Finl
223  }  }
224  /*  /*
225   * $Log$   * $Log$
226     * Revision 1.5  2005/07/08 04:07:46  jgs
227     * Merge of development branch back to main trunk on 2005-07-08
228     *
229   * Revision 1.4  2004/12/15 07:08:32  jgs   * Revision 1.4  2004/12/15 07:08:32  jgs
230   * *** empty log message ***   * *** empty log message ***
231     * Revision 1.1.1.1.2.2  2005/06/29 02:34:47  gross
232     * some changes towards 64 integers in finley
233     *
234     * Revision 1.1.1.1.2.1  2004/11/24 01:37:12  gross
235     * some changes dealing with the integer overflow in memory allocation. Finley solves 4M unknowns now
236   *   *
237   *   *
238   *   *

Legend:
Removed from v.102  
changed lines
  Added in v.123

  ViewVC Help
Powered by ViewVC 1.1.26