Re: [deal.II] ZeroBoundaryCondition when FESystem is used

2017-03-01 Thread hanks0227
Thanks for your reply again. > Yes, of course. Via VectorTools::interpolate_boundary_values and > VectorTools::project_boundary_values. > > Due to your help, I can apply ZeroBoundaryCondition for pressure by adding //test for strong boundary condtion std::map boundary_values_strong;

[deal.II] Re: Announcing the deal.II Code Gallery

2017-03-01 Thread Michael Harmon
doh. thanks! On Wednesday, March 1, 2017 at 12:34:52 AM UTC-5, Jean-Paul Pelteret wrote: > > Wolfgang's answer to Michael's question from a duplicate post: > > >> Just check out the code gallery git repository in a directory parallel to >> the >> examples/ directory and build the deal.II docum

Re: [deal.II] ZeroBoundaryCondition when FESystem is used

2017-03-01 Thread Wolfgang Bangerth
You mention that The difference that in step-3, you are imposing *strong* boundary conditions, whereas in step-20, you are imposing *weak* boundary conditions. Then, is there anyway to use strong boundary condition using FESystem Yes, of course. Via VectorTools::interpolate_bounda

Re: [deal.II] ZeroBoundaryCondition when FESystem is used

2017-03-01 Thread hanks0227
Thank you very much your kind reply. Could I have one more question..? You mention that > > The difference that in step-3, you are imposing *strong* boundary > conditions, > whereas in step-20, you are imposing *weak* boundary conditions. > > Then, is there anyway to use strong boundary cond

Re: [deal.II] Configuring deal.ii error - can't pinpoint what the issue is

2017-03-01 Thread Phil H
Great, thanks! I have started Candi now - thanks to yourself and Uwe. On Wednesday, March 1, 2017 at 1:11:30 PM UTC, Wolfgang Bangerth wrote: > > On 03/01/2017 06:09 AM, Phil H wrote: > > > > Hmmm. So Trilinos may be the issue. Would that explain why LAPACK is set > to > > OFF as well? Do I h

Re: [deal.II] Configuring deal.ii error - can't pinpoint what the issue is

2017-03-01 Thread Wolfgang Bangerth
On 03/01/2017 06:09 AM, Phil H wrote: Hmmm. So Trilinos may be the issue. Would that explain why LAPACK is set to OFF as well? Do I have to specifically ask for LAPACK to be turned on? I think it is an issue with the installation of Trilinos -- some files that are supposed to be there are not

Re: [deal.II] Configuring deal.ii error - can't pinpoint what the issue is

2017-03-01 Thread Phil H
W, I don't end up with a Makefile file in dealii build. Just - -rw-r--r-- 1 klmz32 dgl 71 Mar 1 12:34 revision.log -rw-r--r-- 1 klmz32 dgl 1985 Mar 1 12:34 summary.log -rw-r--r-- 1 klmz32 dgl 16936 Mar 1 12:34 detailed.log drwxr-xr-x 3 klmz32 dgl 4096 Mar 1 12:34 share drwxr-xr-x

Re: [deal.II] ZeroBoundaryCondition when FESystem is used

2017-03-01 Thread Wolfgang Bangerth
Kysusik, I'm studying step-20 tutorial example. But, I can't understand it, because it seems to have different things compared to step-3 in terms of boundary condition. In the introduction of step-3, it is said that ∫Ω∇φ⋅∇u−∫∂Ωφn⋅∇u=∫Ωφf. The test function φ has to satisfy the same kind of

Re: [deal.II] Configuring deal.ii error - can't pinpoint what the issue is

2017-03-01 Thread Wolfgang Bangerth
On 03/01/2017 05:38 AM, Phil H wrote: Here is the full compilation, with the cmake. As I say, I can't see immediately where the issue is. Phil -- I suspect it has to do with the errors during Trilinos configuration, where it tries to find some load file but can't find it. That said, do you

Re: [deal.II] Configuring deal.ii error - can't pinpoint what the issue is

2017-03-01 Thread Phil H
Hi Wolfgang, Here is the full compilation, with the cmake. As I say, I can't see immediately where the issue is. CC=mpicc CXX=mpicxx CFLAGS="-fPIC -O2" cmake -DCMAKE_INSTALL_PREFIX=/local/earthsci/klmz32/Documents/LAB/dealii -DDEAL_II_WITH_TRILINOS=ON -DTRILINOS_DIR=/local/earthsci/klmz32

[deal.II] Re: ZeroBoundaryCondition when FESystem is used

2017-03-01 Thread hanks0227
I'm sorry, I should have include "part for boundary of pressure" What I mean by "part for boundary pressure in assemble_system" is (in step-20) for (unsigned int face_no=0; face_no::faces_per_cell; ++face_no) if (cell->at_boundary(face_no))

[deal.II] ZeroBoundaryCondition when FESystem is used

2017-03-01 Thread hanks0227
Hi all, I'm studying step-20 tutorial example. But, I can't understand it, because it seems to have different things compared to step-3 in terms of boundary condition. In the introduction of step-3, it is said that ∫Ω∇φ⋅∇u−∫∂Ωφn⋅∇u=∫Ωφf. The test function φ has to satisfy the same kind of bou

Re: [deal.II] Extract some blocks from a sparse block matrix

2017-03-01 Thread Daniel Jodlbauer
SparseDirectUMFPACK creates a copy of the matrix using iterators, which are not implemented for the BlockMatrixArray, so these would have to be added (and some minor other functions used within the factorize method). Alternatively, one could implement the factorize method for BlockMatrixArray se

[deal.II] Re: Configuring deal.ii error - can't pinpoint what the issue is

2017-03-01 Thread 'Uwe Köcher' via deal . II User Group
Dear Phil, to compile and use aspect, you need more dependencies (like Trilinos, p4est, MPI) for the deal.II library. To simplify your life, you can use the automated build script as announced on https://www.dealii.org/download.html or the direct link: https://github.com/dealii/candi It will h