Newer
Older
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
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
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
#ifndef POWERTEST_H_
#define POWERTEST_H_
#include <cxxtest/TestSuite.h>
#include <cmath>
#include "WorkspaceCreationHelper.hh"
#include "MantidAlgorithms/Power.h"
#include "MantidAPI/AnalysisDataService.h"
#include "MantidAPI/WorkspaceFactory.h"
#include "MantidDataObjects/Workspace2D.h"
#include "MantidDataObjects/Workspace1D.h"
#include "MantidAPI/WorkspaceProperty.h"
#include "MantidAPI/WorkspaceOpOverloads.h"
using namespace Mantid::API;
using namespace Mantid::Kernel;
using namespace Mantid::Algorithms;
using namespace Mantid::DataObjects;
class MultiplyTest : public CxxTest::TestSuite
{
public:
void testInit()
{
Power alg;
TS_ASSERT_THROWS_NOTHING(alg.initialize());
TS_ASSERT(alg.isInitialized());
//Setting properties to input workspaces that don't exist throws
TSM_ASSERT_THROWS("Base must be numeric", alg.setPropertyValue("LHSWorkspace","n"), std::invalid_argument )
TSM_ASSERT_THROWS("Base must be numeric", alg.setPropertyValue("RHSWorkspace","n"), std::invalid_argument )
}
void testPowerCalculation()
{
WorkspaceSingleValue_sptr baseWs = WorkspaceCreationHelper::CreateWorkspaceSingleValue(2);
WorkspaceSingleValue_sptr exponentWs = WorkspaceCreationHelper::CreateWorkspaceSingleValue(2);
AnalysisDataService::Instance().add("baseWsp",baseWs);
AnalysisDataService::Instance().add("exponentWsp",exponentWs);
Power alg;
alg.initialize();
alg.setPropertyValue("OutputWorkspace","OutputWsp");
alg.setPropertyValue("RHSWorkspace", "exponentWsp");
alg.setPropertyValue("LHSWorkspace", "baseWsp");
alg.execute();
WorkspaceSingleValue_sptr output = boost::dynamic_pointer_cast<WorkspaceSingleValue>(AnalysisDataService::Instance().retrieve("OutputWsp"));
Mantid::MantidVec expectedValue(1,4);
TSM_ASSERT_EQUALS("Power has not been determined correctly", expectedValue, output->dataY());
AnalysisDataService::Instance().remove("exponentWs");
AnalysisDataService::Instance().remove("baseWs");
AnalysisDataService::Instance().remove("OutputWsp");
}
void testSignedExponent()
{
WorkspaceSingleValue_sptr baseWs = WorkspaceCreationHelper::CreateWorkspaceSingleValue(2);
WorkspaceSingleValue_sptr exponentWs = WorkspaceCreationHelper::CreateWorkspaceSingleValue(-2);
AnalysisDataService::Instance().add("baseWs",baseWs);
AnalysisDataService::Instance().add("exponentWs",exponentWs);
Power alg;
alg.initialize();
alg.setPropertyValue("OutputWorkspace","OutputWs");
alg.setPropertyValue("RHSWorkspace", "exponentWs");
alg.setPropertyValue("LHSWorkspace", "baseWs");
alg.execute();
TSM_ASSERT_EQUALS("Algorithm should not run with a negative exponent", alg.isExecuted(),false);
AnalysisDataService::Instance().remove("exponentWs");
AnalysisDataService::Instance().remove("baseWs");
AnalysisDataService::Instance().remove("OutputWs");
}
void testMultivaluedWorkspaceAsExponent()
{
WorkspaceSingleValue_sptr baseWs = WorkspaceCreationHelper::CreateWorkspaceSingleValue(2);
Workspace1D_sptr exponentWs = WorkspaceCreationHelper::Create1DWorkspaceRand(10);
AnalysisDataService::Instance().add("baseWs",baseWs);
AnalysisDataService::Instance().add("exponentWs",exponentWs);
Power alg;
alg.initialize();
alg.setPropertyValue("OutputWorkspace","OutputWs");
alg.setPropertyValue("RHSWorkspace", "exponentWs");
alg.setPropertyValue("LHSWorkspace", "baseWs");
alg.execute();
TSM_ASSERT_EQUALS("Exponent must be a single valued workspace", alg.isExecuted(),false);
AnalysisDataService::Instance().remove("exponentWs");
AnalysisDataService::Instance().remove("baseWs");
AnalysisDataService::Instance().remove("OutputWs");
}
void testPowerErrorCalculation()
{
//Workspace creation helper creates input error as sqrt of input value. So input error = 2.
//if x = p ^ y, then err_x = y * x * err_p / p
WorkspaceSingleValue_sptr baseWs = WorkspaceCreationHelper::CreateWorkspaceSingleValue(4);
WorkspaceSingleValue_sptr exponentWs = WorkspaceCreationHelper::CreateWorkspaceSingleValue(2);
AnalysisDataService::Instance().add("baseWs",baseWs);
AnalysisDataService::Instance().add("exponentWs",exponentWs);
Power alg;
alg.initialize();
alg.setPropertyValue("OutputWorkspace","OutputWs");
alg.setPropertyValue("RHSWorkspace", "exponentWs");
alg.setPropertyValue("LHSWorkspace", "baseWs");
alg.execute();
WorkspaceSingleValue_sptr output = boost::dynamic_pointer_cast<WorkspaceSingleValue>(AnalysisDataService::Instance().retrieve("OutputWs"));
Mantid::MantidVec expectedErrorValue(1,16);
TSM_ASSERT_EQUALS("Power algorithm error has not been determined properly.", expectedErrorValue, output->dataE());
AnalysisDataService::Instance().remove("exponentWs");
AnalysisDataService::Instance().remove("baseWs");
AnalysisDataService::Instance().remove("OutputWs");
}
};
#endif