LaGriT V1.100 Release Notes 1999

This text is converted from old pdf files and may have translation errors. See original pdf for clarification.

LaGriT V1.0 May 1999 PDF Version

A summary of the major Changes found in this release are listed below.

triangulate

triangulate a 2D mesh assuming the ordered nodes in the 2D mesh define the perimeter of a polygon

ung2avs

convert Arclnfo (GIS) Ungenerate files to AVS ung2avs/avs_file_out/ung_file_in/[z_va|ue]

define

allows a number to be associated with a character string, such that the character string can be used in input decks in place of the number.

define/nx/3
define/ny/4
define/nz/5
define/bottom/O.1/
define/top/4.6
define/left/ -4.7
define/right/9.8
surface/s1/reflect/box/0.0,left.bottom/1.0,right,top

colormap

This command builds the colormap. In reality it only builds the material adjacency graph, from which the colormap can be quickly generated when needed. Three actions are possible:

colormap/[add|create|delete]/[cmo_name]
  • add - The material adjacency characteristics of the specified mesh object is added to the existing material adjacency graph, which is created if it didn‘t exist. This is the default action.

  • create - The existing material adjacency graph is deleted and a new one created from the specified mesh object.

  • delete - The material adjacency graph is deleted if it exists. Any specified mesh object is ignored.

Examples:

colormap/create/mesh1
colormapl/mesh2
colormap/delete

massage

added a smoothing operation to the optimization which can be turned off with the nosmooth option

massage/creation/annihilation/toldamage//[ifirst,ilast,istride]/[nosmooth]

smooth

new option aspect will smooth to improve aspect ratio by moving a node toward the neighbor that provides the greatest improvement. New option lpfilter will smooth surface networks (i.e 2D mesh objects or the interface network of a 3D mesh) using a polynomial filter. (filtdeg default 30; k_pb default 0.1)

smooth/position/aspectl[ifirst,ilast,istride/toldamage]
smooth/position//pfilter/[ifirst,ilast,istride/filtdeg/k_pb]

pset

new option surface will identify nodes on the specified surface.

Keyword surface names have the following meaning:

-all- will identify nodes on any surface.
-interface- will identify nodes on any interface surface.
-boundary- will identify nodes on exterior surfaces.

pset/psetname/surface/surface_name/[ifirst,ilast,istride]

refine

new option roughness will refine based on the distance of the endpoint of an edge to the plane determined by the synthetic normal with respect to a specified surface at the other endpoint of the edge.

refine/roughness///edge/ifirst,ilast,istride/
distancelsurface_namelexclusivelinclusive
refinelroughnessllledge/1,0,0l.28lptoplinclusive

new option edge_list will bisect a set of edges specified by the node numbers of the endpoints of the edges.

refine/edge_list///edge/edge_listl

refine/edgeilist///edge/1 2 23 47/ will refine the edge with
endpoints 1 and 2 also the edge with endpoints 23 and 47.

new option interface will bisect a set of non-interface edges of tets all of whose vertices are interface nodes.

refine/interface///edge/pset,get,psetname//// [inclusivelexclusive]l

extract

new option network will extract the network of interfaces (consisting of parent nodes) from a mesh.

extractlnetworklifirst,ilast,istride/cmoout/cmoin

dump

dump/recolor/file_name This command writes the existing colormap to the specified file. (See colormap command.)

dump/fehm/file_name / [cmo_name] / [binary/ ascii | asciic | binaryc] /[scalar, vector, both] / [delatt, keepatt]

The [delatt, keepatt] option gives the user the ability to delete or keep the boundary attributes, top, bottom, left_w, right_e, back_n, front_s, which are created by dump/fehm. The default is delatt. dump/fehm/file_name / [cmo_name] / [binaryc | asciic] produces compressed matrices

dump/gmv/file_name/[cmo_name]/[binary, ascii] 

specify binary or ascii format of GMV file on command line

dump/lagrit/file,name/[cmo_name]/ 

will write an ascii restart file that contains geometry and mesh object information. cmo_name can be ‘-all-‘ in which case all mesh objects are written to the file or it can specify a list of mesh objects to be written.

read

read/lagrit/file_name/[cmo_name]/ 

will read an ascii restart file written by dump/lagrit. All mesh object data is preserved in the file including the cmo_name.

connect

connect will triangulate a 2d planar set of nodes generating a triangular Delaunay grid.

Bug fixes Jan 98 to May 99

multi_material - fixed error for node added that was on both an interface and an exterior boundary might get the wrong itp1 value.

connect - refresh pointers alter call to remove bigtet

ceL_chain - fix bug with memory allocation for mpary array.

massage,getmpary - correctly access pset for massage

try2tob - get pointer to icontab correctly

cel_chain - Check for psetnames = blank

gctbit,sctbit Change declaration of ISHFT to intrinsic

flip2t03,flp2t03b, flp2to3i - update itettyp for new element

recon2d - use cmo.h (icmoget) to pass to testdamage so it knows If it must refresh pointers

dumpavs - close file always before leaving subroutine

refine_edge_list_lg - correct pointer statement

tangent_plane, cer_chain - fix refine on roughness

refine_fix_add - correctly set ier values for added nodes on constrained interfaces

