Class: GSL::Matrix::Complex

Inherits:
Object
  • Object
show all
Defined in:
lib/gs2crmod/astrogk/gsl_tools.rb,
lib/gs2crmod/gsl_tools.rb

Class Method Summary collapse

Instance Method Summary collapse

Class Method Details

.re_im(re, im) ⇒ Object 



48
49
50
51
52
53
54
55
56
57
58
 
# File 'lib/gs2crmod/astrogk/gsl_tools.rb', line 48

def self.re_im(re, im)
	raise "Shape of real and imaginary matrices must match" unless re.shape == im.shape
	rows, cols = re.shape
	mat = alloc(rows, cols)
	for i in 0...rows
		for j in 0...cols
			mat[i,j] = GSL::Complex.alloc([re[i,j], im[i,j]])
		end
	end
	return mat
end

Instance Method Details

#backward_cols_c2c(normalise = false) ⇒ Object 



61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
 
# File 'lib/gs2crmod/astrogk/gsl_tools.rb', line 61

def backward_cols_c2c(normalise = false)
	gm = self.dup
	rows, cols = gm.shape
	table = GSL.cache[[:fft_table, :complex, rows]] ||= GSL::FFT::ComplexWavetable.alloc(rows)
	work = GSL.cache[[:fft_work, :complex, rows]] ||= GSL::FFT::ComplexWorkspace.alloc(rows)
	for i in 0...cols
		vec = gm.col(i)
		vec.backward!(table, work)
		for j in 0...rows
			gm[j,i] =  vec[j]
		end
	end
	gm = gm / rows if normalise
	gm
end

#backward_rows_cc2r(normalise = false) ⇒ Object 



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
 
# File 'lib/gs2crmod/astrogk/gsl_tools.rb', line 92

def backward_rows_cc2r(normalise = false)
	gm = self.dup
	rows, cols = gm.shape
# 	if cols%2 == 0
# 		newcols = cols*2
# 	else
# 		newcols = cols*2 - 1		
# 	end
	was_even = rows.times.inject(true) do |bool, i|
		bool and (gm[i, cols - 1].imag == 0.0)
	end
# 	ep was_even
		
	if was_even
		newcols = cols * 2 - 2
	else
		newcols = cols * 2 - 1
	end
	gm_re = GSL::Matrix.alloc(rows, newcols)
		
	table = GSL.cache[[:fft_table, :real, newcols]] ||= GSL::FFT::RealWavetable.alloc(newcols)
	work = GSL.cache[[:fft_work, :real, newcols]] ||= GSL::FFT::RealWorkspace.alloc(newcols)
	row = GSL::Vector::Complex.alloc(cols)
# 	p rows
	for i in 0...rows
# 		p i
# 		row = gm.row(i)
		(0...cols).each{|j| row[j] = gm[i,j]}
		if was_even
			vec = row.concat(row.subvector(1, row.size - 2).reverse.conjugate) if cols > 2
		else
			vec = row.concat(row.subvector(1, row.size - 1).reverse.conjugate) if cols > 1
		end
		vec.backward!(table, work)
		for j in 0...newcols
			gm_re[i,j] =  vec[j].real
		end
	end
	gm_re = gm_re / newcols.to_f if normalise
	gm_re
end

#forward_cols_c2c ⇒ Object 



77
78
79
80
81
82
83
84
85
86
87
88
89
90
 
# File 'lib/gs2crmod/astrogk/gsl_tools.rb', line 77

def forward_cols_c2c
	gm = self.dup
	rows, cols = gm.shape
	table = GSL.cache[[:fft_table, :complex, rows]] ||= GSL::FFT::ComplexWavetable.alloc(rows)
	work = GSL.cache[[:fft_work, :complex, rows]] ||= GSL::FFT::ComplexWorkspace.alloc(rows)
	for i in 0...cols
		vec = gm.col(i)
		vec.forward!(table, work)
		for j in 0...rows
			gm[j,i] =  vec[j]
		end
	end
	gm
end