![(please configure the [header_logo] section in trac.ini)](/trac/issm/chrome/common/trac_banner.png)
Changeset 27120
- Timestamp:
- 06/29/22 06:34:12 (3 years ago)
- File:
-
- 1 edited
-
issm/trunk-jpl/src/m/mesh/bamg.py (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/mesh/bamg.py
r26358 r27120 23 23 Available options (for more details see ISSM website http://issm.jpl.nasa.gov/): 24 24 25 - domain : followed by an ARGUS file that 25 - domain : followed by an ARGUS file that 26 26 prescribes the domain outline 27 - holes : followed by an ARGUS file that 27 - holes : followed by an ARGUS file that 28 28 prescribes the holes 29 - subdomains : followed by an ARGUS file that 30 prescribes the list of subdomains 29 - subdomains : followed by an ARGUS file that 30 prescribes the list of subdomains 31 31 (that need to be inside domain) 32 32 33 - hmin : minimum edge length (default is 33 - hmin : minimum edge length (default is 34 34 1.0e-100) 35 - hmax : maximum edge length (default is 35 - hmax : maximum edge length (default is 36 36 1.0e100) 37 - hVertices : imposed edge length for each vertex 37 - hVertices : imposed edge length for each vertex 38 38 (geometry or mesh) 39 - hminVertices : minimum edge length for each vertex 39 - hminVertices : minimum edge length for each vertex 40 40 (mesh) 41 - hmaxVertices : maximum edge length for each vertex 41 - hmaxVertices : maximum edge length for each vertex 42 42 (mesh) 43 43 44 - anisomax : maximum ratio between the smallest 45 and largest edges (default is 44 - anisomax : maximum ratio between the smallest 45 and largest edges (default is 46 46 1.0e30) 47 - coeff : coefficient applied to the metric 48 (2 -> twice as many elements, 47 - coeff : coefficient applied to the metric 48 (2 -> twice as many elements, 49 49 default is 1) 50 - cutoff : scalar used to compute the metric 50 - cutoff : scalar used to compute the metric 51 51 when metric type 2 or 3 are applied 52 - err : error used to generate the metric 52 - err : error used to generate the metric 53 53 from a field 54 54 - errg : geometric error (default is 0.1) 55 - field : field of the model that will be 56 used to compute the metric to apply 57 several fields, use one column per 55 - field : field of the model that will be 56 used to compute the metric to apply 57 several fields, use one column per 58 58 field 59 - gradation : maximum ratio between two adjacent 59 - gradation : maximum ratio between two adjacent 60 60 edges 61 - Hessiantype : 0 -> use double P2 projection 61 - Hessiantype : 0 -> use double P2 projection 62 62 (default) 63 63 1 -> use Green formula 64 - KeepVertices : try to keep initial vertices when 65 adaptation is done on an existing 64 - KeepVertices : try to keep initial vertices when 65 adaptation is done on an existing 66 66 mesh (default 1) 67 - NoBoundaryRefinement : do not refine boundary, only follow 68 contour provided (default 0). Allow 69 subdomain boundary refinement 67 - NoBoundaryRefinement : do not refine boundary, only follow 68 contour provided (default 0). Allow 69 subdomain boundary refinement 70 70 though 71 - NoBoundaryRefinementAllBoundaries : do not refine boundary, only follow 71 - NoBoundaryRefinementAllBoundaries : do not refine boundary, only follow 72 72 contour provided (default 0) 73 - maxnbv : maximum number of vertices used to 73 - maxnbv : maximum number of vertices used to 74 74 allocate memory (default is 1.0e6) 75 - maxsubdiv : maximum subdivision of exisiting 75 - maxsubdiv : maximum subdivision of exisiting 76 76 elements (default is 10) 77 - metric : matrix (numberofnodes x 3) used as 77 - metric : matrix (numberofnodes x 3) used as 78 78 a metric 79 - Metrictype : 1 -> absolute error 79 - Metrictype : 1 -> absolute error 80 80 c/(err coeff^2) * Abs(H) (default) 81 2 -> relative error 82 c / (err coeff^2) * Abs(H) / 81 2 -> relative error 82 c / (err coeff^2) * Abs(H) / 83 83 max(s, cutoff * max(s)) 84 3 -> rescaled absolute error 85 c / (err coeff^2) * Abs(H) / 84 3 -> rescaled absolute error 85 c / (err coeff^2) * Abs(H) / 86 86 (smax - smin) 87 - nbjacoby : correction used by Hessiantype = 1 87 - nbjacoby : correction used by Hessiantype = 1 88 88 (default is 1) 89 - nbsmooth : number of metric smoothing 89 - nbsmooth : number of metric smoothing 90 90 procedure (default is 3) 91 - omega : relaxation parameter of the 92 smoothing procedure (default is 91 - omega : relaxation parameter of the 92 smoothing procedure (default is 93 93 1.8) 94 - power : power applied to the metric 94 - power : power applied to the metric 95 95 (default is 1) 96 - splitcorners : split triangles which have 3 97 vertices on the outline (default is 96 - splitcorners : split triangles which have 3 97 vertices on the outline (default is 98 98 1) 99 99 - verbose : level of verbosity (default is 1) 100 100 101 - rifts : followed by an ARGUS file that 101 - rifts : followed by an ARGUS file that 102 102 prescribes the rifts 103 - toltip : tolerance to move tip on an 104 existing point of the domain 103 - toltip : tolerance to move tip on an 104 existing point of the domain 105 105 outline 106 - tracks : followed by an ARGUS file that 107 prescribes the tracks that the mesh 106 - tracks : followed by an ARGUS file that 107 prescribes the tracks that the mesh 108 108 will stick to 109 - RequiredVertices : mesh vertices that are required. 109 - RequiredVertices : mesh vertices that are required. 110 110 [x, y, ref]; ref is optional 111 - tol : if the distance between 2 points of 112 the domain outline is less than 111 - tol : if the distance between 2 points of 112 the domain outline is less than 113 113 tol, they will be merged 114 114 … … 250 250 251 251 # Flag how many edges we have now, that way we know which edges 252 # belong to the subdomain. Will be used later fo required edges 252 # belong to the subdomain. Will be used later fo required edges 253 253 # if NoBoundaryRefinement equals 1. 254 254 new_edge_length = len(bamg_geometry.Edges) 255 edges_required = range((edge_length + 1), (new_edge_length + 1)) # NOTE: Upper bound of range is non-inclusive (compare to src/m/mesh/bamg.m) 256 255 edges_required = np.asarray(range((edge_length + 1), (new_edge_length + 1))) # NOTE: Upper bound of range is non-inclusive (compare to src/m/mesh/bamg.m) 257 256 if i: # NOTE: same as `if i > 0` (MATLAB is `if i > 1`) 258 257 bamg_geometry.SubDomains = np.vstack((bamg_geometry.SubDomains, [2, count + 1, 1, -subdomain_ref])) … … 445 444 np.array([[count, count + 1, 2]]), 446 445 np.hstack(( 447 np.arange(count + 1, count + nods).reshape(-1, ), 446 np.arange(count + 1, count + nods).reshape(-1, ), 448 447 np.arange(count + 2, count + nods + 1).reshape(-1, ), 449 448 (1 + i) * np.ones((nods - 1, 1)) … … 455 454 nods = rift[i]['nods'] - 1 456 455 bamg_geometry.Vertices = np.vstack(( 457 bamg_geometry.Vertices, 456 bamg_geometry.Vertices, 458 457 np.hstack(( 459 458 rift[i]['x'][:],
Note:
See TracChangeset
for help on using the changeset viewer.
