[deal.II] Re: Error with the content "dealii::numbers::is_finite(sum)"

Dear All, I am running in one process. What should i do for solving this error? Does it mean that my vector has complex component? On Tuesday, October 25, 2016 at 2:29:19 PM UTC-5, Daniel Arndt wrote: > > Benhour, > > are you running in serial or with more than one process? In the first > case,

[deal.II] Re: Error in GridIn::read_abaqus for a mesh with two blocks

Awesome! I'm glad that worked for you. I'll update the wiki again and make a note that it works with newer versions of Cubit. If you're up for the challenge, there's the elusive bug in the read_abaqus function that still needs to be found ;-) J-P On Tuesday, October 25, 2016 at 9:16:44 PM UTC+

[deal.II] Re: boundary conditions

It seems to me that the specification of the flux out is not correct. In principle, instead of the Qout, one should specifiy a fixed head or temperature and the flux should adjust itself acoordingly. On Tuesday, October 25, 2016 at 9:14:08 PM UTC+2, Daniel Arndt wrote: > > OK. So do I understan

[deal.II] Re: Error with the content "dealii::numbers::is_finite(sum)"

Benhour, are you running in serial or with more than one process? In the first case, I would suspect that your matrix is singular. Best, Daniel -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- Y

[deal.II] Re: boundary conditions

OK. So do I understand you correctly, that you are imposing Neumann boundary conditions on all *four* boundaries? In that case the solution is only unique up to a constant. It seems to my that in the problem I posed, the boundary condition is not > correct. > Can you be a bit more specific? How

[deal.II] Re: Error with the content "dealii::numbers::is_finite(sum)"

Benhour, On Tuesday, October 25, 2016 at 12:59:29 PM UTC-4, benhour.amiria...@gmail.com wrote: > > I am solving a thermoelastic problem in cylindrical coordinate system. I > used CG solver for thermal problem and Newton-Raphson method for elasticity > section. For solving linear system, I used

[deal.II] Error with the content "dealii::numbers::is_finite(sum)"

Dear All, I am solving a thermoelastic problem in cylindrical coordinate system. I used CG solver for thermal problem and Newton-Raphson method for elasticity section. For solving linear system, I used direct solver SparseDirectUMFPACK with SSOR as a preconditioner. Unfortunately, I got some er

[deal.II] Re: Accessing solution in a vertex given by global index found with find_closest_vertex

Thanks, I will try GridTools::vertex_to_cell_map . I want to set the initial value for the next problem in every point of the mesh (10-100k dofs) and the new mesh point can be outside the

[deal.II] Re: Accessing solution in a vertex given by global index found with find_closest_vertex

Ah, I see. I agree that there's a lack of efficiency here. To complete the previous thought, there's the GridTools:: find_cells_adjacent_to_vertex function, which does the second half of

[deal.II] Re: Accessing solution in a vertex given by global index found with find_closest_vertex

Hi Jean-Paul, thanks for the reply. With cell->vertex_dof_index(v, d) I need to somehow find the cell, which the vertex is a part of. I guess I could loop over every cell to find this, but this doesn't seem very efficient. In principle, the solution vector contains the solutions at every vertex

[deal.II] Re: Accessing solution in a vertex given by global index found with find_closest_vertex

Hi Krei, This has been a popular topic as of late. Here's a list of posts that will help you achieve what you're looking to do. Regards, J-P On Tuesday, October 25, 2016 at 4:08:43 PM UTC+2, krei wrote: > > Hi, I

[deal.II] Accessing solution in a vertex given by global index found with find_closest_vertex

Hi, I have a vector valued solution (2 components), corresponding triangulation and the dofhandler: Triangulation triangulation; DoFHandler dof_handler; Vector solution; I want to find the nearest vertex to an arbitrary point Point &p and then find the values of both of the components of the s

[deal.II] Re: Direct Solvers in Parallel

Hamed, On Monday, October 24, 2016 at 9:17:39 PM UTC-4, Hamed Babaei wrote: > I have installed Dealii and all of its dependent libraries (Petsc, > Trilinos, P4est ...) via Candi (https://github.com/koecher/candi). I was > wondering which direct solver I should use which works the same as > Spa

[deal.II] Re: boundary conditions

In step 7, there are two types of boundary condition. one on T1 and the other on the T2. The T1 boundary is a fixed temperature. In my problem, both boundries are of the heat flux. the step 7 boundary condition is proper. It seems to my that in the problem I posed, the boundary condition is no

Re: [deal.II] Re: How to apply boundary values for a particular point on the boundary instead of the whole boundary surface

Hi You can construct an interpolation of your data and use it for boundary condition, see http://www.dealii.org/developer/doxygen/deal.II/classFunctions_1_1InterpolatedUniformGridData.html http://www.dealii.org/developer/doxygen/deal.II/classFunctions_1_1InterpolatedTensorProductGridData.html B

Re: [deal.II] Re: How to apply boundary values for a particular point on the boundary instead of the whole boundary surface

The problem with using VectorTools::interpolate_boundary_values is that it takes a function f(x,y,..), so if my input boundary condition is discrete, as one value for each original node, i would have to test the nodes positions in order to apply the correct boundary values inside this function. Als

[deal.II] Re: boundary conditions

Retired Replicant, H can be hydraulic head and K is the hydraulic conductivity. > Or H can be temperature with k being diffusion coefficient. > Have a look at step-7, if you just consider a diffusion problem. There Neumann boundary conditions are used on part of the boundary which is also mathe

[deal.II] Re: Direct Solvers in Parallel

Dear Hamed, I think you can only use the TrilinosWrappers::SolverDirect classes for TrilinosWrappers::MPI::Vector 's and of course TrilinosWrappers::Vector 's. Unfortunately, I've no experience with PETScWrappers::SparseDirectMUMPS

[deal.II] Re: Error in GridIn::read_abaqus for a mesh with two blocks

Hi Oded, Thanks for your comment. I've always found it strange that this script (thank you to whoever first developed it) left the output format as .inp. I think it must have been a common suffix for the format, and the fact that it clashes with the ABAQUS file format is just unfortunate. Anyw