Changeset 22319

Timestamp:

01/03/18 10:07:00 (7 years ago)

Author:

Mathieu Morlighem

Message:

CHG: fixed roundmesh, need to use linspace rather than colon operator

Location:

issm/trunk-jpl

Files:

• src/m/mesh/roundmesh.m (modified) (2 diffs)
• test/Archives/Archive1205.arch (modified) ( previous)
• test/Archives/Archive423.arch (modified) ( previous)
• test/Archives/Archive433.arch (modified) ( previous)

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

@@ -10 +10 @@
 
 %First we have to create the domain outline 
+expname = [tempname() '.exp'];
 
 %Get number of points on the circle
-pointsonedge=floor((2.*pi*radius) / resolution);
+pointsonedge=floor((2.*pi*radius) / resolution)+1; %+1 to close the outline
 
 %Calculate the cartesians coordinates of the points
-x_list=ones(pointsonedge,1); y_list=ones(pointsonedge,1);
-theta=(0.:2.*pi/pointsonedge:2.*pi*(1.-1./pointsonedge))';
-x_list=roundsigfig(radius*x_list.*cos(theta),12);
-y_list=roundsigfig(radius*y_list.*sin(theta),12);
+theta=linspace(0,2*pi,pointsonedge)';
+x_list=roundsigfig(radius*cos(theta),12);
+y_list=roundsigfig(radius*sin(theta),12);
 A=struct('x',x_list,'y',y_list,'density',1.);
-expwrite(A,'RoundDomainOutline.exp');
+expwrite(A,expname);
 
 %Call Bamg
-md=triangle(md,'RoundDomainOutline.exp',resolution);
+md=triangle(md,expname,resolution);
 %md=bamg(md,'domain','RoundDomainOutline.exp','hmin',resolution);
 
@@ -32 +32 @@
 
 %delete domain
-delete('RoundDomainOutline.exp')
+delete(expname);
 end