lanl / LaGriT

Los Alamos Grid Toolbox (LaGriT) is a library of user callable tools that provide mesh generation, mesh optimization and dynamic mesh maintenance in two and three dimensions.
https://lanl.github.io/LaGriT/
Other
116 stars 48 forks source link

hextotet / 14 and grid2grid/prismtotet14 seg fault #213

Closed millerta closed 3 years ago

millerta commented 3 years ago

hextotet prism to 3 and 18 work, but prism to 14 tet crashes. Most likely something in the early setup portion of the code. The following hextotet is an expanded version of the hextotet test that should be added to the test suite.

From Carl:

*****

Test hextotet 5, 6, and 24

*****

cmo / create / mohex /// hex createpts/brick/xyz/ 3 3 3 /0. 0. 0. / 1. 1. 1. / 1 1 1 hextotet / 24 / motet / mohex cmo/set_id/ motet/ element / itetclr cmo/select/motet; quality; cmo printatt motet -all- minmax cmo status motet cmo status mohex cmo delete mohex cmo delete motet cmo / create / mohex / / / hex createpts/brick/xyz/ 3 3 3 /0. 0. 0. / 1. 1. 1. / 1 1 1 hextotet / 6 / motet / mohex cmo/set_id/ motet/ element / itetclr cmo/select/motet; quality; cmo printatt motet -all- minmax cmo status motet cmo status mohex cmo delete mohex cmo delete motet cmo / create / mohex / / / hex createpts/brick/xyz/ 3 3 3 /0. 0. 0. / 1. 1. 1. / 1 1 1 hextotet / 5 / motet / mohex cmo/set_id/ motet/ element / itetclr cmo/select/motet; quality; cmo printatt motet -all- minmax cmo status motet cmo status mohex cmo delete mohex cmo delete motet #

Test prism to tet

#

Build a triangle mesh

Use two version of the triangle mesh to stack surfaces and make a prism mesh

Convert prism mesh to tet mesh

# cmo / create / motri /// triplane createpts/xyz/ 3 3 3 / 0. 0. 0. / 1. 1. 0. / 1 1 1 connect cmo / copy / motri2 / motri cmo / select / motri2 trans / 1 0 0 / 0. 0. 0. / 0. 0. 1. dump / avs / tmp_s1_tri.inp / motri dump / avs / tmp_s2_tri.inp / motri2 cmo / create / mostack stack / layers / avs / tmp_s1_tri.inp 1 / tmp_s2_tri.inp stack / fill / moprism / mostack dump / avs / prism.inp / moprism

cmo / status / moprism

hextotet / 3 / mopri3 / moprism hextotet / 18 / mopri18 / moprism

hextotet / 14 / mopri14 / moprism

cmo / status / brief

dump / avs / pri3.inp / mopri3 dump / avs / pri18.inp / mopri18

finish

millerta commented 3 years ago

duplicate issue