Is it possible to solve a variable stiffness problem?

[BraveDrXuTF]

Is it possible to solve a variable stiffness prolem, where the stiffness would change according to the spatial coordinates? For example, consider a Darcy type equation $\nabla k(\mathbf{x}) \nabla T=q$, where $\mathbf{x}$ is the spatial coordinates.

[gdmcbain]

Yes indeed; e.g.

scikit-fem/docs/examples/ex16.py

Line 3 in 6159aa3

This example demonstrates a spatially varying coefficient.

scikit-fem/docs/examples/ex16.py

Lines 32 to 34 in 6159aa3

	@BilinearForm
	def stiffness(u, v, w):
	return (1 - w.x[0]**2) * dot(d(u), d(v))

After integration by parts yours will look like this; k (x) being the factor in parentheses.

[BraveDrXuTF]

Thank you!