Julian Andrej

Yes I configured with --enable-mixed-int. I ran the test on quartz. I did not try to run with the bigint option only, but I can do that if that helps.

Using your suggested options ``` $ srun -n1728 -ppbatch -A ceed ./test/ij -P 12 12 12 -n 1400 1400 1400 -agg_nl 1 ``` works fine. Thanks!

I tried another option since in mfem a simple Laplace problem works fine with mixed-int. Elasticity with the systems option fails. When I run ``` $ srun -n1728 -ppbatch -A...

The test also fails with a much smaller allocation ``` $ srun -n1728 -ppbatch -A ceed ./test/ij -P 12 12 12 -n 700 700 700 -sysL 2 -agg_nl 1 -interptype...

Running ``` $ srun -n1728 -ppbatch -A ceed ./test/ij -P 12 12 12 -n 700 700 700 -sysL 2 -agg_nl 1 -interptype 10 -nf 2 -nodal 1 ``` segfaults without...

Kind of. Derived fields belong to a physics module. Imagine a field that you create on the same function space of the main grid function or similar.

Is this ready for review? It's a small and non-intrusive change that we can review pretty quickly.

If your function is `f(x) = x^2-2` then `df(x)/dx = 2.0*x`. Your Jacobian returns `df(x)/dx = 2.0`.

#2613 aren't reviewed changes. There was no thorough process to check if those changes actually work like intended. The approach that you took seems reasonable. The "discontinuous look" in your...

Sounds great! If you want you can share your effort with others and submit a PR with this as an example?