Re: [deal.II] Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range
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 and the
problem on my computer remained the same as before :-(
Best,
Ehsan
On Monday, July 4, 2016 at 5:00:51 PM UTC-5, Daniel Arndt wrote:
>
> Ehsan,
>
> All I can say: After switching the order of arguments in
> SparseMatrix::add, your code runs for me with a recent developer version
> and Trilinos at least.
>
> Best,
> Daniel
>
> Am Montag, 4. Juli 2016 18:59:05 UTC+2 schrieb Ehsan Esfahani:
>>
>> Dear Professor Bangerth,
>>
>> Thanks for your response. Yes, I did. As I mentioned, I got a backtrace
>> in the debugger (eclipse) and I find out that the problem is in the line I
>> have mentioned but I couldn't find out what's the problem in that line of
>> the code which causes segmentation violation.
>>
>> Best,
>> Ehsan
>>
>>
>> On Sunday, July 3, 2016 at 4:32:16 PM UTC-5, bangerth wrote:
>>>
>>> On 07/03/2016 03:50 PM, Ehsan Esfahani wrote:
>>> > Dear All,
>>> >
>>> > Greetings. I changed step-25 (minor changes) in order to solve my
>>> problem. Now
>>> > I want to change this code for parallel computation based on the
>>> method
>>> > mentioned in step-40. I got several errors and solved them one by one
>>> but the
>>> > following error:
>>> >
>>> > /Number of active cells: 1024
>>> > // Total number of cells: 1365
>>> > //{[0,4224]}//
>>> > //Time step #1; advancing to t = 0.1.
>>> > [...]
>>> > //[0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation
>>> Violation,
>>> > probably memory access out of range
>>> > [...]
>>> >
>>> > /
>>> > /
>>> > Eclipse gives me a backtrace in the following line of the code:
>>> > /solver.solve (system_matrix,
>>> completely_distributed_solution_update,
>>> > system_rhs,/
>>> > / preconditioner);/
>>> > I have no idea why I got this error. The code is running properly for
>>> /fe(1)
>>> > /and /n_global_refinements (4)/ but when I change them to /fe(2)/ and
>>> > n_global_refinments (4) I got that error related to /Segmentation
>>> Violation.
>>> > /Do you know what's going on? Also, I have attached the code here .
>>> Thanks in
>>> > advance for your help.
>>>
>>> Ehsan,
>>>
>>> segmentation violations (SEGV) are typically easy to debug because you
>>> can get
>>> a backtrave in the debugger of the exact place where it happens, and you
>>> can
>>> then look at the local variables to see why this may have happened. Have
>>> you
>>> tried to run the program in a debugger and see what is going on?
>>>
>>> Best
>>> W.
>>>
>>> --
>>>
>>> Wolfgang Bangerth email:bang...@math.tamu.edu
>>> www:
>>> http://www.math.tamu.edu/~bangerth/
>>>
>>>
--
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
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
Re: [deal.II] Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range
Dear Professor Bangerth,
Thanks for your response. Yes, I did. As I mentioned, I got a backtrace in
the debugger (eclipse) and I find out that the problem is in the line I
have mentioned but I couldn't find out what's the problem in that line of
the code which causes segmentation violation.
Best,
Ehsan
On Sunday, July 3, 2016 at 4:32:16 PM UTC-5, bangerth wrote:
>
> On 07/03/2016 03:50 PM, Ehsan Esfahani wrote:
> > Dear All,
> >
> > Greetings. I changed step-25 (minor changes) in order to solve my
> problem. Now
> > I want to change this code for parallel computation based on the method
> > mentioned in step-40. I got several errors and solved them one by one
> but the
> > following error:
> >
> > /Number of active cells: 1024
> > // Total number of cells: 1365
> > //{[0,4224]}//
> > //Time step #1; advancing to t = 0.1.
> > [...]
> > //[0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation
> Violation,
> > probably memory access out of range
> > [...]
> >
> > /
> > /
> > Eclipse gives me a backtrace in the following line of the code:
> > /solver.solve (system_matrix,
> completely_distributed_solution_update,
> > system_rhs,/
> > / preconditioner);/
> > I have no idea why I got this error. The code is running properly for
> /fe(1)
> > /and /n_global_refinements (4)/ but when I change them to /fe(2)/ and
> > n_global_refinments (4) I got that error related to /Segmentation
> Violation.
> > /Do you know what's going on? Also, I have attached the code here .
> Thanks in
> > advance for your help.
>
> Ehsan,
>
> segmentation violations (SEGV) are typically easy to debug because you can
> get
> a backtrave in the debugger of the exact place where it happens, and you
> can
> then look at the local variables to see why this may have happened. Have
> you
> tried to run the program in a debugger and see what is going on?
>
> Best
> W.
>
> --
>
> Wolfgang Bangerth email:bang...@math.tamu.edu
>
> www: http://www.math.tamu.edu/~bangerth/
>
>
--
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
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
[deal.II] Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range
Dear All,
Greetings. I changed step-25 (minor changes) in order to solve my problem.
Now I want to change this code for parallel computation based on the method
mentioned in step-40. I got several errors and solved them one by one but
the following error:
>
> *Number of active cells: 1024*
> * Total number of cells: 1365**{[0,4224]}*
>
> *Time step #1; advancing to t = 0.1.*
>
> *--*
> *MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD **with
> errorcode 59.*
>
> *NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.*
> *You may or may not see output from other processes, depending on*
> *exactly when Open MPI kills them.*
>
> *--*
> *[0]PETSC ERROR:
> *
> *[0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation,
> probably memory access out of range*
> *[0]PETSC ERROR: Try option -start_in_debugger or
> -on_error_attach_debugger*
> *[0]PETSC ERROR: or
> see http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind
> <http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind>*
> *[0]PETSC ERROR: or try http://valgrind.org <http://valgrind.org/> on
> GNU/linux and Apple Mac OS X to find memory corruption errors*
> *[0]PETSC ERROR: configure using --with-debugging=yes, recompile, link,
> and run *
> *[0]PETSC ERROR: to get more information on the crash.*
> *[0]PETSC ERROR: - Error Message
> --*
> *[0]PETSC ERROR: Signal received*
> *[0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html
> <http://www.mcs.anl.gov/petsc/documentation/faq.html> for trouble shooting.*
> *[0]PETSC ERROR: Petsc Release Version 3.6.3, Dec, 03, 2015 *
> *[0]PETSC ERROR:
> /home/ehsan/apps/candi/deal.II-toolchain/deal.II-v8.4.0/examples/step-25
> (mine)/step-25 on a arch-linux2-c-opt named levitasgrad01.me.iastate.edu
> <http://levitasgrad01.me.iastate.edu/> by ehsan Sun Jul 3 13:46:52 2016*
> *[0]PETSC ERROR: Configure options
> --prefix=/home/ehsan/apps/candi/deal.II-toolchain/petsc-3.6.3
> --with-debugging=0 --with-shared-libraries=1 --with-mpi=1 --with-x=0
> --download-hypre=1 CC=mpicc CXX=mpicxx FC=mpif90**[0]PETSC ERROR: #1 User
> provided function() line 0 in unknown file*
Eclipse gives me a backtrace in the following line of the code:
* solver.solve (system_matrix, completely_distributed_solution_update,
system_rhs,*
* preconditioner); *
I have no idea why I got this error. The code is running properly for
*fe(1) *and *n_global_refinements (4)* but when I change them to *fe(2)* and
n_global_refinments (4) I got that error related to *Segmentation
Violation. *Do you know what's going on? Also, I have attached the code
here . Thanks in advance for your help.
Best Regards,
Ehsan Esfahani
--
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
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
/* -
*
* Copyright (C) 2006 - 2015 by the deal.II authors
*
* This file is part of the deal.II library.
*
* The deal.II library is free software; you can use it, redistribute
* it, and/or modify it under the terms of the GNU Lesser General
* Public License as published by the Free Software Foundation; either
* version 2.1 of the License, or (at your option) any later version.
* The full text of the license can be found in the file LICENSE at
* the top level of the deal.II distribution.
*
* -
*
* Author: Ivan Christov, Wolfgang Bangerth, Texas A University, 2006
*/
#include //
#include //
#include //
#include
#include
#include
#include
#include //this and the following are dependent
namespace LA//
{
#if defined(DEAL_II_WITH_PETSC) && !(defined(DEAL_II_WITH_TRILINOS) && defined(FORCE_USE_OF_TRILINOS))
using namespace dealii::LinearAlgebraPETSc;
# define USE_PETSC_LA
#elif defined(DEAL_II_WITH_TRILINOS)
using namespace dealii::LinearAlgebraTrilinos;
#else
# error DEAL_II_WITH_PETSC or DEAL_II_WITH_TRILINOS required
#endif
}
#include
#include
#include
#include //
#include
#include
#include
#include
#include //
#include //
#include //
#include //
#include
#in
[deal.II] Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range
Dear All,
Greetings. I changed step-25 (minor changes) in order to solve my problem.
Now I want to change this code for parallel computation based on the method
mentioned in step-40. I got several errors and solved them one by one but
the following error:
>
> *Number of active cells: 1024*
> * Total number of cells: 1365**{[0,4224]}*
>
> *Time step #1; advancing to t = 0.1.*
>
> *--*
> *MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD **with
> errorcode 59.*
>
> *NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.*
> *You may or may not see output from other processes, depending on*
> *exactly when Open MPI kills them.*
>
> *--*
> *[0]PETSC ERROR:
> *
> *[0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation,
> probably memory access out of range*
> *[0]PETSC ERROR: Try option -start_in_debugger or
> -on_error_attach_debugger*
> *[0]PETSC ERROR: or
> see http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind
> <http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind>*
> *[0]PETSC ERROR: or try http://valgrind.org <http://valgrind.org/> on
> GNU/linux and Apple Mac OS X to find memory corruption errors*
> *[0]PETSC ERROR: configure using --with-debugging=yes, recompile, link,
> and run *
> *[0]PETSC ERROR: to get more information on the crash.*
> *[0]PETSC ERROR: - Error Message
> --*
> *[0]PETSC ERROR: Signal received*
> *[0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html
> <http://www.mcs.anl.gov/petsc/documentation/faq.html> for trouble shooting.*
> *[0]PETSC ERROR: Petsc Release Version 3.6.3, Dec, 03, 2015 *
> *[0]PETSC ERROR:
> /home/ehsan/apps/candi/deal.II-toolchain/deal.II-v8.4.0/examples/step-25
> (mine)/step-25 on a arch-linux2-c-opt named levitasgrad01.me.iastate.edu
> <http://levitasgrad01.me.iastate.edu/> by ehsan Sun Jul 3 13:46:52 2016*
> *[0]PETSC ERROR: Configure options
> --prefix=/home/ehsan/apps/candi/deal.II-toolchain/petsc-3.6.3
> --with-debugging=0 --with-shared-libraries=1 --with-mpi=1 --with-x=0
> --download-hypre=1 CC=mpicc CXX=mpicxx FC=mpif90**[0]PETSC ERROR: #1 User
> provided function() line 0 in unknown file*
Eclipse gives me a backtrace in the following line of the code:
* solver.solve (system_matrix, completely_distributed_solution_update,
system_rhs,*
* preconditioner); *
I have no idea why I got this error. The code is running properly for
*fe(1) *and *n_global_refinements (4)* but when I change them to *fe(2)* and
n_global_refinments (4) I got that error related to *Segmentation
Violation. *Do you know what's going on? Also, I have attached the code
here . Thanks in advance for your help.
Best Regards,
Ehsan Esfahani
--
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
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
/* -
*
* Copyright (C) 2006 - 2015 by the deal.II authors
*
* This file is part of the deal.II library.
*
* The deal.II library is free software; you can use it, redistribute
* it, and/or modify it under the terms of the GNU Lesser General
* Public License as published by the Free Software Foundation; either
* version 2.1 of the License, or (at your option) any later version.
* The full text of the license can be found in the file LICENSE at
* the top level of the deal.II distribution.
*
* -
*
* Author: Ivan Christov, Wolfgang Bangerth, Texas A University, 2006
*/
#include //
#include //
#include //
#include
#include
#include
#include
#include //this and the following are dependent
namespace LA//
{
#if defined(DEAL_II_WITH_PETSC) && !(defined(DEAL_II_WITH_TRILINOS) && defined(FORCE_USE_OF_TRILINOS))
using namespace dealii::LinearAlgebraPETSc;
# define USE_PETSC_LA
#elif defined(DEAL_II_WITH_TRILINOS)
using namespace dealii::LinearAlgebraTrilinos;
#else
# error DEAL_II_WITH_PETSC or DEAL_II_WITH_TRILINOS required
#endif
}
#include
#include
#include
#include //
#include
#include
#include
#include
#include //
#include //
#include //
#include //
#include
#in
[deal.II] Re: changing step-25 based on step-40
I have run this code by Eclipse in debug mode. It's been terminated so I cannot track the error, The error printed on Console is the same as the one I have mentioned here. On Tuesday, June 7, 2016 at 2:47:34 PM UTC-5, Jean-Paul Pelteret wrote: > > Does this code work on one core (mpirun -np 1)? Are there any extra errors > when running in debug mode, or are you already running it in debug mode? > > On Tuesday, June 7, 2016 at 7:02:42 PM UTC+2, Ehsan Esfahani wrote: >> >> Thank for your response. I'm not sure it's related because, previously, >> without distributed triangulation, I modified step-25 in order to solve my >> problem (ginzburg landau eq.) and in that code, I didn't use those lines of >> the code, and it's running without errors. Also, I don't need to implement >> boundary conditions because of GL eq. >> >> On Tuesday, June 7, 2016 at 11:26:02 AM UTC-5, Jean-Paul Pelteret wrote: >>> >>> The one thing that I've noticed is that you're not using a constraint >>> matrix, as is done in step-40: >>> constraints.clear >>> <https://dealii.org/8.4.0/doxygen/deal.II/classConstraintMatrix.html#a24120d0331183f9a63cbe41493a19f6b> >>> >>> (); >>> constraints.reinit >>> <https://dealii.org/8.4.0/doxygen/deal.II/classConstraintMatrix.html#ac2726821354883ac97fe7e6181de9792> >>> >>> (locally_relevant_dofs); >>> DoFTools::make_hanging_node_constraints >>> <https://dealii.org/8.4.0/doxygen/deal.II/group__constraints.html#ga3eaa31a679484e80c193e74e8a967dc8> >>> >>> (dof_handler, constraints); >>> VectorTools::interpolate_boundary_values >>> <https://dealii.org/8.4.0/doxygen/deal.II/namespaceVectorTools.html#af6f700f193e9d5b52e9efe55e9b872d5> >>> >>> (dof_handler, >>> 0, >>> ZeroFunction >>> <https://dealii.org/8.4.0/doxygen/deal.II/classZeroFunction.html>(), >>> constraints); >>> constraints.close >>> <https://dealii.org/8.4.0/doxygen/deal.II/classConstraintMatrix.html#a8056d07faa2a7ed3f158c1b42d56abc8> >>> >>> (); >>> constraints.distribute_local_to_global >>> <https://dealii.org/8.4.0/doxygen/deal.II/classConstraintMatrix.html#aa9f3612a8fc51eafa34252bb436e8ae4> >>> >>> (cell_matrix, >>>cell_rhs, >>>local_dof_indices, >>>system_matrix, >>>system_rhs); >>> Perhaps this has something to do with it? I'm not sufficiently familiar >>> with distributed triangulations to say for sure that this is the problem. >>> It looks like you may be trying to access an entry in the sparsity pattern >>> that doesn't exist. >>> >>> >>> On Tuesday, June 7, 2016 at 6:16:40 PM UTC+2, Jean-Paul Pelteret wrote: >>>> >>>> For those who may be able to offer some advice, this is the full error >>>> message extracted from the discussion here >>>> <https://github.com/dealii/dealii/issues/2670>: >>>> >>>> **> Time step #1 <https://github.com/dealii/dealii/pull/1>; advancing >>>> to t = 0.2. >>>> >>>> >>>> >>>> [0]PETSC ERROR: - Error Message >>>> -- >>>> >>>> [0]PETSC ERROR: Argument out of range >>>> >>>> [0]PETSC ERROR: Inserting a new nonzero at global row/column (0, 0) >>>> into matrix >>>> >>>> [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html >>>> for trouble shooting. >>>> >>>> [0]PETSC ERROR: Petsc Release Version 3.6.3, Dec, 03, 2015 >>>> >>>> [0]PETSC ERROR: >>>> /home/ehsan/apps/candi/deal.II-toolchain/deal.II-v8.4.0/examples/step- >>>> 25 >>>> (mine)/step-25 on a arch-linux2-c-opt named >>>> levitasgrad01.me.iastate.edu by ehsan Mon Jun 6 13:06:02 2016 >>>> >>>> [0]PETSC ERROR: Configure options >>>> --prefix=/home/ehsan/apps/candi/deal.II-toolchain/petsc-3.6.3 >>>> --with-debugging=0 --with-shared-libraries=1 --with-mpi=1 --with-x=0 >>>> --download-hypre=1 CC=mpicc
[deal.II] Re: changing step-25 based on step-40
53s |48% | >> >> +-+---+++ >> >> >> ERROR: Uncaught exception in MPI_InitFinalize on proc 0. Skipping >> MPI_Finalize() to avoid a deadlock. >> >> >> -- >> >> >> Exception on processing: >> >> >> -- >> >> >> An error occurred in line <1424> of file in function >> >> void >> dealii::PETScWrappers::MatrixBase::add(dealii::PETScWrappers::MatrixBase >> ::size_type, >> dealii::PETScWrappers::MatrixBase::size_type, const size_type*, const >> PetscScalar*, bool, bool) >> >> The violated condition was: >> >> ierr == 0 >> >> The name and call sequence of the exception was: >> >> ExcPETScError(ierr) >> >> Additional Information: >> >> >> An error with error number 63 occurred while calling a PETSc function >> >> Aborting! >> >> ** >> >> >> >> >> On Monday, June 6, 2016 at 4:14:12 PM UTC+2, Ehsan Esfahani wrote: >>> >>> Dear Deal-ii Users, >>> >>> Greetings. I am trying to change step-25 for MPI based on the method >>> described in step-40. Unfortunately, I got stuck with an error :( >>> >>>> >>>> >>>> >>>> >>>> >>>> >>>> >>>> >>>> *An error occurred in line <1424> of file >>>> >>>> >>>> in functionvoid >>>> dealii::PETScWrappers::MatrixBase::add(dealii::PETScWrappers::MatrixBase::size_type, >>>> >>>> dealii::PETScWrappers::MatrixBase::size_type, const size_type*, const >>>> PetscScalar*, bool, >>>> bool)+-+---+++The >>>> >>>> violated condition was: ierr == 0The name and call sequence of the >>>> exception was:ExcPETScError(ierr)Additional Information: An error with >>>> error number 63 occurred while calling a PETSc function* >>> >>> *I've attached the code* if you like to see a full description of the >>> error. I cannot understand why I'm getting this error because I'm just >>> trying to change step-25 based on step-40. It was thinking it would be easy >>> but it seems I was mistaken. Do you have any suggestion about this error? I >>> find out that I am inserting a new nonzero component in a sparse matrix and >>> it violates the sparsitypattern but I cannot find out where I'm doing this >>> because the code is aborting and I cannot debug it correctly in >>> Eclipse-Mars II. >>> Also, is it really possible to change step-25 for MPI by using step-40 >>> or I'm completely on a wrong road? >>> Thanks for your help in advance. >>> >>> Best Regards, >>> Ehsan >>> >> -- 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 --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
[deal.II] "run configurations" problem in Eclipse (mars)
Dear Dealii users,
Greetings. As a person who is trying to learn Deal-ii, I'm trying to change
step-25 for parallel computation based on step-40, I have changed this step
for my problem (I solved my problem with modifying this step previously
without MPI). Unfortunately, I got several errors. The problem is: when I
want to use "run configurations" in eclipse for debugging the code and
understanding what's going wrong in my code, I cannot create a C++
application. I can create that for all other steps, but for mine, it
doesn't work and I cannot use debug in Eclipse. Also, attached I am sending
my code, and I wish there would be someone who faced the same problem.
Thanks in advance for your help.
Best,
--
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
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
/* -
*
* Copyright (C) 2006 - 2015 by the deal.II authors
*
* This file is part of the deal.II library.
*
* The deal.II library is free software; you can use it, redistribute
* it, and/or modify it under the terms of the GNU Lesser General
* Public License as published by the Free Software Foundation; either
* version 2.1 of the License, or (at your option) any later version.
* The full text of the license can be found in the file LICENSE at
* the top level of the deal.II distribution.
*
* -
*
* Author: Ivan Christov, Wolfgang Bangerth, Texas A University, 2006
*/
#include //
#include //
#include //
#include
#include
#include
#include
#include //this and the following are dependent
namespace LA//
{
#if defined(DEAL_II_WITH_PETSC) && !(defined(DEAL_II_WITH_TRILINOS) && defined(FORCE_USE_OF_TRILINOS))
using namespace dealii::LinearAlgebraPETSc;
# define USE_PETSC_LA
#elif defined(DEAL_II_WITH_TRILINOS)
using namespace dealii::LinearAlgebraTrilinos;
#else
# error DEAL_II_WITH_PETSC or DEAL_II_WITH_TRILINOS required
#endif
}
#include
#include
#include
#include //
#include
#include
#include
#include
#include //
#include //
#include //
#include //
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include //
#include
#include
#include
#include //
#include //
#include
#include
namespace Step25
{
using namespace dealii;
template
class SineGordonProblem
{
public:
SineGordonProblem ();
void run ();
private:
void make_grid_and_dofs ();
void assemble_system ();
void right_hand_side ();
void compute_nl_term (const LA::MPI::Vector _data,
const LA::MPI::Vector _data,
LA::MPI::Vector _term) const;//
void compute_nl_matrix (const LA::MPI::Vector _data,
const LA::MPI::Vector _data,
LA::MPI::SparseMatrix _matrix) const;//
unsigned int solve ();
void output_results (const unsigned int timestep_number) const;
MPI_Comm mpi_communicator;//
parallel::distributed::Triangulation triangulation;//
FE_Qfe;
DoFHandler dof_handler;
IndexSet locally_owned_dofs;
IndexSet locally_relevant_dofs;
//SparsityPattern sparsity_pattern;
LA::MPI::SparseMatrix system_matrix, laplace_matrix, mass_matrix;//system_matrix is the matrix that we want to invert it.
const unsigned int n_global_refinements;
double time;
const double final_time, time_step, theta, DeltaGbar;
//solution_update is equal to delta u in NR
LA::MPI::Vector locally_relevant_solution, locally_relevant_solution_update, locally_relevant_old_solution, system_rhs, M_x_velocity;//
const unsigned int output_timestep_skip;
ConditionalOStreampcout;//
TimerOutput computing_timer;//
};
template
class ExactSolution : public Function
{
public:
ExactSolution (const unsigned int n_components = 1,
const double time = 0.) : Function(n_components, time) {}
virtual double value (const Point ,
const unsigned int component = 0) const;
};
template
double ExactSolution::value (const Point ,
const unsigned int /*component*/) const
{
double t = this->get_time ();
switch (dim)
{
case 1:
{
const double m = 0.5;
const double c1 = 0.;
const double c2 = 0.;
return -4.*std::atan (m /
