[deal.II] problem adding two blocks of BlockSparseMatrix

Dear All, I need to subtract two blocks of a BlockSparseMatrix system_matrix. I am trying to achieve this by creating SpraseMatrix tmp1, then tmp1.copy_from(system_matrix.block(0,0)); tmp1.add(-1.0, system_matrix.block(0,1)); However, I am getting the following error message: "When copying

Re: [deal.II] Re: Status of support for complex equations?

On Tue, Jul 25, 2017, at 14:24 CDT, Denis Davydov wrote: > which Matthias wanted to do over the weekend in July 2016 :-) <_< Long weekend! *cough* Best, Matthias -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see

[deal.II] Re: Status of support for complex equations?

Have a look at this issue which aggregates efforts in this direction: https://github.com/dealii/dealii/issues/2033 Essentially, the main missing piece is ConstraintMatrix https://github.com/dealii/dealii/issues/1760 which Matthias wanted to do over the weekend in July 2016 :-) I think there is

[deal.II] Re: website petsc tutorial

Hi Yi-Chung, We typically update that documentation for the versions that we specifically test against, so unless the documentation was accidentally not updated then it is best to assume that the version mentioned on the website is the most recent version that is, to the best of our knowledge,

[deal.II] Re: website petsc tutorial

Hi Jean-Paul, Thank you for explaining it. I have the question because the page shows deal.ii is compatible with petsc 3.6.0. As I may recall, 8.4.1 is compatible with 3.6.4, right? All best, YC Chen On Tuesday, July 25, 2017 at 3:50:41 PM UTC+1, Jean-Paul Pelteret wrote: > > Dear

[deal.II] Re: website petsc tutorial

Dear Yi-Chung, Thanks for the information, but this is done by design (so there is no problem). At the time of the release of deal.II 8.5.0, the content linked to on http://www.dealii.org/8.5.0/external-libs/petsc.html

[deal.II] Re: Parallel GridIn of large file - PETSC Error

Hi Mathias Would you mind to check the memory usage? I guess you have about 3M DOFs right (based on the file size of msh)? My experience told me that 32GB can handle 10M DOFs(roughly). Your case is about 10GB memory based on the calculation.If you are using distributed mesh, which means you

[deal.II] website petsc tutorial

Hi all, one link errors of the website regarding PETSC installation (I assume) if visiting http://www.dealii.org/8.5.0/readme.html and click More information on configuring and building PETSc can be found here .<= this link it will

Re: [deal.II] Mapping points of mesh generated by gmsh

Dear Jean-Paul, (I tried also version 8.4.1 and same problem). Ok, thanks a lot! Best, Giovanni On Tuesday, July 25, 2017 at 9:08:53 AM UTC+2, Jean-Paul Pelteret wrote: > > Dear Giovanni, > > Thank you very much for the clear and simple test case! I've played around > with it and I don't see

Re: [deal.II] Mapping points of mesh generated by gmsh

Dear Giovanni, Thank you very much for the clear and simple test case! I've played around with it and I don't see any fault in the way that you're trying to use the mapping functionality, so I think that its quite clear that there's something wrong with Mapping::transform_real_to_unit_cell.