[deal.II] Re: Assembling sparse matrix from Matrix-free vmult with constrains

2020-03-09 Thread 'peterrum' via deal.II User Group
Hi Michal, any chance that you post or send me a small runnable test program. By the way, there is an open PR (https://github.com/dealii/dealii/pull/9343) for the computation of the diagonal in a matrix-free manner. Once this is merged, I will work the matrix-free assembly of sparse matrices.

Re: [deal.II] Compiler linking

2020-03-09 Thread Matthias Maier
If you want to use binary packages I highly recommend to use the PPA here in order to get the newest version: https://launchpad.net/~ginggs/+archive/ubuntu/deal.ii-backports After that make sure that you actually install the developer version. All in one: sudo add-apt-repository ppa:ginggs/d

[deal.II] Compiler linking

2020-03-09 Thread Robert Kopp
I downloaded the source code for deal.II in Ubuntu, which contains the example files. (I'm a beginner.) I noted, however, that a binary is available for Ubuntu (18.04) and installed that, instead of compiling the library. Then I attempted to run step-1.cc from the source code package, after det

[deal.II] deal.II Newsletter #113

2020-03-09 Thread Rene Gassmoeller
Hello everyone! This is deal.II newsletter #113. It automatically reports recently merged features and discussions about the deal.II finite element library. ## Below you find a list of recently proposed or merged features: #9642: Add missing include file for step-69 (proposed by kronbichler; m

[deal.II] Assembling sparse matrix from Matrix-free vmult with constrains

2020-03-09 Thread Michał Wichrowski
Dear all, I've got matrix-free multigrid solver for Stokes problem. The main bottleneck is solution of coarse problem, so I tried to assemble the regular sparse matrix and use direct solver. Since the coarse problem is (relatively) small, I used vmults by unit vector to obtain columns of the m

[deal.II] Matrix-free: Continuous elements and boundary (face) integrals

2020-03-09 Thread Michał Wichrowski
Deal all, I want to weakly impose boundary conditions on, let's say, Laplace problem (in fact I'm working on Stokes). I'm considering [FE_Q<2>(2)^2] ( that is FESystem<2>[FE_Q<2>(2)^2] + FE_Q<2>(1) for Stokes). As a MatrixFree loop I'm using: this->data->loop(&SymGradMass::local_apply, &Sy