1751 |
|
|
1752 |
// check the results |
// check the results |
1753 |
tmp = std::sin((double)p); |
tmp = std::sin((double)p); |
1754 |
assert(dataView()==tmp); |
assert(std::abs(dataView()-tmp)<1.e-10); |
1755 |
|
|
1756 |
if (p<npoints-1) { |
if (p<npoints-1) { |
1757 |
dataView.incrOffset(); |
dataView.incrOffset(); |
1845 |
for (int j=0;j<shape[1];j++) { |
for (int j=0;j<shape[1];j++) { |
1846 |
for (int k=0;k<shape[2];k++) { |
for (int k=0;k<shape[2];k++) { |
1847 |
for (int l=0;l<shape[3];l++) { |
for (int l=0;l<shape[3];l++) { |
1848 |
assert(dataView(i,j,k,l)==std::log((double)dataView.index(i,j,k,l))); |
assert(std::abs(dataView(i,j,k,l)-std::log((double)dataView.index(i,j,k,l)))<1.e-10); |
1849 |
} |
} |
1850 |
} |
} |
1851 |
} |
} |