[deal.II] Choosing between adaptive BDF and IRK integrators: adaptive meshing effect?

2016-11-15 Thread kosakian
Dear deal.ii team, I am currently at the stage of choosing a transient solver for our deal.ii-based software OpenFCST. The solver will be applied to solving nonlinear transport and reaction problems. After days of research, I found two options: 1. Variable order, variable step BDF schemes (

Re: [deal.II] condition number

2016-11-15 Thread Wolfgang Bangerth
Anup, I wanted to know the condition number of my system_matrix and hence I have used the line in blue (below) in my CG solve function. It shows the condition number. But in the documentation the definition of this is not mentioned. You mean the documentation is missing a definition of what t

Re: [deal.II] Re: Publications based on deal.II

2016-11-15 Thread Michael Harmon
Thank you! -- 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

Re: [deal.II] Re: mesh quality during deformations using MappingQEulerian

2016-11-15 Thread Wolfgang Bangerth
Thomas, My initial mesh is not great - here's where I begin: | |GridGenerator::hyper_sphere(triangulation,Point(0,0,0),1.0); triangulation.set_all_manifold_ids(0); triangulation.refine_global(initial_global_refinement); | |/* interpolate the map from triangulation to desired geometry for the

[deal.II] Re: mesh quality during deformations using MappingQEulerian

2016-11-15 Thread thomas stephens
> > I can't say whether that's the same problem that is commonly reported in > the > literature, but this sort of mesh entanglement is the motivation for > people to > look at Arbitrary Lagrangian Eulerian (ALE) methods, instead of just > Eulerian > methods. > > Can you show a picture of yo

Re: [deal.II] Re: Publications based on deal.II

2016-11-15 Thread Wolfgang Bangerth
On 11/15/2016 10:10 AM, Michael Harmon wrote: Not sure if its too late for this, but here are some publications thanks to deal.ii! Never too late! @phdthesis{Harmon2016, Title= Numerical algorithms based on Galerkin methods for the modeling of reactive interfaces in pho

[deal.II] Re: Error when plotting curved interior cells with vtk

2016-11-15 Thread Wolfgang Bangerth
Just to follow up on this, many moons later: The fix for this problem turns out to be only one line long. That didn't make it easier to find, though :-) It can be found at https://github.com/dealii/dealii/pull/3589 and should become part of the next release. Best Wolfgang -- The deal.II pr

[deal.II] meshworker getting gradient values

2016-11-15 Thread Sudarshan Kumar
dim=2 template void Step12::integrate_cell_term (DoFInfo& dinfo, CellInfo& info) { const FEValuesBase& fe_v = info.fe_values(); const std::vector>& grad = info.gradients[0][0]; } When I call grad[0][0] it gives segmentation fault. Could any one tell how to access the gradie

[deal.II] Re: Periodic Boundary Condition

2016-11-15 Thread Daniel Arndt
Hamed, Am Dienstag, 15. November 2016 19:09:18 UTC+1 schrieb Hamed Babaei: > > Hi all, > > Following step-45, I am going to implement periodic boundary condition on > my code for a Thermo-elastic problem, in which the thermal and elastic > parts are solved separately (I have two different dof_ha

[deal.II] Re: Periodic Boundary Condition

2016-11-15 Thread Hamed Babaei
In addition I just tried to apply periodic boundary condition constraints at the end of assembly using, void ConstraintMatrix::condense. However, it seems this function doesn't work with TrilinosWrappers::sparsematrix but dealii ones. Thanks -- The deal.II project is located at http://www.dea

[deal.II] Periodic Boundary Condition

2016-11-15 Thread Hamed Babaei
Hi all, Following step-45, I am going to implement periodic boundary condition on my code for a Thermo-elastic problem, in which the thermal and elastic parts are solved separately (I have two different dof_handelers for temprature and displacement fields without any block structure). I am won

Re: [deal.II] question on doxygen

2016-11-15 Thread Timo Heister
> Anyone else have this issue? If you inspect the page on chrome you will get: index.html:1 Mixed Content: The page at 'https://urldefense.proofpoint.com/v2/url?u=https-3A__mdh266.github.io_PECS_index.html&d=CwIBaQ&c=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4&r=4k7iKXbjGC8LfYxVJJXiaYVu6FRWmEjX38S

[deal.II] Re: question on doxygen

2016-11-15 Thread 'Uwe Köcher' via deal . II User Group
Hej Mike, sometimes I also use latex for code documentation with doxygen, for an equation line: \f[ \partial_t u -\nabla \cdot \epsilon \nabla u + \boldsymbol b \cdot \nabla u = f\,,\quad \text{in } \Omega \times I\,, \f] and \f$ \Omega = (0,1)^d \f$ for latex within the text line. I'

[deal.II] Re: Publications based on deal.II

2016-11-15 Thread Michael Harmon
Not sure if its too late for this, but here are some publications thanks to deal.ii! @Article{HaGaRe-JCP2016, Title= Numerical algorithms based on Galerkin methods for the modeling of reactive interfaces in photoelectrochemical (PEC) solar cells}, Author

[deal.II] question on doxygen

2016-11-15 Thread Michael Harmon
Hi, I used deal for my phd work, loved using it and wanted to put my code up on github. I spent a lot of time documenting it with doxygen, but latex equations only work on my local machine and not on github: https://mdh266.github.io/PECS/index.html Anyone else have this issue? Thanks! Mike

[deal.II] condition number

2016-11-15 Thread Anup Basak
Hello all, I wanted to know the condition number of my system_matrix and hence I have used the line in blue (below) in my CG solve function. It shows the condition number. But in the documentation the definition of this is not mentioned. I mean when I store my system matrix and calculate its max

Re: [deal.II] reuse FEEvaluation within a cell for RHS-like calculation

2016-11-15 Thread Denis Davydov
Hi Martin, Cool, thanks for the prompt reply. Kind regards, Denis > On 15 Nov 2016, at 14:59, Martin Kronbichler > wrote: > > Hi Denis, > > > >> I wonder if it possible to re-use FEEvaluation within each cell to calculate >> integrals of multiple right-hand-side like vectors. > Yes, you

Re: [deal.II] reuse FEEvaluation within a cell for RHS-like calculation

2016-11-15 Thread Martin Kronbichler
Hi Denis, I wonder if it possible to re-use FEEvaluation within each cell to calculate integrals of multiple right-hand-side like vectors. Yes, you can re-use an FEEvaluation object. As soon as you call submit_value(), you will write into the data field that is used by integrate(), which in t

[deal.II] reuse FEEvaluation within a cell for RHS-like calculation

2016-11-15 Thread Denis Davydov
I wonder if it possible to re-use FEEvaluation within each cell to calculate integrals of multiple right-hand-side like vectors. Something along these lines: for (unsigned int cell=0; cellhttp://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/deali