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

commentary

Browse Source
This commit is contained in:
philippun 2023-06-08 15:58:30 +02:00
parent 1df7f82553
commit ecd3e95477
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
src/BTCS.cpp
View File

@ -138,7 +138,7 @@ auto fillMatrixFromRow(const DVectorRowMajor &alpha,
"BC_TYPE_CONSTANT are currently not supported."); "BC_TYPE_CONSTANT are currently not supported.");
} }
A_matrix.insert(1, 1) = 1; A_matrix.insert(1, 1) = 1;
} else { } else { // if BC_UNSET
sx = (alpha[0] * time_step) / (dx * dx); sx = (alpha[0] * time_step) / (dx * dx);
A_matrix.insert(1, 1) = -1. - 3. * sx; A_matrix.insert(1, 1) = -1. - 3. * sx;
A_matrix.insert(1, 0) = 2. * sx; A_matrix.insert(1, 0) = 2. * sx;
@ -153,7 +153,7 @@ auto fillMatrixFromRow(const DVectorRowMajor &alpha,
} }
A_matrix.insert(j, j) = 1; A_matrix.insert(j, j) = 1;
continue; continue;
} } // if BC_UNSET
sx = (alpha[k] * time_step) / (dx * dx); sx = (alpha[k] * time_step) / (dx * dx);
A_matrix.insert(j, j) = -1. - 2. * sx; A_matrix.insert(j, j) = -1. - 2. * sx;
@ -167,7 +167,7 @@ auto fillMatrixFromRow(const DVectorRowMajor &alpha,
"BC_TYPE_CONSTANT are currently not supported."); "BC_TYPE_CONSTANT are currently not supported.");
} }
A_matrix.insert(A_size - 2, A_size - 2) = 1; A_matrix.insert(A_size - 2, A_size - 2) = 1;
} else { } else { // if BC_UNSET
sx = (alpha[size - 1] * time_step) / (dx * dx); sx = (alpha[size - 1] * time_step) / (dx * dx);
A_matrix.insert(A_size - 2, A_size - 2) = -1. - 3. * sx; A_matrix.insert(A_size - 2, A_size - 2) = -1. - 3. * sx;
A_matrix.insert(A_size - 2, A_size - 3) = sx; A_matrix.insert(A_size - 2, A_size - 3) = sx;
@ -203,7 +203,7 @@ auto fillVectorFromRow(const DVectorRowMajor &c, const DVectorRowMajor &alpha,
} }
b_vector[j + 1] = bc_inner[j].value; b_vector[j + 1] = bc_inner[j].value;
continue; continue;
} } // if BC_UNSET
b_vector[j + 1] = -c[j] + d_ortho[j]; b_vector[j + 1] = -c[j] + d_ortho[j];
} }