Yup, I use petsc. Thanks a lot.
Miguel
On Fri, Sep 5, 2014 at 10:40 AM, Roy Stogner
wrote:
>
>
> On Fri, 5 Sep 2014, Miguel Angel Salazar de Troya wrote:
>
> Now it's working. I added the vectors I mentioned before to the system
>> using the routine System::add_vector() and the cg works after
On Fri, 5 Sep 2014, Miguel Angel Salazar de Troya wrote:
Now it's working. I added the vectors I mentioned before to the system using
the routine System::add_vector() and the cg works after the
refinement. It's weird because it was working before after refinement when
using gmres.
Could you
Now it's working. I added the vectors I mentioned before to the system
using the routine System::add_vector() and the cg works after the
refinement. It's weird because it was working before after refinement when
using gmres.
Could you help me out with the preconditioners please? Just to mention ho
On Thu, 4 Sep 2014, Miguel Angel Salazar de Troya wrote:
> I'm using a FEMSystem with a NewtonSolver. I realized that in the FEMSystem,
> in the assembly routines, DofMap::constrain_element_matrix_and_vector
> is called with asymmetric_constraint_rows = false. Also, because it's a
> nonlinear s
I'm using a FEMSystem with a NewtonSolver. I realized that in the
FEMSystem, in the assembly routines,
DofMap::constrain_element_matrix_and_vector is called with
asymmetric_constraint_rows = false. Also, because it's a nonlinear solver.
DofMap::enforce_constraints_exactly(system) is called as well.
Could you please tell me how to set up the conditioner like you mention?
AMG on each processor wins on problems where it works, ILU (usually
with ~4 levels) wins otherwise, tweaking Block Jacobi vs Additive
Schwarz to tie it together can speed things up slightly.
What are the command line options
On Thu, 4 Sep 2014, Miguel Angel Salazar de Troya wrote:
> Thanks. Could you answer the last question please? Basically I want to know
> if I still can use cg once I use mesh refinement? I did quick
> calculations with equation (7) and I see the matrix is not symmetric so I
> think I cannot us
Thanks. Could you answer the last question please? Basically I want to know
if I still can use cg once I use mesh refinement? I did quick calculations
with equation (7) and I see the matrix is not symmetric so I think I cannot
use it. I would like to double check with you. Which preconditioner woul
On Wed, 3 Sep 2014, Miguel Angel Salazar de Troya wrote:
> In this paper
>
> http://arxiv.org/pdf/1009.4975.pdf
>
> In page 9, at the bottom, they mention this:
>
> As indicated by the circled nodes in Figure 3(a), level-one mesh
> incompatibility results in hanging nodes. The libMesh package ha
Hello all
In this paper
http://arxiv.org/pdf/1009.4975.pdf
In page 9, at the bottom, they mention this:
As indicated by the circled nodes in Figure 3(a), level-one mesh
incompatibility results in hanging nodes. The libMesh package handles those
hanging nodes by using the projection method to en
10 matches
Mail list logo
