144{
145 if constexpr (Rows == 1 || Cols == 1)
146 {
147 std::vector<Scalar>
values(matrix.data(), matrix.data() + matrix.rows() * matrix.cols());
148 nlohmann::to_json(json,
values);
149 }
150 else
151 {
152 const auto nrows = matrix.rows();
153 const auto ncols = matrix.cols();
154 std::vector<std::vector<Scalar>>
values(nrows, std::vector<Scalar>(ncols));
155 for (const auto i : make_range(nrows))
156 for (const auto j : make_range(ncols))
157 values[i][j] = matrix(i, j);
158 nlohmann::to_json(json,
values);
159 }
160}
std::array< Real, 2 > values
for(PetscInt i=0;i< nvars;++i)