Class: GSL::Tensor

Inherits:
Object
  • Object
show all
Defined in:
lib/gs2crmod/gsl_data_3d.rb

Direct Known Subclasses

TensorComplex

Instance Attribute Summary collapse

Class Method Summary collapse

Instance Method Summary collapse

Constructor Details

#initialize(narray) ⇒ Tensor

Returns a new instance of Tensor.



32
33
34
 
# File 'lib/gs2crmod/gsl_data_3d.rb', line 32

def initialize(narray)
	@narray = narray
end

Dynamic Method Handling

This class handles dynamic methods through the method_missing method

#method_missing(meth, *args) ⇒ Object 



70
71
72
73
74
75
76
77
78
79
 
# File 'lib/gs2crmod/gsl_data_3d.rb', line 70

def method_missing(meth, *args)
	result = @narray.send(meth, *args.reverse)
	if result.kind_of? NArray
		self.class.new(result)
	else
		result
	end
rescue NoMethodError
	self.class.new(NMath.send(meth, @narray))
end

Instance Attribute Details

#narray ⇒ Object (readonly)

Returns the value of attribute narray.



28
29
30
 
# File 'lib/gs2crmod/gsl_data_3d.rb', line 28

def narray
  @narray
end

Class Method Details

.alloc(*args) ⇒ Object 



29
30
31
 
# File 'lib/gs2crmod/gsl_data_3d.rb', line 29

def self.alloc(*args)
	new(NArray.float(*args.reverse))
end

.method_missing(meth, *args) ⇒ Object 



21
22
23
24
25
26
 
# File 'lib/gs2crmod/gsl_data_3d.rb', line 21

def method_missing(meth, *args)
	#ep 'calling... ', meth, args
	ans = new(NArray.send(meth, *args.reverse))
	#ep 'got', ans
	ans
end

Instance Method Details

#[](*args) ⇒ Object 



38
39
40
41
42
43
44
45
46
47
48
49
50
51
 
# File 'lib/gs2crmod/gsl_data_3d.rb', line 38

def [](*args)
	#if args.inject(true){|b,i| b and i.kind_of? Integer}
		#@narray[*args.reverse]
	#else
		#self.class.new(@narray[*args.reverse])
	#end
 	case 	ans =	@narray[*args.reverse]
	when Numeric
		ans
	else
		self.class.new(@narray[*args.reverse])
	end
		
end

#[]=(*args, value) ⇒ Object 



52
53
54
55
56
 
# File 'lib/gs2crmod/gsl_data_3d.rb', line 52

def []=(*args, value)
#def []=(*args)
	#ep 'args', args, value
	@narray[*args.reverse] = value
end

#inspect ⇒ Object 



35
36
37
 
# File 'lib/gs2crmod/gsl_data_3d.rb', line 35

def inspect
	"GSL::Tensor.new(#{@narray.inspect})"
end

#iterate(&block) ⇒ Object 



80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
 
# File 'lib/gs2crmod/gsl_data_3d.rb', line 80

def iterate(&block)
	shp = shape
	cumul = 1
	cumulshp = []
	for i in 1..shape.size
		cumulshp[shp.size-i] = cumul
		cumul *= shp[shp.size-i]
	end
	 #= shape.reverse.map{|dim| cumul*=(dim); cumul.to_i}.reverse
	#ep cumulshp; gets
	(cumulshp[0]*shp[0]).times do |n|
		#indexes = cumulshp.reverse.map{|cumul| rem = n%cumul; n-=rem; rem}.reverse
		indexes = cumulshp.map{|cumul| idx = (n/cumul).floor; n -= idx*cumul; idx}
		yield(*indexes)
	end
end

#iterate_row_maj(&block) ⇒ Object 



96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
 
# File 'lib/gs2crmod/gsl_data_3d.rb', line 96

def iterate_row_maj(&block)
	shp = shape
	cumul = 1
	cumulshp = []
	for i in 0...shape.size
		cumulshp[i] = cumul
		cumul *= shp[i]
	end
	 #= shape.reverse.map{|dim| cumul*=(dim); cumul.to_i}.reverse
	#ep cumulshp; gets
	(cumulshp[-1]*shp[-1]).times do |n|
		#indexes = cumulshp.reverse.map{|cumul| rem = n%cumul; n-=rem; rem}.reverse
		indexes = cumulshp.reverse.map{|cumul| idx = (n/cumul).floor; n -= idx*cumul; idx}.reverse
		yield(*indexes)
	end
end

#reshape!(*args) ⇒ Object 



60
61
62
63
 
# File 'lib/gs2crmod/gsl_data_3d.rb', line 60

def reshape!(*args)
	#ep 'rags', args
	@narray.reshape!(*args.reverse)
end

#shape ⇒ Object 



57
58
59
 
# File 'lib/gs2crmod/gsl_data_3d.rb', line 57

def shape
	@narray.shape.reverse
end

#to_a ⇒ Object 



64
65
66
 
# File 'lib/gs2crmod/gsl_data_3d.rb', line 64

def to_a
	@narray.transpose(*(0...@narray.shape.size).to_a.reverse).to_a
end

#transpose(*args) ⇒ Object 



67
68
69
 
# File 'lib/gs2crmod/gsl_data_3d.rb', line 67

def transpose(*args)
	self.class.new(@narray.transpose(*args))
end