[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