On 6/23/21 1:57 AM, Chen R wrote:
>
> I am mainly trying to integrate the dealii system with our own geometric
> modeling program to do FE calculations. Right now created a basic elasticity
> solver based on the tutorials.
> The thing is I want to create a triangulation from our geometry and calculate
> on it then present the results back to in our program, but not entirely sure
> in what order the results are outputted using the SolverCG.
>
> My idea was to use some object like DataOut which can ouput the solution in a
> predetermined way but I see that this is not the way to do it.
Right. DataOut produces data in some internal format that is useful for later
outputting it in graphical form. What that internal format is is unimportant,
but it is not what you want anyway.
If you want to know what the solution is on an individual cell, ask the cell via
cell->get_dof_indices(...)
what its degrees of freedom are, and then query the solution vector for these
components. Or just call
cell->get_dof_values(...)
right away for the values of degrees of freedom located on a cell.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email:
bang...@colostate.edu
www:
http://www.math.colostate.edu/~bangerth/