[SciPy-user] do complex numbers default to double precision?
Ryan Krauss
ryanfedora at comcast.net
Thu Jun 30 08:08:17 CDT 2005
Thanks to Nils and Robert for their quick responses. This is definitely
one thing to love about SciPy. I posted my null space question
yesterday and Robert responded in 20 minutes and Nils in 30. I posted
my optimization question this morning and Nils reponded in 8 minutes!
It is almost like chatting with tech support. You make SciPy a pleasure
to use! Thanks.
Coincidentally, these two problems are linked and I don't know if my
numerical error problems are from fmin or the null space/svd stuff.
Once I find the input that drives my matrix to have a null space, I find
the vector that corresponds to the null space (assuming the sub-matrix
is only rank deficient by 1). Then I combined the null space vector
with zeros that correspond with the boundary condition on one end of the
problem. So, a 4x1 null space vector would give me an 8x1 full vector.
I then take the 8x8 full matrix which has a 4x4 submatrix whose
determinant is roughly 0 and multiply it by the 8x1 vector of the
boudnary conditions I just solved for. This should then give me an 8x1
vector of the boundary conditions on the other end. The problem is that
there are some elements of this second vector that should be 0 because
of the boundary conditions and they are actually of order 1e-10, if the
vector is normalized so that its magnitude is 1. This physically means
that I have a cantilever beam with a free end that has just a little bit
of force and moment at the free tip.
Ryan
Ryan Krauss wrote:
> The matrix is currently 4x4 but will grow to probably 6x6. It is
> definitely nonlinear. The matrix contains sinh, cosh, sin, and cos.
> I am using the transfer matrix method to analyze structures. When you
> say two-parameter, do you mean the real and imaginary part of the
> independent variable? I guess you are right that I don't necessarily
> need to use the determinant. In order to satisfy the boundary
> conditions of the problem this 4x4 or 6x6 matrix (which is really a
> submatrix of an 8x8 or 12x12) must have a null space. So, what would
> be the better thing to look for? An eignevalue that approaches 0?
>
> Ryan
>
> Nils Wagner wrote:
>
>> Ryan Krauss wrote:
>>
>>> I have a matrix that is a function of a complex valued input. I am
>>> trying to find that value of that input that drives the determinant
>>> of the matrix to zero. I am searching for this value using fmin.
>>> The error I am trying to minimize is the abs(det(complex matrix)).
>>
>>
>>
>> It's not a good idea to use the determinant directly since det(A) is
>> a rapidly varying function. As far as I understand your problem,
>> you are interested in the solution of a two-parameter nonlinear
>> eigenvalue problem. Is that correct ? How about the size of your
>> complex matrix A ?
>>
>> Nils
>>
>>> I don't seem to be able to drive this error lower that roughly
>>> 9e-17, regardless of the values for ftol and xtol I use.
>>>
>>> Am I hitting some internal limitation? Are complex values by
>>> default single or double precision?
>>>
>>> Thanks,
>>>
>>> Ryan
>>>
>>> _______________________________________________
>>> SciPy-user mailing list
>>> SciPy-user at scipy.net
>>> http://www.scipy.net/mailman/listinfo/scipy-user
>>
>>
>>
>>
>>
>>
>> _______________________________________________
>> SciPy-user mailing list
>> SciPy-user at scipy.net
>> http://www.scipy.net/mailman/listinfo/scipy-user
>>
>
> _______________________________________________
> SciPy-user mailing list
> SciPy-user at scipy.net
> http://www.scipy.net/mailman/listinfo/scipy-user
>
More information about the SciPy-user
mailing list