Class: GSL::MultiFit::MultidLMNumDiff

Inherits:
MultidLM
  • Object
show all
Defined in:
lib/gsl_extras.rb

Instance Attribute Summary

Attributes inherited from MultidLM

#chi2, #covar, #dof, #my_chi2, #position

Class Method Summary collapse

Instance Method Summary collapse

Methods inherited from MultidLM

#eval, #solve

Constructor Details

#initialize(yproc, ndata, ndims, nparams) ⇒ MultidLMNumDiff

Returns a new instance of MultidLMNumDiff.



934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
 
# File 'lib/gsl_extras.rb', line 934

def initialize(yproc, ndata, ndims, nparams)
	fproc = Proc.new do |x, *gridpoints, y, sigma, f|
# 				gridpoints = (0...@ndims).to_a.map do |i| 
# 					@gridpoints.col(i)
# 				end
		for i in 0...ndata.size do
			f[i] = (@yproc.call(x, *gridpoints.map{|vec| vec[i]}) - y[i])/sigma[i]
		end
	end
	dfproc = Proc.new do |x, *gridpoints, y, sigma, jac|
			for j in 0...nparams do
				xj = x[j]
				xplus = x.dup
				xplus[j] = xj + @delt[j]
				xminus = x.dup
				xminus[j] = xj - @delt[j]

				for i in 0...ndata do
					gp = gridpoints.map{|vec| vec[i]}
					yplus = @yproc.call(xplus, *gp)
					yminus = @yproc.call(xminus, *gp)
					#p "delt", @delt, "sigma", @sigma, "jac", jac, "jac.shape", jac.shape, "result", (yplus - yminus)/2*@delt[j]/sigma[i]
					jac.set(i, j, (yplus - yminus)/2*@delt[j]/sigma[i])
				end
			end
	end	
	super(yproc, fproc, dfproc, ndata, ndims, nparams)


end

Class Method Details

.alloc(*params) ⇒ Object 



930
931
932
 
# File 'lib/gsl_extras.rb', line 930

def self.alloc(*params)
	new(*params)
end

Instance Method Details

#set_data(xstart, delt, *gridpoints, y, sigma) ⇒ Object 



964
965
966
967
 
# File 'lib/gsl_extras.rb', line 964

def set_data(xstart, delt, *gridpoints, y, sigma)
	@delt = (delt || GSL::Vector.alloc([1e-5]*xstart.size))
	super(xstart, *gridpoints, y, sigma)
end