On Fri, Aug 24, 2012 at 3:48 PM, Blaise Bourdin <span dir="ltr"><<a href="mailto:bourdin@lsu.edu" target="_blank">bourdin@lsu.edu</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word">I am attaching a debugged version.</div></blockquote><div><br></div><div>Checked in as DMComplex ex1f90.F</div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word"><span class="HOEnZb"><font color="#888888"><div>Blaise</div><div><br></div><div></div></font></span></div><br><div style="word-wrap:break-word"><div><br><div><div>On Aug 24, 2012, at 3:30 PM, Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>> wrote:</div>
<br><blockquote type="cite">On Fri, Aug 24, 2012 at 2:16 PM, Chris Eldred <span dir="ltr"><<a href="mailto:chris.eldred@gmail.com" target="_blank">chris.eldred@gmail.com</a>></span> wrote:<br><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
The following code snippet does not appear to be working:<br></blockquote><div><br></div><div>1) I hate Fortran. A lot. It is not debuggable.</div><div><br></div><div>2) Here is my code that should work, but the do loop somehow does not work:</div>
<div><br></div><div> program main</div><div> implicit none</div><div>!</div><div>#include <finclude/petsc.h90></div><div>#include <finclude/petscdmcomplex.h90></div><div>!</div><div> DM dm</div>
<div> PetscInt, target, dimension(4) :: EC</div><div> PetscInt, pointer :: pEC(:)</div><div> PetscInt c, firstCell, numCells</div><div> PetscErrorCode ierr</div><div><br></div><div> call PetscInitialize(PETSC_NULL_CHARACTER,ierr)</div>
<div> call DMComplexCreate(PETSC_COMM_WORLD, dm, ierr)</div><div> firstCell = 0</div><div> numCells = 5</div><div> call DMComplexSetChart(dm, 0, numCells, ierr)</div><div> do c=firstCell,numCells</div>
<div> write(*,*) 'c',c</div><div> call DMComplexSetConeSize(dm, c, 4)</div><div> end do</div><div> call DMSetUp(dm, ierr)</div><div><br></div><div> EC(1) = 1</div><div> EC(2) = 2</div>
<div> EC(3) = 3</div><div> EC(4) = 4</div><div> pEC => EC</div><div> write(*,*) 'cell',c,pEC</div><div> c = 0</div><div> call DMComplexSetCone(dm, c , pEC, ierr)</div><div> CHKERRQ(ierr)</div>
<div> call DMComplexGetCone(dm, c , pEC, ierr)</div><div> CHKERRQ(ierr)</div><div> write(*,*) 'cell',c,pEC</div><div><br></div><div> call DMDestroy(dm,ierr)</div><div> call PetscFinalize(ierr)</div>
<div> end</div><div><br></div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
PetscInt, dimension(4) :: EC<br>
<br>
EC(1) = 1<br>
EC(2) = 2<br>
EC(3) = 3<br>
EC(4) = 4<br>
write(*,*) 'cell',EC<br>
call DMComplexSetCone(model_mesh, lcell , EC, ierr)<br>
CHKERRQ(ierr)<br>
call DMComplexGetCone(model_mesh, lcell , EC, ierr)<br>
CHKERRQ(ierr)<br>
write(*,*) 'cell', EC<br>
<br>
I get the following as output:<br>
<br>
cell 1 2 3 4<br>
cell 18673792 0 -1 -1<br>
<br>
DMComplexGetConeSize returns 4 as expected. I am using the latest<br>
version of petsc-dev.<br>
<br>
Later on there is a segmentation fault violation, probably related to<br>
this issue.<br>
<br>
Any ideas?<br>
<span><font color="#888888"><br>
--<br>
Chris Eldred<br>
DOE Computational Science Graduate Fellow<br>
Graduate Student, Atmospheric Science, Colorado State University<br>
B.S. Applied Computational Physics, Carnegie Mellon University, 2009<br>
<a href="mailto:chris.eldred@gmail.com" target="_blank">chris.eldred@gmail.com</a><br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<br>
</blockquote></div><br><div>
<span style="border-collapse:separate;font-size:12px;border-spacing:0px"><span style="text-indent:0px;letter-spacing:normal;font-variant:normal;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:12px;white-space:normal;font-family:Helvetica;word-spacing:0px"><div style="word-wrap:break-word">
<span style="text-indent:0px;letter-spacing:normal;font-variant:normal;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:12px;white-space:normal;font-family:Helvetica;word-spacing:0px"><div style="word-wrap:break-word">
<span style="text-indent:0px;letter-spacing:normal;font-variant:normal;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:12px;white-space:normal;font-family:Helvetica;word-spacing:0px"><div style="word-wrap:break-word">
<div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px">-- </div><div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px">Department of Mathematics and Center for Computation & Technology</div>
<div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px">Louisiana State University, Baton Rouge, LA 70803, USA</div><div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px">Tel. <a href="tel:%2B1%20%28225%29%20578%201612" value="+12255781612" target="_blank">+1 (225) 578 1612</a>, Fax <a href="tel:%2B1%20%28225%29%20578%204276" value="+12255784276" target="_blank">+1 (225) 578 4276</a> <a href="http://www.math.lsu.edu/~bourdin" target="_blank">http://www.math.lsu.edu/~bourdin</a></div>
<div><br></div><div><br></div><br></div></span></div></span><br></div></span><br></span><br>
</div>
<br></div></div><br></blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<br>