mirror of
https://git.gfz-potsdam.de/naaice/tug.git
synced 2025-12-13 17:38:23 +01:00
cleanup code
This commit is contained in:
parent
96a2d1cc5b
commit
85278bcaff
@ -80,12 +80,9 @@ void BTCSDiffusion::simulate(std::vector<double> &c, std::vector<double> &alpha,
|
||||
else
|
||||
b[A_line] = this->bc[1];
|
||||
|
||||
// std::cout << b << std::endl;
|
||||
|
||||
Eigen::SparseMatrix<double> A(size, size);
|
||||
A.setFromTriplets(tripletList.begin(), tripletList.end());
|
||||
|
||||
// std::cout << A << std::endl;
|
||||
Eigen::SparseQR<Eigen::SparseMatrix<double>, Eigen::COLAMDOrdering<int>>
|
||||
solver;
|
||||
|
||||
|
||||
@ -15,8 +15,6 @@ int main(int argc, char *argv[]) {
|
||||
|
||||
bc_left.push_back(5. * std::pow(10, -6));
|
||||
bc_right.push_back(-1);
|
||||
// input[x + 2] = 5.5556554 * std::pow(10, -6);
|
||||
// input[x + 3] = 5.234564213 * std::pow(10, -6);
|
||||
|
||||
BTCSDiffusion diffu(x);
|
||||
|
||||
@ -25,7 +23,6 @@ int main(int argc, char *argv[]) {
|
||||
|
||||
for (int i = 0; i < 100; i++) {
|
||||
diffu.simulate(input, alpha, 1.);
|
||||
// BTCS1D(x, input, alpha, 1., bc);
|
||||
}
|
||||
|
||||
return 0;
|
||||
|
||||
Loading…
x
Reference in New Issue
Block a user