1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
|
#!/usr/bin/env python3
# encoding: utf-8
"""
_Diagnostics.py
Created by Graham Dennis on 2009-06-15.
Copyright (c) 2009-2012, Graham Dennis
This program is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 2 of the License, or
(at your option) any later version.
This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
You should have received a copy of the GNU General Public License
along with this program. If not, see <http://www.gnu.org/licenses/>.
"""
from xpdeint.Features._Feature import _Feature
from xpdeint.Geometry.UniformDimensionRepresentation import UniformDimensionRepresentation
from xpdeint.Function import Function
class _Diagnostics (_Feature):
def nonlocalAccess(self, dict):
"""
The purpose of this function is to safety-check nonlocal dimension access. The only place where this is potentially
unsafe is with integer-valued dimensions as we permit users to use any code to access a dimension. For speed reasons
we don't run safety checks. But in the diagnostics feature, we can sacrifice speed for more safety.
For any nonlocally-accessed dimReps, we replace the usual #define with a function call that includes the line number
from which it was called from the original .xmds file. The function then does a bounds check on all potentially-unsafe
accesses and stops the simulation if any are found.
"""
availableDimReps = dict['availableDimReps']
dimRepsNeeded = dict['dimRepsNeeded']
if not any(dimRep.type == 'long' and isinstance(dimRep, UniformDimensionRepresentation) for dimRep in dimRepsNeeded):
return
vector = dict['vector']
nonlocalAccessVariableName = dict['nonlocalAccessVariableName']
nonlocalAccessString = dict['nonlocalAccessString']
componentName = dict['componentName']
args = [('long', dimRep.loopIndex) for dimRep in availableDimReps]
args.extend([('const char*', 'filename'), ('int', 'line_number')])
implementation = lambda func: self.nonlocalAccessValidationFunctionContents(dimRepsNeeded, nonlocalAccessString, componentName, func)
self.functions[nonlocalAccessVariableName] = Function(
'_' + nonlocalAccessVariableName,
args,
implementation,
returnType = 'inline ' + vector.type + '&',
)
defineArgumentsString = ', '.join(dimRep.loopIndex for dimRep in dimRepsNeeded)
functionArgumentsString = ', '.join(dimRep.loopIndex for dimRep in availableDimReps)
defineString = '#define %(nonlocalAccessVariableName)s(%(defineArgumentsString)s) _%(nonlocalAccessVariableName)s(%(functionArgumentsString)s, __FILE__, __LINE__)\n' % locals()
dict['defineString'] = defineString
|
