It is not clear to me how to access an object in the structure 
“Postprocessor”. The idea is that I compute the solution, but I am also 
computing secondary variables using the postprocessor structure in this way:

dealii::DataOut<dim> data_out;

data_out.attach_dof_handler(dof_handler);

Postprocessor postprocessor;

data_out.add_data_vector(solution, postprocessor);

data_out.build_patches();

When declaring the postprocessor class, I declare a function 
evaluate_vector_field where I compute all the secondary variables that are 
written on a vtk file later:

void evaluate_vector_field(const 
dealii::DataPostprocessorInputs::Vector<dim>& input_data, 
std::vector<dealii::Vector<double>>& computed_quantities) const override;

                …..

    for (int q = 0; q < n_q_points; ++q) {

        // Displacements

        computed_quantities[q](0) = input_data.solution_values[q](0);

        computed_quantities[q](1) = input_data.solution_values[q](1);

        // Stress

        computed_quantities[q](2)  = sigma_x

        computed_quantities[q](3)  = sigma_y

    }

These variables are saved on the data_out object, so I use 
“data_out.write_vtk(output);” and I get them on the vtk file. However, I do 
not know how to access to these new variables (computed_quantities) from 
another file different from the one where I declare the functions. 

The function evaluate_vector_field is a virtual function, and it is 
declared as const override. I do not know how to export the variables 
computed_quantities or modify them, because it is only allowed to read them 
since the function is declared as constant. Is there any way to access the 
structure postprocessor? Something like:

postprocessor.getVariables(xxxxxxxx)

Or an example where I can see how to implement it. Thank you very much in 
advance

Raul

-- 
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.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/d9a60280-0181-427b-a372-2eee29ab3257n%40googlegroups.com.

Reply via email to