[deal.II] Re: cmake-gui mpi

Dear Peng, This part of the installation documentation suggests that you can explicitly indicate which MPI compilers to use with the following flags passed to CMake: cmake -DDEAL_II_WITH_MPI:BOOL=ON -DCMAKE_C_COMPILER="" > -DCMAKE_CXX_COMPIL

[deal.II] cmake-gui mpi

Dear professor: I use cmake-gui to configure deall.II. the package mpi has installed to /usr/local, but cmake gui can not found it. i tried to run cmake withou gui, it also give the same problems how to resolve this problem. -- MPI not found but DEAL_II_WITH_MPI is set to TRUE. Try again w

[deal.II] Re: Adding periodicity to a thin sample with one element in thickness direction

Ah, Daniel beat me to it. I took too much time putting together another demo code :-) Aligned with Daniel suggested, I'd bet you're running in debug mode. Here's the output of the attached code running in release mode for 1000x1000x1 repetitions: $ make release > [100%] Switch CMAKE_BUILD_TYP

[deal.II] Re: Adding periodicity to a thin sample with one element in thickness direction

Dear Hamed, However, I think it is not the producing of the final refined mesh that > hangs the code but the add_periodicity when repetition parameter is > relatively large(100 for instance). Since, it has been mentioned in the > documentation that refinement of the course mesh should be applie

[deal.II] Re: Adding periodicity to a thin sample with one element in thickness direction

Hello all, I am so sorry, Wolfgang, if I failed to express my thoughts in a clear and proper manner. Let me restate my reply to J-P's comment, Hope it is better this time ! > I've attached a slightly amended version of your code that I think will > highlight what the perceived problem is. Wh

Re: [deal.II] Re: Adding periodicity to a thin sample with one element in thickness direction

On 05/16/2017 12:33 PM, Hamed Babaei wrote: Thank you for your help. For small number of repetition argument, for example less than 50, the code reaches the end. However, I just checked 100 for repetition and it has been for about an hour that the code is still stuck at the "triangulation.add_pe

[deal.II] Re: Adding periodicity to a thin sample with one element in thickness direction

Dear Jean-Paul, Thank you for your help. For small number of repetition argument, for example less than 50, the code reaches the end. However, I just checked 100 for repetition and it has been for about an hour that the code is still stuck at the "triangulation.add_periodicity(periodicity_vect

[deal.II] Re: Adding periodicity to a thin sample with one element in thickness direction

Dear Hamed, I could not find a bug in the code you provided. It worked perfectly fine for me using deal.II 8.5 and with the number of MPI processes = {1,2}. I've attached a slightly amended version of your code that I think will highlight what the perceived problem is. When you use GridGenerat

[deal.II] Re: Adding periodicity to a thin sample with one element in thickness direction

Dear Jean-Paul and Daniel, I provided a minimal code which demonstrates the issue. There is no problem when I use GridGenerator::hyper_rectangle to create a course mesh and after adding periodicity, triangulation.refine_global is called. However, when I use subdivided_hyper_rectangle with repet

Re: [deal.II] Re: Step 13 & 14 -- Different BCs on different boundary indicators

Elena, I am not sure I understand exactly what you want to do. Can you write a simple pseudo code to show the problem? The way I understand what you want to do, you need to replace interpolate_boundary_values function from step 13 & 14 (which allows you to use only one BC) and instead used the one

[deal.II] Re: Step 13 & 14 -- Different BCs on different boundary indicators

Dear Bruno Thanks for your reply. I know how to use interpolate_boundary_values but in steps 13 & 14, if I want to keep Solver::assemble_linear_system in the general form, I don't know how to implement new boundary conditions. In steps 13 & 14, Solver::assemble_linear_system assembles the line

[deal.II] Re: Step 13 & 14 -- Different BCs on different boundary indicators

Elena, On Tuesday, May 16, 2017 at 7:19:08 AM UTC-4, Elena Greco wrote: > > May be it would be better to edit my question and remove the "cylinder" > part, but I do not have access to modify the question. > Please consider my question for a general case of a geometry with multiple > boundary ind

[deal.II] Re: Getting error during cmake configuration step and later during make run

Daniel, Thanks for your reply. I am attaching 2 files. I do not see detailed.log file anywhere in the build directory. Vikram On Monday, May 15, 2017 at 6:25:38 AM UTC-5, Daniel Arndt wrote: > > Vikram, > >> I am new to dealii. I installed deal.II 8.5.0 using the dmg package from >> dealii web

Re: [deal.II] Error in the second derivative of solution...?

Hi Kyusik, What makes me confused here is that, according to your explanation, two > different DoFHandler could be used simultaneously. > Yes, that it correct. You can convert between cell iterators. So I do the following: typename Triangulation::active_cell_iterator cell =

Re: [deal.II] Error in the second derivative of solution...?

Dear Jean-Paul, Thank you again for a detailed explanation. However, please excuse my ignorance... > > 1. Create a vector-valued finite element system (maybe an > FESystem(FE_Q(1),dim) ) and use it to distribute DoFs onto a new > DoFHandler. > 2. Create a mass matrix "M" with the new DoFHandl

[deal.II] Re: Step 13 & 14 -- Different BCs on different boundary indicators

May be it would be better to edit my question and remove the "cylinder" part, but I do not have access to modify the question. Please consider my question for a general case of a geometry with multiple boundary indicators and different BCs on them. Thanks. Elena On Tuesday, May 16, 2017 at 1:10

[deal.II] Re: Step 13 & 14 -- Different BCs on different boundary indicators

Dear Uwe Thanks for you prompt reply. But if I import a mesh from gmesh and for a not regular geometry I would face a problem. Do you have any suggestion for this case? Thanks. Elena On Tuesday, May 16, 2017 at 12:50:36 PM UTC+2, Elena Greco wrote: > > Hello, > > In Steps 13 and 14, for primal

[deal.II] Re: Step 13 & 14 -- Different BCs on different boundary indicators

Elena, in your case you do not really need a second boundary indicator. You have the position argument "Point x" in the value function and from that you can compute the radius and decide if the first or second boundary condition holds and return the respective value. Then you have not to change a

[deal.II] Step 13 & 14 -- Different BCs on different boundary indicators

Hello, In Steps 13 and 14, for primal problem, the boundary condition is calculate based on an analytical solution, therefore it is easy to set different boundary indicators and apply the BC on different boundary indicators. In my problem I have a cylinder (boundary indicator of R_in is 0 and