Hi Jean-Paul,
I can already run my code in the cluster with the lapack and blas module.
Thank you for your time.
Joaquin
2017-03-29 2:52 GMT-04:00 Jean-Paul Pelteret :
> Hi Joaquin,
>
> Its been a really long time since I compiled Lapack manually and linked it
> against deal.II (at that time,
Always thanks for you.
I have read step-30, and now i understand how it is done,
I was trying to implement on my problem as like
Frist, refine_mesh function , I have designated which cell to be refined,
template
void
*StokesProblem::refine_mesh (**const unsigned int
refinem
Either that or create an inverse mapping that maps back from the current to
the reference configuration. This method could potentially introduce to
some numerical error though.
On Thursday, March 30, 2017 at 6:57:12 PM UTC+2, Wolfgang Bangerth wrote:
>
>
> Rajat,
>
> > I am updating my mesh aft
Rajat,
I am updating my mesh after every time-increment. So, cell->vertex(v)
contains the current coordinates of the mesh. I know the coordinates of
the cell (and whole mesh) at t = 0, I wanted
to know how I can initialize the fe_values object wrt to the coordinates
of the cell at t = 0.
If po
There is currently no automated way to parse this, but I just created
an example that might help you:
https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_dealii_dealii_pull_4143&d=DwIFaQ&c=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4&r=4k7iKXbjGC8LfYxVJJXiaYVu6FRWmEjX38S7JmlS9Vw&m=dgi6rq
On 03/30/2017 03:38 PM, 'Seyed Ali Mohseni' via deal.II User Group wrote:
Dear Thomas Wick,
I see. So I should change it manually everywhere. I thought the degree
variable will do this automatically, but why did you then implement it
like this? Is it for another purpose?
It is just copy&pas
> To be honest, I don't even see the point of having 9 integration points for
> a 4-noded quadrilateral or am I overseeing something?
Yes, you are missing some background information:
Quadrature formulas are exact for polynomials up to a certain degree.
For Gauss quadrature you need degree n to in
Dear Thomas Wick,
I see. So I should change it manually everywhere. I thought the degree
variable will do this automatically, but why did you then implement it like
this? Is it for another purpose?
Kind regards,
S. A. Mohseni
--
The deal.II project is located at http://www.dealii.org/
For mai
Dear S. A. Mohseni,
the Gauss points are assigned in that program as always in deal.II:
QGauss quadrature_formula(degree + 2);
and in the argument, you can change to whatever number you want:
QGauss quadrature_formula(n);
with n=1 or n=2 or n = 3 for instance.
Best regards,
Thomas W.
--
+
Dear Timo Heister and Thomas Wick,
I am trying to run your phase-field crack propagation example, namely the
single edge notched tension test, with less than 9 Gauss points.
Hence, the fe variable is initialized using FE_Q(degree) where degree
cannot be 0. This means degree should always be >1
Hi
I'm trying to create a file with all the parameters I need for my program.
I'm using the ParameterHandler class for this. One parameter could be
defined as a map and for this I try to use Patterns::Map
std::map my_map;
however, when I try to parse the parameter as:
my_map = prm.get ("map st
11 matches
Mail list logo
