Real antiperiodic space leads to an error

NGSolve
,
1

Dear NGsolve developers,

I would like to report a small issue. I noticed that one of my old codes is not working with the latest version 6.2.2606.post40.dev0, while it was running fine on version 6.2.2506. Below is a minimal example

from ngsolve import *
from netgen.occ import *

# Geometry & mesh with bottom-top periodicity

box = Rectangle(1,1).Face()
box.edges[0].name = "bottom"
box.edges[2].name = "top"
box.edges[0].Identify(box.edges[2], "top")
mesh = Mesh(OCCGeometry(box, dim=2).GenerateMesh(maxh=0.1))

# Anti-periodic space

fes = Periodic(H1(mesh), phase = [-1]) # runs in older version, not in the latest 

...

In the latest version, I get the error message

NgException: Shouldn't get here: QuasiPeriodicFESpace::TransformMR, space with complex factors should be complexin Assemble BilinearForm 'biform_from_py'

NGSolve does not report an error with complex = True when calling Periodic, but I would prefer not to use a complex finite element space when the “phase” and other coefficients of the physical problem are purely real.

2

Thank you for the report, we will fix that. A quick fix for you is to give the phase as float:


 phase = [-1.])