solve via nullspace

[pvazteixeira]

ApproxManifoldProducts.jl/src/ApproxManifoldProducts.jl

Line 26 in 963a36e

r = NLsolve.nlsolve(gg, [0.0])

[dehann]

we were able to make this work up to a scale factor. We think the scale should be easy to find, however, did not get to it when we were doing PoC. Either way, dont think we ever copied the working code anywhere. Hopefully we can remember how to do this with nullspace and copy MWE in here...

[dehann]

Ha, it's in the tests:

ApproxManifoldProducts.jl/test/ex_1D.jl

Lines 35 to 47 in f82da7c

	## These work, but we are not sure about how to get the scale for product of more than 2 elements
	# preS = nullspace(Lambdas*(Lambdas'))
	#
	# sc = sum(abs.(preS)) / abs(μ1-μ2)
	#
	# mu = vec(preS / sc)
	# Λ = getΛ(Lambdas)
	#
	#
	# u,s,v = svd(Lambdas*(Lambdas'))
	#
	# s[1,1]*u[:,2]