[deal.II] Transfer local cell data to global cell vector

2017-03-17 Thread Sumedh Yadav
Hello, I need to do some post-processing on computed solution data (say Vector Solution(dof_handler.n_dofs()) is the solution vector). For post-processing field (say Vector Curvature(triangulation.n_dofs()) is the scalar field to be computed) I decided to have one entry per cell. Once I calcula

Re: [deal.II] Solver convergence failure with L2 residual diverging for SolverBicgstab

2017-03-17 Thread Vinayak Gholap
Hello, Thank you for your suggestions. In essence, what happens here is that your solver has just diverged. > This may be for three reasons: > * The matrix or right hand side already contain INF or NaN values. You >can test this by printing the norm of the matrix and rhs before you >se

Re: [deal.II] Re: AMR , how to pass solution vector to refined mesh

2017-03-17 Thread Jaekwang Kim
After the first cycle, I am doing template void nonlinear::refine_grid () { Vector estimated_error_per_cell (triangulation.n_active_cells()); KellyErrorEstimator::estimate (dof_handler, QGauss(3), typename F

[deal.II] How to compute Nédélec gradient

2017-03-17 Thread justtry
Hi, I am a new deal.II user. I'm solving a 2D problem on quadrilateral cartesian mesh using FE_Q(1) elements. I would like to compute the gradient of the bilinear solution, known to belong to the FE_Nedelec(0) space: how can I manage to do it? In particular I'm only interested in computing the

Re: [deal.II] Evaluate solution of face elements

2017-03-17 Thread Wolfgang Bangerth
On 03/14/2017 10:54 AM, l.j.corbijnvanwillenswa...@student.utwente.nl wrote: I'm trying to evaluate the solution of face elements (FE_FaceQ) at specific points on a face. For uniformly refined meshes the following works perfectly (as far as I can tell) std::vector> result_values (1, Vector(1));

Re: [deal.II] Re: AMR , how to pass solution vector to refined mesh

2017-03-17 Thread Wolfgang Bangerth
On 03/16/2017 08:38 PM, Jaekwang Kim wrote: I think this two function do everything I needed... but when I run my code, I run into error message when I assemble system in second time. what might be the problem.? As the error message says, you are trying to write into a matrix entry that d

[deal.II] Re: Installation using Spack fails during 'ncurses' installation

2017-03-17 Thread Stephen DeWitt
Bruno, Jean-Paul, and Denis, Thank you all for your responses. As suggested, I'll send a message to the Spack forum and/or open an issue on their GitHub site. In the meantime, I'll continue debugging the manual PETSc installation to see if I can get that to work. I think the filesystem I'm tryin

[deal.II] Re: Installation using Spack fails during 'ncurses' installation

2017-03-17 Thread Denis Davydov
p.p.s. Just to be sure, i wiped my installations of Spack on Ubuntu16.04+gcc5.4.0 PC and Centos7+gcc4.8.5 cluster and reinstalled everything from scratch using commit 1124bdc99ee84c26201c40536d9b04dac74d7f6a (this is the current HEAD in Spack, I updated dealii's wiki to mention it). I still t

[deal.II] Re: Installation using Spack fails during 'ncurses' installation

2017-03-17 Thread Denis Davydov
On Thursday, March 16, 2017 at 6:43:13 PM UTC+1, Bruno Turcksin wrote: > > Steve, > > I have tried to use spack several times on two clusters and it never > worked for me (but it works fine on my own machine). I usually have to > patch a bunch of things and at the > p.s. maybe a good idea is to

[deal.II] Re: Installation using Spack fails during 'ncurses' installation

2017-03-17 Thread Denis Davydov
@Stephen DeWitt I think the best place to report such problems is https://github.com/LLNL/spack/issues Looking at your error, i think it happens within your cluster, as it's ldconfig which complains /sbin/ldconfig: Can't create temporary cache file /etc/ld.so.cache~: > Permission denied Bu