2025-06-13 09:31:05 -04:00
|
|
|
from serif import config
|
|
|
|
|
from serif.polytrope import PolySolver
|
2025-06-16 12:01:03 -04:00
|
|
|
from serif.mfem import Matrix
|
|
|
|
|
from serif.mfem import Operator
|
|
|
|
|
from serif.mfem.forms import BilinearForm
|
2025-06-13 09:31:05 -04:00
|
|
|
|
|
|
|
|
config.loadConfig('../../testsConfig.yaml')
|
|
|
|
|
n = config.get("Tests:Poly:Index", 0.0)
|
|
|
|
|
|
|
|
|
|
polytrope = PolySolver(n, 1)
|
2025-06-16 12:01:03 -04:00
|
|
|
polytrope.solve()
|
|
|
|
|
theta = polytrope.get_theta()
|
|
|
|
|
print(theta)
|
|
|
|
|
|
|
|
|
|
FESpace = theta.FESpace()
|
|
|
|
|
print(FESpace)
|
|
|
|
|
|
|
|
|
|
print(theta.as_numpy())
|