User Portlet
Featured Contributor
Discussions |
---|
Unfortunately I cannot replicate those spikes using version 13.2 or 14.0 on Ubuntu Linux. It might be platform dependent. |
We'll look into it. |
We take a derivative after that. Will the logs be off by something other than a constant in the variable `phi`? |
Interesting method I guess. It is different from the one I'll show, but at least the similarities of being iterative and geometrically convergent. The (reasonably) well-known "power method" can be used to find the largest root. It amounts to... |
A web search for “geometric algebra Mathematica package” turns up several hits. At least one was presented at a Wolfram Technology Conference, author being Terje Vold (I recall it being a good talk but I’ve not used the package). |
That appears to be the series for `1/(2+x^2)`. In[5]:= Series[1/(2+x^2),{x,0,8}]//InputForm Out[5]//InputForm= SeriesData[x, 0, {1/2, 0, -1/4, 0, 1/8, 0, -1/16, 0, 1/32}, 0, 9, 1] |
Ah. My mistake, I misread where the error message was located. |
You are setting precision higher than that of the input. This means some bits are effectively junk. As for assessing the condition number, one way is like so. {vv, ss, ww} = SingularValueDecomposition[g]; sv = Diagonal[ss]; ... |
A reference implementation of a recent variant can be found in [this prior Community thread](https://community.wolfram.com/groups/-/m/t/2344199?p_p_auth=UZ0BHnEu). |
Depending on the input, maybe you would want to do a Taylor expansion instead. |