I'm assembling a vector F that represent f(u); the action of f(u) on a
vector field v with FE expansion coefficients V can be represented as
F*V.
OK, so F_i would then be (f(u),\varphi_i), right?
I think the question is whether you do or do not apply
ConstraintMatrix::condense to F. I never
>
> ...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 vecto
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 the
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 -- i
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<
> Dear Rajat,
>
>
>> ... which I then write in a file on
> 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:
https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_dealii_dealii_blob_master_examples_step-2D40_step-2D40.cc-23L610&d=DwIBaQ&c=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLx
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 s
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 AlignedVector:
> On 23 Feb 2017, at 17:08, Daniel Shapero wrote:
>
> 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 con
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 myfile;
>
> 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 yo
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 w
Dear Wolfgang,
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 some simplifications.
Best, Aycil
13 matches
Mail list logo
