[deal.II] Issue with MKL 2023.1.0

Hello! I have installed deal.II with MKL 2023.1.0: -DDEAL_II_WITH_LAPACK=ON -DBLAS_DIR=${MKLROOT} -DLAPACK_DIR=${MKLROOT} However, when I run the step-1.release progoram, I encountered an error: ./step-1.release: symbol lookup error: /opt/intel/oneapi/mkl/2023.1.0/lib/intel64/libmkl_gf_lp64.so.2:

[deal.II] Missing documentation for Tensor class?

Hi Everyone, I am not 100% certain of this, but it looks like there is some documentation missing from the Tensor class for v 9.5.0 - the first entry in the contraction section is for the l1-norm, although some of the text appea

Re: [deal.II] construct TrilinosWrappers::SparseMatrix from TrilinosWrappers::BlockSparseMatrix to use Amesos direct solver

I have some experience with the direct solvers of Trilinos and it requires very few modifications to change the assembly to a SparseMatrix instead of a BlockSparseMatrix for vector valued problems. Especially the assembly with distribute_local_to_global does not actually care apart from the matr

Re: [deal.II] Issue with PETSc vector

On 12/18/23 13:23, Giselle Sosa Jones wrote: // Initialize vector dof_handler_dg0.distribute_dofs(fe_dg0); const std::vector locally_owned_dofs_per_proc_dg0 = DoFTools::locally_owned_dofs_per_subdomain(dof_handler_dg0); locally_owned_dofs_dg0 = locally_owned_dofs_per_proc_dg0[this_mpi_process];

Re: [deal.II] construct TrilinosWrappers::SparseMatrix from TrilinosWrappers::BlockSparseMatrix to use Amesos direct solver

On 12/18/23 12:44, Simon wrote: I use a BlockSolver for the solution and right hand side of the linear system to compute the various norms (non-linear solver) of the the individual blocks. For a direct solver, I would now use a SparseMatrix as matrix type and a BlockVector as vector type. At e

[deal.II] Issue with PETSc vector

Hello, I have an issue that seems very simple but I cannot figure out. I have a piecewise constant parallel vector that I defined using a DG-0 space. I have an auxiliary function that calculates the value according to the cell, and then in my main file, I would like to assign the value to the v

Re: [deal.II] construct TrilinosWrappers::SparseMatrix from TrilinosWrappers::BlockSparseMatrix to use Amesos direct solver

"This. I don't see a different way of doing things. Though I might just always create both objects and then call distribute_local_to_global() with one or the other matrix object, depending on your run-time flags. " Your approach sounds more reasonable than copying data back and forth at ever

Re: [deal.II] construct TrilinosWrappers::SparseMatrix from TrilinosWrappers::BlockSparseMatrix to use Amesos direct solver

On 12/18/23 04:47, Simon wrote: Can you think of an elegant solution or must the non-zero values be copied via a loop over the individual blocks? This. I don't see a different way of doing things. Though I might just always create both objects and then call distribute_local_to_global() with

[deal.II] construct TrilinosWrappers::SparseMatrix from TrilinosWrappers::BlockSparseMatrix to use Amesos direct solver

Dear all, there are several discussions, e.g. https://groups.google.com/g/dealii/c/KU7g_YCNVRE/m/h7ab-2gRCAAJ?hl=en or https://groups.google.com/g/dealii/c/DdMiTjuPx6Q/m/UeH4fmmBAwAJ?hl=en about the usage of TrilinosWrappers::BlockSparseMatrix in the context of direct solvers in a distributed