... | @@ -320,11 +320,7 @@ void SWE::SWE_ADERDG::boundaryValues(const double* const x,const double t,const |
... | @@ -320,11 +320,7 @@ void SWE::SWE_ADERDG::boundaryValues(const double* const x,const double t,const |
|
out.p(0) = in.p(0);
|
|
out.p(0) = in.p(0);
|
|
out.p(1) = in.p(1);
|
|
out.p(1) = in.p(1);
|
|
out.p(direction)= -in.p(direction);
|
|
out.p(direction)= -in.p(direction);
|
|
|
|
\\....
|
|
fluxout.h() = fluxin.h();
|
|
|
|
fluxout.p(0) = fluxin.p(0);
|
|
|
|
fluxout.p(1) = fluxin.p(1);
|
|
|
|
}
|
|
|
|
```
|
|
```
|
|
|
|
|
|
To get the new flux we apply the inverted velocity vector to the flux term (in this example for normalNonZero=0):
|
|
To get the new flux we apply the inverted velocity vector to the flux term (in this example for normalNonZero=0):
|
... | @@ -345,9 +341,9 @@ which results in the implementation |
... | @@ -345,9 +341,9 @@ which results in the implementation |
|
ReadOnlyVariables f_in (fluxIn) ;
|
|
ReadOnlyVariables f_in (fluxIn) ;
|
|
Variables f_out(fluxOut);
|
|
Variables f_out(fluxOut);
|
|
|
|
|
|
f_out.h() = -f_in.h();
|
|
f_out.h() = f_in.h();
|
|
f_out.p(normalNonZero) = f_in.p(normalNonZero);
|
|
f_out.p(0) = f_in.p(0);
|
|
f_out.p(normalZero) = -f_in.p(normalZero);
|
|
f_out.p(1) = f_in.p(1);
|
|
}
|
|
}
|
|
```
|
|
```
|
|
|
|
|
... | | ... | |