You may be interested in the "Checkpoint" mesh format:
https://github.com/libMesh/libmesh/blob/master/include/mesh/checkpoint_io.h

It is still XDR based... but basically just writes out the raw mesh as it
sits on each processor and then reads it back in during a restart.  Because
of this you have to restart using the same number of processors as the
initial run.  It's the format we use for restarting in MOOSE.

Derek

On Thu, Jan 14, 2016 at 6:10 PM Salazar De Troya, Miguel <
salazardet...@llnl.gov> wrote:

> I didn’t know the numbering could change. Now it makes sense. Thanks.
>
> Sincerely,
>
> Miguel A. Salazar de Troya
> Graduate Scholar, Lawrence Livermore National Laboratory
> Ph.D. Candidate, University of Illinois at Urbana-Champaign
> T1878
> Rm: 112
> Ext: 2-6411
> L-code: 288
>
>
>
>
> On 1/14/16, 11:13 AM, "Roy Stogner" <royst...@ices.utexas.edu> wrote:
>
> >
> >On Thu, 14 Jan 2016, Salazar De Troya, Miguel wrote:
> >
> >> Sorry the last statement is wrong. Leaving partition_agnostic=true with
> >>a
> >> SerialMesh doesn’t return the same assertion error, but the vector
> >>written
> >> and the vector read are not the same when I print them. The values seem
> >>to
> >> be in different order.
> >
> >Yes, that's as expected.  They're written out in a Hilbert
> >space-filling curve order, to make initial repartitioning easier when
> >they're read in again.
> >
> >>>> In that case, what would be the compatible way to restart the mesh if
> >>>>I
> >>>> still leave partition_agnostic=false?
> >
> >I'm not positive, since I don't use this myself.  In general it's
> >*not* safe to have your code depend on indexing.  We keep dof indices
> >contiguous on processors, which requires shifting every index with
> >every refinement/coarsening step, for example.
> >
> >Probably set skip_renumbering on before reading the Mesh?
> >
> >>>> Leaving partition_agnostic=true doesn¹t avoid the assertion either
> >>>> for me.
> >
> >You know what?  It doesn't for me either, my previous reports to the
> >contrary.  But that use case *should* work, and does work on lots of
> >other meshes I've tried.  I'll keep looking into it.
> >
> >Thanks,
> >---
> >Roy
>
>
> ------------------------------------------------------------------------------
> Site24x7 APM Insight: Get Deep Visibility into Application Performance
> APM + Mobile APM + RUM: Monitor 3 App instances at just $35/Month
> Monitor end-to-end web transactions and take corrective actions now
> Troubleshoot faster and improve end-user experience. Signup Now!
> http://pubads.g.doubleclick.net/gampad/clk?id=267308311&iu=/4140
> _______________________________________________
> Libmesh-users mailing list
> Libmesh-users@lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/libmesh-users
>
------------------------------------------------------------------------------
Site24x7 APM Insight: Get Deep Visibility into Application Performance
APM + Mobile APM + RUM: Monitor 3 App instances at just $35/Month
Monitor end-to-end web transactions and take corrective actions now
Troubleshoot faster and improve end-user experience. Signup Now!
http://pubads.g.doubleclick.net/gampad/clk?id=267308311&iu=/4140
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to