> I have tried a simple testcase for my problem.
>> I don't know if it's suitable as a testcase but if I use X=0 in my
>> equations above, I know the solution.
>>
>> The solutions is u = v = 0, \phi = (1-x^2 -y^2) and the pressure is then
>> P = 8 \phi_lin^2 (1-x²-y²)
>> It fits with my bound
>
> If the reference solution is 0, you can't see if you are missing a
> constant factor in your equation even when the rate of convergence looks
> correct.
>
I understand thank you !
> Yes, but then you also have to be careful with the boundary conditions. If
> you have an analytical solut
One final remark for today, if I remove the two asymmetrical terms in my
equation so it becomes :
*Then all three codes gives back the same solution...*
So I thought it would come from my assembly of this part of the system but
I have checked and checked again... It should not come from my wea
iting it as 1.0/2.0 or just 0.5.
>
> I’ve got no idea if this is the only source of error in your code, but its
> certainly not helping you :-) Good luck in sorting this out!
>
> Best,
> Jean-Paul
>
> On 25 Jul 2019, at 16:10, Félix Bunel >
> wrote:
>
> One final
Hello everyone.
I'm having some trouble to understand the convergence rate i'm observing in
my code.
Here is what i'm solving :
- I'm in 2D on a round mesh.
- I'm solving a simple Poisson equation on this mesh for a variable named
Phi the solution is known for this and is 1-x^2-y^2
- With this
RESET
<< std::endl;
Le jeudi 10 octobre 2019 14:00:23 UTC+2, Bruno Blais a écrit :
>
> A quick question, since you are working on a sphere, are you specifying a
> mapping of the same order as your phi?
>
> On Wednesday, 9 October 2019 08:57:45 UTC-4,
achieve the correct convergence rates. In your case,
> you are using a MappingQ1 (the default mapping, if you don’t specify
> anything). With that, you won’t go beyond 2nd order in L2. If you use
> MappingQ2, you should see again the optimal convergence rates.
>
> Best
I found why it's not going further down ! It was just that the GPL file
saved is with a precision of 6 digits which corresponds with the error per
dofs when I norm my error plot by the solution.
Thanks for the help everyone !
On Friday, October 11, 2019 at 9:51:40 AM UTC+2, Félix Bunel
Hi,
I have two simple question :
The first one is :
Is it possible to integrate source terms to make them easier to calculate ?
Basically, I can write this : the function xi is my test function, and
sigma is just an independant function.
But will I be able to calculate it in my right hand side