Changeset 13524

Timestamp:

10/03/12 15:05:41 (12 years ago)

Author:

jschierm

Message:

NEW: Preliminary python versions of some rift functions (not useful until TriMeshProcessRifts module is converted).

Location:

issm/trunk-jpl/src/m/mesh

Files:

• ElementsFromEdge.m (modified) (1 diff)
• ElementsFromEdge.py (added)
• rifts/meshprocessoutsiderifts.m (modified) (2 diffs)
• rifts/meshprocessoutsiderifts.py (added)
• rifts/meshprocessrifts.py (added)

issm/trunk-jpl/src/m/mesh/ElementsFromEdge.m

@@ -1 @@
-function edgeelements=ElementsFromEdge(elements,A,B) 
+function edgeelements=ElementsFromEdge(elements,A,B)
 %ELEMENTSFROMEDGE: find elements connected to one edge defined by nodes A and B
 %
-% Usage: edgeelements=ElementsFromEdge(elements,A,B) 
+%   Usage: edgeelements=ElementsFromEdge(elements,A,B) 
 %
-% Eg:    edgeelements=ElementsFromEdge(md.mesh.elements,tip1,tip2)
+%   Eg:    edgeelements=ElementsFromEdge(md.mesh.elements,tip1,tip2)
 %
 %
 edgeelements=find(...
-(elements(:,1)==A & elements(:,2)==B )| ...
-(elements(:,1)==A & elements(:,3)==B )| ...
-(elements(:,2)==A & elements(:,3)==B )| ...
-(elements(:,2)==A & elements(:,1)==B )| ...
-(elements(:,3)==A & elements(:,1)==B )| ...
-(elements(:,3)==A & elements(:,2)==B ));
+        (elements(:,1)==A & elements(:,2)==B )| ...
+        (elements(:,1)==A & elements(:,3)==B )| ...
+        (elements(:,2)==A & elements(:,3)==B )| ...
+        (elements(:,2)==A & elements(:,1)==B )| ...
+        (elements(:,3)==A & elements(:,1)==B )| ...
+        (elements(:,3)==A & elements(:,2)==B ));

issm/trunk-jpl/src/m/mesh/rifts/meshprocessoutsiderifts.m

@@ -41 +41 @@
                 elements=[];
 
-                while  flags(B), %as long as B does not belong to the domain outline, keep looking.
+                while flags(B), %as long as B does not belong to the domain outline, keep looking.
                         %detect elements on edge A,B:
                         edgeelements=ElementsFromEdge(md.mesh.elements,A,B);
@@ -69 +69 @@
                 %deal with segments
                 tipsegments=find((md.mesh.segments(:,1)==tip) | (md.mesh.segments(:,2)==tip));
-                for  k=1:length(tipsegments),
+                for k=1:length(tipsegments),
                         segment_index=tipsegments(k);
                         pos=find(md.mesh.segments(segment_index,1:2)~=tip);