Re: [deal.II] FullMatrix to SparseMatrix

On 10/11/20 6:09 PM, Nikki Holtzer wrote: I was hoping to take the outer product which is formed in a FullMatrix, called Kronecker_matrix below, multiply it by a coefficient, and add a matrix to it called nonlinear_part_full. This entire operation I have tried to put into a Sparse matrix, calle

Re: [deal.II] Re: Fluid-Structure interaction

Thanks Bruno, Jean-Paul, and Thomas. I going to study each answer that you suggest me in detail, looks like it is exactly what I need. El dom., 11 oct. 2020 a las 10:41, blais...@gmail.com (< blais.br...@gmail.com>) escribió: > Dear Ramiro, > You could also take a look at step-70 by Luca and I.

Re: [deal.II] Re: APPLY BOUNDARY CONDTION ON INTERNAL FACE

Xiang, Adding to Bruno's response: If the interface you want to apply Dirichlet boundary conditions to stays the same through the whole simulation, then treating the corresponding faces as disjoint while creating your mesh might be a better option. Setting constraints yourself requires you to know

Re: [deal.II] Issue with MeshWorker::InegrationInfo.values

You will need to tell the MeshWorker to compute the values with something like info.cell_selector.add("solution", true, ..., ...); On Mon, Oct 12, 2020 at 8:18 AM Gregor Pfau wrote: > > Dear Timo, > > first of all thank you for your reply and advice. > > As you suspected, info.values[0][0] and

[deal.II] Re: APPLY BOUNDARY CONDTION ON INTERNAL FACE

Xiang, You cannot apply a boundary condition on an internal face. You need to use the AffineConstraint class to impose the constraints yourself. This is what we use internally to apply Dirichlet boundary condition. Best,

Re: [deal.II] Application of the non-homogeneous boundary conditions in Finite Elasticity (Compressible case)

Hi Jean, Thanks a lot for your detailed answer to my questions. Now I am able to run my program and the newton method is converging. I completely misunderstood the documentation. Also, thank you for opening up the issue. It would really be helpful if the complications involved in making constra

Re: [deal.II] Issue with MeshWorker::InegrationInfo.values

Dear Timo, first of all thank you for your reply and advice. As you suspected, info.values[0][0] and info.values[0] are both of size zero. I have tried this with the example given in my previous question, as well as with the original code. Since the original problem is time dependent, the

Re: [deal.II] Tools for parallel debugging

Thank you, Wolfgang and Daniel. Seems like I will go with the command line then. I was just wondering if people here use Eclipse's PTP which sounded like a good graphical tool. In my case it frequently crashes or simply gets stuck. Best, Konrad On Sunday, October 11, 2020 at 11:35:34 PM UTC+2

Re: [deal.II] Tools for parallel debugging

Oh, seems like I missed that. Thanks you! Konrad On Sunday, October 11, 2020 at 11:35:34 PM UTC+2 Wolfgang Bangerth wrote: > On 10/11/20 3:26 PM, Daniel Arndt wrote: > > > > have a look at > > > https://github.com/dealii/dealii/wiki/Frequently-Asked-Questions#how-do-i-debug-mpi-programs > >