Todo List

Subprogram bc::bc_finalize_base (this, only_facets)

add 2D case

Subprogram coefs::coef_generate_mass (c)

cleanup once we have device math in place

Subprogram coefs::coef_init_all (this, gs_h)

Be clever and try to avoid allocating zeroed geom. factors

Subprogram dofmap::dofmap_init (this, msh, xh)

implement for 2d elements

Subprogram dofmap::dofmap_number_faces (this)

don't assume lx = ly = lz

Subprogram fld_file::fld_file_write (this, data, t)

fix support for single precision output?

Subprogram fld_file::fld_file_write_masked (this, data, mask, t)

fix support for single precision output?

Subprogram fluid_plan1::fluid_plan1_t::init (this, msh, lx, param)

Remaning plan1 related data, ax, precon etc

Subprogram fluid_scheme_incompressible::fluid_scheme_bc_apply_prs (this, time)

Why can't we call the interface here?

Subprogram fluid_scheme_incompressible::fluid_scheme_bc_apply_vel (this, time, strong)

Why can't we call the interface here?

Subprogram glb_intrp_comm::glb_intrp_comm_sendrecv (this, send, recv, n_send, n_recv)

Check size etc against status

Subprogram gs_cpu::gs_scatter_cpu (this, v, m, dg, u, n, gd, nb, b, bo, shrd, event)

Make the kernel abstract

Subprogram gs_mpi::gs_nbwait_mpi (this, u, n, op, strm)

Check size etc against status

Subprogram gs_sx::gs_scatter_sx (this, v, m, dg, u, n, gd, nb, b, bo, shrd, event)

Make the kernel abstract

Subprogram htable::htable_iter_data (this, data)

Remove once we figure out how to do this with value() for arbitrary data types

Subprogram intersection_detector::intersect_detector_free (this)

cleanup the aabb tree

Subprogram intersection_detector::intersect_detector_init (this, msh, padding)

Rework this part once the new mesh structure is in place

Subprogram krylov::krylov_free (this)

add calls to destroy precon. if necessary

Module map

Figure out a better name for this module

Subprogram map_file::map_file_read (this, data)

Check if neli matches mapnel

Subprogram mesh::mesh_generate_conn (this)

These loop nests needs a lot of love...

Subprogram mesh::mesh_generate_edge_conn (this)

move this into distdata

Subprogram mesh::mesh_generate_facet_numbering (this)

Since we now the neigh. we can actually do p2p here...

Move this into distdata as a method...

move this into distdata

Subprogram mesh::mesh_get_local_point (this, p)

why do we still need to do this?

Subprogram mesh::mesh_have_point_glb_idx (this, index)

Consider moving this to distdata

Subprogram mesh::mesh_init_common (this)

resize onces final size is known Only init if we generate connectivity

resize onces final size is known

Subprogram mesh::mesh_t::get_global (this, e)

Add similar mappings for element ids

Subprogram mesh::mesh_t::get_local (this, p)

Add similar mappings for element ids

Struct mesh_field::mesh_fld_t

Add support for different data types

Subprogram quad::quad_facet_id (this, t, side)

sort this

Subprogram re2_file::re2_file_read (this, data)

Add support for curved side data Skip curved side data

Subprogram re2_file::re2_file_read_bcs (msh, nbcs, dist, fh, mpi_offset, v2_format)

Use element offset in parallel

Subprogram re2_file::re2_file_write (this, data, t)

Add support for curved side data

Struct rea::rea_t

add missing data fields

Subprogram speclib::endw1 (n, alpha, beta)

document ENDW1

Subprogram speclib::endw2 (n, alpha, beta)

document ENDW2

Subprogram speclib::gammaf (x)

document GAMMAF

Subprogram speclib::pnormj (n, alpha, beta)

document PNORMJ

Subprogram tensor::tensr3 (v, nv, u, nu, a, bt, ct, w)

Add 2d case

Subprogram tri::tri_facet_id (this, t, side)

sort this