[deal.II] Re: Periodic boundary conditions seem to be not applied

2016-06-08 Thread Daniel Arndt
Bastien, Unfortunately, you won't be able to achieve what you want directly with make_periodicity_constraints. The offset vector would help you with something like u(0,y)=u(1,y+lambda). For what you want to do you would first create the constraints for u(0,y)=u(1,y) and modify the corresponding

[deal.II] Re: test failing with message /usr/lib/openmpi/lib/libmpi_cxx.so: error: undefined reference to 'opal_list_item_t_class'

2016-06-08 Thread Bruno Turcksin
Giovanna, can you run a simple "hello world" program with MPI? Best, Bruno On Wednesday, June 8, 2016 at 4:48:24 PM UTC-4, Giovanna Bucci wrote: > > Hello, > > I have compiled deal II without problems (see attached detailed.log file > for settings) on a new machine with Ubuntu 16.04. However a

[deal.II] test failing with message /usr/lib/openmpi/lib/libmpi_cxx.so: error: undefined reference to 'opal_list_item_t_class'

2016-06-08 Thread Giovanna Bucci
Hello, I have compiled deal II without problems (see attached detailed.log file for settings) on a new machine with Ubuntu 16.04. However all the tests fail. It looks like there is some problem with linking libraries and it is related to openmpi (see quicktest.log file attached), but I can't fi

Re: [deal.II] Re: deal.ii and Trilinos

2016-06-08 Thread Wolfgang Bangerth
On 06/08/2016 09:26 AM, Marco Delchini wrote: Our objective is to solve for two-phase flows in multi-D using a fully implicit temporal scheme. In the short term, we would need the Trilinos solver/perconditioner capabilities and also are very interested in the Sacado package that does automatic d

[deal.II] Re: Periodic boundary conditions seem to be not applied

2016-06-08 Thread Bastien Lauras
Hi, Thanks a lot, it works well now! I don't know if I misread step-45, but I don't think it is written in it that boundary_ids should be setup before any refinement. And if you take something like step-44, boundry_ids are setup after refining the mesh. So maybe it would be wise to detail it in

[deal.II] Re: Output incorrect along parallel partition

2016-06-08 Thread Artur Safin
Jonathan, Sorry, you're right -- the distribute() function is not the culprit here. I remember having the exact same issue a while ago, but don't quite remember what caused it. Can you post your program here? Artur -- The deal.II project is located at http://www.dealii.org/ For mailing list/f

[deal.II] Defining line elements at the interface of planner domain

2016-06-08 Thread asim khan
Respected user, I am new to dealii and just started solving test problems using tutorials. My objective is to solve the mathematical model having integral over the area and integral over the internal boundary. this internal boundary can be considered as the interface between two material. My pl

[deal.II] Re: Using the solution nodal values at previous time steps

2016-06-08 Thread Mohammad Sabawi
Dear Daniel, I am sorry It is a typo error, actually, I have old_s = old_solution_values[q_point](1) which represents the solution at previous time step i.e. U(t_{n-1}) and old_s0 = old_solution_values0 [q_point](1) which represents the solution at time (t_{n-2}) i.e. U(t_{n-2}) and old_s1 =

[deal.II] Re: Thermoelastic Problem

2016-06-08 Thread Hamed Babaei
Thank you for your helpful information. I was wondering if there would be any difficulties for parallelization of the code when using block matrices and vectors. would you please let me know if you have paralleled your code. Thanks, Hamed On Wednesday, June 8, 2016 at 4:48:49 AM UTC-5, ha.ba..

Re: [deal.II] Re: deal.ii and Trilinos

2016-06-08 Thread Bruno Turcksin
Marco, 2016-06-08 10:26 GMT-04:00 Marco Delchini : > What would be the main difference between Epetra and Tpetra? What would be > the benefits of using Tpetra over Epetra? With Tpetra, you can use Kokkos and you need Kokkos to run Trilinos on GPUs. > Our objective is to solve for two-phase flows

Re: [deal.II] Re: deal.ii and Trilinos

2016-06-08 Thread Marco Delchini
Jean-Paul and Bruno, thanks for your answer, this is helpful. What would be the main difference between Epetra and Tpetra? What would be the benefits of using Tpetra over Epetra? Our objective is to solve for two-phase flows in multi-D using a fully implicit temporal scheme. In the short term, w

[deal.II] Re: Using the solution nodal values at previous time steps

2016-06-08 Thread Daniel Arndt
Mohammad, first of all const double old_s = old_solution_values[q_point](1); const double old_s = old_solution_values0[q_point](1); const double old_s = old_solution_values1[q_point](1); looks weird. Are you redefining old_s all the time? Apart from that it looks so far good to me. Can you

[deal.II] Re: Using the solution nodal values at previous time steps

2016-06-08 Thread Mohammad Sabawi
Dear all, Any advice will be appreciated. Best regards Mohammad Sabawi On Tuesday, June 7, 2016 at 10:29:45 AM UTC+1, Mohammad Sabawi wrote: > > Dear all, > > > > I am solving a 3x3 block system arising form DG time discretisation of the > semilinear parabolic equations in the form u_t - \De

[deal.II] Re: deal.ii and Trilinos

2016-06-08 Thread Bruno Turcksin
Marco, like Jean-Paul said the Trilinos wrappers are mainly for the linear algebra. This doesn't mean that you can't use a package that we don't have wrapper for. You can always get the underlying Trilinos object from the wrapper and use it when you need to call another package. When you choose

[deal.II] Re: Thermoelastic Problem

2016-06-08 Thread ha . badnava
I have solve a thermoelastic problem using block matrices and block vectors and a same dofhandler for the displacement and temperature fields. Fully coupled thermal-stress analysis is needed when the stress analysis is dependent on the temperature distribution and the temperature distributio