"README.md" did not exist on "cab39a770e0d49e8a571f8e82c41be35b6021b2d"
Newer
Older
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
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
#include <cxxtest/TestSuite.h>
#include <cmath>
#include "WorkspaceCreationHelper.hh"
#include "MantidAlgorithms/Plus.h"
#include "MantidAPI/AnalysisDataService.h"
#include "MantidAPI/WorkspaceFactory.h"
#include "MantidDataObjects/Workspace2D.h"
#include "MantidDataObjects/Workspace1D.h"
#include "MantidAPI/WorkspaceProperty.h"
using namespace Mantid::API;
using namespace Mantid::Kernel;
using namespace Mantid::Algorithms;
using namespace Mantid::DataObjects;
class ComplexOpTest : public Algorithm
{
public:
ComplexOpTest() : Algorithm() {}
virtual ~ComplexOpTest() {}
void init()
{
declareProperty(new WorkspaceProperty<Workspace>("InputWorkspace_1","",Direction::Input));
declareProperty(new WorkspaceProperty<Workspace>("InputWorkspace_2","",Direction::Input));
declareProperty(new WorkspaceProperty<Workspace>("OutputWorkspace","",Direction::Output));
}
void exec()
{
Workspace_sptr in_work1 = getProperty("InputWorkspace_1");
Workspace_sptr in_work2 = getProperty("InputWorkspace_2");
Workspace_sptr out_work = (in_work1 + in_work2)/(in_work1 - in_work2);
setProperty("OutputWorkspace",out_work);
}
};
class ChainedOperatorTest : public CxxTest::TestSuite
{
public:
void testChainedOperator()
{
int sizex = 10,sizey=20;
// Register the workspace in the data service
AnalysisDataService* ADS = AnalysisDataService::Instance();
Workspace_sptr work_in1 = WorkspaceCreationHelper::Create2DWorkspace123(sizex,sizey);
Workspace_sptr work_in2 = WorkspaceCreationHelper::Create2DWorkspace154(sizex,sizey);
ComplexOpTest alg;
std::string wsNameIn1 = "testChainedOperator_in21";
std::string wsNameIn2 = "testChainedOperator_in22";
std::string wsNameOut = "testChainedOperator_out";
ADS->add(wsNameIn1, work_in1);
ADS->add(wsNameIn2, work_in2);
alg.initialize();
alg.setPropertyValue("InputWorkspace_1",wsNameIn1);
alg.setPropertyValue("InputWorkspace_2",wsNameIn2);
alg.setPropertyValue("OutputWorkspace",wsNameOut);
TS_ASSERT_THROWS_NOTHING(alg.execute());
TS_ASSERT( alg.isExecuted() );
Workspace_sptr work_out1;
TS_ASSERT_THROWS_NOTHING(work_out1 = ADS->retrieve(wsNameOut));
checkData(work_in1, work_in2, work_out1);
ADS->remove(wsNameIn1);
ADS->remove(wsNameIn2);
ADS->remove(wsNameOut);
}
private:
void checkData( Workspace_sptr work_in1, Workspace_sptr work_in2, Workspace_sptr work_out1)
{
int ws2LoopCount;
if (work_in2->size() > 0)
{
ws2LoopCount = work_in1->size()/work_in2->size();
}
ws2LoopCount = (ws2LoopCount==0) ? 1 : ws2LoopCount;
for (int i = 0; i < work_out1->size(); i++)
{
checkDataItem(work_in1,work_in2,work_out1,i,i/ws2LoopCount);
}
}
void checkDataItem (Workspace_sptr work_in1, Workspace_sptr work_in2, Workspace_sptr work_out1, int i, int ws2Index)
{
double sig1 = work_in1->dataY(i/work_in1->blocksize())[i%work_in1->blocksize()];
double sig2 = work_in2->dataY(ws2Index/work_in1->blocksize())[ws2Index%work_in1->blocksize()];
double sig3 = work_out1->dataY(i/work_in1->blocksize())[i%work_in1->blocksize()];
TS_ASSERT_DELTA((sig1 + sig2)/(sig1-sig2), sig3, 0.0001);
//Note err calculation not checked due to complexity.
}
};