Use of ParameterHandler in the FSI code by T. Wick

29 views
Skip to first unread message

Joaquin

unread,
Jul 15, 2017, 10:33:37 PM7/15/17
to deal.II User Group
Dear all:

I am trying to use ParameterHandler for runtime parameters in the ALE FSI code developed by Thomas Wick. I get the following error:


jomivalen@Nalia ~/ale/jvale-1 $ ./jvale-1
ParameterHandler::read_input: could not open file <ale_fem.prm> for reading.
Trying to make file <ale_fem.prm> with default values for you.

--------------------------------------------------------
An error occurred in line <367> of file </home/jomivalen/thesis/dealii-8.4.1/include/deal.II/base/smartpointer.h> in function
    T* dealii::SmartPointer<T, P>::operator->() const [with T = const dealii::Triangulation<2, 2>; P = dealii::DoFHandler<2>]
The violated condition was:
    t != 0
The name and call sequence of the exception was:
    ExcNotInitialized()
Additional Information:
(none)

Stacktrace:
-----------
#0  /home/jomivalen/cfem/dealii/lib/libdeal_II.g.so.8.4.1: dealii::SmartPointer<dealii::Triangulation<2, 2> const, dealii::DoFHandler<2, 2> >::operator->() const
#1  /home/jomivalen/cfem/dealii/lib/libdeal_II.g.so.8.4.1: void dealii::internal::DoFHandler::Implementation::reserve_space<2>(dealii::DoFHandler<2, 2>&)
#2  /home/jomivalen/cfem/dealii/lib/libdeal_II.g.so.8.4.1: dealii::DoFHandler<2, 2>::distribute_dofs(dealii::FiniteElement<2, 2> const&)
#3  ./jvale-1: ALEFEM<2>::setup_system()
#4  ./jvale-1: ALEFEM<2>::run()
#5  ./jvale-1: main
--------------------------------------------------------

Aborted


What does this error mean?

Thanks to anyone who may help me.

Joaquín

Thomas Wick

unread,
Jul 16, 2017, 3:08:58 AM7/16/17
to dea...@googlegroups.com
Dear Joaquín !

I never implemented a parameter file for this code. For this reason it is difficult to say in view of the original code what
exactly you have changed.

Three remarks:

1. It seems that the parameter file cannot be found and for this reason some (maybe strange!?) values
    are given that cause the error in the program. Please try to link your code to the parameter
    file and see if it runs again.

2. Can you use gdb to find out more about the error message?

3. The original ANS code was developed for deal.II version 7.3.0 and will not run
    on deal.II version 8.5.0. You know that I have updates on my webpage.
    (But honestly I do not think that this is the error using the old FSI code with the
     newest deal.II version.)

Best regards,

Thomas

-- 
++--------------------------------------------++
Dr. Thomas Wick
Maitre de conferences / Assistant Professor

Centre de Mathematiques Appliquees (CMAP)
Ecole Polytechnique
91128 Palaiseau cedex, France

Email:  thoma...@cmap.polytechnique.fr
Phone:  0033 1 69 33 4579
www:    http://www.cmap.polytechnique.fr/~wick/
++--------------------------------------------++
--
--
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+un...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

-- 
++--------------------------------------------++
Dr. Thomas Wick
Maitre de conferences / Assistant Professor

Centre de Mathematiques Appliquees (CMAP)
Ecole Polytechnique
91128 Palaiseau cedex, France

Email:  thoma...@cmap.polytechnique.fr
Phone:  0033 1 69 33 4579
www:    http://www.cmap.polytechnique.fr/~wick/
++--------------------------------------------++
--

Wolfgang Bangerth

unread,
Jul 16, 2017, 11:19:33 AM7/16/17
to 'Joaquin' via deal.II User Group

Joaquin,

> I am trying to use ParameterHandler for runtime parameters in the ALE FSI code
> developed by Thomas Wick. I get the following error:
>
>
> *jomivalen@Nalia ~/ale/jvale-1 $ ./jvale-1
The first and second errors are unrelated. The ParameterHandler class issues
the message about creating the file for you, and continues. Then, at a later
point and unrelatedly, you get the second error.

I suspect that the error happens because you created a DoFHandler and try to
distribute DoFs on it, but have forgotten to previously associate a
triangulation with the DoFHandler. Could that be?

Best
W.

--
------------------------------------------------------------------------
Wolfgang Bangerth email: bang...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/

Reply all
Reply to author
Forward
0 new messages