115     create_h1_basis(c::Ceed, topo::Topology, ncomp, nnodes, nqpts, interp, grad, qref, qweight)
124 - `nqpts`:   Total number of quadrature points.
125 - `interp`:  Matrix of size `(nqpts, nnodes)` expressing the values of nodal basis functions
127 - `grad`:    Array of size `(dim, nqpts, nnodes)` expressing derivatives of nodal basis
129 - `qref`:    Matrix of size `(dim, nqpts)` holding the locations of quadrature points on the
131 - `qweight`: Array of length `nqpts` holding the quadrature weights on the reference
139     nqpts,  argument
146     @assert size(interp) == (nqpts, nnodes)
147     @assert size(grad) == (dim, nqpts, nnodes)
148     @assert size(qref) == (dim, nqpts)
149     @assert length(qweight) == nqpts
162         nqpts,
173     create_hdiv_basis(c::Ceed, topo::Topology, ncomp, nnodes, nqpts, interp, div, qref, qweight)
182 - `nqpts`:   Total number of quadrature points.
183 - `interp`:  Matrix of size `(dim, nqpts, nnodes)` expressing the values of basis functions
185 - `div`:     Array of size `(nqpts, nnodes)` expressing divergence of basis functions at
187 - `qref`:    Matrix of size `(dim, nqpts)` holding the locations of quadrature points on the
189 - `qweight`: Array of length `nqpts` holding the quadrature weights on the reference
197     nqpts,  argument
204     @assert size(interp) == (dim, nqpts, nnodes)
205     @assert size(div) == (nqpts, nnodes)
206     @assert size(qref) == (dim, nqpts)
207     @assert length(qweight) == nqpts
220         nqpts,
231     create_hcurl_basis(c::Ceed, topo::Topology, ncomp, nnodes, nqpts, interp, curl, qref, qweight)
240 - `nqpts`:   Total number of quadrature points.
241 - `interp`:  Matrix of size `(dim, nqpts, nnodes)` expressing the values of basis functions
243 - `curl`:    Matrix of size `(curlcomp, nqpts, nnodes)`, `curlcomp = 1 if dim < 3 else dim`)
245 - `qref`:    Matrix of size `(dim, nqpts)` holding the locations of quadrature points on the
247 - `qweight`: Array of length `nqpts` holding the quadrature weights on the reference
255     nqpts,  argument
263     @assert size(interp) == (dim, nqpts, nnodes)
264     @assert size(curl) == (curlcomp, nqpts, nnodes)
265     @assert size(qref) == (dim, nqpts)
266     @assert length(qweight) == nqpts
279         nqpts,
414     nqpts = Ref{CeedInt}()
415     C.CeedBasisGetNumQuadraturePoints(b[], nqpts)
416     nqpts[]