Viewing file: user_array.py (7.3 KB) -rw-r--r-- Select action/file-type: (+) | (+) | (+) | Code (+) | Session (+) | (+) | SDB (+) | (+) | (+) | (+) | (+) | (+) |
""" Standard container-class for easy multiple-inheritance. Try to inherit from the ndarray instead of using this class as this is not complete. """
from numpy.core import array, asarray, absolute, add, subtract, multiply, \ divide, remainder, power, left_shift, right_shift, bitwise_and, \ bitwise_or, bitwise_xor, invert, less, less_equal, not_equal, equal, \ greater, greater_equal, shape, reshape, arange, sin, sqrt, transpose
class container(object): def __init__(self, data, dtype=None, copy=True): self.array = array(data, dtype, copy=copy)
def __repr__(self): if len(self.shape) > 0: return self.__class__.__name__+repr(self.array)[len("array"):] else: return self.__class__.__name__+"("+repr(self.array)+")"
def __array__(self,t=None): if t: return self.array.astype(t) return self.array
# Array as sequence def __len__(self): return len(self.array)
def __getitem__(self, index): return self._rc(self.array[index])
def __getslice__(self, i, j): return self._rc(self.array[i:j])
def __setitem__(self, index, value): self.array[index] = asarray(value,self.dtype) def __setslice__(self, i, j, value): self.array[i:j] = asarray(value,self.dtype)
def __abs__(self): return self._rc(absolute(self.array)) def __neg__(self): return self._rc(-self.array)
def __add__(self, other): return self._rc(self.array+asarray(other)) __radd__ = __add__
def __iadd__(self, other): add(self.array, other, self.array) return self
def __sub__(self, other): return self._rc(self.array-asarray(other)) def __rsub__(self, other): return self._rc(asarray(other)-self.array) def __isub__(self, other): subtract(self.array, other, self.array) return self
def __mul__(self, other): return self._rc(multiply(self.array,asarray(other))) __rmul__ = __mul__ def __imul__(self, other): multiply(self.array, other, self.array) return self
def __div__(self, other): return self._rc(divide(self.array,asarray(other))) def __rdiv__(self, other): return self._rc(divide(asarray(other),self.array)) def __idiv__(self, other): divide(self.array, other, self.array) return self
def __mod__(self, other): return self._rc(remainder(self.array, other)) def __rmod__(self, other): return self._rc(remainder(other, self.array)) def __imod__(self, other): remainder(self.array, other, self.array) return self
def __divmod__(self, other): return (self._rc(divide(self.array,other)), self._rc(remainder(self.array, other))) def __rdivmod__(self, other): return (self._rc(divide(other, self.array)), self._rc(remainder(other, self.array)))
def __pow__(self,other): return self._rc(power(self.array,asarray(other))) def __rpow__(self,other): return self._rc(power(asarray(other),self.array)) def __ipow__(self,other): power(self.array, other, self.array) return self
def __lshift__(self,other): return self._rc(left_shift(self.array, other)) def __rshift__(self,other): return self._rc(right_shift(self.array, other)) def __rlshift__(self,other): return self._rc(left_shift(other, self.array)) def __rrshift__(self,other): return self._rc(right_shift(other, self.array)) def __ilshift__(self,other): left_shift(self.array, other, self.array) return self def __irshift__(self,other): right_shift(self.array, other, self.array) return self
def __and__(self, other): return self._rc(bitwise_and(self.array, other)) def __rand__(self, other): return self._rc(bitwise_and(other, self.array)) def __iand__(self, other): bitwise_and(self.array, other, self.array) return self
def __xor__(self, other): return self._rc(bitwise_xor(self.array, other)) def __rxor__(self, other): return self._rc(bitwise_xor(other, self.array)) def __ixor__(self, other): bitwise_xor(self.array, other, self.array) return self
def __or__(self, other): return self._rc(bitwise_or(self.array, other)) def __ror__(self, other): return self._rc(bitwise_or(other, self.array)) def __ior__(self, other): bitwise_or(self.array, other, self.array) return self
def __neg__(self): return self._rc(-self.array) def __pos__(self): return self._rc(self.array) def __abs__(self): return self._rc(abs(self.array)) def __invert__(self): return self._rc(invert(self.array))
def _scalarfunc(self, func): if len(self.shape) == 0: return func(self[0]) else: raise TypeError, "only rank-0 arrays can be converted to Python scalars."
def __complex__(self): return self._scalarfunc(complex) def __float__(self): return self._scalarfunc(float) def __int__(self): return self._scalarfunc(int) def __long__(self): return self._scalarfunc(long) def __hex__(self): return self._scalarfunc(hex) def __oct__(self): return self._scalarfunc(oct)
def __lt__(self,other): return self._rc(less(self.array,other)) def __le__(self,other): return self._rc(less_equal(self.array,other)) def __eq__(self,other): return self._rc(equal(self.array,other)) def __ne__(self,other): return self._rc(not_equal(self.array,other)) def __gt__(self,other): return self._rc(greater(self.array,other)) def __ge__(self,other): return self._rc(greater_equal(self.array,other))
def copy(self): return self._rc(self.array.copy())
def tostring(self): return self.array.tostring()
def byteswap(self): return self._rc(self.array.byteswap())
def astype(self, typecode): return self._rc(self.array.astype(typecode))
def _rc(self, a): if len(shape(a)) == 0: return a else: return self.__class__(a)
def __array_wrap__(self, *args): return self.__class__(args[0])
def __setattr__(self,attr,value): if attr == 'array': object.__setattr__(self, attr, value) return try: self.array.__setattr__(attr, value) except AttributeError: object.__setattr__(self, attr, value)
# Only called after other approaches fail. def __getattr__(self,attr): if (attr == 'array'): return object.__getattribute__(self, attr) return self.array.__getattribute__(attr)
############################################################# # Test of class container ############################################################# if __name__ == '__main__': temp=reshape(arange(10000),(100,100))
ua=container(temp) # new object created begin test print dir(ua) print shape(ua),ua.shape # I have changed Numeric.py
ua_small=ua[:3,:5] print ua_small ua_small[0,0]=10 # this did not change ua[0,0], which is not normal behavior print ua_small[0,0],ua[0,0] print sin(ua_small)/3.*6.+sqrt(ua_small**2) print less(ua_small,103),type(less(ua_small,103)) print type(ua_small*reshape(arange(15),shape(ua_small))) print reshape(ua_small,(5,3)) print transpose(ua_small)
|