<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);">
Hello,</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 am trying to use the HDF5 viewer to store the solution of my system at every iteration.
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
I am storing the solutions in order to compare two different results, one from an old sequential program and one from a new parallel program that is parallelized using the DMDA structure.</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
I am using matlab to compare the results.<br>
</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);">
In the sequential program, the vector for the solution is just a simple one dimensional vector made from VecCreate. In the parallel program, the vector is made from DMCreateGlobalVector.</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);">
In both versions of the program, I create the viewer as:</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);">
<div style="color: #000000; background-color: #ffffff; font-family: Consolas, 'Courier New', monospace; font-weight: normal; font-size: 14px; line-height: 19px">
<div><span style="color: #af00db">call</span><span> </span><span style="color: #795e26">PetscViewerHDF5Open</span><span>(PETSC_COMM_WORLD,</span><span style="color: #a31515">"filename.hdf5"</span><span>,&</span></div>
<div><span> FILE_MODE_WRITE,hdf5viewer,ierr)</span></div>
<div><span><br>
</span></div>
<div><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; line-height: normal; color: rgb(0, 0, 0); background-color: rgba(0, 0, 0, 0);">When I load the files into MATLAB after words, the sequential program creates a new dataset
each iteration, while the parallel version using DMCreateGlobalVector only stores a single dataset.</span><span><br>
</span></div>
<div><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; line-height: normal; color: rgb(0, 0, 0); background-color: rgba(0, 0, 0, 0);">Is there a way I can force both viewers to create a new dataset in the same way each time vecview
is called? Or is there someway I can manually create a new dataset?</span><span><br>
</span></div>
<div><span><br>
</span></div>
<div><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; line-height: normal; color: rgb(0, 0, 0); background-color: rgba(0, 0, 0, 0);">I am using Open MPI 3.1.4, Petsc 3.12.4, Intel Fortran 2019.5, and hdf5/1.10.5.</span></div>
<div><span><br>
</span></div>
<div><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; line-height: normal; color: rgb(0, 0, 0); background-color: rgba(0, 0, 0, 0);">Thanks,</span></div>
<div><span><br>
</span></div>
<div><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; line-height: normal; color: rgb(0, 0, 0); background-color: rgba(0, 0, 0, 0);">Lucas Banting</span><span><br>
</span></div>
</div>
<br>
</div>
</body>
</html>