Я пытаюсь выразить структуру разреженной матрицы (и структуры, возникающей в результате операций редкой матрицы) во время компиляции с параметрами шаблона и функциями contexpr. < /p>
Я определил эту структуру, чтобы представить структуру матрицы CSC. < /p>
Код: Выделить всё
template
struct CscStructure {
// ...
};
< /code>
Тогда я хотел бы иметь возможность иметь операции contexpr, такие как умножение матрицы-вектора. Вычисление элементов NNZ, а затем индексы строки, например ,::
constexpr size_t nnz = get_mat_vec_nnz();
constexpr std::array row_indices = get_mat_vec_row_indices();
...
CscStructure mat_vec_structure;
Однако существует значительное дублирование между get_mat_vec_nnz и get_mat_vec_row_indices Потому что первая функция выполняет все одинаковые работы, что и вторая, чтобы вычислить NNZ, просто без записи индексов. Например, кажется, нельзя использовать std :: vector, потому что nnz = vector_row_indices.size () не удовлетворит требования для шаблона.
Подробнее здесь:
https://stackoverflow.com/questions/796 ... a-constexp