<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Hi all,</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
I might be missing something obvious, but I can't tell from the documentation if 1D coordinate DMSWARMs (immersed in interval DMPlex meshes) are supported. Does anyone know?</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
I've successfully implemented 2D and 3D coordinate DMSwarms in 2D and 3D DMPlexes using DMSwarmSetPointCoordinates in Firedrake using petsc4py (via swarm.setPointCoordinates. The petsc4py function forces the input numpy array to have 2 dimensions with 1 column
per dimension. In the 1D case, where each row of the coordinates array therefore has one column, I get the following error:</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span>E petsc4py.PETSc.Error: error code 63<br>
</span>
<div>E [0] DMSwarmSetPointCoordinates() line 305 in /Users/rwh10/firedrake/src/petsc/src/dm/impls/swarm/swarmpic.c<br>
</div>
<div>E [0] DMLocatePoints() line 6499 in /Users/rwh10/firedrake/src/petsc/src/dm/interface/dm.c<br>
</div>
<div>E [0] DMLocatePoints_Plex() line 744 in /Users/rwh10/firedrake/src/petsc/src/dm/impls/plex/plexgeometry.c<br>
</div>
<div>E [0] DMPlexLocatePoint_Internal() line 462 in /Users/rwh10/firedrake/src/petsc/src/dm/impls/plex/plexgeometry.c<br>
</div>
<div>E [0] Argument out of range<br>
</div>
<span>E [0] No point location for cell 0 with type segment</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span><br>
</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Thanks</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Reuben Hill</div>
</body>
</html>