Changeset 8260 for trunk/anuga_core/source/anuga_parallel/parallel_api.py
- Timestamp:
- Nov 30, 2011, 7:57:30 PM (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/anuga_core/source/anuga_parallel/parallel_api.py
r8114 r8260 20 20 # Mesh partitioning using Metis 21 21 from anuga_parallel.distribute_mesh import build_submesh 22 from anuga_parallel.distribute_mesh import pmesh_divide_metis 22 from anuga_parallel.distribute_mesh import pmesh_divide_metis_with_map 23 23 24 24 from anuga_parallel.parallel_shallow_water import Parallel_domain … … 97 97 points, vertices, boundary, quantities,\ 98 98 ghost_recv_dict, full_send_dict,\ 99 number_of_full_nodes, number_of_full_triangles =\ 99 number_of_full_nodes, number_of_full_triangles,\ 100 s2p_map, p2s_map =\ 100 101 distribute_mesh(domain, verbose=verbose) 101 102 103 # Send serial to parallel (s2p) and parallel to serial (p2s) triangle mapping to proc 1 .. numprocs 104 for p in range(1, numprocs): 105 send(s2p_map, p) 106 send(p2s_map, p) 102 107 103 108 if verbose: print 'Communication done' … … 112 117 rec_submesh(0, verbose) 113 118 119 # Recieve serial to parallel (s2p) and parallel to serial (p2s) triangle mapping 120 s2p_map = receive(0) 121 p2s_map = receive(0) 122 114 123 115 124 #------------------------------------------------------------------------ … … 125 134 number_of_full_nodes=number_of_full_nodes, 126 135 number_of_full_triangles=number_of_full_triangles, 127 geo_reference=georef) ## jj added this 136 geo_reference=georef, 137 tri_map = s2p_map, 138 inv_tri_map = p2s_map) ## jj added this 128 139 129 140 #------------------------------------------------------------------------ … … 165 176 # Subdivide the mesh 166 177 if verbose: print 'Subdivide mesh' 167 nodes, triangles, boundary, triangles_per_proc, quantities = \ 168 pmesh_divide_metis(domain, numprocs) 169 170 178 nodes, triangles, boundary, triangles_per_proc, quantities, s2p_map, p2s_map = \ 179 pmesh_divide_metis_with_map(domain, numprocs) 180 181 #PETE: s2p_map (maps serial domain triangles to parallel domain triangles) 182 #p2_map (maps parallel domain triangles to domain triangles) 171 183 172 184 … … 227 239 return points, vertices, boundary, quantities,\ 228 240 ghost_recv_dict, full_send_dict,\ 229 number_of_full_nodes, number_of_full_triangles 230 231 232 233 241 number_of_full_nodes, number_of_full_triangles, s2p_map, p2s_map 242 243 244 245
Note: See TracChangeset
for help on using the changeset viewer.