escript/src/DataAbstract.cpp

// $Id$
/*
 ************************************************************
 *          Copyright 2006 by ACcESS MNRF                   *
 *                                                          *
 *              http://www.access.edu.au                    *
 *       Primary Business: Queensland, Australia            *
 *  Licensed under the Open Software License version 3.0    *
 *     http://www.opensource.org/licenses/osl-3.0.php       *
 *                                                          *
 ************************************************************
*/

#include "DataAbstract.h"
#include "DataException.h"

using namespace std;

namespace escript {

DataAbstract::DataAbstract(const FunctionSpace& what):
    m_noDataPointsPerSample(what.getNumDPPSample()),
    m_noSamples(what.getNumSamples()),
    m_functionSpace(what)
{
}

DataAbstract::~DataAbstract() 
{
}

void
DataAbstract::setPointDataView(const DataArrayView& input)
{
    m_pointDataView.reset(new DataArrayView(input.getData(),input.getShape(),input.getOffset()));
}

void
DataAbstract::resetPointDataView()
{
    m_pointDataView.reset(new DataArrayView());
}

void
DataAbstract::operandCheck(const DataAbstract& right) const
{
    if ((right.getNumDPPSample()!=getNumDPPSample()) ||
    (right.getNumSamples()!=getNumSamples()) ||
    (right.getFunctionSpace()!=getFunctionSpace())) {
      stringstream temp;
      temp << "Error - Right hand argument sample shape or function space "
       << "incompatible with left." << endl 
       << "LHS: (" << getNumSamples() << ","
       << getNumDPPSample() << ") " << getFunctionSpace().toString() 
       << endl
       << "RHS: (" << right.getNumSamples() << "," 
       << right.getNumDPPSample() << ") " 
       << right.getFunctionSpace().toString();
      throw DataException(temp.str());
    }

    //
    // Check the shape of the point data, a rank of 0(scalar) is okay
    if (!((right.getPointDataView().getRank()==0) || 
      (right.getPointDataView().getShape()==getPointDataView().getShape())))
      {
        stringstream temp;
    temp << "Error - Right hand argument point data shape: " 
         << DataArrayView::shapeToString(right.getPointDataView().getShape())
         << " doesn't match left: " 
         << DataArrayView::shapeToString(getPointDataView().getShape());
    throw DataException(temp.str());
      }
}

void
DataAbstract::dump(const std::string fileName) const
{
    throw DataException("Error - DataAbstract:: dump: not implemented.");
}


DataAbstract::ValueType::value_type*
DataAbstract::getSampleDataByTag(int tag)
{
    throw DataException("Error - DataAbstract::getSampleDataByTag: Data type does not have tag values.");
}

void
DataAbstract::setTaggedValue(int tagKey,
                             const DataArrayView& value)
{
    throw DataException("Error - DataAbstract::setTaggedValue: Data type does not have tag values.");
}

int
DataAbstract::getTagNumber(int dpno)
{
    throw DataException("Error - DataAbstract::getTagNumber: Data type cannot be accessed by tag values.");
    return (0);
}

int
DataAbstract::archiveData(ofstream& archiveFile,
                          const ValueType::size_type noValues) const
{
  return 0;
}

int
DataAbstract::extractData(ifstream& archiveFile,
                          const ValueType::size_type noValues)
{
  return 0;
}

void
DataAbstract::copyAll(const boost::python::numeric::array& value)
{
    throw DataException("Error - DataAbstract::copying data from numarray objects is not supported.");
}
void
DataAbstract::copyToDataPoint(const int sampleNo, const int dataPointNo, const double value)
{
    throw DataException("Error - DataAbstract::copying data from double value to a single data point is not supported.");
}
void
DataAbstract::copyToDataPoint(const int sampleNo, const int dataPointNo, const boost::python::numeric::array& value)
{
    throw DataException("Error - DataAbstract::copying data from numarray objects to a single data point is not supported.");
}

void
DataAbstract::symmetric(DataAbstract* ev) 
{
    throw DataException("Error - DataAbstract::symmetric is not supported.");
}

void
DataAbstract::nonsymmetric(DataAbstract* ev) 
{
    throw DataException("Error - DataAbstract::nonsymmetric is not supported.");
}

void
DataAbstract::trace(DataAbstract* ev, int axis_offset) 
{
    throw DataException("Error - DataAbstract::trace is not supported.");
}

void
DataAbstract::swapaxes(DataAbstract* ev, int axis0, int axis1)
{
    throw DataException("Error - DataAbstract::component swapaxes is not supported.");
}
void
DataAbstract::transpose(DataAbstract* ev, int axis_offset) 
{
    throw DataException("Error - DataAbstract::transpose is not supported.");
}

void
DataAbstract::eigenvalues(DataAbstract* ev) 
{
    throw DataException("Error - DataAbstract::eigenvalues is not supported.");
}
void
DataAbstract::eigenvalues_and_eigenvectors(DataAbstract* ev,DataAbstract* V,const double tol)
{
    throw DataException("Error - DataAbstract::eigenvalues_and_eigenvectors is not supported.");

}
void
DataAbstract::setToZero() 
{
    throw DataException("Error - DataAbstract:: cannot set values to zero.");
}


}  // end of namespace
1	// $Id$
2	/*
3	************************************************************
4	* Copyright 2006 by ACcESS MNRF *
5	* *
6	* http://www.access.edu.au *
7	* Primary Business: Queensland, Australia *
8	* Licensed under the Open Software License version 3.0 *
9	* http://www.opensource.org/licenses/osl-3.0.php *
10	* *
11	************************************************************
12	*/
13
14	#include "DataAbstract.h"
15	#include "DataException.h"
16
17	using namespace std;
18
19	namespace escript {
20
21	DataAbstract::DataAbstract(const FunctionSpace& what):
22	m_noDataPointsPerSample(what.getNumDPPSample()),
23	m_noSamples(what.getNumSamples()),
24	m_functionSpace(what)
25	{
26	}
27
28	DataAbstract::~DataAbstract()
29	{
30	}
31
32	void
33	DataAbstract::setPointDataView(const DataArrayView& input)
34	{
35	m_pointDataView.reset(new DataArrayView(input.getData(),input.getShape(),input.getOffset()));
36	}
37
38	void
39	DataAbstract::resetPointDataView()
40	{
41	m_pointDataView.reset(new DataArrayView());
42	}
43
44	void
45	DataAbstract::operandCheck(const DataAbstract& right) const
46	{
47	if ((right.getNumDPPSample()!=getNumDPPSample()) \|\|
48	(right.getNumSamples()!=getNumSamples()) \|\|
49	(right.getFunctionSpace()!=getFunctionSpace())) {
50	stringstream temp;
51	temp << "Error - Right hand argument sample shape or function space "
52	<< "incompatible with left." << endl
53	<< "LHS: (" << getNumSamples() << ","
54	<< getNumDPPSample() << ") " << getFunctionSpace().toString()
55	<< endl
56	<< "RHS: (" << right.getNumSamples() << ","
57	<< right.getNumDPPSample() << ") "
58	<< right.getFunctionSpace().toString();
59	throw DataException(temp.str());
60	}
61
62	//
63	// Check the shape of the point data, a rank of 0(scalar) is okay
64	if (!((right.getPointDataView().getRank()==0) \|\|
65	(right.getPointDataView().getShape()==getPointDataView().getShape())))
66	{
67	stringstream temp;
68	temp << "Error - Right hand argument point data shape: "
69	<< DataArrayView::shapeToString(right.getPointDataView().getShape())
70	<< " doesn't match left: "
71	<< DataArrayView::shapeToString(getPointDataView().getShape());
72	throw DataException(temp.str());
73	}
74	}
75
76	void
77	DataAbstract::dump(const std::string fileName) const
78	{
79	throw DataException("Error - DataAbstract:: dump: not implemented.");
80	}
81
82
83
84	DataAbstract::ValueType::value_type*
85	DataAbstract::getSampleDataByTag(int tag)
86	{
87	throw DataException("Error - DataAbstract::getSampleDataByTag: Data type does not have tag values.");
88	}
89
90	void
91	DataAbstract::setTaggedValue(int tagKey,
92	const DataArrayView& value)
93	{
94	throw DataException("Error - DataAbstract::setTaggedValue: Data type does not have tag values.");
95	}
96
97	int
98	DataAbstract::getTagNumber(int dpno)
99	{
100	throw DataException("Error - DataAbstract::getTagNumber: Data type cannot be accessed by tag values.");
101	return (0);
102	}
103
104	int
105	DataAbstract::archiveData(ofstream& archiveFile,
106	const ValueType::size_type noValues) const
107	{
108	return 0;
109	}
110
111	int
112	DataAbstract::extractData(ifstream& archiveFile,
113	const ValueType::size_type noValues)
114	{
115	return 0;
116	}
117
118	void
119	DataAbstract::copyAll(const boost::python::numeric::array& value)
120	{
121	throw DataException("Error - DataAbstract::copying data from numarray objects is not supported.");
122	}
123	void
124	DataAbstract::copyToDataPoint(const int sampleNo, const int dataPointNo, const double value)
125	{
126	throw DataException("Error - DataAbstract::copying data from double value to a single data point is not supported.");
127	}
128	void
129	DataAbstract::copyToDataPoint(const int sampleNo, const int dataPointNo, const boost::python::numeric::array& value)
130	{
131	throw DataException("Error - DataAbstract::copying data from numarray objects to a single data point is not supported.");
132	}
133
134	void
135	DataAbstract::symmetric(DataAbstract* ev)
136	{
137	throw DataException("Error - DataAbstract::symmetric is not supported.");
138	}
139
140	void
141	DataAbstract::nonsymmetric(DataAbstract* ev)
142	{
143	throw DataException("Error - DataAbstract::nonsymmetric is not supported.");
144	}
145
146	void
147	DataAbstract::trace(DataAbstract* ev, int axis_offset)
148	{
149	throw DataException("Error - DataAbstract::trace is not supported.");
150	}
151
152	void
153	DataAbstract::swapaxes(DataAbstract* ev, int axis0, int axis1)
154	{
155	throw DataException("Error - DataAbstract::component swapaxes is not supported.");
156	}
157	void
158	DataAbstract::transpose(DataAbstract* ev, int axis_offset)
159	{
160	throw DataException("Error - DataAbstract::transpose is not supported.");
161	}
162
163	void
164	DataAbstract::eigenvalues(DataAbstract* ev)
165	{
166	throw DataException("Error - DataAbstract::eigenvalues is not supported.");
167	}
168	void
169	DataAbstract::eigenvalues_and_eigenvectors(DataAbstract* ev,DataAbstract* V,const double tol)
170	{
171	throw DataException("Error - DataAbstract::eigenvalues_and_eigenvectors is not supported.");
172
173	}
174	void
175	DataAbstract::setToZero()
176	{
177	throw DataException("Error - DataAbstract:: cannot set values to zero.");
178	}
179
180
181
182	} // end of namespace
Name	Value
svn:eol-style	native
svn:keywords	Author Date Id Revision