(x)
| 309 | e_Q_local[dofs] = e_eval.flatten() |
| 310 | |
| 311 | def e_exact(x): |
| 312 | T = x[0] + 2.0 * x[1] |
| 313 | K = 1.0 / (A.value + B.value * T) |
| 314 | |
| 315 | grad_T = np.zeros((2, x.shape[1])) |
| 316 | grad_T[0, :] = 1.0 |
| 317 | grad_T[1, :] = 2.0 |
| 318 | |
| 319 | e = B.value * K**2 * grad_T |
| 320 | return e |
| 321 | |
| 322 | # # FIXME: Below is only for testing purposes, |
| 323 | # # never to be used in user code! |
no outgoing calls
no test coverage detected