Re: [deal.II] Query using hdf5 with deal.ii

2017-02-14 Thread Wolfgang Bangerth
Rajat, I am using deal.ii to solve a 3D solid mechanics problem. I am using Petsc and P4est along with deal.ii There is a previous post related to my query but it does not talk in detail

[deal.II] Query using hdf5 with deal.ii

2017-02-14 Thread RAJAT ARORA
Hello all, I am using deal.ii to solve a 3D solid mechanics problem. I am using Petsc and P4est along with deal.ii There is a previous post related to my query but it does not talk in det

Re: [deal.II] Parameter map default value

2017-02-14 Thread Wolfgang Bangerth
On 02/13/2017 06:15 AM, Franco Milicchio wrote: Hello, I know this is a weird question, but is there a way of telling a *Map* that there is a default value? For instance I'd like to express something like "*stress:4, power*" instead of what I have to do right now, such as "*stress:4, power:8*"

Re: [deal.II] Re: Recommendation: BCs in config

2017-02-14 Thread Wolfgang Bangerth
On 02/13/2017 07:09 AM, Bruno Turcksin wrote: you could write everything into one string and then split it. So it would like this: subsection bc set value = 0.5*x+y; x*y^2 end Take a look here http://stackoverflow.com/questions/236129/split-a-string-in-c for some ideas how to split a string

Re: [deal.II] Problem with MUMPS Solver

2017-02-14 Thread RAJAT ARORA
Hello Timo, Thanks for your detailed reply giving many possible reasons. There are no constraints on the equation. It is a first order time dependent problem, so there are just initial conditions that are taken care of separately. The code is for 3D brick elements with 9 dofs per node. I just u

Re: [deal.II] Re: Who is owning the instance in an exception raised by Subscriptor?

2017-02-14 Thread Bruno Turcksin
2017-02-14 11:30 GMT-05:00 'Uwe Köcher' via deal.II User Group : > I'm wondering, if it is impossible to call dof.clear(), > which should remove the link to triangulation object, > in python or the c++ code beyond? > > I've included such a line in a destructor within my c++ codes, to get rid of > a

Re: [deal.II] Re: Who is owning the instance in an exception raised by Subscriptor?

2017-02-14 Thread 'Uwe Köcher' via deal . II User Group
I'm wondering, if it is impossible to call dof.clear(), which should remove the link to triangulation object, in python or the c++ code beyond? I've included such a line in a destructor within my c++ codes, to get rid of a similar message as posted above. On Tuesday, February 14, 2017 at 4:21:00

Re: [deal.II] Re: Who is owning the instance in an exception raised by Subscriptor?

2017-02-14 Thread Bruno Turcksin
2017-02-14 10:13 GMT-05:00 Sebastian.Gonzalez-Pintor : > I am doing some reading about this now. I do not have explicit destructors > neither in c++ nor in python, so I want to be carefull that I can define the > destructors to add the cout without changing the what the destructor would > do the wo

Re: [deal.II] Re: Who is owning the instance in an exception raised by Subscriptor?

2017-02-14 Thread Sebastian.Gonzalez-Pintor
Bruno, I am doing some reading about this now. I do not have explicit destructors neither in c++ nor in python, so I want to be carefull that I can define the destructors to add the cout without changing the what the destructor would do the work. Thanks for the hint. Best, Sebas. > -- The

Re: [deal.II] Re: Who is owning the instance in an exception raised by Subscriptor?

2017-02-14 Thread Bruno Turcksin
Sebas, 2017-02-14 9:29 GMT-05:00 Sebastian.Gonzalez-Pintor : > I do not know why this happens, because I do not wrap the interface for > Subscriptor (it is called internally inside the c++ classes, whithout python > knowing about its existence). I will try to ask for this behaviour somewhere > els

[deal.II] Re: Who is owning the instance in an exception raised by Subscriptor?

2017-02-14 Thread Sebastian.Gonzalez-Pintor
Hi @Bruno. In this case the Triangulation is created inside a class called ProbGeom1D, and the DoFHandler is in another class called State1D, created later and that uses an instance of ProbGeom1D as an input. geo = pygeom.ProbGeom1D(inp) state = pystate.State1D(inp, geo); If I finish my scri

[deal.II] Re: Who is owning the instance in an exception raised by Subscriptor?

2017-02-14 Thread Bruno Turcksin
Hi, most of time this problem happens because the DoFHandler was declared before the Triangulation. In the destructor, the Triangulation will be destroyed before the DoFHandler. Best, Bruno On Tuesday, February 14, 2017 at 8:35:10 AM UTC-5, Sebastian.Gonzalez-Pintor wrote: > > > > Hello! I a

[deal.II] Who is owning the instance in an exception raised by Subscriptor?

2017-02-14 Thread Sebastian.Gonzalez-Pintor
Hello! I am using swig to parse some code that I have wrote, and this code uses dealii. When I run the code with python, I get an error raised by Subscriptor as follows An error occurred in line <114> of file in function virtual dea

Re: [deal.II] Parallelization of step-2: DynamicMatrix to SparsityMatrix

2017-02-14 Thread Kartik Jujare
Thank you Timo for your reply. On Sunday, February 12, 2017 at 3:54:19 PM UTC+1, Timo Heister wrote: > > Yes, it looks like you can not copy a DSP with an IndexSet into a > SparsityPattern. We could make that work but it also is not a very > useful operation. Can you try replacing > DynamicSp