The next release, expected in
early 2011, will have 64 bit memory pointers. This will allow access to
much more real and virtual memory so that the present limitations on
memory and mesh size that can be built and manipulated will be greatly
expanded.
Modifications and New Capabilities Index:
dump / zone_outside / ...
Changed FEHM outside area calculation to default to Voronoi area
associated with nodes of a 3D tetrahedral mesh instead of computing
Median area. FEHM file file_root_outside.area changed to
file_root_outside_vor.area For dump/zone, added keywords keepatt_area
or keepatt_voronoi which will compute and keep voronoi vector areas
xn_varea, yn_varea, zn_varea and keepatt_median will compute area/num
nodes on face and keep attributes xn_area, yn_area, zn_area The written
file file_root_outside_vor.area or file_name_outside_med.area is a list
of 2D area vectors (Ax_i,Ay_i,Az_i) associated with each node. https://lagrit.lanl.gov/docs/DUMP3.html
dump/zone_outside
cmo / addatt /
voronoi_varea / ...
This module will do the same voronoi calculation on triangles as is
done with the outside area for a 3D tetrahderal mesh. The call will
create vector components for each node and fill the node attributes
xn_varea, yn_varea, zn_varea. https://lagrit.lanl.gov/docs/cmo_addatt.html
extract / surfmesh / ...
Now creates attributes to hold element local face numbers of 3D input
mesh that occur on either side of output mesh face, idface0 and
idface1. Now copies user-created node-based attributes from source mesh
into extracted sink mesh. https://lagrit.lanl.gov/docs/EXTRACT1.html
addmesh / excavate / ...
excavate - remove nodes and elements if they fall with the circumsphere
of triangles on the input mesh.
mesh1 must be a 3D mesh (of any geometry) and mesh2 must be a 2D
triangular mesh. This command then excavates a volume in mesh1 around
mesh2, such that the surface could then be inserted into the 3D mesh
(such as to insert a fault into a background terrain mesh). The
background mesh, minus the excavated/removed nodes, is put into mesh3.
If the optional [bfs] argument is given, the routine will use a
breadth-first search algorithm to find nodes to remove, as opposed to
the default KD-tree algorithm. If the optional [connect] argument is
given, the program will, after excavation, execute an addmesh/append,
and then a connect, to produce a fully connected mesh with the surface
(mesh2) inserted into the background (mesh1). https://lagrit.lanl.gov/docs/ADDMESH.html
interpolate / ...
Changed interpolate to "find" more points on edges this will permit
nodes to find a nearest edge or point and be "inside" the triangle for
extreme small or large numbers where epsilon values are difficult to
evaluate correctly. Note, this changed test results for interpolate,
test/level01 results were updated for these improvements. https://lagrit.lanl.gov/docs/main_interpolate.html
intrp_gtg.f
A bug was fixed in interpolation that would sometimes save
a node id in pt_gtg or el_gtg attributes that was not related
to the found candidate and value. This could occur where there are
multiple candidates for the source and if epsilon values are near
machine limits. The test in level01/intrp_2D_sizes
was changed to capture and evaluate these issues.
intrp_gtg.f, inside_lg.f
There are changes to interpolate using tests for finding points
that are inside or on edges or vertices of an element. The epsilon
tests have been relaxed to allow points that are "near" to be
found on edge - if within the chosen epsilon. The interpolation
has been changed to evaluate candidate points based on the
confidence of being inside the associated triangle. A result
indicating the point is inside will "win" over a candidate result
that is on edge or vertice. If idebug attribute is set to a
number of 5 or greater, there will be many more statements
written that are related to the inside triangle and epsilon
tests.
Substantial changes and additions to DUMP2.html which describes
all the dump file_types and include more descriptions of FEHM files.
DUMP3.html for dump/fehm now includes descriptions of FEHM files.
https://lagrit.lanl.gov/docs/DUMP2.html
https://lagrit.lanl.gov/docs/DUMP3.html
read /...
Three token read implemented. Files are recognized based on their
suffix (AVS, GMV, ...) rather than requiring that the second token
specify the file type.
https://lagrit.lanl.gov/docs/READ.html
dump /...
Two and three token dump implemented. Files types (AVS, GMV, ...) are
recognized based on their
suffix rather than requiring that the second token specify the file
type. Two token write does not require a MO name. The default MO is
used.
https://lagrit.lanl.gov/docs/DUMP2.html
cmo / attribute_union /
...
Change two meshes so they both share the same set of attributes (taking
the union of their sets of attributes)
https://lagrit.lanl.gov/docs/cmo_att_derive.html
compute /
linear_extrapolate / ...
linear_extrapolate - keyword for an extrapolation from an attribute
value in a surface onto every node of a 3D mesh. Given a 3D mesh and a
2D surface, this command will extrapolate a scalar value from that
surface onto every point of the mesh. This can be used to (for example):
* Propogate head values from a surface onto all
nodes of a mesh.
* Expand a mesh to fit a surface, by propogating the
appropriate spatial coordinate.
* Compute the depth relative to a topographic
surface to each node of a mesh. https://lagrit.lanl.gov/docs/COMPUTE.html
grid2grid /...
grid2grid wrapper for hextotet. Use to convert:
- quadtotri2 quad to 2 triangles, no new
points.
- prismtotet3 prism to 3 tets, no new points.
- quattotri4 quad to 4 triangles, with one
new point.
- pyrtotet4 pyramid to 4 tets, with
one new point.
- hextotet5 hex to 5 tets, no new
points.
- hextotet6 hex to 6 tets, no new
points.
- prismtotet14 prism to 14 tets, four new points (1 + 3
faces).
- prismtotet18 prism to 18 tets, six new points (1 + 5
faces).
- hextotet24 hex to 24 tets, seven new
points (1 + 6 faces).
- tree_to_fe quadtree or octree grid to
grid with no parent-type elements.
https://lagrit.lanl.gov/docs/GRID2GRID.html
dump / stor / ...
anothermatbld3d_wrapper.f
Create two new node vectors, ccoef, ij_ccoef
Put the negative ij coefficient value into the two nodes connected to
the ij edge.
The vector ij_coef will assign the j index value to node i so that one
can determine which edge is associated with the neative coefficient
that is assigned to nodes. https://lagrit.lanl.gov/docs/DUMP2.html
https://lagrit.lanl.gov/docs/DUMP3.html
Changes to include TranslateTetToZero for geometric calculations
(not sure it is really helping with some of our neg ccoef issues)
anothermatbld3d_wrapper.f
Extensive chages to error handling and messages, but not to the logic
of program
This code has same logic as matbld3d - but uses linked lists instead of
mmgetblk calls
Use io_type to toggle creation of attribute for voronoi volumes or to
write to stor file
added extensive error checking to eliminate segmentation faults
added error check and message for every mmgetblk and mmrelblk
added calls to mmprint when mm calls fail
cleaned up variable declarations and added comments
added istatus to check for errors and completion of matrix
changed all routine messages to start with AMatbld3d_stor to
distinguish from matrix built with Matbld3d_stor
added idebug options
added status report at end of routine
matbld3d_stor.f
Extensive chages to error handling and messages, but not to the logic
of program
This code uses many mmgetblk calls and about 40 percent more memory
than linked list version
added extensive error checking to eliminate segmentation faults
added error check and message for every mmgetblk and mmrelblk
added calls to mmprint when mm calls fail
cleaned up variable declarations and added comments
added istatus to check for errors and completion of matrix
added idebug options
added status report at end of routine
matbld3d_stor.f
add warning for newlen call with uncertain effect
sparseMatrix.c
initialize list pointers to null
assign null to pointers after free
add warning messages for failure to free
dumpfehm.f
Add compress_opt to dumpfehm arguments
add comments and error checking to clarify code logic
check options and set for 2D or 3D calls to matbld
use matbld3d_stor for compress options none and coefs
use anothermatbld3d_wrapper for compress options all and graph
Note anothermatbld3d_wrapper can write only scalar coef values
pset /...
Add option to pset/ / zone for user specified zone id number.
https://lagrit.lanl.gov/docs/PSET.html
dump
/
writedump.f
declare implicit none and initialize variables
add comments to clarify the case switches
add more error checking and messages
change syntax for dump/ fehm and dump/ stor
old keywords not needed include alternate_scalar, binaryc, asciic
compression keywords are now none, coefs, graph, or all old syntax
still works, but now code checks
for keywords after filename and cmo and sets
options for the fehm and stor routine calls The man pages are updated
and corrected. https://lagrit.lanl.gov/docs/DUMP2.html
memory
/
...
New options to print and check memory manager and report memory usage.
This superseeds old utilities mmprint, mmcheck, etc. https://lagrit.lanl.gov/docs/memory.html
Test Cases Modified:
test/level01
Change reference files to results from Linux
Updated zone_outside and zone_outside_minmax to current voronoi version
|