Chris Richardson

Gauss-Legendre-Lobatto points used for the definition of the Lagrange elements

1

**[Original report](https://bitbucket.org/fenics-project/FIAT/issue/22) by Ben Crestel (Bitbucket: [bcrestel](https://bitbucket.org/bcrestel), GitHub: [bcrestel](https://github.com/bcrestel)).** ---------------------------------------- There is a general agreement that these interpolation points are better suited for high-order elements. It seems like this should...

redundant apts[:,i]-apts[:,i] expressions

3

**[Original report](https://bitbucket.org/fenics-project/FIAT/issue/12) by Nico Schlömer (Bitbucket: [nschloe](https://bitbucket.org/nschloe), GitHub: [nschloe](https://github.com/nschloe)).** ---------------------------------------- Lines like https://bitbucket.org/fenics-project/fiat/src/550063505522529f9a6c3875007a3c0277d819f0/FIAT/newdubiner.py?at=master#cl-145 ``` results[0,:] = 1.0 + apts[:,0]-apts[:,0]+apts[:,1]-apts[:,1]+apts[:,2]-apts[:,2] ``` look buggy ``` 0 == apts[:,0]-apts[:,0] + apts[:,1]-apts[:,1] + apts[:,2]-apts[:,2]...

Order-of-magnitude slowdown in FIAT

7

**[Original report](https://bitbucket.org/fenics-project/FIAT/issue/13) by Andrew McRae (Bitbucket: [amcrae](https://bitbucket.org/amcrae), GitHub: [amcrae](https://github.com/amcrae)).** ---------------------------------------- Hello all, We at Firedrake recently incorporated the changes between 5500635 and ec68699 into our fork of FIAT (i.e. all...

**[Original report](https://bitbucket.org/fenics-project/FIAT/issue/16) by Garth N. Wells (Bitbucket: [garth-wells](https://bitbucket.org/garth-wells), GitHub: [garth-wells](https://github.com/garth-wells)).** ---------------------------------------- Using SymPy in FIAT is like using a sledgehammer to crack a peanut. It's a heavy dependency that imports...

Reimplement PointSource

10

PointSource has been removed and needs to be reimplemented.

After checking a coverage report, it seems like `dolfinx.nls.petsc.NewtonSolver` is not used or tested in the test suite.

Optional documentation seems hard to achieve, see: https://github.com/FEniCS/dolfinx/blob/chris/petsc-optional/python/doc/source/api.rst How can we keep PETSc documentation, but build correctly when it is not installed?

Split pointers for diagonal and off-diagonal parts of MatrixCSR

6

For better compatibility with 3rd party libraries, it would be better to have the diagonal and off-diagonal blocks of the sparse matrix in `MatrixCSR` stored separately. The data structures would...