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

2016-10-16 Thread Jean-Paul Pelteret
For later reference, this conversation has now been moved to this thread . On Friday, October 14, 2016 at 6:01:22 PM UTC+2, Hamed Babaei wrote: > > Hi Daniel, > > I hope you are doing well. I've been struggling to run step-40 using > Tr

Re: [deal.II] Problem with PETScWrappers::SolverCG

2016-10-16 Thread Jean-Paul Pelteret
Dear Hamed, You are duplicating your posts in two threads. This is not productive, so I'm going to ask you to please focus this topic within your other thread . You have been given some useful advice there that you should follow up on.

Re: [deal.II] Re: How to use Trilinos instead of Petsc in step-40

2016-10-16 Thread Wolfgang Bangerth
On 10/16/2016 03:46 PM, Hamed Babaei wrote: Another point I need to mention is that I had the same problem when as my first try of parallelization, I paralleled a much simpler code, the step-25, based on the step-40 . I received the same Segmentation Violation error from Petsc. At that time, I re

[deal.II] Easy way to calculate second-deviatoric tensor

2016-10-16 Thread JAEKWANG KIM
Hi, all. I am really receiving a lot of help from this group whenever I visit here. I have a question today. I am developing FEM code for analyzing non-newtonian fluids using deal.ii library, which is amazing tool. In many cases, I have to calculate 'The second invariant of rate-of-str

[deal.II] Re: How to use Trilinos instead of Petsc in step-40

2016-10-16 Thread Hamed Babaei
Another point I need to mention is that I had the same problem when as my first try of parallelization, I paralleled a much simpler code, the step-25, based on the step-40 . I received the same Segmentation Violation error from Petsc. At that time, I replaced PreconditionerAMG with Precondition

Re: [deal.II] Problem with PETScWrappers::SolverCG

2016-10-16 Thread Hamed Babaei
Another point I need to mention is that I had the same problem when as my first try of parallelization, I paralleled a much simpler code, the step-25, based on the step-40 . I received the same Segmentation Violation error from Petsc. At that time, I replaced PreconditionerAMG with Precondition

Re: [deal.II] Problem with PETScWrappers::SolverCG

2016-10-16 Thread Hamed Babaei
Timo, Step-40 works on my machine perfectly. The code and its dependent parameter file are attached. I appreciate it if you could run the code on your machine and see what's the problem with it. Thanks On Sunday, October 16, 2016 at 7:55:35 AM UTC-5, Timo Heister wrote: > > > I have paralleli

[deal.II] Re: How to use Trilinos instead of Petsc in step-40

2016-10-16 Thread Hamed Babaei
Even if I implement Trilinos It terminates with the following Petsc error: "PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range" On Sunday, October 16, 2016 at 10:39:49 AM UTC-5, Daniel Arndt wrote: > > In case there is an error with Trilinos,

[deal.II] Re: How to use Trilinos instead of Petsc in step-40

2016-10-16 Thread Hamed Babaei
Daniel, Let me explain the problem more. I heve parallelized a code based on step-40. The non-parallel code works for all the initial values and boundary conditions But the parallel one fails at the first call to solverCG with this error: "[0]PETSC ERROR: Caught signal number 11 SEGV: Segment

Re: [deal.II] How to addapt MappingQEulerian class for Parallel codes using the PETSc wrapper classes

2016-10-16 Thread Hamed Babaei
Timo, Thank you for your prompt reply which resolved my problem. On Sunday, October 16, 2016 at 7:58:40 AM UTC-5, Timo Heister wrote: > > > "To enable the use of the MappingQ1Eulerian class also in the context of > > parallel codes using the PETSc wrapper classes, the type of the vector > can

[deal.II] Re: How to use Trilinos instead of Petsc in step-40

2016-10-16 Thread Daniel Arndt
In case there is an error with Trilinos, what is it? Am Sonntag, 16. Oktober 2016 17:27:26 UTC+2 schrieb Daniel Arndt: > > Hamed, > > It still works with Petsc and Trilinos is not implemented at all. I appreciate it if you could let me know what changes should be made in step-40 t

[deal.II] Re: How to use Trilinos instead of Petsc in step-40

2016-10-16 Thread Daniel Arndt
Hamed, >>> It still works with Petsc and Trilinos is not implemented at all. I >>> appreciate it if you could let me know what changes should be made in >>> step-40 to make it work with Trilinos. >>> >> What exactly, do you mean by that? So you have both PETSc and Trilinos installed and step-4

Re: [deal.II] How to addapt MappingQEulerian class for Parallel codes using the PETSc wrapper classes

2016-10-16 Thread Timo Heister
> "To enable the use of the MappingQ1Eulerian class also in the context of > parallel codes using the PETSc wrapper classes, the type of the vector can > be specified as template parameter EulerVectorType Not specifying this > template argument in applications using the PETSc vector classes leads t

Re: [deal.II] Problem with PETScWrappers::SolverCG

2016-10-16 Thread Timo Heister
> I have parallelized a code for a Thermoelastic problem based on step-40 . Does step-40 work on your system? > It seems from debugger backtracethat that error comes from the > PETScWrappers::SolverCG. Can you share the backtrace with us? -- Timo Heister http://www.math.clemson.edu/~heister/

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

2016-10-16 Thread Timo Heister
> Timo -- I think what confused the reporter was that the comment says > "uncomment the following #define", but there is just an #if and a > complicated one to boot. Maybe ought to update the comment -- I think I > recall the #if having become significantly more complicated recently. I don't under