basix
basix copied to clipboard
FEniCSx finite element basis evaluation library
Transposing the last two indices will mean we don't have to transpose the data when assigning it. This will need changes in DOLFINx and FFCx.
Judging by the build errors at https://buildd.debian.org/status/package.php?p=basix, looks like Nedelec 2nd kind H(curl) is not compatible with 32 bit systems. e.g. i386 https://buildd.debian.org/status/fetch.php?pkg=basix&arch=i386&ver=0.0.1%7Egit20210122.4f10ef2-1&stamp=1612391676&raw=0 armel https://buildd.debian.org/status/fetch.php?pkg=basix&arch=armel&ver=0.0.1%7Egit20210122.4f10ef2-1&stamp=1612398463&raw=0 armhf https://buildd.debian.org/status/fetch.php?pkg=basix&arch=armhf&ver=0.0.1%7Egit20210122.4f10ef2-1&stamp=1612392128&raw=0 Is this problem...
- Currently basix only provides the interpolation matrix which has shape ` shape = [Npoints, ndofs * value_size]`. The interpolation data for any v in value_size is not contiguous, and...
These are a low-order, H1-conforming pair that discretise the stokes complex https://arxiv.org/abs/1710.08044 The pressure space is easy, but the velocity space is rather fiddly, it's P1 + normal components of...
The number of floating point operations per cell is currently `O((P0+1)^d *(P1+1)^d) =~ O(P^2*d)` where P0 is the degree of interpolation from and P1 the degree for interpolation to. If...
Currently all the tests use the Python interface (via wrapper.cpp). It would be good to add tests that directly interface through cpp so we can test this without having to...
Curently there is `tabulate` that returns a matrix, and `tabulate_m` that directly acts on memory https://github.com/FEniCS/basix/blob/e2137c874f8640988cb71b2aa347829a102dbe8b/cpp/basix/finite-element.h#L267-L316 The tabulate returning a matrix should be removed, then the memory version can be...
## GLL warped These points on a Pyramid are currently causing tests to intermittently fail (see #269) so have been disabled. ## GLL Isaac In #267, Isaac's GLL points on...