Skip to content
Snippets Groups Projects
Commit 1ed86f65 authored by Dmitri Naumov's avatar Dmitri Naumov
Browse files

[MaL] Abort if LisMatrix is not assembled.

parent a1fb2aee
No related branches found
No related tags found
No related merge requests found
......@@ -15,7 +15,9 @@
#include "LisMatrix.h"
#include <cmath>
#include <cassert>
#include <cstdlib>
#include <logog/include/logog.hpp>
#include "LisVector.h"
#include "LisCheck.h"
......@@ -79,7 +81,11 @@ int LisMatrix::add(std::size_t rowId, std::size_t colId, double v)
void LisMatrix::write(const std::string &filename) const
{
assert(_is_assembled);
if (!_is_assembled)
{
ERR("LisMatrix::write(): matrix not assembled.");
std::abort();
}
lis_output_matrix(_AA, LIS_FMT_MM, const_cast<char*>(filename.c_str()));
}
......@@ -103,7 +109,11 @@ double LisMatrix::getMaxDiagCoeff()
void LisMatrix::multiply(const LisVector &x, LisVector &y) const
{
assert(_is_assembled);
if (!_is_assembled)
{
ERR("LisMatrix::multiply(): matrix not assembled.");
std::abort();
}
int ierr = lis_matvec(_AA, const_cast<LisVector*>(&x)->getRawVector(), y.getRawVector());
checkLisError(ierr);
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment