max/TugJulia
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

change values of bc @ A and b

Browse Source
This commit is contained in:
Max Luebke 2021-12-06 20:18:33 +01:00
parent 76640da6cb
commit 0a5c271a4b
1 changed files with 5 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
src/BTCSDiffusion.cpp
View File

@ -106,8 +106,11 @@ void BTCSDiffusion::simulate1D(std::vector<double> &c, double bc_left,
A_matrix.resize(size, size); A_matrix.resize(size, size);
A_matrix.reserve(Eigen::VectorXi::Constant(size, 3)); A_matrix.reserve(Eigen::VectorXi::Constant(size, 3));
A_matrix.insert(0, 0) = bc_left; A_matrix.insert(0, 0) = 1;
A_matrix.insert(size - 1, size - 1) = bc_right; A_matrix.insert(size - 1, size - 1) = 1;
b_vector[0] = bc_left;
b_vector[size-1] = bc_right;
for (int i = 1; i < this->dim_x + 1; i++) { for (int i = 1; i < this->dim_x + 1; i++) {
double sx = (alpha[i - 1] * time_step) / (dx * dx); double sx = (alpha[i - 1] * time_step) / (dx * dx);