[deal.II] Adding inhomogeneity to Periodic BC

2021-10-22 Thread Raghunandan Pratoori
Hello, I have a cube with periodicity in all 3 directions. I want to add inhomogeneity in addition to Periodic BC on one face. Below is what I have done - { IndexSet selected_dofs_z; std::set boundary_ids_z = std::set(); boundary_ids_x.insert(4);

Re: [deal.II] Re: Trilinos AMG ML .reinit() ... how to?

2021-10-22 Thread Wolfgang Bangerth
On 10/22/21 2:53 AM, richard@gmail.com wrote: I think I know what it is: the ML guide on p 16 says in footnote 11: "Note that the hierarchy produced by ReComputePreconditioner() can differ from that produced by ComputePreconditioner() for non-zero threshold values." So when the matrix chang

Re: [deal.II] Efficient Matrix-based Calculation

2021-10-22 Thread Jean-Paul Pelteret
Dear Masoud, > As you mentioned, we expect this to be slower than the equivalent code in > step-8; however, it is very slower than the same code that I wrote before in > Matlab. A quick question: are you sure that you’re running your problem in release mode, rather than debug mode? I had a col

Re: [deal.II] Missing (block_sparse_matrix)

2021-10-22 Thread Daniel Arndt
Ahmed, The file and the respective classes were deleted in https://github.com/dealii/dealii/pull/9761 and deprecated in https://github.com/dealii/dealii/pull/5352. If you need to replace BlockMatrixArray, you should have a look at the documentation in the latter pull request to use BlockLinearOper

[deal.II] Re: Trilinos AMG ML .reinit() ... how to?

2021-10-22 Thread richard....@gmail.com
HA! I think I know what it is: the ML guide on p 16 says in footnote 11: "Note that the hierarchy produced by ReComputePreconditioner() can differ from that produced by ComputePreconditioner() for non-zero threshold values." So when the matrix changes too much, the const double aggregation_thres

[deal.II] Missing (block_sparse_matrix)

2021-10-22 Thread Ahmed Galal
I am trying to compile a deal.ii code but it failed because of missing . I do not know how to include it in the recent dealii version. I found it available deal.ii 8.5.0. Thanks in advance! -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https