Commit 654b0cb0 authored by GILLES Sebastien's avatar GILLES Sebastien

#1422 Add comments about the fact ascii files can't be used to reload a PETSc vector.

parent 7bb02445
......@@ -507,6 +507,10 @@ namespace MoReFEM
* \copydoc doxygen_hide_invoking_file_and_line
* \copydetails doxygen_hide_mpi_param
* \param[in] output_file File into which the vector content will be written.
*
* \attention To my knowledge there are no way to reload a PETSc vector from the file; if you need
* to do so rather use binary format.
*
*/
void View(const Mpi& mpi, const std::string& output_file, const char* invoking_file, int invoking_line,
PetscViewerFormat format = PETSC_VIEWER_DEFAULT) const;
......@@ -541,6 +545,14 @@ namespace MoReFEM
* \param[in] output_file File into which the vector content will be written.
* \param[in] binary_or_ascii_choice Whether the vector should be printed as binary or ascii. Default
* value takes its cue from the choice written in the input data file.
*
* \attention If the purpose is to be able to reload the file from disk later on, format MUST be
* binary (PETSc [VecLoad](https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Vec/VecLoad.html)
* works only with binary or HDF5). The file is not enough to rebuild the vector: you will also need
* to know:
* - The number of elements processor-wise.
* - The number of elements program-wise.
* - The program-wise indexes on the ghosted elements.
*/
template<MpiScale MpiScaleT>
void Print(const Mpi& mpi,
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment