Commit 1ccc9848 authored by Robert K's avatar Robert K

[!316] [bugfix][PetscLinearOperator] There seems to be a bug with preallocation.

Merge branch 'bugfix/petsc-linear-operator-setup' into 'master'

ref:dune-fem/dune-fem This is a temporary quick-fix of the matrix memory and
vector valued spaces but block size 1 for the matrix.

See merge request [!316]

  [!316]: gitlab.dune-project.org/dune-fem/dune-fem/merge_requests/316
parents d2f715a7 fd9fc3e4
Pipeline #18423 passed with stage
in 43 minutes and 59 seconds
......@@ -319,7 +319,7 @@ namespace Dune
::Dune::Petsc::MatSetType( petscMatrix_, MATAIJ );
}
if( isSimpleStencil || std::is_same< StencilType,SimpleStencil<DomainSpaceType,RangeSpaceType> >::value )
if( ! blockedMode_ || isSimpleStencil || std::is_same< StencilType,SimpleStencil<DomainSpaceType,RangeSpaceType> >::value )
{
::Dune::Petsc::MatSetUp( petscMatrix_, bs, domainLocalBlockSize * stencil.maxNonZerosEstimate() );
}
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment