33 coorOffset=-coorOffset;
35 else if (
Size(coorOffset) >=
dim_[i]) {
36 coorOffset = 2*(
dim_[i]-1) - coorOffset;
38 return myIndex + coorOffset*
spacing_[i];
50 if (coorOffset1 < 0) {
51 coorOffset1=-coorOffset1;
53 else if (
Size(coorOffset1) >=
dim_[i1]) {
54 coorOffset1 = 2*(
dim_[i1]-1) - coorOffset1;
58 if (coorOffset2 < 0) {
59 coorOffset2=-coorOffset2;
61 else if (
Size(coorOffset2) >=
dim_[i2]) {
62 coorOffset2 = 2*(
dim_[i2]-1) - coorOffset2;
72 std::vector<Size> coordinates = iterator.
coordinates();
76 coorOffset=-coorOffset;
78 else if (
Size(coorOffset) >=
dim_[i]) {
79 coorOffset = 2*(
dim_[i]-1) - coorOffset;
81 coordinates[i] =
Size(coorOffset);
const std::vector< Size > & coordinates() const
FdmLinearOpIterator iter_neighbourhood(const FdmLinearOpIterator &iterator, Size i, Integer offset) const
std::vector< Size > spacing_
Size neighbourhood(const FdmLinearOpIterator &iterator, Size i, Integer offset) const
Size index(const std::vector< Size > &coordinates) const
memory layout of a fdm linear operator
QL_INTEGER Integer
integer number
std::size_t Size
size of a container