max/tug
1
0
Fork 0
mirror of https://git.gfz-potsdam.de/naaice/tug.git synced 2025-12-13 17:38:23 +01:00
Code Issues Packages Projects Releases Wiki Activity

Apply new scheme to model (only 1D)

Browse Source
This commit is contained in:
Max Luebke 2022-04-20 09:55:38 +02:00
parent 57cfa682c5
commit e9a1d06784
1 changed files with 11 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

23
src/BTCSDiffusion.cpp
View File

@ -239,8 +239,8 @@ void Diffusion::BTCSDiffusion::fillMatrixFromRow(
A_matrix.insert(1, 1) = 1; A_matrix.insert(1, 1) = 1;
} else { } else {
double sx = (alpha[0] * time_step) / (dx * dx); double sx = (alpha[0] * time_step) / (dx * dx);
A_matrix.insert(1, 1) = -1. - 2. * sx; A_matrix.insert(1, 1) = -1. - 3. * sx;
A_matrix.insert(1, 0) = sx; A_matrix.insert(1, 0) = 2. * sx;
A_matrix.insert(1, 2) = sx; A_matrix.insert(1, 2) = sx;
} }
@ -261,9 +261,9 @@ void Diffusion::BTCSDiffusion::fillMatrixFromRow(
A_matrix.insert(A_size - 2, A_size - 2) = 1; A_matrix.insert(A_size - 2, A_size - 2) = 1;
} else { } else {
double sx = (alpha[size - 1] * time_step) / (dx * dx); double sx = (alpha[size - 1] * time_step) / (dx * dx);
A_matrix.insert(A_size - 2, A_size - 2) = -1. - 2. * 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;
A_matrix.insert(A_size - 2, A_size - 1) = sx; A_matrix.insert(A_size - 2, A_size - 1) = 2. * sx;
} }
A_matrix.insert(A_size - 1, A_size - 1) = 1; A_matrix.insert(A_size - 1, A_size - 1) = 1;
@ -294,15 +294,14 @@ void Diffusion::BTCSDiffusion::fillVectorFromRow(
b_vector[j + 1] = -c[j] - t0_c_j; b_vector[j + 1] = -c[j] - t0_c_j;
} }
if (!left_constant) { // this is not correct currently.We will fix this when we are able to define
// this is not correct currently.We will fix this when we are able to define // FLUX boundary conditions
// FLUX boundary conditions b_vector[0] =
b_vector[0] = getBCFromFlux(left, c[0], alpha[0]); (left_constant ? left.value : getBCFromFlux(left, c[0], alpha[0]));
}
if (!right_constant) { b_vector[b_size - 1] =
b_vector[b_size - 1] = getBCFromFlux(right, c[size - 1], alpha[size - 1]); (right_constant ? right.value
} : getBCFromFlux(right, c[size - 1], alpha[size - 1]));
} }
void Diffusion::BTCSDiffusion::setTimestep(double time_step) { void Diffusion::BTCSDiffusion::setTimestep(double time_step) {