mulimoen/SummationByParts
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

divide instead of multiply h

Browse Source
This commit is contained in:
Magnus Ulimoen
2020-04-10 01:07:57 +02:00
parent a51ee210fe
commit 147f3040c3
1 changed files with 4 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
sbp/src/euler.rs
View File

@@ -630,7 +630,7 @@ fn SAT_characteristics<SBP: SbpOperator>(
) {
// North boundary
{
let hi = (k.ny() - 1) as Float * SBP::h()[0];
let hi = (k.ny() - 1) as Float / SBP::h()[0];
let sign = -1.0;
let tau = 1.0;
let slice = s![y.ny() - 1, ..];
@@ -648,7 +648,7 @@ fn SAT_characteristics<SBP: SbpOperator>(
}
// South boundary
{
let hi = (k.ny() - 1) as Float * SBP::h()[0];
let hi = (k.ny() - 1) as Float / SBP::h()[0];
let sign = 1.0;
let tau = -1.0;
let slice = s![0, ..];
@@ -666,7 +666,7 @@ fn SAT_characteristics<SBP: SbpOperator>(
}
// West Boundary
{
let hi = (k.nx() - 1) as Float * SBP::h()[0];
let hi = (k.nx() - 1) as Float / SBP::h()[0];
let sign = 1.0;
let tau = -1.0;
let slice = s![.., 0];
@@ -684,7 +684,7 @@ fn SAT_characteristics<SBP: SbpOperator>(
}
// East Boundary
{
let hi = (k.nx() - 1) as Float * SBP::h()[0];
let hi = (k.nx() - 1) as Float / SBP::h()[0];
let sign = -1.0;
let tau = 1.0;
let slice = s![.., y.nx() - 1];