| 226 | } |
| 227 | |
| 228 | void build_layer_map_edge(const Mesh &Th2, map< int, int > &mapemil, map< int, int > &mapezmax, |
| 229 | map< int, int > &mapezmin) { |
| 230 | int numero_label = 0; |
| 231 | |
| 232 | for (int ii = 0; ii < Th2.neb; ii++) { |
| 233 | const Mesh::BorderElement &K(Th2.be(ii)); |
| 234 | map< int, int >::const_iterator imap1 = mapezmax.find(K.lab); |
| 235 | map< int, int >::const_iterator imap2 = mapemil.find(K.lab); |
| 236 | map< int, int >::const_iterator imap3 = mapezmin.find(K.lab); |
| 237 | |
| 238 | if (imap1 == mapezmax.end( )) { |
| 239 | mapezmax[K.lab] = numero_label; |
| 240 | numero_label = numero_label + 1; |
| 241 | } |
| 242 | |
| 243 | if (imap2 == mapemil.end( )) { |
| 244 | mapemil[K.lab] = numero_label; |
| 245 | numero_label = numero_label + 1; |
| 246 | } |
| 247 | |
| 248 | if (imap3 == mapezmin.end( )) { |
| 249 | mapezmin[K.lab] = numero_label; |
| 250 | numero_label = numero_label + 1; |
| 251 | } |
| 252 | } |
| 253 | } |
| 254 | |
| 255 | Mesh3 *build_layer(const Mesh &Th2, const int Nmax, const int *tab_Ni, const double *tab_zmin, |
| 256 | const double *tab_zmax, const map< int, int > &maptet, |
no test coverage detected