Re: [deal.II] Issue encountered while solving Step-40 in 1 dimension

2022-08-23 Thread syed ansari
Thanks, Daniel. This is very helpful.

Best Regards,
Syed Ansari S.

On Mon, Aug 22, 2022 at 7:28 PM Daniel Arndt  wrote:

> Syed,
>
> Yes, you should be able to use parallel::shared::Triangulation instead.
>
> Best,
> Daniel
>
> On Sat, Aug 20, 2022 at 5:25 AM syed ansari  wrote:
>
>> Thanks Daniel for your quick reply. Is it possible to solve the
>> same problem with parallel::shared::Triangulation for dim ==1?
>>
>> On Fri, 19 Aug 2022, 8:18 pm Daniel Arndt, 
>> wrote:
>>
>>> Syed,
>>>
>>> parallel::distributed::Triangulation is just not implemented for dim==1
>>> so you can't run step-40 for the one-dimensional case.
>>>
>>> Best,
>>> Daniel
>>>
>>> On Fri, Aug 19, 2022 at 7:07 AM syed ansari  wrote:
>>>
 Dear all,
 I was trying to run step-40 in 1 dimension and encountered
 the error corresponding to MeshSmoothing in the constructor. The details of
 the error are as follows:
 
 An error occurred in line <3455> of file
 
 in function
 dealii::parallel::distributed::Triangulation<1,
 spacedim>::Triangulation(ompi_communicator_t* const&, typename
 dealii::Triangulation<1, spacedim>::MeshSmoothing,
 dealii::parallel::distributed::Triangulation<1, spacedim>::Settings) [with
 int spacedim = 1; MPI_Comm = ompi_communicator_t*; typename
 dealii::Triangulation<1, spacedim>::MeshSmoothing =
 dealii::Triangulation<1, 1>::MeshSmoothing]
 The violated condition was:
 false
 Additional information:
 You are trying to use functionality in deal.II that is currently not
 implemented. In many cases, this indicates that there simply didn't
 appear much of a need for it, or that the author of the original
 code
 did not have the time to implement a particular case. If you hit
 this
 exception, it is therefore worth the time to look into the code to
 find out whether you may be able to implement the missing
 functionality. If you do, please consider providing a patch to the
 deal.II development sources (see the deal.II website on how to
 contribute).

 Stacktrace:
 ---
 #0  /home/syed/dealii-candi/deal.II-v9.3.2/lib/libdeal_II.g.so.9.3.2:
 dealii::parallel::distributed::Triangulation<1,
 1>::Triangulation(ompi_communicator_t* const&, dealii::Triangulation<1,
 1>::MeshSmoothing, dealii::parallel::distributed::Triangulation<1,
 1>::Settings)
 #1  ./step-40: Step40::LaplaceProblem<1>::LaplaceProblem()
 #2  ./step-40: main
 

 Calling MPI_Abort now.
 To break execution in a GDB session, execute 'break MPI_Abort' before
 running. You can also put the following into your ~/.gdbinit:
   set breakpoint pending on
   break MPI_Abort
   set breakpoint pending auto

 
 An error occurred in line <3455> of file
 
 in function
 dealii::parallel::distributed::Triangulation<1,
 spacedim>::Triangulation(ompi_communicator_t* const&, typename
 dealii::Triangulation<1, spacedim>::MeshSmoothing,
 dealii::parallel::distributed::Triangulation<1, spacedim>::Settings) [with
 int spacedim = 1; MPI_Comm = ompi_communicator_t*; typename
 dealii::Triangulation<1, spacedim>::MeshSmoothing =
 dealii::Triangulation<1, 1>::MeshSmoothing]
 The violated condition was:
 false
 Additional information:
 You are trying to use functionality in deal.II that is currently not
 implemented. In many cases, this indicates that there simply didn't
 appear much of a need for it, or that the author of the original
 code
 did not have the time to implement a particular case. If you hit
 this
 exception, it is therefore worth the time to look into the code to
 find out whether you may be able to implement the missing
 functionality. If you do, please consider providing a patch to the
 deal.II development sources (see the deal.II website on how to
 contribute).

 Stacktrace:
 ---
 #0  /home/syed/dealii-candi/deal.II-v9.3.2/lib/libdeal_II.g.so.9.3.2:
 dealii::parallel::distributed::Triangulation<1,
 1>::Triangulation(ompi_communicator_t* const&, dealii::Triangulation<1,
 1>::MeshSmoothing, dealii::parallel::distributed::Triangulation<1,
 1>::Settings)
 #1  ./step-40: Step40::LaplaceProblem<1>::LaplaceProblem()
 #2  ./step-40: main
 

 Calling MPI_Abort now.
 To break execution in a GDB session, execute 'break MPI_Abort' before
 running. You can also put the following into your ~/.gdbinit:
   set breakpoint pending on
   break MPI_Abort
   set breakpoint pending auto

 

