Skip to content
GitLab
Menu
Projects
Groups
Snippets
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Sign in / Register
Toggle navigation
Menu
Open sidebar
Julia
Nektar
Commits
cce1dc18
Commit
cce1dc18
authored
Jan 27, 2014
by
Chris Cantwell
Browse files
Fixed code formatting.
parent
7eadb1d2
Changes
5
Hide whitespace changes
Inline
Side-by-side
library/LocalRegions/HexExp.cpp
View file @
cce1dc18
...
...
@@ -1598,9 +1598,10 @@ namespace Nektar
break
;
case
StdRegions
::
eLaplacian
:
{
if
(
m_metricinfo
->
GetGtype
()
==
SpatialDomains
::
eDeformed
||
mkey
.
GetNVarCoeff
()
||
(
mkey
.
ConstFactorExists
(
StdRegions
::
eFactorSVVCutoffRatio
)))
if
(
m_metricinfo
->
GetGtype
()
==
SpatialDomains
::
eDeformed
||
mkey
.
GetNVarCoeff
()
||
mkey
.
ConstFactorExists
(
StdRegions
::
eFactorSVVCutoffRatio
))
{
NekDouble
one
=
1.0
;
DNekMatSharedPtr
mat
=
GenMatrix
(
mkey
);
...
...
library/LocalRegions/PrismExp.cpp
View file @
cce1dc18
...
...
@@ -1341,9 +1341,10 @@ namespace Nektar
case
StdRegions
::
eLaplacian
:
{
if
((
m_metricinfo
->
GetGtype
()
==
SpatialDomains
::
eDeformed
)
||
(
mkey
.
GetNVarCoeff
()
>
0
)
||
(
mkey
.
ConstFactorExists
(
StdRegions
::
eFactorSVVCutoffRatio
)))
if
(
m_metricinfo
->
GetGtype
()
==
SpatialDomains
::
eDeformed
||
mkey
.
GetNVarCoeff
()
>
0
||
mkey
.
ConstFactorExists
(
StdRegions
::
eFactorSVVCutoffRatio
))
{
NekDouble
one
=
1.0
;
DNekMatSharedPtr
mat
=
GenMatrix
(
mkey
);
...
...
library/LocalRegions/PyrExp.cpp
View file @
cce1dc18
...
...
@@ -451,9 +451,10 @@ namespace Nektar
break
;
case
StdRegions
::
eLaplacian
:
{
if
((
m_metricinfo
->
GetGtype
()
==
SpatialDomains
::
eDeformed
)
||
(
mkey
.
GetNVarCoeff
()
>
0
)
||
(
mkey
.
ConstFactorExists
(
StdRegions
::
eFactorSVVCutoffRatio
)))
if
(
m_metricinfo
->
GetGtype
()
==
SpatialDomains
::
eDeformed
||
mkey
.
GetNVarCoeff
()
>
0
||
mkey
.
ConstFactorExists
(
StdRegions
::
eFactorSVVCutoffRatio
))
{
NekDouble
one
=
1.0
;
DNekMatSharedPtr
mat
=
GenMatrix
(
mkey
);
...
...
library/StdRegions/StdExpansion3D.cpp
View file @
cce1dc18
...
...
@@ -197,8 +197,8 @@ namespace Nektar
Array
<
OneD
,
NekDouble
>
&
outarray
,
const
StdRegions
::
StdMatrixKey
&
mkey
)
{
if
((
mkey
.
GetNVarCoeff
()
==
0
)
&&
!
mkey
.
ConstFactorExists
(
eFactorSVVCutoffRatio
))
if
(
mkey
.
GetNVarCoeff
()
==
0
&&
!
mkey
.
ConstFactorExists
(
eFactorSVVCutoffRatio
))
{
// This implementation is only valid when there are no
// coefficients associated to the Laplacian operator
...
...
library/StdRegions/StdQuadExp.cpp
View file @
cce1dc18
...
...
@@ -1459,10 +1459,12 @@ namespace Nektar
{
if
(
j
+
k
>=
cutoff
)
//to filter out only the "high-modes"
{
//orthocoeffs[cnt] *= (1.0+SvvDiffCoeff*exp(-(j-nmodes)*(j-nmodes)/((NekDouble)((j-cutoff+epsilon)*(j-cutoff+epsilon))))*exp(-(k-nmodes)*(k-nmodes)/((NekDouble)((k-cutoff+epsilon)*(k-cutoff+epsilon)))));
orthocoeffs
[
j
*
nmodes_b
+
k
]
*=
(
1.0
+
SvvDiffCoeff
*
exp
(
-
(
j
+
k
-
nmodes
)
*
(
j
+
k
-
nmodes
)
/
((
NekDouble
)((
j
+
k
-
cutoff
+
1
)
*
(
j
+
k
-
cutoff
+
1
)))));
orthocoeffs
[
j
*
nmodes_b
+
k
]
*=
(
1.0
+
SvvDiffCoeff
*
exp
(
-
(
j
+
k
-
nmodes
)
*
(
j
+
k
-
nmodes
)
/
((
NekDouble
)((
j
+
k
-
cutoff
+
1
)
*
(
j
+
k
-
cutoff
+
1
)))));
}
cnt
++
;
cnt
++
;
}
}
...
...
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment