Re: [deal.II] inverse operation to AffineConstraints.condense() (recovering full system matrix)

2023-10-28 Thread Wolfgang Bangerth
On 10/28/23 13:41, Mathieu wrote: What do you think is an appropriate data structure for the storage of the entries associated with the constrained dofs? A FullMatrix or a separate AffineConstraints? Later, I have to somehow merge this data structure with the condensed system matrix since I

Re: [deal.II] Mixed codimensionality

2023-10-28 Thread Wolfgang Bangerth
On 10/27/23 12:21, Daniel Arndt wrote: If you can formulate your problem/approach in a way (iteratively?) that you could solve separately on the three triangulations (using input from the solution on the other triangulations), that might be easier. That corresponds to an operator splitting s

Re: [deal.II] inverse operation to AffineConstraints.condense() (recovering full system matrix)

2023-10-28 Thread Mathieu
This makes sense. What do you think is an appropriate data structure for the storage of the entries associated with the constrained dofs? A FullMatrix or a separate AffineConstraints? Later, I have to somehow merge this data structure with the condensed system matrix since I have to solve a linea

Re: [deal.II] inverse operation to AffineConstraints.condense() (recovering full system matrix)

2023-10-28 Thread Wolfgang Bangerth
On 10/28/23 12:46, Simon wrote: After assembling and solving the linear system, I need the uncondensed system matrix and right hand side to perform a sensitivity analysis. In particular, I need the rows and entries of the system matrix and RHS corresponding to constrained dofs. However, the co

[deal.II] inverse operation to AffineConstraints.condense() (recovering full system matrix)

2023-10-28 Thread Simon
Dear all, the "Constraints on Degrees of Freedom" module introduces several approaches to deal with constraints. I implemented the "second" approach, that is, distributing the constraints on the fly when assembling the linear system using the constraints.distribute_local_to_global() function.