unsigned int ElasticProblem<dim>::solve()
{
SolverControl solver_control(solution.size(), 1e-8 * system_rhs.l2_norm());
- PETScWrappers::SolverCG cg(solver_control, mpi_communicator);
+ PETScWrappers::SolverCG cg(solver_control);
PETScWrappers::PreconditionBlockJacobi preconditioner(system_matrix);
SolverControl solver_control(dof_handler.n_dofs(),
1e-16 * system_rhs.l2_norm());
- PETScWrappers::SolverCG cg(solver_control, mpi_communicator);
+ PETScWrappers::SolverCG cg(solver_control);
PETScWrappers::PreconditionBlockJacobi preconditioner(system_matrix);
PETScInverse(const dealii::PETScWrappers::MatrixBase &A,
dealii::SolverControl & cn,
const MPI_Comm &mpi_communicator = PETSC_COMM_SELF)
- : solver(cn, mpi_communicator)
+ : solver(cn)
, matrix(A)
, preconditioner(matrix)
{}
PETScInverse(const dealii::PETScWrappers::MatrixBase &A,
dealii::SolverControl & cn,
const MPI_Comm &mpi_communicator = PETSC_COMM_SELF)
- : solver(cn, mpi_communicator)
+ : solver(cn)
, matrix(A)
, preconditioner(matrix)
{}
{
pcout << "...solving" << std::endl;
SolverControl solver_control(prm.max_iterations, prm.tolerance, false, false);
- PETScWrappers::SolverCG cg(solver_control, mpi_communicator);
+ PETScWrappers::SolverCG cg(solver_control);
PETScWrappers::PreconditionSOR preconditioner(system_matrix);
SolverControl solver_control(dof_handler.n_dofs(), 1e-12);
- PETScWrappers::SolverCG solver(solver_control, mpi_communicator);
+ PETScWrappers::SolverCG solver(solver_control);
#ifndef PETSC_USE_COMPLEX
PETScWrappers::PreconditionBoomerAMG preconditioner(
SolverControl solver_control(dof_handler.n_dofs(), 1e-12);
- PETScWrappers::SolverCG solver(solver_control, mpi_communicator);
+ PETScWrappers::SolverCG solver(solver_control);
#ifndef PETSC_USE_COMPLEX
PETScWrappers::PreconditionBoomerAMG preconditioner(
SolverControl solver_control(dof_handler.n_dofs(), 1e-12, false, false);
- PETScWrappers::SolverCG solver(solver_control, mpi_communicator);
+ PETScWrappers::SolverCG solver(solver_control);
#ifndef PETSC_USE_COMPLEX
// Ask for a symmetric preconditioner by setting the first parameter in
SolverControl solver_control(dof_handler.n_dofs(), 1e-12);
- PETScWrappers::SolverCG solver(solver_control, mpi_communicator);
+ PETScWrappers::SolverCG solver(solver_control);
#ifndef PETSC_USE_COMPLEX
PETScWrappers::PreconditionBoomerAMG preconditioner(
SolverControl solver_control(dof_handler.n_dofs(), 1e-12);
- PETScWrappers::SolverCG solver(solver_control, mpi_communicator);
+ PETScWrappers::SolverCG solver(solver_control);
#ifndef PETSC_USE_COMPLEX
PETScWrappers::PreconditionBoomerAMG preconditioner(
SolverControl solver_control(dof_handler.n_dofs(), 1e-12);
- PETScWrappers::SolverCG solver(solver_control, mpi_communicator);
+ PETScWrappers::SolverCG solver(solver_control);
#ifndef PETSC_USE_COMPLEX
PETScWrappers::PreconditionBoomerAMG preconditioner(
SolverControl solver_control(dof_handler.n_dofs(), 1e-12);
- PETScWrappers::SolverCG solver(solver_control, MPI_COMM_WORLD);
- PETScWrappers::PreconditionBDDC<2> preconditioner;
+ PETScWrappers::SolverCG solver(solver_control);
+ PETScWrappers::PreconditionBDDC<2> preconditioner;
PETScWrappers::PreconditionBDDC<2>::AdditionalData data;
// Now we setup the dof coordinates if a sufficiently new PETSc is used
1e-16 * system_rhs.l2_norm(),
false,
false);
- PETScWrappers::SolverCG cg(solver_control, mpi_communicator);
+ PETScWrappers::SolverCG cg(solver_control);
PETScWrappers::PreconditionBlockJacobi preconditioner(system_matrix);
cg.solve(system_matrix,
distributed_incremental_displacement,
1e-16 * system_rhs.l2_norm(),
false,
false);
- PETScWrappers::SolverCG cg(solver_control, mpi_communicator);
+ PETScWrappers::SolverCG cg(solver_control);
PETScWrappers::PreconditionBlockJacobi preconditioner(system_matrix);
cg.solve(system_matrix,
distributed_incremental_displacement,
ElasticProblem<dim>::solve()
{
SolverControl solver_control(solution.size(), 1e-8 * system_rhs.l2_norm());
- PETScWrappers::SolverCG cg(solver_control, mpi_communicator);
+ PETScWrappers::SolverCG cg(solver_control);
// PreconditionBlockJacobi depends on the processor count
PETScWrappers::PreconditionBlockJacobi preconditioner(system_matrix);
1e-16 * system_rhs.l2_norm());
#ifdef DEAL_II_WITH_PETSC
- PETScWrappers::SolverCG cg(solver_control, mpi_communicator);
+ PETScWrappers::SolverCG cg(solver_control);
PETScWrappers::PreconditionBlockJacobi preconditioner(system_matrix);
1e-15,
/*log_history*/ false,
/*log_results*/ false);
- PETScWrappers::SolverCG linear_solver(linear_solver_control,
- mpi_communicator);
+ PETScWrappers::SolverCG linear_solver(linear_solver_control);
linear_solver.initialize(*preconditioner);
dealii::SolverControl solver_control(100,
1e-15,
/*log_history*/ false,
/*log_results*/ false);
- PETScWrappers::SolverCG linear_solver(linear_solver_control,
- mpi_communicator);
+ PETScWrappers::SolverCG linear_solver(linear_solver_control);
linear_solver.initialize(*preconditioner);
dealii::SolverControl solver_control(100,
1e-15,
/*log_history*/ false,
/*log_results*/ false);
- PETScWrappers::SolverCG linear_solver(linear_solver_control,
- mpi_communicator);
+ PETScWrappers::SolverCG linear_solver(linear_solver_control);
linear_solver.initialize(*preconditioner);
dealii::SolverControl solver_control(100,