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

fixed a bug resulting in too high results

Browse Source
This commit is contained in:
philippun 2023-08-09 17:33:58 +02:00
parent 18abb0fd0e
commit 3cce7fb357
2 changed files with 132 additions and 160158 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
examples/second_example.cpp
View File

@ -17,8 +17,9 @@ int main(int argc, char *argv[]) {
vector<double> alpha(n * m, 1); vector<double> alpha(n * m, 1);
vector<double> field(n * m, 0); vector<double> field(n * m, 0);
field[n * 19] = 2000; field[0] = 2000;
field[n * 19 + 19] = 2000; // field[n * 19] = 2000;
// field[n * 19 + 19] = 2000;
// for (int i = 1; i<20; i++) { // for (int i = 1; i<20; i++) {
// for (int j = 0; j<20; j++ ) { // for (int j = 0; j<20; j++ ) {
@ -42,7 +43,7 @@ int main(int argc, char *argv[]) {
input_param.setDomainSize(n, m); input_param.setDomainSize(n, m);
BoundaryCondition bc(n, m); BoundaryCondition bc(n, m);
boundary_condition bc_constant; boundary_condition bc_constant;
bc_constant.type = BC_TYPE_CLOSED; bc_constant.type = BC_TYPE_CONSTANT;
bc_constant.value = 0; bc_constant.value = 0;
bc.setSide(BC_SIDE_LEFT, bc_constant); bc.setSide(BC_SIDE_LEFT, bc_constant);
bc.setSide(BC_SIDE_RIGHT, bc_constant); bc.setSide(BC_SIDE_RIGHT, bc_constant);
@ -73,7 +74,7 @@ int main(int argc, char *argv[]) {
exit(1); exit(1);
} }
int iterations = 1000; int iterations = 100;
for (int t = 0; t < iterations; t++) { for (int t = 0; t < iterations; t++) {
double result = ADI_2D(input_param, &field[0], &alpha[0]); double result = ADI_2D(input_param, &field[0], &alpha[0]);

160281
proto/BTCS.ipynb
View File

File diff suppressed because one or more lines are too long