Changeset 3437


Ignore:
Timestamp:
04/07/10 16:52:49 (15 years ago)
Author:
Mathieu Morlighem
Message:

Prognostic2 is all up-to-date

Location:
issm/trunk/src/c/ModelProcessorx/Prognostic2
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/c/ModelProcessorx/Prognostic2/CreateLoadsPrognostic2.cpp

    r3378 r3437  
    1313void    CreateLoadsPrognostic2(DataSet** ploads, IoModel* iomodel,ConstDataHandle iomodel_handle){
    1414
     15        /*Intermediary*/
    1516        int i,j;
    1617        int i1,i2;
     
    1819        double e1,e2;
    1920
    20         extern int my_rank;
    21         extern int num_procs;
    22 
     21        /*Output*/
    2322        DataSet*       loads        = NULL;
    24         Numericalflux* numericalflux= NULL;
    2523
    2624        /*numericalflux intermediary data: */
     
    4442                e2=iomodel->edges[4*i+3]-1; //edges are [node1 node2 elem1 elem2]
    4543
    46                 #ifdef _PARALLEL_
    47                 if (iomodel->epart[(int)e1]!=my_rank){
    48                         /*This load does not belong to this cluster node, as it references an element
    49                          *that does not belong to this node's partition. Drop this 'i':*/
    50                         continue;
    51                 }
    52                 #endif
     44                /*Now, if this element is not in the partition, pass: */
     45                if(!iomodel->my_elements[i])continue;
    5346
    5447                /*Create load*/
     
    10396                }
    10497
    105                 numericalflux = new Numericalflux(numericalflux_id,numericalflux_type,numericalflux_node_ids,numericalflux_elem_id);
    106 
    107                 loads->AddObject(numericalflux);
     98                loads->AddObject(new Numericalflux(numericalflux_id,numericalflux_type,numericalflux_node_ids,numericalflux_elem_id));
    10899        }
    109100
  • issm/trunk/src/c/ModelProcessorx/Prognostic2/CreateParametersPrognostic2.cpp

    r3382 r3437  
    4444                for(j=0;j<3;j++){
    4545                        part[3*i+j]=(int)elements[3*i+j]-1; //Matlab to C indexing
    46                         ISSMASSERT(part[3*i+j]<iomodel->numberofnodes);
     46                        ISSMASSERT(part[3*i+j]<iomodel->numberofvertices);
    4747                }
    4848        }
Note: See TracChangeset for help on using the changeset viewer.