Thanks Prof. Bangerth and Bruno for your valuable suggestions. It is much
more clear to me now.
On Wednesday, July 19, 2023 at 5:44:09 AM UTC+5:30 blais...@gmail.com wrote:
> In addition to Wolfgang's answer. Here are some suggestions.
>
> 1. You system is not symmetric. Consequently, it is no
In addition to Wolfgang's answer. Here are some suggestions.
1. You system is not symmetric. Consequently, it is not a viable idea to
use the CG linear solver
(see https://en.wikipedia.org/wiki/Conjugate_gradient_method). As
suggested, using something like GMRES would be a better idea for your
On 7/18/23 00:44, sabyasachi chatterjee wrote:
Note that my stiffness matrix is not symmetric ( system_matrix.is_symmetric()
is returning 0 ).
Then there is your mistake. You can only use the CG solver for symmetric and
positive definite matrices. You will want to use GMRES for other matrices.
Hello,
I have come up with a solution, although there may be other ones using
iterative solvers. I changed the solver to SparseDirectMUMPS (from tutorial
62). It seems to have stopped generating the error. Note that my stiffness
matrix is not symmetric ( system_matrix.is_symmetric() is retur
I just wanted to correct one small thing. What I mean by relative error
before is the residual norm, which is a relative number.
On Thursday, July 13, 2023 at 12:57:23 PM UTC+5:30 sabyasachi chatterjee
wrote:
> Hello,
>
> I am trying to solve a nonlinear contact problem using penalty method.
