FEniCS / ffcx
89%

DEFAULT BRANCH: main
Repo Added
Files 60
Badge
Embed ▾
README BADGES
x

If you need to use a raster PNG badge, change the '.svg' to '.png' in the link

Markdown

Textile

RDoc

HTML

Rst

LAST BUILD ON BRANCH main
branch: SELECT
CHANGE BRANCH
x

coverage: 88.921% (+0.2%) from 88.75%
41411093991

push

github

web-flow 
Enable sum factorisation by default and improve O2 kernel generation (#906)

* Fix sum factorisation with several or one point quadrature rules

- 1D tensor factor tables were named FE_TF0, FE_TF1, ... from zero for
  each quadrature rule, so the tables of integrals with different
  quadrature degrees clashed and one rule read the other's tables, giving
  a wrong element tensor. Include the quadrature rule in their name.
- A piecewise table with tensor factors, such as a coefficient on a one
  point rule, is accessed outside the quadrature loop, where the
  quadrature index can't be split by direction, which failed an
  assertion. Its factors are constant over the points, so read point 0.

Co-Authored-By: Claude Opus 5.5 (1M context) <noreply@anthropic.com>

* Make quadrature weights tensor products

* Add opt-in sum-factorised contraction for tensor product cell integrals

With `sum_factorization_contraction`, cell integrals on tensor product
rules cache the integrand on the whole rule and contract one direction at
a time, instead of updating every element tensor entry at every quadrature
point: O(N^2 Q^2 + N^4 Q) rather than O(Q^2 N^4) in 2D, and O(N^6 Q)
rather than O(Q^3 N^6) in 3D.

- The 1D quadrature weights are folded into the first argument's 1D
  tables at code generation time.
- Terms with the same dof loops share one loop nest; stage arrays and
  contributions shared by several terms are computed once, and each
  element tensor entry is updated once.
- Integrals whose buffers exceed a 256 KiB stack budget, or with an
  argument that doesn't factorise, keep the fused loop.
- Without the option, kernels use the fused loop with the full weights
  table, as before 59e849d6.
- Tensor factor tables are indexed by their own single permutation and
  entity.
- Unread element tables are pruned when the contraction is used.

Co-Authored-By: Claude Opus 5.5 (1M context) <noreply@anthropic.com>

* Plan sum-factorised contractions once and declare only ... (continued)

425 of 426 new or added lines in 11 files covered. (99.77%)

13 existing lines in 4 files now uncovered.

5209 of 5858 relevant lines covered (88.92%)

0.89 hits per line

Relevant lines Covered
5858 RELEVANT LINES 5209 COVERED LINES
0.89 HITS PER LINE
Source Files on main

Recent builds

Builds Branch Commit Type Ran Committer Via Coverage
41411093991 main Enable sum factorisation by default and improve O2 kernel generation (#906) * Fix sum factorisation with several or one point quadrature rules - 1D tensor factor tables were named FE_TF0, FE_TF1, ... from zero for each quadrature rule, so the ... push web-flow github
88.92
40291595106 main Enable sum factorisation by default and improve O2 kernel generation (#906) * Fix sum factorisation with several or one point quadrature rules - 1D tensor factor tables were named FE_TF0, FE_TF1, ... from zero for each quadrature rule, so the ... push web-flow github
88.92
40228706075 main Enable sum factorisation by default and improve O2 kernel generation (#906) * Fix sum factorisation with several or one point quadrature rules - 1D tensor factor tables were named FE_TF0, FE_TF1, ... from zero for each quadrature rule, so the ... push web-flow github
88.92
40227735998 gh-readonly-queue/main/pr-906-efb27548c8a2453afcf7894ee196936d9e10f7c8 Enable sum factorisation by default and improve O2 kernel generation (#906) * Fix sum factorisation with several or one point quadrature rules - 1D tensor factor tables were named FE_TF0, FE_TF1, ... from zero for each quadrature rule, so the ... push web-flow github
88.92
40210089223 main Sum-factorise coefficient evaluation in sum-factorised contractions (#904) * Sum-factorise coefficient evaluation in sum-factorised contractions With `sum_factorization_contraction`, the integrand was cached on the whole tensor product rule, but... push web-flow github
88.75
40208823333 gh-readonly-queue/main/pr-904-4be7614d0089ac5173720e7b3e1670a6679872fc Sum-factorise coefficient evaluation in sum-factorised contractions (#904) * Sum-factorise coefficient evaluation in sum-factorised contractions With `sum_factorization_contraction`, the integrand was cached on the whole tensor product rule, but... push web-flow github
88.75
40207796927 gh-readonly-queue/main/pr-901-b3a4d493855ee7ea713969e94f044fd6cd56ed3d Sum-factorised contraction for tensor product cell integrals (#901) * Fix sum factorisation with several or one point quadrature rules - 1D tensor factor tables were named FE_TF0, FE_TF1, ... from zero for each quadrature rule, so the tables ... push web-flow github
88.61
40136004794 main Add documentation version switcher (#931) * Add documentation version switcher Enable the pydata-sphinx-theme version dropdown. The version list is generated by FEniCS/docs from the published documentation, and CI sets DOCS_VERSION to the branch... push web-flow github
88.1
40135077700 gh-readonly-queue/main/pr-931-044a03e9643057cdd49427f22ee526572c710257 Add documentation version switcher (#931) * Add documentation version switcher Enable the pydata-sphinx-theme version dropdown. The version list is generated by FEniCS/docs from the published documentation, and CI sets DOCS_VERSION to the branch... push web-flow github
88.1
40087880283 main Include the Basix and UFL versions in the JIT signature (#926) Generated kernels embed Basix's tabulated values and depend on UFL's form processing, but the cache key only had the FFCx version, so a Basix or UFL patch release could reuse stale ke... push web-flow github
88.1
See All Builds (3299)