Re: [deal.II] bug in program only with adaptive mesh refinement

2017-02-23 Thread Daniel Shapero
> > ...which you compute via quadrature? Or do you compute a vector F that > corresponds to f(u) somehow? > > If you go the route via vectors you have to pay attention to *what kind > of vector you have*, namely one that does or does not incorporate > constraints. Dealing with dual space

[deal.II] Re: FE_Q constructor

2017-02-23 Thread Aycil Cesmelioglu
Thanks Wolfgang. I'll give it a try assuming everything depends on x and y. Best, Aycil. -- 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

Re: [deal.II] bug in program only with adaptive mesh refinement

2017-02-23 Thread Wolfgang Bangerth
Hi Dan, * In the formula above, P(.) is a functional, I assume, i.e., it takes a function and returns a number, right? * If so, what exactly does f(u) . v actually mean? How do you compute this? * Same for the second derivatives? Sorry if that was unclear --

[deal.II] Re: How to output a single scalar in a parallel code

2017-02-23 Thread RAJAT ARORA
Hello, Please look at the solution by Timo. In your case, you can do something like. if (Utilities::MPI::this_mpi_process(mpi_communicator) == 0) { std::ofstream myfile; myfile.open ("resultant_strass.txt"); myfile<< resultant_stress<

Re: [deal.II] Re: How to output a single scalar in a parallel code

2017-02-23 Thread Timo Heister
> Do I need to make any change in it? You need to only do it on processor with rank 0, for example how we do it in step-40:

Re: [deal.II] Re: FE_Q constructor

2017-02-23 Thread Wolfgang Bangerth
On 02/23/2017 08:47 AM, Aycil Cesmelioglu wrote: I have an elasticity type problem; the horizontal displacement depends on (x,y) but the vertical component of the displacement doesn't depend on y. It's like it is 1.5 dimensional instead of 2. Because of this the variational formulation also has

Re: [deal.II] Segfault in EvaluatorTensorProduct::apply when AVX is enabled

2017-02-23 Thread Stephen DeWitt
Hi Martin, Thanks! I tried with the updated deal.II version and it worked like a charm. Thanks for the help and the bugfix! Best, Steve On Thursday, February 23, 2017 at 1:55:30 AM UTC-5, Martin Kronbichler wrote: > > Dear Stephen, > > The first problem you are seeing is a bug in

[deal.II] Re: How to output a single scalar in a parallel code

2017-02-23 Thread Hamed Babaei
Dear Rajat, > ... which I then write in a file on a master process. > I was wondering how to write in a file on a master process so that there would be just one output file not as many as processor exist. I use the following commands to write in my file in serial code : std::ofstream

[deal.II] Re: bug in program only with adaptive mesh refinement

2017-02-23 Thread Daniel Shapero
> > The vectors you have in this equation P(u+h*v) =..., which of those have > constrains distributed and which zeroed? > If you assembly matrices with ConstraintMatrix.distribute_local_to_global() > the diagonal elements corresponding to constrained DoFs > ... > Along the same lines: does

Re: [deal.II] bug in program only with adaptive mesh refinement

2017-02-23 Thread Daniel Shapero
Hi Wolfgang, thanks for looking! > * In the formula above, P(.) is a functional, I assume, i.e., it takes a > function and returns a number, right? > * If so, what exactly does > f(u) . v > actually mean? How do you compute this? > * Same for the second derivatives? > Sorry if that