- Thank you received: 0
VectorL2
7 years 6 days ago #258
by pschroe
Hi guys,
I have a problem with the new VectorL2 FE space...
Attached you will find a simple SIP-DG method for a vector-valued Poisson problem.
Unfortunately, when I run the script with the VectorL2 space, I get the following error:
RuntimeError: undefined fespace 'VectorL2'
When the VectorH1 space is used, everything works just fine.
What am I doing wrong?
Thank you for your help!
Best,
Philipp
I have a problem with the new VectorL2 FE space...
Attached you will find a simple SIP-DG method for a vector-valued Poisson problem.
Unfortunately, when I run the script with the VectorL2 space, I get the following error:
RuntimeError: undefined fespace 'VectorL2'
When the VectorH1 space is used, everything works just fine.
What am I doing wrong?
Thank you for your help!
Best,
Philipp
Attachments:
7 years 6 days ago #259
by cwinters
Hi Philipp,
I get the same error when using the VectorL2. I'll have a closer look at the VectorL2 tomorrow.
You could also use a space with more components. This can be done with the "dim" argument.
Using your formulation it is still not working. When I switch to a formulation on the element boundaries, everything is fine. I'll check if there is an issue with the skeleton BFI's.
Attached is a file which should do what you want, using and L2 space with as many components as your mesh dimension and an element boundary formulation.
Best,
Christoph
I get the same error when using the VectorL2. I'll have a closer look at the VectorL2 tomorrow.
You could also use a space with more components. This can be done with the "dim" argument.
Code:
V = L2(mesh, order=order,dim=mesh.dim, flags={"dgjumps":True})
Using your formulation it is still not working. When I switch to a formulation on the element boundaries, everything is fine. I'll check if there is an issue with the skeleton BFI's.
Attached is a file which should do what you want, using and L2 space with as many components as your mesh dimension and an element boundary formulation.
Best,
Christoph
Attachments:
The following user(s) said Thank You: pschroe
7 years 6 days ago #263
by pschroe
Hi again,
is it possible that now somehow the dimensions of the div operator are not clean?
When I add the div-div term
I get:
RuntimeError: SymblicBFI needs scalar-valued CoefficientFunction
This does not happen with the VectorH1.
Best,
Philipp
is it possible that now somehow the dimensions of the div operator are not clean?
When I add the div-div term
Code:
a += SymbolicBFI( div(u)*div(v) )
I get:
RuntimeError: SymblicBFI needs scalar-valued CoefficientFunction
This does not happen with the VectorH1.
Best,
Philipp
Time to create page: 0.111 seconds