Thank you for the hints. It works now.
Best regards,
Paras
On Wed, Mar 6, 2024 at 12:36 AM Wolfgang Bangerth
wrote:
> On 3/5/24 09:59, Paras Kumar wrote:
> >
> > Thank you for your response. That is exactly what I am trying to do via
> the
> > deserialize fun
the call to
dealii::Utilities::pack().
Any clues on how to achieve this would be really helpful.
Best regards,
Paras
On Tue, Mar 5, 2024 at 5:40 PM Wolfgang Bangerth
wrote:
> On 3/5/24 07:46, Paras Kumar wrote:
> >
> > Within the context of an MPI parallel code, I ne
.
Could someone please help in understanding how I could solve this issue?
Thanks and best regards,
Paras Kumar
--
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
Wasim,
In case, what you do is close to the phase-field fracture model (which is
essentially a gradient damage model, and in this case a one-way coupled
one), you might find the following master thesis and the associated code
helpful: https://github.com/Jatandeep/Thesis
Best regards,
Paras
On
Dear Wolfgang,
Thank you for the response.
> Paras:
> I'm not sure any of us have experience with Amesos:SuperLU, so I'm not
> sure
> anyone will know right away what the problem may be.
>
> I was wondering if, while writing the wrappers and testing them out,
someone managed to figure out the
installed on the
cluster. The gcc compiler and corresponding openmpi@4.1.2 are available
form the cluster.
Any ideas on solving this issue would be of great help.
Kind regards,
Paras Kumar
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
HI Amit,
Focusing on the question of mesh generation, it is possible to generate an
automated HEX mesh using Gmsh, but since that mesh is based on Delaunay
based approaches, quality may not be very good. See discussions in this
thread:
a second such object that has
> the
> opposite sign.
>
> Best
> W.
>
>
> On 4/14/22 11:28, Paras Kumar wrote:
> > *** Caution: EXTERNAL Sender ***
> >
> > Dear dealii Community,
> >
> > I am working on a solid mechanics related problem
Dear dealii Community,
I am working on a solid mechanics related problem involving anti-periodic
constraints, i.e.
u_leftFace = - u_RightFace,
where u is the unknown field (a component of the dim-dimensional
displacement field).
One possible solution, I am aware of, involves marking the
>
> > 2. The ExceptionBase::print_stack_trace() function does not print
> > anything if I catch an exception thrown using AssetThrow(cond,
> > dealii::ExcMessage()) or even my own MyException class derived from
> > ExceptionBase. How do I ensure that the stack trace is "populated". I
> > tried
Dear Wolfgang,
Thank you for the quick response.
This is actually very difficult to achieve. I made that work for the 9.3
> release when you use `new_task()`, and I think one could probably make
> that
> work for `new_thread()` as well based on features C++11/14 provides. But I
> don't think
Dear deal.ii Community,
I am trying to work on the exception handling aspect of my code, which
employs dealii::Workstream::run() for thread based parallelization of FE
system assembly. The code employs the Assert and AssertThrow macros, in
conjunction with dealii::ExcMessage(), at several
Hi Corbin,
The cell-wise solution corresponding to the p field could be determined as
follows:
```
dealii::FEValuesExtractorScalar pressureExtractor(2); //assuming you have
three scalar fields and pressure is the last one
dealii::Vector solVecPerCell(nDOFsPerCell);
Hi Simon,
This could be helpful:
https://github.com/dealii/dealii/wiki/Frequently-Asked-Questions#how-to-get-the-mapped-position-of-support-points-of-my-element
Here's a code snippet for Option-1, which I usually use.
const auto & feSystem = dofHandler.get_fe();
Dear Wolfgang,
Thank you so much for the suggestion. Just to avoid confusion for future
readers, the to_value() function needs a unique_ptr, and hence the
aforementioned line needs to be replaced by the following lines:
std::unique_ptr nItersListPattern(new
Dear All,
I am trying to read a list of integers from the parameter file using the
dealii::Patterns::List class.
However, I get a runtime error which can be reproduced using the attached
MWE.
I also tried the approach described (for tensor) in
some other boost functions work fine. Also, the
same code compiles on my local machine having Ubuntu 18.04
Could someone please help me understand where I go wrong?
Best regards,
Paras Kumar
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
htt
shan...@gmail.com> wrote:
>>>
>>>> Hello Paras,
>>>>
>>>> Thanks for your response. I have been using the GUI of Gmsh to generate
>>>> the mesh. Do you have any suggestions about what to do in Gmsh to avoid the
>>>> negative v
I use the .inp format of Abaqus, since my code relies on material-IDs which
(in my limited experience) are easy to handle with .inp format. I generate
meshes using both GMSH (here as you mentioned, one needs to be careful to
add curves/surface parts in the correct order to avoid negative volumes)
Dear Prof. Bangerth,
Thank you for the appreciation.
In the quest of alternative solutions, I came across the Mesquite library
and was curious to know if that could be helpful to improve the quality of
the Delaunay based Hex mesh so as to avoid the issue with stress
irregularities.
It also
Bangerth wrote:
>
> On 8/7/20 11:20 AM, Paras Kumar wrote:
> >
> > I have a query regarding writing the solution to a VTK file for
> visualization
> > with PARAVIEW. If one wanted to do this via a mapping based on the
> > displacement field, then one could simpl
Dear Deal.II Community,
I am trying to model the phase-field fracture problem, which, essentially
is a multi field problem involving a vector valued displacement field
coupled with a scalar crack phase-field. Both the variables are
approximated using FE_Q type ansatz but the polynomial degree
ed? I could not find it in the
source code.
Best regards,
Paras
On Tuesday, August 4, 2020 at 6:15:47 PM UTC+2, Paras Kumar wrote:
>
> Hi,
>
> Could you please provide an MWE to describe how the
> Patterns::Convert::to_value() function would work in this case.
>
> Is it must to use pr
Hi,
Could you please provide an MWE to describe how the
Patterns::Convert::to_value() function would work in this case.
Is it must to use prm.add_parameter() to be able to do so? I usually use
prm.declare_entry() and prm.get().
Best regards,
Paras
On Wednesday, April 15, 2020 at 6:01:53 PM
Hi Timo,
The video links don't seem to work. Could you please have a look.
Thanks and best regards,
Paras
On Wednesday, May 20, 2020 at 2:40:31 AM UTC+2, Timo Heister wrote:
>
> Hi all,
>
> we would like to announce a one day-deal.II workshop on May 26, 2020 with
> several talks about recent
Dear Wolfgang,
Thank you for the reply. I will try to add it soon.
Best,
Paras
On Thursday, May 21, 2020 at 5:04:19 PM UTC+2, Wolfgang Bangerth wrote:
>
> On 5/21/20 4:55 AM, Paras Kumar wrote:
> >
> > Could someone please confirm if deal.II offers an analogue of
>
Dear Deal.II community,
Could someone please confirm if deal.II offers an analogue of MPI_Bcast(),
i.e. a function to broadcast arbitrary types, since there are functions
like dealii::Utilities::MPI::gather().
An alternative solution could be to pack the type into a stream using
Ahmad,
On Tuesday, April 14, 2020 at 5:39:06 PM UTC+2, Ahmad Shahba wrote:
>
> I don't know if it is the optimal way but I would use the following
> approach
>
>1. Read tensor components as Patterns::List
>2. Use get
>
>
w, if there is a smarter way (avoiding this manual
copying by looping over the tensor entries) to do so using some other type
of dealii::Patterns object or some other way.
Thanks in advance and best regards,
Paras Kumar
--
The deal.II project is located at http://www.dealii.org/
For mailing list/f
Dear Bruno,
Thank you for your response.
On Friday, February 14, 2020 at 2:34:41 PM UTC+1, Bruno Turcksin wrote:
>
> Paras,
>
> You could try to use SuperLU_MT (see
> https://portal.nersc.gov/project/sparse/superlu/) but we don't have
> wrapper for it. The Krylov solver in deal.II are
Dear deal.ii Community,
I am working on finite deformation hyperelasticity problem which is
essentially a nonlinear-vector-valued problem with displacement as the
unknown at each support point(dim=2,3). *With regards to parallelism, we
currently restrict ourselves to shared memory parallel
Dear Wolfgang,
Thank you so much for your help. The function DataOut::set_cell_selection()
seems to be quite useful.
Best regards,
Paras
On Monday, February 3, 2020 at 12:01:39 AM UTC+1, Wolfgang Bangerth wrote:
>
>
> Finally the fix: https://github.com/dealii/dealii/pull/9466
>
> Paras:
Hi Amy,
I post below one possible solution which worked for me.
On Tuesday, October 1, 2019 at 4:44:25 AM UTC+2, Amy Kaczmarowski wrote:
>
> Hi,
>
> I'm attempting to apply a displacement boundary condition on a surface in
> my problem. However, I would like the direction of the displacement
33 matches
Mail list logo
