We need to implement the ability to read the matrix from the format that I have been developing, for outputting the Jacobian from SOD2D.
This consists of a hdf containing (for now) 4 arrays, for a matrix in CSC format:
- col_ptr (number of non-zero entries in each column)
- row_idxs (the row indices of the non-zero entries of a given column)
- vals (values of those nonzero entries)
- work_node_coords (x, y, z coordinates of the mesh's work nodes)
The first three arrays allow the matrix itself to be read, while the last allows the eigenmodes of the Jacobian to be mapped to the original problem mesh.
We need to implement the ability to read the matrix from the format that I have been developing, for outputting the Jacobian from SOD2D.
This consists of a hdf containing (for now) 4 arrays, for a matrix in CSC format:
The first three arrays allow the matrix itself to be read, while the last allows the eigenmodes of the Jacobian to be mapped to the original problem mesh.