Re: [deal.II] MatrixTools::apply_boundary_values for BlockSparseMatrix> and BlockVector> entries

2021-09-20 Thread Мария Бронзова
Dear Daniel, Thank you many time, I changed the source files, as suggested, and recompiled the whole thing again. Now it seems to be working! Thank you again! With kind regards, Mariia Bronzova пятница, 17 сентября 2021 г. в 18:23:42 UTC+2, d.arnd...@gmail.com: > Mariia, > > I also tried

Re: [deal.II] MatrixTools::apply_boundary_values for BlockSparseMatrix> and BlockVector> entries

2021-09-17 Thread Daniel Arndt
Mariia, I also tried fixing this within deal.II and discovered that some more changes are necessary, see https://github.com/dealii/dealii/pull/12761. Would you mind checking if the changes proposed there work for you? Best, Daniel Am Fr., 17. Sept. 2021 um 06:25 Uhr schrieb Мария Бронзова <

Re: [deal.II] MatrixTools::apply_boundary_values for BlockSparseMatrix> and BlockVector> entries

2021-09-17 Thread Мария Бронзова
Dear Daniel, Thank you a lot for the solution ways! I was trying to modify the /dealii-9.3.1/source/numerics/matrix_tools.inst.in as follows: // - // // Copyright (C) 2010 - 2018 by the deal.II authors // // This file is part

Re: [deal.II] MatrixTools::apply_boundary_values for BlockSparseMatrix> and BlockVector> entries

2021-09-15 Thread Daniel Arndt
Mariia, That's a problem a little more complicated to solve. Basically, we are not instantiating the template for complex and BlockSparseMatrix. If you are OK with recompiling deal.II, you should modify source/numerics/ matrix_tools.inst.in so that all overloads are also compiled for complex

Re: [deal.II] MatrixTools::apply_boundary_values for BlockSparseMatrix> and BlockVector> entries

2021-09-15 Thread Мария Бронзова
Dear Daniel, Thank you a lot for the hint, it makes perfect sense! I also correctted the* interpolate_boundary_values* call as follows: *templateclass DirichletBoundaryValues : public Function> { public: DirichletBoundaryValues(): Function>(1) {} virtual

Re: [deal.II] MatrixTools::apply_boundary_values for BlockSparseMatrix> and BlockVector> entries

2021-09-15 Thread Мария Бронзова
Dear Daniel, Thank you a lot for the hint, makes perfect sense! The MatrixTools::apply_boundary_values works well now! But I believe i also need to adjust the VectorTools::interpolate_boundary_values accordingly. Right now i have defined a template for pressure component boundary values:

Re: [deal.II] MatrixTools::apply_boundary_values for BlockSparseMatrix> and BlockVector> entries

2021-09-14 Thread Daniel Arndt
Mariia, You need to make sure that the std::map also uses values of type std::complex. The compile error you shared shows that you use doubles. Best, Daniel Am Di., 14. Sept. 2021 um 07:30 Uhr schrieb Мария Бронзова < masianic...@gmail.com>: > Hi everyone, > > I ran into the following

[deal.II] MatrixTools::apply_boundary_values for BlockSparseMatrix> and BlockVector> entries

2021-09-14 Thread Мария Бронзова
Hi everyone, I ran into the following difficulty: I have a vector-valued problem with complex numbers as entries for my BlockSparseMatrix system_matrix and BlockVector system_rhs and solution vectors and I am trying to apply Dirichlet boundary conditions, using the