Re: [deal.II] log of a symmetric positive definite tensor

2019-10-15 Thread Wolfgang Bangerth
On 10/15/19 5:17 PM, Maien Hamed wrote: > > But this is not reliable, as I've run into situations where the returned > eigenvalues are fewer than `dim` (I'm not sure how > `compute_eigenvalues_symmetric` treats repeated eigenvalues). I've looked at > using `compute_svd`, but it seems to me that

[deal.II] log of a symmetric positive definite tensor

2019-10-15 Thread Maien Hamed
Hello everyone, Is there a reliable way to get the natural logarithm of a symmetric positive definite tensor in deal.ii? I'm trying `LAPACKFullMatrix.compute_eigenvalues_symmetric`, taking the `std::log` of the eigenvalues and then multiplying by the eigenvectors appropriately: const Symmetric

[deal.II] Re: Particle contact detection

2019-10-15 Thread Rene Gassmöller
Hi Shahab, Adding to Luca's reply you probably want to look into the `ParticleHandler` class, in particular the `particles_in_cell` function that returns you all particles in a given cell (https://dealii.org/developer/doxygen/deal.II/classParticles_1_1ParticleHandler.html#acaf1232ffce0746baa641

Re: [deal.II] Error installing with Candi

2019-10-15 Thread David Ryan
Hi, That's unfortunate. Thank you for the reply Bye! David On Tue, 15 Oct 2019, 01:12 'Uwe Köcher' via deal.II User Group, < dealii@googlegroups.com> wrote: > Dear David, > > sorry, macOS Mojave (and further Catalina) do not have candi support - > there is something strange happening > in the

Re: [deal.II] Error installing with Candi

2019-10-15 Thread 'Uwe Köcher' via deal . II User Group
Dear David, sorry, macOS Mojave (and further Catalina) do not have candi support - there is something strange happening in the recent past with the base compilers and homebrew installed mpi compilers. (I had a similar problem as you reported) Currently I recommend to install the deal.II toolch

[deal.II] assemble matrix after mesh refinement

2019-10-15 Thread Mauro Murer
Hello everyone! I'm trying to rewrite Step-35 (Navier-Stokes tutorial) without splitting the implementation for component (In Step 35 one thread resolve the x component and the other one the y component). Proceeding in this way I can use "AffineConstraints" class to manipulate constraints in