/usr/share/doc/python-ufl-doc/demo/StiffnessAD.ufl is in python-ufl-doc 1.6.0-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 | #
# Author: Martin Sandve Alnes
# Date: 2008-10-30
#
element = FiniteElement("Lagrange", triangle, 1)
w = Coefficient(element)
# H1 semi-norm
f = inner(grad(w), grad(w))/2*dx
# grad(w) : grad(v)
b = derivative(f, w)
# stiffness matrix, grad(u) : grad(v)
a = derivative(b, w)
# adjoint, grad(v) : grad(u)
astar = adjoint(a)
# action of adjoint, grad(v) : grad(w)
astaraction = action(astar)
forms = [f, b, a, astar, astaraction]
|