Re: [deal.II] Question step5 in the lectures(about Coefficient...)

2016-07-05 Thread Praveen C
Hi In value_list function, you can get (x,y) like this double x = points[i][0]; double y = points[i][1]; Best praveen On Wed, Jul 6, 2016 at 12:13 PM, wrote: > Hi, > > I am trying to study step-5.cc. But I have some problem... > > I would like to make Coefficient 1/x (As you know, In step-5.c

[deal.II] Question step5 in the lectures(about Coefficient...)

2016-07-05 Thread hanks0227
Hi, I am trying to study step-5.cc. But I have some problem... I would like to make Coefficient 1/x (As you know, In step-5.cc, The Coefficient is originally 20(if r<0.5) or 1(if r>0.5)) of course not on the circle grid(circle-grid.inp), but on the other grid not including the point that x=0

Re: [deal.II] Dealing with conflicting constraints

2016-07-05 Thread Ce Qin
Dear Wolfgang, I think the first one does the correct thing whereas the second may not. I think so. But in my case, when following the instructions in the first one, the linear solver does not converge at all. The second one does. Since @Markus implemented the preject_boundary_values_curl_confor

Re: [deal.II] Re: introduce random initial condition

2016-07-05 Thread Bruno Turcksin
Anup, On 07/05/2016 08:17 PM, Anup Basak wrote: But my question is that is there any function in dealii through which I can generate random numbers within a range at the nodal points. No, you will need to do it yourself using rand (http://www.cplusplus.com/reference/cstdlib/rand/). Best, Br

Re: [deal.II] Re: introduce random initial condition

2016-07-05 Thread Anup Basak
Hello Bruno, Thank you very much for your reply. But my question is that is there any function in dealii through which I can generate random numbers within a range at the nodal points. Thanks and regards, Anup. On Tue, Jul 5, 2016 at 7:13 PM, Bruno Turcksin wrote: > Anup, > > I assume that you

[deal.II] Re: introduce random initial condition

2016-07-05 Thread Bruno Turcksin
Anup, I assume that your solution at time t_1 depends on the solution at time t_0. If you use Lagrangian FE, you just need to fill your vector solution at t_0 with random values between T_0 and T_1. Best, Bruno On Tuesday, July 5, 2016 at 6:02:32 PM UTC-4, Anup Basak wrote: > > Dear All, > >

[deal.II] Re: Restoring Gauss Point History during checkpoint/restart

2016-07-05 Thread Bruno Turcksin
Rajat, On Sunday, July 3, 2016 at 1:19:26 PM UTC-4, RAJAT ARORA wrote: > > The Quadrature point history is a std::vector<> of a struct PointHistory > just like in step 18 > . > > Is there any way to do this ? > PointHistory is just a Symmetr

[deal.II] introduce random initial condition

2016-07-05 Thread Anup Basak
Dear All, I am solving a heat conduction equation and want to impose random distribution of temperature, say between T_0 to T_1, as an initial condition over the entire domain. I shall be thankful if someone can tell me how to introduce in dealii. Thanks and regards, Anup. -- The deal.II p

Re: [deal.II] (Multigrid) preconditioner for the PML method

2016-07-05 Thread Wolfgang Bangerth
On 07/05/2016 12:28 PM, Artur Safin wrote: I have convergence to the correct solution for very coarse meshes (I have an analytic free-space solution that I can compare to), and I have also been able to solve the problem with a direct solver. The issue of non-convergence crops up when I use a fin

Re: [deal.II] (Multigrid) preconditioner for the PML method

2016-07-05 Thread Artur Safin
Wolfgang, I have convergence to the correct solution for very coarse meshes (I have an analytic free-space solution that I can compare to), and I have also been able to solve the problem with a direct solver. The issue of non-convergence crops up when I use a finer mesh, and as a result, have a

[deal.II] New: step-55 and step-56

2016-07-05 Thread Timo Heister
I am happy to announce that we have two new tutorial programs in deal.II: step-55: "Solving the Stokes problem in parallel" https://www.dealii.org/developer/doxygen/deal.II/step_55.html contributed by Timo Heister The tutorial shows how to solve a Stokes problem with a block diagonal precondition

Re: [deal.II] (Multigrid) preconditioner for the PML method

2016-07-05 Thread Wolfgang Bangerth
On 07/05/2016 10:46 AM, Artur Safin wrote: I am working on solving a Helmholtz-type problem, where I am using the PML method to simulate absorbing boundary conditions. Originally (before adding PML) I was using the complex GAMG preconditioner from the PETSc package, which was working fine. But n

Re: [deal.II] Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range

2016-07-05 Thread Wolfgang Bangerth
On 07/05/2016 10:08 AM, Ehsan Esfahani wrote: thanks for your response. the problem is in preconditioner. I find out that my code is running on my friend's machine but I get the error related to SEGV. I have changed preconditioner to Jacobi and it's running but with more iterations than AMG. Als

[deal.II] (Multigrid) preconditioner for the PML method

2016-07-05 Thread Artur Safin
Hello, I am working on solving a Helmholtz-type problem, where I am using the PML method to simulate absorbing boundary conditions. Originally (before adding PML) I was using the complex GAMG preconditioner from the PETSc package, which was working fine. But now I am having a difficulty where t

Re: [deal.II] Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range

2016-07-05 Thread Ehsan Esfahani
Daniel, thanks for your response. the problem is in preconditioner. I find out that my code is running on my friend's machine but I get the error related to SEGV. I have changed preconditioner to Jacobi and it's running but with more iterations than AMG. Also, it's not running on my computer an

Re: [deal.II] Dealing with conflicting constraints

2016-07-05 Thread Wolfgang Bangerth
On 07/04/2016 08:52 AM, Ce Qin wrote: Dear all, I am now solving a problem using adaptive finite element method. The boundary condition is Dirichlet condition. When the cells on boundary are locally refined, the hanging node constraints are conflicting with inhomogeneous constraints. In Const

[deal.II] Re: problems with installation ; Unknown CMake command "DEAL_II_ADD_LIBRARY"

2016-07-05 Thread Pablo Perez del Castillo
Yes, i tried without opencascade, and i doewnloaded tar, and then zip, and then from git pull. El sábado, 2 de julio de 2016, 19:21:57 (UTC+2), Pablo Perez del Castillo escribió: > > I can not install dealii, it was downloaded with Git > > pablo@pablo-HP-Spectre-13-Ultrabook:~/dealii/build$ cmak