sheet - explicitly specify -def— for mesh object name

rzbrick - fix ratio fiag

control_command_lg - correctly remove unnecessary blanks from command lines

cmo_create - make interpolation type be ‘and’ for isetwd and xtetwd

cmo_interpolate - fix interpolation for isetwd and xtetwd

pset - idebug delared as integer

rmmat - fix error return flag

resetpts - fix error return flag

surfset - fix memory management error

gctmpary - sct defaults correctly by testing nwds

closed_surfaces - fix arguments to getregv2 call

refine_edge_add - modify pset membership for new nodes.

cmo_select, cmo_get_name - remove null character from end of name

recon2d - set itetoff

lpfilter,LowPassFilterModule - avoid overwiting data

Code Improvements Jan 98 to Mar 99

smooth - new option smooth/position/aspect will smooth to improve aspect ratios.

smooth, extract - new option smooth/position/lpfilter will smooth surface networks. New extract option extract/network will an interface network from a 3D mesh.

pset - New options for surface (surface names: -all-, -interface-, -boundary- have the obvious special meanings)

delaunay - Insert nodes in mesh in random order. Replace n**2 a1gorithm to find matching faces with a linked list approach

reeon2d - changed test to use consistent volume calculation.

refine, tangent_plane, cer_chain, refine_edge_list_lg, lpfilter, LowPassFilterModule, GmphModule - new command options

triangulate_lg, msgtty - add triangulate command

pntlimc - check for pset named ‘-def—‘ or empty string

corrected warning that showed up on the DEC compile in the following routines:

addmesh, addmesh delete, addlnesh pyramid, boundary components, chkreg, chkregv, closed surfaces, cmo_delatt_def, cmo_interpolate, cmo_setatt, cmo_release, connect, correctpc, derefine, dopmat, dumpchad, filholes, geniee, get_mregions, get_regions, get_surfaces, getreg, getregv, grid_to_grid, hextotet_att, l1n1en1adjb, hsb2seta, ifacept, initx3d, math, occonv, pstatus, readgmv_binary, refine_coupling_coef, refine_edge_add, refine_face, refine_face_add, rmregmn, rmsurf, rwdpmw, r2, search2d, sortbins, taylor_error, translate, volume_tet, voron2d, writedump, refine_edge_list_lg

recon2, mega_error - restrict existence of ‘mega’ related attributes to recon loop. change IO disposition to not write to GMV files

refine, refine_interface_elements_lg - new refine option to refine non—intcrface edges of tets, all of whose vertices are interface nodes.

dump/fehm, writedump, matbld3d_stor - generate compressed matrix for geometric coefficients .stor file

cel_chain, cer_chain, refine_edge_add - set pset membership of child nodes in refine_edge_add_tet

refine_edge_add - pset is inherited from ‘anding’ the pset of the endpoints of refined edgeh

Code Changes Nov 98 to Apr 99

agd3d massage - add smoothing operation to optimization loop in massage.

sgd, primestep - smoothing now automatic in massage, turn it off with ‘nosmooth’

cel_chain - remove call to recon from inside refine/rivara loop.

agd3d - allow more merges of nodes that do not have unique successors and predecessors

dumpavs - allow for ranksfl and limit coordinate range to (—1 16—30, 11e+30)

aratio_tet - handle extreme aspect ratio tets correctly

agd3d, aratio_tet, aratio_tri - remove assumption that fp errors would not be trapped

massage - set ipointi to 1 and ipointj to nnodes

intradd - use a more memory efficient al gorithm to create child nodes

agd3d - change error to warning when material match in question (skip merge)

all common blocks - moved common statements after declarations added ‘save’ statement

dump_recolor_lg, neighbor_recolor_lg, writedump - add dump/recolor command (see above)

dump - fehm option to keep/delete boundary attributes on fehm files dump_outside_list

ung2avs - option to convert Ungenerate files to AVS files

dumpgmv_hybrid - read binary/ascii from command line

llip3t02, llip4to4, llip2t00, llip3t021 [lip4to4i llip2t00b In[lip recon IceonZ fiiplt00, flip2t02, - remove calls to fluxing routines and clean up associated memory usage

control_eommand_lg - new method of command processing

writloga, writinit, dotask, dotaskx3d, initlagrit, msgtty, control lg.h, lagrith dumpgmvihybrid cmo attribute 7def- is modified so that it Will not be written to gmv files. writcdump,rcaddump dump/lagrit and rcad/lagrit - now write and read ascii geometry files dumpilagrit, eventually this command will also dump the mesh objects dump_lagrit_geom,

read_lagrit - read_lagrit_geom

cmo_dump_cmo dump/lagrit and read/lagrit - now write and read ascii re start files

cmo_read_dump_cmo - that contain geometry and mesh object intbnnation

matbld2dstor - add max connections to output, make consistent with matbld3d_stor

eset - don’t print element number of member of set

quality - print if idebug set to 1

connect2d_lg - new code to connect 2d planar node distributions into 2d grids

delaunay2d_lg, delaunay2d_connect_lg, multi_material2d_lg, fix_small_triangles_lg, make_big_triangle_lg scale_lg, msgtty change subroutine name ‘scale’ to ‘scale_lg’ to avoid conflicts with other libraries