Changeset 4495
- Timestamp:
- May 28, 2007, 10:07:56 AM (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/source/anuga/abstract_2d_finite_volumes/util.py
r4487 r4495 1582 1582 1583 1583 1584 def remove_lone_verts(verts, triangles): 1584 def remove_lone_verts(verts, triangles, number_of_full_nodes=None): 1585 """ 1586 Removes vertices that are not associated with any triangles. 1587 1588 verts is a list/array of points 1589 triangles is a list of 3 element tuples. Each tuple represents a triangle. 1590 1591 number_of_full_nodes relate to parallelism when a mesh has an 1592 extra layer of ghost points. 1593 1594 """ 1585 1595 verts = ensure_numeric(verts) 1586 1596 triangles = ensure_numeric(triangles) 1597 1587 1598 N = len(verts) 1599 1588 1600 # initialise the array to easily find the index of the first loner 1589 1601 loners=arange(2*N, N, -1) # if N=3 [6,5,4] 1590 1602 1591 for i,t in enumerate(triangles):1603 for t in triangles: 1592 1604 for vert in t: 1593 1605 loners[vert]= vert # all non-loners will have loners[i]=i … … 1603 1615 verts = verts[0:lone_start] 1604 1616 else: 1605 # change the loners list so it can be used to modify triangle 1617 # change the loners list so it can be used to modify triangles 1606 1618 # Remove the loners from verts 1607 1619 # Could've used X=compress(less(loners,N),loners)
Note: See TracChangeset
for help on using the changeset viewer.