Ignore:
Timestamp:
01/23/19 10:31:44 (6 years ago)
Author:
yushan
Message:

dev on ADA. add flag switch _usingEP/_usingMPI

File:
1 edited

Legend:

Unmodified
Added
Removed
  • XIOS/dev/branch_openmp/extern/remap/src/libmapper.cpp

    r1335 r1642  
    1515#include "cputime.hpp" // cputime 
    1616 
    17 using namespace ep_lib; 
    18  
    1917using namespace sphereRemap ; 
    2018 
     
    2220   and deallocated during the second step (computing the weights) */ 
    2321Mapper *mapper; 
    24 #pragma omp threadprivate(mapper) 
     22 
    2523 
    2624/** xxx_bounds_yyy is of length n_vert_per_cell_xxx*n_cell_xxx 
     
    4341        assert(n_cell_dst >= 4); 
    4442        assert(1 <= order && order <= 2); 
    45  
    46   mapper = new Mapper(MPI_COMM_WORLD); 
     43  double* src_area=NULL ; 
     44  double* dst_area=NULL ; 
     45  mapper = new Mapper(EP_COMM_WORLD); 
    4746  mapper->setVerbosity(PROGRESS) ; 
    48   mapper->setSourceMesh(src_bounds_lon, src_bounds_lat, n_vert_per_cell_src, n_cell_src, src_pole ) ; 
    49   mapper->setTargetMesh(dst_bounds_lon, dst_bounds_lat, n_vert_per_cell_dst, n_cell_dst, dst_pole ) ; 
     47  mapper->setSourceMesh(src_bounds_lon, src_bounds_lat, src_area, n_vert_per_cell_src, n_cell_src, src_pole ) ; 
     48  mapper->setTargetMesh(dst_bounds_lon, dst_bounds_lat, dst_area, n_vert_per_cell_dst, n_cell_dst, dst_pole ) ; 
    5049 
    5150/* 
Note: See TracChangeset for help on using the changeset viewer.