Re: [deal.II] Issue encountered while solving Step-40 in 1 dimension

2022-08-22 Thread Daniel Arndt
Syed,

Yes, you should be able to use parallel::shared::Triangulation instead.

Best,
Daniel

On Sat, Aug 20, 2022 at 5:25 AM syed ansari  wrote:

> Thanks Daniel for your quick reply. Is it possible to solve the
> same problem with parallel::shared::Triangulation for dim ==1?
>
> On Fri, 19 Aug 2022, 8:18 pm Daniel Arndt,  wrote:
>
>> Syed,
>>
>> parallel::distributed::Triangulation is just not implemented for dim==1
>> so you can't run step-40 for the one-dimensional case.
>>
>> Best,
>> Daniel
>>
>> On Fri, Aug 19, 2022 at 7:07 AM syed ansari  wrote:
>>
>>> Dear all,
>>> I was trying to run step-40 in 1 dimension and encountered
>>> the error corresponding to MeshSmoothing in the constructor. The details of
>>> the error are as follows:
>>> 
>>> An error occurred in line <3455> of file
>>> 
>>> in function
>>> dealii::parallel::distributed::Triangulation<1,
>>> spacedim>::Triangulation(ompi_communicator_t* const&, typename
>>> dealii::Triangulation<1, spacedim>::MeshSmoothing,
>>> dealii::parallel::distributed::Triangulation<1, spacedim>::Settings) [with
>>> int spacedim = 1; MPI_Comm = ompi_communicator_t*; typename
>>> dealii::Triangulation<1, spacedim>::MeshSmoothing =
>>> dealii::Triangulation<1, 1>::MeshSmoothing]
>>> The violated condition was:
>>> false
>>> Additional information:
>>> You are trying to use functionality in deal.II that is currently not
>>> implemented. In many cases, this indicates that there simply didn't
>>> appear much of a need for it, or that the author of the original code
>>> did not have the time to implement a particular case. If you hit this
>>> exception, it is therefore worth the time to look into the code to
>>> find out whether you may be able to implement the missing
>>> functionality. If you do, please consider providing a patch to the
>>> deal.II development sources (see the deal.II website on how to
>>> contribute).
>>>
>>> Stacktrace:
>>> ---
>>> #0  /home/syed/dealii-candi/deal.II-v9.3.2/lib/libdeal_II.g.so.9.3.2:
>>> dealii::parallel::distributed::Triangulation<1,
>>> 1>::Triangulation(ompi_communicator_t* const&, dealii::Triangulation<1,
>>> 1>::MeshSmoothing, dealii::parallel::distributed::Triangulation<1,
>>> 1>::Settings)
>>> #1  ./step-40: Step40::LaplaceProblem<1>::LaplaceProblem()
>>> #2  ./step-40: main
>>> 
>>>
>>> Calling MPI_Abort now.
>>> To break execution in a GDB session, execute 'break MPI_Abort' before
>>> running. You can also put the following into your ~/.gdbinit:
>>>   set breakpoint pending on
>>>   break MPI_Abort
>>>   set breakpoint pending auto
>>>
>>> 
>>> An error occurred in line <3455> of file
>>> 
>>> in function
>>> dealii::parallel::distributed::Triangulation<1,
>>> spacedim>::Triangulation(ompi_communicator_t* const&, typename
>>> dealii::Triangulation<1, spacedim>::MeshSmoothing,
>>> dealii::parallel::distributed::Triangulation<1, spacedim>::Settings) [with
>>> int spacedim = 1; MPI_Comm = ompi_communicator_t*; typename
>>> dealii::Triangulation<1, spacedim>::MeshSmoothing =
>>> dealii::Triangulation<1, 1>::MeshSmoothing]
>>> The violated condition was:
>>> false
>>> Additional information:
>>> You are trying to use functionality in deal.II that is currently not
>>> implemented. In many cases, this indicates that there simply didn't
>>> appear much of a need for it, or that the author of the original code
>>> did not have the time to implement a particular case. If you hit this
>>> exception, it is therefore worth the time to look into the code to
>>> find out whether you may be able to implement the missing
>>> functionality. If you do, please consider providing a patch to the
>>> deal.II development sources (see the deal.II website on how to
>>> contribute).
>>>
>>> Stacktrace:
>>> ---
>>> #0  /home/syed/dealii-candi/deal.II-v9.3.2/lib/libdeal_II.g.so.9.3.2:
>>> dealii::parallel::distributed::Triangulation<1,
>>> 1>::Triangulation(ompi_communicator_t* const&, dealii::Triangulation<1,
>>> 1>::MeshSmoothing, dealii::parallel::distributed::Triangulation<1,
>>> 1>::Settings)
>>> #1  ./step-40: Step40::LaplaceProblem<1>::LaplaceProblem()
>>> #2  ./step-40: main
>>> 
>>>
>>> Calling MPI_Abort now.
>>> To break execution in a GDB session, execute 'break MPI_Abort' before
>>> running. You can also put the following into your ~/.gdbinit:
>>>   set breakpoint pending on
>>>   break MPI_Abort
>>>   set breakpoint pending auto
>>>
>>> 
>>> An error occurred in line <3455> of file
>>> 
>>> in function
>>> dealii::parallel::distributed::Triangulation<1,
>>> spacedim>::Triangulation(ompi_communicator_t* const&, typename
>>> 

Re: [deal.II] Issue encountered while solving Step-40 in 1 dimension

2022-08-20 Thread syed ansari
Thanks Daniel for your quick reply. Is it possible to solve the
same problem with parallel::shared::Triangulation for dim ==1?

On Fri, 19 Aug 2022, 8:18 pm Daniel Arndt,  wrote:

> Syed,
>
> parallel::distributed::Triangulation is just not implemented for dim==1 so
> you can't run step-40 for the one-dimensional case.
>
> Best,
> Daniel
>
> On Fri, Aug 19, 2022 at 7:07 AM syed ansari  wrote:
>
>> Dear all,
>> I was trying to run step-40 in 1 dimension and encountered
>> the error corresponding to MeshSmoothing in the constructor. The details of
>> the error are as follows:
>> 
>> An error occurred in line <3455> of file
>> 
>> in function
>> dealii::parallel::distributed::Triangulation<1,
>> spacedim>::Triangulation(ompi_communicator_t* const&, typename
>> dealii::Triangulation<1, spacedim>::MeshSmoothing,
>> dealii::parallel::distributed::Triangulation<1, spacedim>::Settings) [with
>> int spacedim = 1; MPI_Comm = ompi_communicator_t*; typename
>> dealii::Triangulation<1, spacedim>::MeshSmoothing =
>> dealii::Triangulation<1, 1>::MeshSmoothing]
>> The violated condition was:
>> false
>> Additional information:
>> You are trying to use functionality in deal.II that is currently not
>> implemented. In many cases, this indicates that there simply didn't
>> appear much of a need for it, or that the author of the original code
>> did not have the time to implement a particular case. If you hit this
>> exception, it is therefore worth the time to look into the code to
>> find out whether you may be able to implement the missing
>> functionality. If you do, please consider providing a patch to the
>> deal.II development sources (see the deal.II website on how to
>> contribute).
>>
>> Stacktrace:
>> ---
>> #0  /home/syed/dealii-candi/deal.II-v9.3.2/lib/libdeal_II.g.so.9.3.2:
>> dealii::parallel::distributed::Triangulation<1,
>> 1>::Triangulation(ompi_communicator_t* const&, dealii::Triangulation<1,
>> 1>::MeshSmoothing, dealii::parallel::distributed::Triangulation<1,
>> 1>::Settings)
>> #1  ./step-40: Step40::LaplaceProblem<1>::LaplaceProblem()
>> #2  ./step-40: main
>> 
>>
>> Calling MPI_Abort now.
>> To break execution in a GDB session, execute 'break MPI_Abort' before
>> running. You can also put the following into your ~/.gdbinit:
>>   set breakpoint pending on
>>   break MPI_Abort
>>   set breakpoint pending auto
>>
>> 
>> An error occurred in line <3455> of file
>> 
>> in function
>> dealii::parallel::distributed::Triangulation<1,
>> spacedim>::Triangulation(ompi_communicator_t* const&, typename
>> dealii::Triangulation<1, spacedim>::MeshSmoothing,
>> dealii::parallel::distributed::Triangulation<1, spacedim>::Settings) [with
>> int spacedim = 1; MPI_Comm = ompi_communicator_t*; typename
>> dealii::Triangulation<1, spacedim>::MeshSmoothing =
>> dealii::Triangulation<1, 1>::MeshSmoothing]
>> The violated condition was:
>> false
>> Additional information:
>> You are trying to use functionality in deal.II that is currently not
>> implemented. In many cases, this indicates that there simply didn't
>> appear much of a need for it, or that the author of the original code
>> did not have the time to implement a particular case. If you hit this
>> exception, it is therefore worth the time to look into the code to
>> find out whether you may be able to implement the missing
>> functionality. If you do, please consider providing a patch to the
>> deal.II development sources (see the deal.II website on how to
>> contribute).
>>
>> Stacktrace:
>> ---
>> #0  /home/syed/dealii-candi/deal.II-v9.3.2/lib/libdeal_II.g.so.9.3.2:
>> dealii::parallel::distributed::Triangulation<1,
>> 1>::Triangulation(ompi_communicator_t* const&, dealii::Triangulation<1,
>> 1>::MeshSmoothing, dealii::parallel::distributed::Triangulation<1,
>> 1>::Settings)
>> #1  ./step-40: Step40::LaplaceProblem<1>::LaplaceProblem()
>> #2  ./step-40: main
>> 
>>
>> Calling MPI_Abort now.
>> To break execution in a GDB session, execute 'break MPI_Abort' before
>> running. You can also put the following into your ~/.gdbinit:
>>   set breakpoint pending on
>>   break MPI_Abort
>>   set breakpoint pending auto
>>
>> 
>> An error occurred in line <3455> of file
>> 
>> in function
>> dealii::parallel::distributed::Triangulation<1,
>> spacedim>::Triangulation(ompi_communicator_t* const&, typename
>> dealii::Triangulation<1, spacedim>::MeshSmoothing,
>> dealii::parallel::distributed::Triangulation<1, spacedim>::Settings) [with
>> int spacedim = 1; MPI_Comm = ompi_communicator_t*; typename
>> dealii::Triangulation<1, spacedim>::MeshSmoothing =
>> dealii::Triangulation<1, 

Re: [deal.II] Issue encountered while solving Step-40 in 1 dimension

2022-08-19 Thread Daniel Arndt
Syed,

parallel::distributed::Triangulation is just not implemented for dim==1 so
you can't run step-40 for the one-dimensional case.

Best,
Daniel

On Fri, Aug 19, 2022 at 7:07 AM syed ansari  wrote:

> Dear all,
> I was trying to run step-40 in 1 dimension and encountered the
> error corresponding to MeshSmoothing in the constructor. The details of the
> error are as follows:
> 
> An error occurred in line <3455> of file
> 
> in function
> dealii::parallel::distributed::Triangulation<1,
> spacedim>::Triangulation(ompi_communicator_t* const&, typename
> dealii::Triangulation<1, spacedim>::MeshSmoothing,
> dealii::parallel::distributed::Triangulation<1, spacedim>::Settings) [with
> int spacedim = 1; MPI_Comm = ompi_communicator_t*; typename
> dealii::Triangulation<1, spacedim>::MeshSmoothing =
> dealii::Triangulation<1, 1>::MeshSmoothing]
> The violated condition was:
> false
> Additional information:
> You are trying to use functionality in deal.II that is currently not
> implemented. In many cases, this indicates that there simply didn't
> appear much of a need for it, or that the author of the original code
> did not have the time to implement a particular case. If you hit this
> exception, it is therefore worth the time to look into the code to
> find out whether you may be able to implement the missing
> functionality. If you do, please consider providing a patch to the
> deal.II development sources (see the deal.II website on how to
> contribute).
>
> Stacktrace:
> ---
> #0  /home/syed/dealii-candi/deal.II-v9.3.2/lib/libdeal_II.g.so.9.3.2:
> dealii::parallel::distributed::Triangulation<1,
> 1>::Triangulation(ompi_communicator_t* const&, dealii::Triangulation<1,
> 1>::MeshSmoothing, dealii::parallel::distributed::Triangulation<1,
> 1>::Settings)
> #1  ./step-40: Step40::LaplaceProblem<1>::LaplaceProblem()
> #2  ./step-40: main
> 
>
> Calling MPI_Abort now.
> To break execution in a GDB session, execute 'break MPI_Abort' before
> running. You can also put the following into your ~/.gdbinit:
>   set breakpoint pending on
>   break MPI_Abort
>   set breakpoint pending auto
>
> 
> An error occurred in line <3455> of file
> 
> in function
> dealii::parallel::distributed::Triangulation<1,
> spacedim>::Triangulation(ompi_communicator_t* const&, typename
> dealii::Triangulation<1, spacedim>::MeshSmoothing,
> dealii::parallel::distributed::Triangulation<1, spacedim>::Settings) [with
> int spacedim = 1; MPI_Comm = ompi_communicator_t*; typename
> dealii::Triangulation<1, spacedim>::MeshSmoothing =
> dealii::Triangulation<1, 1>::MeshSmoothing]
> The violated condition was:
> false
> Additional information:
> You are trying to use functionality in deal.II that is currently not
> implemented. In many cases, this indicates that there simply didn't
> appear much of a need for it, or that the author of the original code
> did not have the time to implement a particular case. If you hit this
> exception, it is therefore worth the time to look into the code to
> find out whether you may be able to implement the missing
> functionality. If you do, please consider providing a patch to the
> deal.II development sources (see the deal.II website on how to
> contribute).
>
> Stacktrace:
> ---
> #0  /home/syed/dealii-candi/deal.II-v9.3.2/lib/libdeal_II.g.so.9.3.2:
> dealii::parallel::distributed::Triangulation<1,
> 1>::Triangulation(ompi_communicator_t* const&, dealii::Triangulation<1,
> 1>::MeshSmoothing, dealii::parallel::distributed::Triangulation<1,
> 1>::Settings)
> #1  ./step-40: Step40::LaplaceProblem<1>::LaplaceProblem()
> #2  ./step-40: main
> 
>
> Calling MPI_Abort now.
> To break execution in a GDB session, execute 'break MPI_Abort' before
> running. You can also put the following into your ~/.gdbinit:
>   set breakpoint pending on
>   break MPI_Abort
>   set breakpoint pending auto
>
> 
> An error occurred in line <3455> of file
> 
> in function
> dealii::parallel::distributed::Triangulation<1,
> spacedim>::Triangulation(ompi_communicator_t* const&, typename
> dealii::Triangulation<1, spacedim>::MeshSmoothing,
> dealii::parallel::distributed::Triangulation<1, spacedim>::Settings) [with
> int spacedim = 1; MPI_Comm = ompi_communicator_t*; typename
> dealii::Triangulation<1, spacedim>::MeshSmoothing =
> dealii::Triangulation<1, 1>::MeshSmoothing]
> The violated condition was:
> false
> Additional information:
> You are trying to use functionality in deal.II that is currently not
> implemented. In many cases, this indicates that there simply didn't
> appear much of a need for it, or that the author of the 

[deal.II] Issue encountered while solving Step-40 in 1 dimension

2022-08-19 Thread syed ansari
Dear all,
I was trying to run step-40 in 1 dimension and encountered the
error corresponding to MeshSmoothing in the constructor. The details of the
error are as follows:

An error occurred in line <3455> of file

in function
dealii::parallel::distributed::Triangulation<1,
spacedim>::Triangulation(ompi_communicator_t* const&, typename
dealii::Triangulation<1, spacedim>::MeshSmoothing,
dealii::parallel::distributed::Triangulation<1, spacedim>::Settings) [with
int spacedim = 1; MPI_Comm = ompi_communicator_t*; typename
dealii::Triangulation<1, spacedim>::MeshSmoothing =
dealii::Triangulation<1, 1>::MeshSmoothing]
The violated condition was:
false
Additional information:
You are trying to use functionality in deal.II that is currently not
implemented. In many cases, this indicates that there simply didn't
appear much of a need for it, or that the author of the original code
did not have the time to implement a particular case. If you hit this
exception, it is therefore worth the time to look into the code to
find out whether you may be able to implement the missing
functionality. If you do, please consider providing a patch to the
deal.II development sources (see the deal.II website on how to
contribute).

Stacktrace:
---
#0  /home/syed/dealii-candi/deal.II-v9.3.2/lib/libdeal_II.g.so.9.3.2:
dealii::parallel::distributed::Triangulation<1,
1>::Triangulation(ompi_communicator_t* const&, dealii::Triangulation<1,
1>::MeshSmoothing, dealii::parallel::distributed::Triangulation<1,
1>::Settings)
#1  ./step-40: Step40::LaplaceProblem<1>::LaplaceProblem()
#2  ./step-40: main


Calling MPI_Abort now.
To break execution in a GDB session, execute 'break MPI_Abort' before
running. You can also put the following into your ~/.gdbinit:
  set breakpoint pending on
  break MPI_Abort
  set breakpoint pending auto


An error occurred in line <3455> of file

in function
dealii::parallel::distributed::Triangulation<1,
spacedim>::Triangulation(ompi_communicator_t* const&, typename
dealii::Triangulation<1, spacedim>::MeshSmoothing,
dealii::parallel::distributed::Triangulation<1, spacedim>::Settings) [with
int spacedim = 1; MPI_Comm = ompi_communicator_t*; typename
dealii::Triangulation<1, spacedim>::MeshSmoothing =
dealii::Triangulation<1, 1>::MeshSmoothing]
The violated condition was:
false
Additional information:
You are trying to use functionality in deal.II that is currently not
implemented. In many cases, this indicates that there simply didn't
appear much of a need for it, or that the author of the original code
did not have the time to implement a particular case. If you hit this
exception, it is therefore worth the time to look into the code to
find out whether you may be able to implement the missing
functionality. If you do, please consider providing a patch to the
deal.II development sources (see the deal.II website on how to
contribute).

Stacktrace:
---
#0  /home/syed/dealii-candi/deal.II-v9.3.2/lib/libdeal_II.g.so.9.3.2:
dealii::parallel::distributed::Triangulation<1,
1>::Triangulation(ompi_communicator_t* const&, dealii::Triangulation<1,
1>::MeshSmoothing, dealii::parallel::distributed::Triangulation<1,
1>::Settings)
#1  ./step-40: Step40::LaplaceProblem<1>::LaplaceProblem()
#2  ./step-40: main


Calling MPI_Abort now.
To break execution in a GDB session, execute 'break MPI_Abort' before
running. You can also put the following into your ~/.gdbinit:
  set breakpoint pending on
  break MPI_Abort
  set breakpoint pending auto


An error occurred in line <3455> of file

in function
dealii::parallel::distributed::Triangulation<1,
spacedim>::Triangulation(ompi_communicator_t* const&, typename
dealii::Triangulation<1, spacedim>::MeshSmoothing,
dealii::parallel::distributed::Triangulation<1, spacedim>::Settings) [with
int spacedim = 1; MPI_Comm = ompi_communicator_t*; typename
dealii::Triangulation<1, spacedim>::MeshSmoothing =
dealii::Triangulation<1, 1>::MeshSmoothing]
The violated condition was:
false
Additional information:
You are trying to use functionality in deal.II that is currently not
implemented. In many cases, this indicates that there simply didn't
appear much of a need for it, or that the author of the original code
did not have the time to implement a particular case. If you hit this
exception, it is therefore worth the time to look into the code to
find out whether you may be able to implement the missing
functionality. If you do, please consider providing a patch to the
deal.II development sources (see the deal.II website on how to
contribute).

Stacktrace:
---
#0