Re: [deal.II] how to use solution.

2016년 10월 5일 수요일 오후 6시 55분 19초 UTC+9, Daniel Arndt 님의 말: > > Kyusik, > > As Wolfgang already mentioned, you want to use > FEValues::get_function_values() > to extract from a global FE Vector the local values at the quadrature > points in the cell you are on. > These quadrature points normally

Re: [deal.II] Re: Running deal.ii with METIS

On 10/05/2016 08:41 AM, Bruno Turcksin wrote: Do you know whether automatic > partitioning uses METIS? But still, I think the old step-17 and 18 should > work well with version 8.3... Automatic partitioning does not use METIS and does not work across nodes. I believe we only do the automatic p

Re: [deal.II] Using p-refinement with high order elements

On 10/05/2016 09:03 AM, Deepak Gupta wrote: Since it works for the combination of p=1 and 2, I do not think hanging nodes are a problem. I looked at the hp-paper as well and do not see something that might be going wrong related to the hanging nodes. I am still trying to narrow down the problem

[deal.II] Using p-refinement with high order elements

Dear All, I am trying to solve an elasticity problem and am dealing with p-refinement. With global refinement, my solution is as expected (correct). I am using an initial mesh of 4X2 elements. If I keep p=1 for all the elements, and make p = 2 for the 7th element, the solution is still correct

Re: [deal.II] Re: Running deal.ii with METIS

2016-10-05 10:19 GMT-04:00 Chenchen Liu : > For "#PBS -l nodes=1:ppn=x", > when x=1, and use "mpirun -np 1 ./step-17", it works well as a single > processor. > when x=1, and use "mpirun -np 4 ./step-17", memory error happens. > For "#PBS -l nodes=4:ppn=1", > when I use "mpirun -np 1 ./step-17", it

[deal.II] Re: Nonlinear Material

Lisa On Wednesday, October 5, 2016 at 10:17:17 AM UTC-4, Lisa Collins wrote: > I would like to solve a problem with material non-linearity. > Can you please tell me which tutorial steps could be the base of my code. > Take a look at step 15 (http://dealii.org/developer/doxygen/deal.II/step_15.ht

Re: [deal.II] Re: Running deal.ii with METIS

Hi Bruno, Thank you for the reply. For "#PBS -l nodes=1:ppn=x", when x=1, and use "mpirun -np 1 ./step-17", it works well as a single processor. when x=1, and use "mpirun -np 4 ./step-17", memory error happens. For "#PBS -l nodes=4:ppn=1", when I use "mpirun -np 1 ./step-17", it works well as

[deal.II] Nonlinear Material

Hello, I would like to solve a problem with material non-linearity. Can you please tell me which tutorial steps could be the base of my code. Thanks. Lisa -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=

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

Daniel, Thanks again. To summarize, I want to implement the gradients and the laplacian functions for the SymmetricTensor class. The laplacian (trace of hessian) for the fe_values block extractor is straightforward, but the gradients (of the solution field) for the SymmetricTensor has some iss

Re: [deal.II] Re: Running deal.ii with METIS

Chenchen, 2016-10-04 18:15 GMT-04:00 Chenchen Liu : > Second, for "#PBS -l nodes=1:ppn=x", no matter what the value x is, the > error always exists. But this is not running across multiple nodes. You are only using one node, so this the same as running the code on your laptop. What happens if you

Re: [deal.II] Re: Issue with Trilinos installation (issue with blas?)

On 10/04/2016 07:12 PM, Matt Weller wrote: Okay I tried as you suggested, and same issue. It can't find -l/usr/lib64/libblas.so.3.6.1 . There is something off with the library but I have not the foggiest what it is. I don't know what to suggest other than to play with your package manager unt

Re: [deal.II] how to use solution.

Kyusik, As Wolfgang already mentioned, you want to use FEValues::get_function_values() to extract from a global FE Vector the local values at the quadrature points in the cell you are on. These quadrature points normally don't coincide with with the mapped unit support points of the FiniteElem

Re: [deal.II] how to use solution.

2016년 10월 5일 수요일 오전 1시 17분 31초 UTC+9, Wolfgang Bangerth 님의 말: > > On 10/04/2016 09:12 AM, hank...@gmail.com wrote: > > > > template > > void Step6::get_profile () > > { > > const QGauss quadrature_formula(3); > > FEValues fe_values (fe, quadrature_formula, > >

Re: [deal.II] Change Boost version

Hello, I was trying to set the boost directory because it failed while configuring. What I had to do to get it running on windows, was to set the environmental variable BOOST_DIR with four backslashes. That way I had no parse errors, because of the set command in cmake (see here