Compare commits

...

4 Commits

Author SHA1 Message Date
d5fd78d493 Revert "Add restrict keyword to pointers (not sharing the same memory for multiple fields)"
This reverts commit 60b70bb332.

It appears to have minimal performance impact, and fails to compile on
nvidia cards.
2019-07-17 00:55:28 -07:00
fe2f4f4510 style fixes 2019-07-15 00:06:43 -07:00
314e36a3cc fix for 1-sized grids 2019-07-15 00:06:43 -07:00
a37df3a74f expand gitignores 2019-07-15 00:06:43 -07:00
4 changed files with 15 additions and 10 deletions

8
.gitignore vendored
View File

@ -1,5 +1,9 @@
.idea/
__pycache__
*.h5
*.pyc
__pycache__
*.py[cod]
build/
dist/
*.egg-info/

View File

@ -78,7 +78,8 @@ int p{{r}} = + (int) di{{r}};
int wrap_{{r}} = (s{{r}} - 1) * (int) di{{r}};
if ( {{r}} == 0 ) {
m{{r}} = wrap_{{r}};
} else if ( {{r}} == s{{r}} - 1 ) {
}
if ( {{r}} == s{{r}} - 1 ) {
p{{r}} = -wrap_{{r}};
}
{% endfor %}

View File

@ -170,7 +170,7 @@ class Simulation(object):
ctype = type_to_C(self.arg_type)
def ptr(arg: str) -> str:
return ctype + ' * restrict ' + arg
return ctype + ' *' + arg
base_fields = OrderedDict()
base_fields[ptr('E')] = self.E
@ -286,7 +286,7 @@ class Simulation(object):
alpha_max = pml['cfs_alpha']
def par(x):
scaling = ((x / (pml['thickness'])) ** pml['m'])
scaling = (x / pml['thickness']) ** pml['m']
sigma = scaling * sigma_max
kappa = 1 + scaling * (kappa_max - 1)
alpha = ((1 - x / pml['thickness']) ** pml['ma']) * alpha_max