Re: [Libmesh-users] Constrains using penalty

2015-10-06 Thread Вася Бубликов
Works like a charm. Thank you again. Roman 2015-10-06 17:49 GMT+03:00 David Knezevic : > On Tue, Oct 6, 2015 at 10:46 AM, Вася Бубликов wrote: > >> David >> >> >> It pains me to trouble you again but >> dynamic_cast*>(system.matrix); >> returns NULL! >> > > > Make sure you do the cast after equ

Re: [Libmesh-users] Constrains using penalty

2015-10-06 Thread David Knezevic
On Tue, Oct 6, 2015 at 10:46 AM, Вася Бубликов wrote: > David > > > It pains me to trouble you again but > dynamic_cast*>(system.matrix); > returns NULL! > Make sure you do the cast after equation_systems.init(); David > 2015-10-06 17:11 GMT+03:00 David Knezevic : > >> Please reply to the

Re: [Libmesh-users] Constrains using penalty

2015-10-06 Thread Вася Бубликов
David It pains me to trouble you again but dynamic_cast*>(system.matrix); returns NULL! Roman 2015-10-06 17:11 GMT+03:00 David Knezevic : > Please reply to the list as well. > > Instead of creating a new matrix, you should cast system.matrix to a > PetscMatrix and then call MatSetOption on tha

Re: [Libmesh-users] Constrains using penalty

2015-10-06 Thread David Knezevic
Please reply to the list as well. Instead of creating a new matrix, you should cast system.matrix to a PetscMatrix and then call MatSetOption on that. David On Tue, Oct 6, 2015 at 10:06 AM, Вася Бубликов wrote: > David > > Thank you for your answer. Here is a short program to try out your f

Re: [Libmesh-users] Constrains using penalty

2015-10-06 Thread Vasileios Vavourakis
Hi Can I naively ask how "big" the penalty factor introduced here is wrt the mechanical stiffness? Have you confirmed that it is (somehow) at least 5-6 orders of magnitude bigger than the max element of the system matrix? cheers, Vas On 6 October 2015 at 11:13, Вася Бубликов wrote: > Dear Si

Re: [Libmesh-users] Constrains using penalty

2015-10-06 Thread David Knezevic
This is a PETSc error that you're hitting because you're adding new non-zeros in the matrix. PETSc tries to help you by throwing an error if you do this, because allocating new non-zeros (outside the specified sparsity pattern) can slow down your program a lot. The easy fix is indicated in the PET

[Libmesh-users] Constrains using penalty

2015-10-06 Thread Вася Бубликов
Dear Sirs, I'm trying to make Libmesh work with a simple linear elasticity problem using linear systems example 4 as an example. I need to force the top boundary remain flat, i.e. all normal displacements are of the same value to be computed. To obtain this, for example to make the displacements f