[SciPy-user] argmin and optimization
Nils Wagner
nwagner@iam.uni-stuttgart...
Thu Jul 26 07:05:51 CDT 2007
dmitrey wrote:
> Nils Wagner wrote:
>
>> dmitrey wrote:
>>
>>
>>> So your Omega is just all x: (x,x)<1?
>>>
>>>
>>>
>> (x,x) <= 1 ("closed ball constraint")
>>
>>
>>
>>> (Can you attach non-Latex notation elseware?)
>>>
>>>
>>>
>> I can't see what you mean here.
>>
>>
> I meant couldn't you describe what does Omega and P_Omega mean w/o Latex?
> And/or describe howto see Latex strings in Linux?
>
>>> Why can't you just set P_Omega(y) = {
>>> y, for ||y|| <=1
>>> y/||y||, for ||y|| > 1
>>> }
>>> ?
>>>
>>>
>>>
>> I don't see your point.
>>
>> Nils
>>
>>
> I meant is this true that the func P_Omega will yield the result you intend?
> from numpy.linalg import norm
> def P_Omega(y):
> n_y = norm(y)
> if n_y <= 1: return y
> else: return y/ny
> HTH, D
>
>>
>> _______________________________________________
>> SciPy-user mailing list
>> SciPy-user@scipy.org
>> http://projects.scipy.org/mailman/listinfo/scipy-user
>>
>>
>>
>>
>>
>
> _______________________________________________
> SciPy-user mailing list
> SciPy-user@scipy.org
> http://projects.scipy.org/mailman/listinfo/scipy-user
>
Dmitrey,
Thank you very much for your help. It works fine for me.
I have attached the script for the sake of completeness.
Cheers,
Nils
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test_golub.py
Type: text/x-python
Size: 1151 bytes
Desc: not available
Url : http://projects.scipy.org/pipermail/scipy-user/attachments/20070726/341c1f05/attachment.py
More information about the SciPy-user
mailing list