Ignore:
Timestamp:
06/03/10 07:37:08 (15 years ago)
Author:
Mathieu Morlighem
Message:

Updated CreateNodes

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/c/modules/ModelProcessorx/SlopeCompute/CreateNodesSlopeConpute.cpp

    r3999 r4005  
    1717        /*Intermediary*/
    1818        int i;
    19         int totalnodes;
    2019        bool continuous_galerkin=true;
    2120
    2221        /*DataSets: */
    2322        DataSet*    nodes = NULL;
     23
     24        /*Recover pointer: */
     25        nodes=*pnodes;
     26
     27        /*Create nodes if they do not exist yet*/
     28        if(!nodes) nodes = new DataSet(NodesEnum);
    2429       
    25         /*First create nodes*/
    26         nodes = new DataSet(NodesEnum);
    27 
    28         /*Recover number of nodes already created in other analyses: */
    29         totalnodes=iomodel->nodecounter;
    30 
    3130        /*Continuous Galerkin partition of nodes: */
    3231        NodesPartitioning(&iomodel->my_nodes,iomodel->my_elements, iomodel->my_vertices, iomodel->my_bordervertices, iomodel, iomodel_handle,continuous_galerkin);
     
    4746                       
    4847                        /*Add node to nodes dataset: */
    49                         nodes->AddObject(new Node(i+1,i,iomodel));
     48                        nodes->AddObject(new Node(iomodel->nodecounter+i+1,i,iomodel));
    5049
    5150                }
     
    6059        xfree((void**)&iomodel->gridoniceshelf);
    6160       
    62         /*All our datasets are already order by ids. Set presort flag so that later on, when sorting is requested on these
    63          * datasets, it will not be redone: */
    64         nodes->Presort();
    65 
    6661        /*Assign output pointer: */
    6762        *pnodes=nodes;
    68 
    6963}
Note: See TracChangeset for help on using the changeset viewer.