[MATRIX-SIG] umath.logical_blah vs normal python functions.

[Phil Austin]

>>>>> "Rob" == Rob Hooft <R.Hooft@EuroMail.com> writes:

Rob>>        min=Numeric.minimum(image0,image1)
Rob>>        sum=image0+image1
Rob>>        difsquare=(image0-image1)**2
Rob>>        condition=Numeric.greater(difsquare,4*sum)
Rob>>        # if dif^2>4*sum return 2*min, else return sum
Rob>>        result=Numeric.choose(condition,(sum,2*min))

and the obligatory posting from an S-plus expatriate:

if Python adopted Splus-style indexing this could look like:

condition= difsquare > 4*sum
result[condition]=sum
result[!condition]=2*min


Regards, Phil


Phil Austin		INTERNET: phil@geog.ubc.ca
(604) 822-2175		FAX:	  (604) 822-6150

http://www.geog.ubc.ca/~phil
Associate Professor
Atmospheric Sciences Programme
Geography #217
University of British Columbia
1984 W Mall
Vancouver, BC  V6T 1Z2
CANADA



_______________
MATRIX-SIG  - SIG on Matrix Math for Python

send messages to: matrix-sig@python.org
administrivia to: matrix-sig-request@python.org
_______________