Ignore:
Timestamp:
08/23/18 11:09:21 (7 years ago)
Author:
Mathieu Morlighem
Message:

CHG: moving my_edges and my_faces to iomodel to avoid recomputation

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/c/modules/ModelProcessorx/EdgesPartitioning.cpp

    r17692 r23167  
    88#include "./ModelProcessorx.h"
    99
    10 void EdgesPartitioning(bool** pmy_edges,IoModel* iomodel){
     10void EdgesPartitioning(IoModel* iomodel){
    1111
    12         /*Intermediaries*/
    13         int elementnbe;
     12        /*If faces are already present, exit*/
     13        if(iomodel->my_edges) return;
    1414
    1515        /*Get edges and elements*/
     
    1818
    1919        /*Mesh dependent variables*/
     20        int elementnbe;
    2021        switch(iomodel->meshelementtype){
    2122                case TriaEnum:  elementnbe = 3; break;
     
    2627
    2728        /*output: */
    28         bool* my_edges=xNewZeroInit<bool>(iomodel->numberofedges);
     29        iomodel->my_edges=xNewZeroInit<bool>(iomodel->numberofedges);
    2930
    3031        for(int i=0;i<iomodel->numberofelements;i++){
    3132                if(iomodel->my_elements[i]){
    3233                        for(int j=0;j<elementnbe;j++){
    33                                 my_edges[iomodel->elementtoedgeconnectivity[i*elementnbe+j]] = true;
     34                                iomodel->my_edges[iomodel->elementtoedgeconnectivity[i*elementnbe+j]] = true;
    3435                        }
    3536                }
    3637        }
    37 
    38         /*Free data and assign output pointers */
    39         *pmy_edges=my_edges;
    4038}
Note: See TracChangeset for help on using the changeset viewer.