I am using a procedure similar to step-40.
DataOut<dim> data_out;
data_out.attach_dof_handler(dof_handler);
// ========================[ CONFIGURATIONAL FORCES
]==========================
std::vector<std::string> configurational_forces_magnitude(dim,
"config_forces");
std::vector<DataComponentInterpretation::DataComponentInterpretation>
configurational_forces_interpretation(dim,
DataComponentInterpretation::component_is_part_of_vector);
data_out.add_data_vector(configurational_forces,
configurational_forces_magnitude, DataOut<dim>::type_dof_data,
configurational_forces_interpretation);
// ================================[ WRITE VTU
]===============================
// Construct data structures
data_out.build_patches();
// Write output
const std::string filename = ("solution-" +
Utilities::int_to_string(timestep_no, 4) + "." +
Utilities::int_to_string(triangulation.locally_owned_subdomain(), 4));
std::ofstream output(("output/" + filename + ".vtu").c_str()); // Store
file name string as output stream
data_out.write_vtu(output); // Write .vtu output file
// ================================[ WRITE PVTU
]==============================
// Write "master record" (names of the various files that combined
represents the graphical data for the entire domain)
if ( Utilities::MPI::this_mpi_process(mpi_com) == 0 ) // only 1st processor
runs this task!
{
std::vector<std::string> filenames;
for (unsigned int i = 0; i < Utilities::MPI::n_mpi_processes(mpi_com); ++i)
filenames.push_back("solution-" + Utilities::int_to_string(timestep_no, 4)
+ "." + Utilities::int_to_string(i, 4) + ".vtu");
std::ofstream master_output(("output/solution-" +
Utilities::int_to_string(timestep_no, 4) + ".pvtu").c_str());
data_out.write_pvtu_record(master_output, filenames);
}
The above is my complete procedure for output of the material forces only.
Kind regards,
S. A. Mohseni
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
