-
Notifications
You must be signed in to change notification settings - Fork 0
/
test.py
40 lines (31 loc) · 1019 Bytes
/
test.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
# -*- fenics/dependencies: ("ffc"); pythonic-environment: "/docker:fenics@fenics-dev:/usr" -*-
# Recall that fenics-dev is 127.0.0.1 in /etc/hosts at the docker host
from __future__ import print_function
from dolfin import *
import ffc
from ffc.log import add_logfile, set_level, DEBUG
parameters["form_compiler"]["representation"] = "quadrature"
set_level(DEBUG)
add_logfile("/tmp/fenics.log")
mesh = UnitSquareMesh(32, 32)
V = FunctionSpace(mesh, "Hermite", 3)
u = TrialFunction(V)
v = TestFunction(V)
a = inner(div(grad(u)), div(grad(v)))*dx
with open("/tmp/biharmonic.h", "wt") as f:
out = ffc.compile_form(a, prefix="bh")
#print(out[1])
f.write(out[0])
################################################################################
if False:
from dolfin import *
mesh = UnitIntervalMesh(10)
V = FunctionSpace(mesh, 'Hermite', 3)
f = Argument(V, 0)
R = assemble(f*dx)
print("done")
def x2(x=0):
return 2*x
def x3(x=0):
return 3*x
print(x)