Forum Message

 

 

We have moved the forum to https://forum.ngsolve.org . This is an archived version of the topics until 05/05/23. All the topics were moved to the new forum and conversations can be continued there. This forum is just kept as legacy to not invalidate old links. If you want to continue a conversation just look for the topic in the new forum.

Notice

The forum is in read only mode.

VectorL2

More
6 years 5 months ago #258 by pschroe
VectorL2 was created 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
Attachments:
More
6 years 5 months ago #259 by cwinters
Replied by cwinters on topic VectorL2
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.
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
More
6 years 5 months ago #260 by pschroe
Replied by pschroe on topic VectorL2
Thank you, Christoph, for your fast and helpful answer! :)
More
6 years 5 months ago #261 by joachim
Replied by joachim on topic VectorL2
Hi Philipp,

VectorL2 was actually not registered, I just fixed it. Your example is working now (with the next nightly build).

Joachim
The following user(s) said Thank You: pschroe
More
6 years 5 months ago #262 by pschroe
Replied by pschroe on topic VectorL2
Thank you very much, Joachim.
The VectorL2 is very convenient :)
More
6 years 5 months ago #263 by pschroe
Replied by pschroe on topic VectorL2
Hi again,

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.135 seconds

Login