Re: [deal.II] Memory error from utilities.cc

[Raghunandan Pratoori]

Thank you for your suggestion Prof. Bangerth. I was able to pinpoint the 
problem in my code. I am trying to write all the stress and strain 
components to the output *vtu files. Below is how initialize the variables 
needed for that before filling them up with quadrature point history. The 
lines of code in red are where the memory error is arising from. Is there a 
better way of handling this process, or am I making some fundamental 
mistake?

//Output stress and strain componenets
std::vector>> history_field_stress(dim, 
std::vector>(dim)),

local_history_values_at_qpoints(dim, std::vector>(dim)),

local_history_fe_values(dim, std::vector>(dim));

std::vector>> history_field_strain(dim, 
std::vector>(dim));

for (unsigned int i=0; i On 8/19/22 13:29, Raghunandan Pratoori wrote:
> > 
> > I am trying to run a simulation with grid refine factor 7. I know this 
> is 
> > significantly large and any improper code will raise memory issues. I am 
> in 
> > fact getting memory issues after completion of first time step and I am 
> not 
> > able to pin point where I probably am making a mistake.
>
> Raghunandan:
> somewhere in your code, there is a place where something wants to allocate 
> 134 
> MB of memory, but your system does not have this much memory left. You 
> need to 
> find out where this is, either by strategically placing print statements 
> into 
> your code to see how far it runs, or (better!) by running the program in a 
> debugger.
>
> Best
> W.
>
> -- 
> 
> Wolfgang Bangerth email: bang...@colostate.edu
> www: http://www.math.colostate.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.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/9394451e-b300-4348-9d69-db43186dff0bn%40googlegroups.com.

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

[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, 1>::Me