Pb with conditional Dirichlet
Posted: 03 Apr 2018, 19:10
Hello;
It seems that Solver % Matrix % ConstrainedDOF is never reinitailised to FALSE so that when a node get a Dirichlet condition, it will always apply a Dirichlet condition. This could be a problem when using a "Condition" for the dirichlet condition and when the condition change from TRUE to FALSE.
A fix for the default beahaviour is:
Line 4384 of DefUtils.F90 (in DefaultDirichletBCs):
A % ConstrainedDOF = .FALSE.
should come after the ENDIF so that it will be re-initialized at each call of DefaultDirichletBCs;
However this do not solve the issue for solvers that do not use the DefaultDirichletBCs routine.
best regards,
Fabien
It seems that Solver % Matrix % ConstrainedDOF is never reinitailised to FALSE so that when a node get a Dirichlet condition, it will always apply a Dirichlet condition. This could be a problem when using a "Condition" for the dirichlet condition and when the condition change from TRUE to FALSE.
A fix for the default beahaviour is:
Line 4384 of DefUtils.F90 (in DefaultDirichletBCs):
A % ConstrainedDOF = .FALSE.
should come after the ENDIF so that it will be re-initialized at each call of DefaultDirichletBCs;
However this do not solve the issue for solvers that do not use the DefaultDirichletBCs routine.
best regards,
Fabien