[deal.II] Re: Steps 14 in 3D

2017-07-07 Thread Jean-Paul Pelteret
Dear Elena, You have asked a question that is beyond the scope of expectation for this forum . Right now you're essentially asking someone to do your preliminary work for you, which is unfair. Part of the way to learn deal.II and becom

Re: [deal.II] transitioning from parallel::shared::Triangulation to parallel::distributed::Triangulation

2017-07-07 Thread Jean-Paul Pelteret
Dear Alberto, > Given the code snippet above, I wonder whether > locally_accumulated_displacement is a vector that only has the locally > owned > dofs, or whether it is a vector that has at least the locally active > elements > (or maybe also the ghost/locally relevant elements). If you cal

Re: [deal.II] different results with different number of processors

2017-07-07 Thread Weixiong Zheng
Appreciated. I'll try this idea out. Sent from my iPhone > On Jul 7, 2017, at 16:47, Wolfgang Bangerth wrote: > >> On 07/06/2017 06:25 PM, Weixiong Zheng wrote: >> I suspect it would not be about the weak form directly. So, every time I >> give an set of initial cells, I need to make sure tota

Re: [deal.II] transitioning from parallel::shared::Triangulation to parallel::distributed::Triangulation

2017-07-07 Thread Wolfgang Bangerth
Alberto, I am bothering today about a code I have been writing in the shared triangulation framework and its transition to the parallel::distributed::Triangulation framework. The shared version works just fine, and so does the parallel if run with 1 single process. I understand therefore th

Re: [deal.II] different results with different number of processors

2017-07-07 Thread Wolfgang Bangerth
On 07/06/2017 06:25 PM, Weixiong Zheng wrote: I suspect it would not be about the weak form directly. So, every time I give an set of initial cells, I need to make sure total initial cell number N_cell_tot % n_process == 0. If this condition is fulfilled, solution will always be correct. Say,

Re: [deal.II] Re: CMake extra libraries?

2017-07-07 Thread Martin Kronbichler
This code triggers a safety mechanism that I introduced recently (April 26). We also see this on one of the test machines: https://cdash.kyomu.43-1.org/viewBuildError.php?buildid=8758 but I previously thought this is only due to a leftover in the configuration and I did not see it in real life b

Re: [deal.II] Re: CMake extra libraries?

2017-07-07 Thread Bruno Turcksin
2017-07-07 15:42 GMT-04:00 Victor Eijkhout : > It's weird. I just noticed taht I didn't have any -O2 flags in trilinos, so > I added. And now the deal install has compile lines with O2 and the AVX > flags, but also a bunch of compile lines with O0. And only sometimes those > blow up: Can you config

[deal.II] Re: define 2D/3D Point

2017-07-07 Thread Daniel Arndt
Lisa, you can write something like Point p; switch(dim) { case 1: p = Point (1); break; case 2: p = Point (1, 2); break; case 3: p = Point (1, 2, 3); break; default: Assert(false, ExcNotImplemented()); } const Point test_point = p; Best, Daniel -- The dea

Re: [deal.II] Re: CMake extra libraries?

2017-07-07 Thread Victor Eijkhout
On Friday, July 7, 2017 at 7:22:14 AM UTC-5, Bruno Turcksin wrote: > > This is very strange. Do you set the AVX flags yourself? > Yes. The CMakeOutput.log starts with: The system is: Linux - 3.10.0-514.21.2.el7.x86_64 - x86_64 Compiling the C compiler identification source file "CMakeCCompile

[deal.II] Steps 14 in 3D

2017-07-07 Thread Elena Greco
Dear Friends, I want to change step 14 to run a 3D model. I know that in deal ii many parts are dim independent and by changing dim to 3, automatically many parts work for 3D and do not need any modification. I want to know, for 3D model which parts should I change in step 14. In general I know

[deal.II] define 2D/3D Point

2017-07-07 Thread Lisa Collins
Hello, I would like to define a constant 2D/3D point, but I can not do it. Maybe my question is more a C++ question than a Deal ii, but I would appreciate it if someone help me. I want to write something like: switch (dim) { case 2: const Point test_point (value1, value2);

Re: [deal.II] Re: CMake extra libraries?

2017-07-07 Thread Bruno Turcksin
This is very strange. Do you set the AVX flags yourself? @Martin have you ever seen something like that? Bruno 2017-07-07 3:50 GMT-04:00 Victor Eijkhout : > Hm. That second compile line should be: > > > > > > > > > > [ 63%] Building CXX object > source/numerics/CMakeFiles/obj_numerics_debug.dir/

[deal.II] Re: Example for three coupled equations in dealII?

2017-07-07 Thread Bruno Turcksin
Hey, On Friday, July 7, 2017 at 5:03:05 AM UTC-4, Maxi Miller wrote: > > I have the following three PDEs (reduced from a more complex version): > > >

[deal.II] Example for three coupled equations in dealII?

2017-07-07 Thread 'Maxi Miller' via deal.II User Group
I have the following three PDEs (reduced from a more complex version): Now I am stuck. Usually I convert them to a discretized form (as for exam

[deal.II] Re: CMake extra libraries?

2017-07-07 Thread Victor Eijkhout
Hm. That second compile line should be: [ 63%] Building CXX object source/numerics/CMakeFiles/obj_numerics_debug.dir /vector_tools_boundary.cc.o cd /tmp/dealii-build/source/numerics && /opt/intel/ compilers_and_libraries_2017.4.196/linux/mpi/intel64/bin/mpicxx -DDEBUG - DTBB_DO_ASSERT=1

[deal.II] Re: CMake extra libraries?

2017-07-07 Thread Victor Eijkhout
Yay! Finally compiling. Mind if I just continue in this ticket even though the question is moving on? I think I'm correctly specifying a bunch of AVX flags to trilinos. In the deal build they show up in compile lines such as: [ 63%] Building CXX object source/physics/CMakeFiles/obj_physics_rel