Thanks Chris,
I'm using IDL in a supercomputer and it has 8.2 :( (by the way, it's listed as a "graphics package" wtf?) http://nci.org.au/nci-systems/national-facility/peak-system/ raijin/application-software/
Will check what are the chances of upgrading to 8.3 but I doubt it will be possible in the short term.
One of the reasons I need to use this python bridge is to be able to read .vrt files (see https://groups.google.com/d/msg/comp.lang.idl-pvwave/YF2EaUq nP_k/bhMcY_Iqk8MJ ) but now I'm stuck again :(
JP
On Friday, 3 October 2014 01:33:39 UTC+10, Chris Torrence wrote:
> On Thursday, October 2, 2014 5:24:36 AM UTC-6, JP wrote:
>
>> well,
>
>>
>
>> my IDL didn't crash with your code, Have you tried asking the people in www.jacquette.com ?
>
>>
>
>>
>
>>
>
>> JP
>
>>
>
>>
>
>>
>
>>
>
>>
>
>>
>
>>
>
>> [jpg599@raijin2 ~]$ idl
>
>>
>
>> IDL Version 8.2.1 (linux x86_64 m64). (c) 2012, Exelis Visual Information Solutions, Inc.
>
>>
>
>> Installation number: 237570.
>
>>
>
>> Licensed for use by: ANU Supercomputer Facility
>
>>
>
>>
>
>>
>
>> IDL> py = pyimport('__main__')
>
>>
>
>> py->exec, 'import matplotlib.mlab as mlab'
>
>>
>
>> py->exec, 'import numpy as np'
>
>>
>
>> py->exec, 'import random as r'
>
>>
>
>> py->exec, 'n = 50'
>
>>
>
>> py->exec, 'a = np.array([r.random() for i in range(n)]).reshape(1,n)'
>
>>
>
>> py->exec, 'b = np.array([r.random() for i in range(n)]).reshape(1,n)'
>
>>
>
>> py->exec, 'r = np.concatenate((a,b), 0).transpose()'
>
>>
>
>> Slither 3.0, by Jacquette Consulting, Inc. Built on 2014-02-03 11:40:48.314070
>
>>
>
>> Copyright (c) 2011 by Jacquette Consulting, Inc. All rights reserved.
>
>>
>
>> http://www.jacquette.com
>
>>
>
>>
>
>>
>
>> Jacquette Consulting, Inc. grants you a non-exclusive and non-transferable
>
>>
>
>> license to use Slither (hereinafter "the software").
>
>>
>
>>
>
>>
>
>> The granting of this license does not obligate Jacquette Consulting, Inc. to
>
>>
>
>> provide any sort of technical support.
>
>>
>
>>
>
>>
>
>> THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
>
>>
>
>> IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
>
>>
>
>> FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
>
>>
>
>> AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
>
>>
>
>> LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
>
>>
>
>> OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
>
>>
>
>> SOFTWARE.
>
>>
>
>> % Loaded DLM: SLITHER.
>
>>
>
>> IDL> py->exec, 'import matplotlib.mlab as mlab'
>
>>
>
>> IDL> py->exec, 'import numpy as np'
>
>>
>
>> IDL> py->exec, 'import random as r'
>
>>
>
>> IDL> py->exec, 'n = 50'
>
>>
>
>> IDL> py->exec, 'a = np.array([r.random() for i in range(n)]).reshape(1,n)'
>
>>
>
>> IDL> py->exec, 'b = np.array([r.random() for i in range(n)]).reshape(1,n)'
>
>>
>
>> IDL> py->exec, 'r = np.concatenate((a,b), 0).transpose()'
>
>>
>
>> IDL> py->exec, 'p = mlab.PCA(r)' ; it crashes here
>
>>
>
>> IDL>
>
>>
>
>>
>
>>
>
>>
>
>>
>
>>
>
>>
>
>> On Thursday, 2 October 2014 20:04:45 UTC+10, Fabien wrote:
>
>>
>
>>> Hi JP,
>
>>
>
>>>
>
>>
>
>>>
>
>>
>
>>>
>
>>
>
>>> On 02.10.2014 04:01, JP wrote:
>
>>
>
>>>
>
>>
>
>>>> So i get an error when trying to do 'a = np.array([r.random() for n in range(n)]).reshape(1,n)'
>
>>
>
>>>
>
>>
>
>>>> but in any case, no crash when I get to the mlab.PCA.
>
>>
>
>>>
>
>>
>
>>>
>
>>
>
>>>
>
>>
>
>>> yes this was a stupid error in the python code. This should work better:
>
>>
>
>>>
>
>>
>
>>>
>
>>
>
>>>
>
>>
>
>>> py = pyimport('__main__')
>
>>
>
>>>
>
>>
>
>>> py->exec, 'import matplotlib.mlab as mlab'
>
>>
>
>>>
>
>>
>
>>> py->exec, 'import numpy as np'
>
>>
>
>>>
>
>>
>
>>> py->exec, 'import random as r'
>
>>
>
>>>
>
>>
>
>>> py->exec, 'n = 50'
>
>>
>
>>>
>
>>
>
>>> py->exec, 'a = np.array([r.random() for i in range(n)]).reshape(1,n)'
>
>>
>
>>>
>
>>
>
>>> py->exec, 'b = np.array([r.random() for i in range(n)]).reshape(1,n)'
>
>>
>
>>>
>
>>
>
>>> py->exec, 'r = np.concatenate((a,b), 0).transpose()'
>
>>
>
>>>
>
>>
>
>>> py->exec, 'p = mlab.PCA(r)' ; it crashes here
>
>>
>
>>>
>
>>
>
>>>
>
>>
>
>>>
>
>>
>
>>> And this still makes IDL crash on my machine. I gave up on using slither
>
>>
>
>>>
>
>>
>
>>> long ago now, which is a pity since I'd like to play around with python
>
>>
>
>>>
>
>>
>
>>> in IDL much more often...
>
>
>
> Hi Fabien and JP,
>
>
>
> I've played around with Slither a bit, and I have some suggestions. First, I would strongly recommend not bothering to try it in IDL 8.2.x or earlier. It was compiled against IDL 8.3, and if you can somehow get it to work with an earlier version, you are getting extremely lucky. If it crashes against IDL 8.2.x, I wouldn't even try to dig deeper. This is because the symbol tables tend to shift around with each new release of IDL, and your DLM ends up pointing at the wrong symbols (IDL_Deltmp instead of IDL_Gettmp, for example :-).
>
>
>
> With IDL 8.4, we tried very hard to keep the symbol table the same as IDL 8.3, specifically so that Slither would keep working - so you should be able to upgrade to IDL 8.4 (in a few weeks) without any problems.
>
>
>
> Second, I think using "Exec" is not the right approach. Slither was really designed to mimic the Python API from within IDL. For example, I would rewrite your code as:
>
>
>
> mlab = pyimport("matplotlib.mlab")
>
> np = pyimport("numpy")
>
> r = pyimport("numpy.random")
>
> n = 50
>
> a = r.rand(n) ; why not use IDL's mersenne twister??
>
> b = r.rand(n)
>
> c = TRANSPOSE([[a],[b]])
>
> p = mlab.PCA(c)
>
>
>
> Now, obviously if you want to use things like Python generators or Lambda functions, you will need to use Exec. But in general, it will make your code cleaner and easier to debug if you simply treat the Python objects just like IDL objects and call methods directly on the imported libraries.
>
>
>
> Hope this helps a bit.
>
>
>
> Cheers,
>
> Chris
>
> ExelisVIS
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
