Implement Blockstructuring in a sum factorization frame of mind
As talked about in Bonn. Here is a quick run-through:
Part 1, CG:
-
Implement theta matrices with hat functions for CG -
Optimize innermost mat-mat-mult loop to honor sparsity -
Determine total cost of unoptimized boundary stuff -
Add configuration and compiler code path for combined approach (no new interface necessary!)
Part 2, DG:
-
Implement patchwise DG1 basis -
Calculate inner interior_facet integrals in alpha_volume (this the hard one, details in my notes) -
Calculate outer interior_facet integrals in alpha_skeleton (same considerations as above!) -
Adapt driver macro element degree calculation