Re: [deal.II] Derivatives (gradients) of Symmetric Tensors

2016-09-29 Thread Metin C.
Dear Wolfgang, I would like to implement the gradients (or symmetric_gradient) and laplacians functions for the symmetric_tensor class. I need help to figure out some functions in fe_values.h and fe_values.cc sources. Could you please help me understand what's happening in the divergence

[deal.II] Re: Vector-valued gradient of solution vector

2016-09-29 Thread Daniel Arndt
Joel, Yes, this looks about right. I would really replace "std::vector double x,y,z" by "Vector x,y,z". In the end, entries in these vectors only make sense with the corresponding DoFHandler. This would also allow you to use these Vectors for postprocessing in deal.II. Best, Daniel -- The

[deal.II] Relation between Solution Error Behavior and Polynomial Approximation Degree

2016-09-29 Thread JAEKWANG KIM
Hi all, I have question on error behavior of FEM. I thought that the order of error is O(h^p) where h is a mesh-size and p is polynomial degree we use in

Re: [deal.II] Derivatives (gradients) of Symmetric Tensors

2016-09-29 Thread Wolfgang Bangerth
Metin, I need help to figure out some functions in fe_values.h and fe_values.cc sources. Could you please help me understand what's happening in the divergence function (SymmetricTensor<2, dim, spacedim>::divergence) and the symmetric_gradient function for the vectors

Re: [deal.II] Neumann bc in step-20

2016-09-29 Thread Erik Svensson
Ok, thanks Wolfgang and Daniel! I will follow your advise and post updates on this issue. /Erik 2016-09-30 5:54 GMT+02:00 Wolfgang Bangerth : > On 09/28/2016 11:32 PM, Erik Svensson wrote: > >> >> >> I’m trying to implement homogeneous Neumann bc on part of the boundary

Re: [deal.II] Neumann bc in step-20

2016-09-29 Thread Wolfgang Bangerth
On 09/28/2016 11:32 PM, Erik Svensson wrote: I’m trying to implement homogeneous Neumann bc on part of the boundary in tutorial example step-20. The Neumann bc should be imposed strongly (as for Dirichlet bc in normal, non-mixed, formulation). Is this correct? In any case, I’m trying to

[deal.II] Re: Trilinos SparseMatrix mmult bug

2016-09-29 Thread Ben Shields
So this turned out to be a bug in my own code after I made the original change in trilinos_sparse_matrix.cc. When I was constructing the sparsity patterns for the matrices using DoFTools::make_sparsity_pattern, there is a flag I had set keep_constrained_dofs = false. Setting this flag to true

[deal.II] Re: macOS Sierra dynamic linker limit

2016-09-29 Thread Denis Davydov
after digging more i found out that it has to do with extra not needed LC_RPATH load commands added by Spack https://github.com/LLNL/spack/issues/1874 On Wednesday, September 28, 2016 at 5:58:07 PM UTC+2, Denis Davydov wrote: > > I am hitting run-time error of the form: > > 36: dyld: Library not

[deal.II] HR event for developers seeking new projects in Berlin

2016-09-29 Thread ngimeno
Hi everyone! we are launching our first event Talents@Unicorn.Berlin and we are currently looking for developers who are seeking jobs in the IT industry in Berlin. You can register and attend for free our event on the 1st of November and meet with selected top companies (Babbel

[deal.II] a few questions about matrix-free

2016-09-29 Thread Denis Davydov
Dear all, I have accumulated a few matrix-free related questions which i would like to clarify ask. (1) FEEvaluation has read_dof_values (const std::vector< VectorType > , const unsigned int first_index=0), but i don’t see a function to get to values of each DoF vector at quadrature points.

[dealii-developers] HR event for developers seeking new projects in Berlin

2016-09-29 Thread ngimeno
Hi everyone! we are launching our first event Talents@Unicorn.Berlin and we are currently looking for developers who are seeking jobs in the IT industry in Berlin. You can register and attend for free our event on the 1st of November and meet with selected top companies (Babbel