could you give me a hint how to dump out information about the Jacobian? I would like to get a better feeling for one of the problems that I try to solve. For that I would like to get a feeling how the singular values look like. The problem is a steady state problem that only depends on the boundary values. This means I do not have to give initial conditions or so. My idea was roughly the following.
- Set up the simulation to have two “time steps” in preCICE, i.e. setting
<timestep-length value="1" />in the preCICE config file. The first preCICE “time step” is used to solve the problem. The second “time step” is purely used to mess around with preCICE.
- Use the IMVJ postprocessing.
- Solve the problem until reaching the steady state in the first “time step”.
- Set the options of the IMVJ postprocessing such that drops all singular values at the beginning of the “second” time step. Ideally, it would tell me what the values where.
- Abort simulation after dumping the information such that “second time step” is never carried out.
Does this approach make sense/work? Can I find a bit more information about the settings for the IMVJ postprocessing than in the XML reference and the IMVJ wiki page?
Thank you in advance!