Re: [Numpy-discussion] Fortran was dead ... [was Re: rewriting NumPy code in C or C++ or similar]

2011-03-23 Thread Sturla Molden
Den 16.03.2011 14:46, skrev Neal Becker: Also: * can it adopt external memory? Yes. Using CRAY pointers with libc malloc/free is e.g. a common way to get dynamic memory in Fortran 77. * can it interwork with numpy? (kinda required for this audience) Yes, that is why NumPy has f2py :-) -

Re: [Numpy-discussion] Fortran was dead ... [was Re: rewriting NumPy code in C or C++ or similar]

2011-03-23 Thread Sturla Molden
I've been busy lately, so I haven't had time to answer. The idea is of course to use Fortran with Python, not to build a framework in Fortran. I also find C useful to interface with the operating system, and Cython to write C extensions for Python. Also, if the computationally demaning parts

Re: [Numpy-discussion] Fortran was dead ... [was Re: rewriting NumPy code in C or C++ or similar]

2011-03-16 Thread Sturla Molden
Den 16.03.2011 00:01, skrev Neal Becker: Here is how Fortran compares: * 1-d, 2-d only or N-d?? Any of those. * support for slice views? What exactly kind of support? Fortran 90 pointers create a view. real*8, target :: array(n,m) real*8, pointer :: view view = array(::2, ::2) Slicing

Re: [Numpy-discussion] Fortran was dead ... [was Re: rewriting NumPy code in C or C++ or similar]

2011-03-16 Thread Neal Becker
Sturla Molden wrote: Den 16.03.2011 13:25, skrev Sturla Molden: Fortran 90 pointers create a view. real*8, target :: array(n,m) real*8, pointer :: view view = array(::2, ::2) Pardon, the second line should be real*8, pointer :: view(:,:) Sturla Also: * can it adopt external

Re: [Numpy-discussion] Fortran was dead ... [was Re: rewriting NumPy code in C or C++ or similar]

2011-03-16 Thread Yung-Yu Chen
On Wed, Mar 16, 2011 at 09:46, Neal Becker ndbeck...@gmail.com wrote: Sturla Molden wrote: Den 16.03.2011 13:25, skrev Sturla Molden: Fortran 90 pointers create a view. real*8, target :: array(n,m) real*8, pointer :: view view = array(::2, ::2) Pardon, the second line

Re: [Numpy-discussion] Fortran was dead ... [was Re: rewriting NumPy code in C or C++ or similar]

2011-03-16 Thread Ravi
On Monday 14 March 2011 15:02:32 Sebastian Haase wrote: Sturla has been writing so much about Fortran recently, and Ondrej now says he has done the move from C/C++ to Fortran -- I thought Fortran was dead ... !? ;-) What am I missing here ? Comparing Fortran with C++ is like comparing

Re: [Numpy-discussion] Fortran was dead ... [was Re: rewriting NumPy code in C or C++ or similar]

2011-03-16 Thread Dag Sverre Seljebotn
On 03/16/2011 08:10 PM, Ravi wrote: On Monday 14 March 2011 15:02:32 Sebastian Haase wrote: Sturla has been writing so much about Fortran recently, and Ondrej now says he has done the move from C/C++ to Fortran -- I thought Fortran was dead ... !? ;-) What am I missing here ? Comparing

Re: [Numpy-discussion] Fortran was dead ... [was Re: rewriting NumPy code in C or C++ or similar]

2011-03-15 Thread Sebastien Binet
On Mon, 14 Mar 2011 23:10:13 +0100, Matthieu Brucher matthieu.bruc...@gmail.com wrote: Hi, Intel Fortran is an excellent Fortran compiler. Why is Fortran still better than C and C++? - some rules are different, like arrays passed to functions are ALWAYS supposed to be independent in

Re: [Numpy-discussion] Fortran was dead ... [was Re: rewriting NumPy code in C or C++ or similar]

2011-03-15 Thread Sturla Molden
Den 14.03.2011 23:02, skrev Sebastian Haase: Sturla has been writing so much about Fortran recently, and Ondrej now says he has done the move from C/C++ to Fortran -- I thought Fortran was dead ... !? ;-) What am I missing here No, it is just that Fortran receives less hype. If Fortran was

Re: [Numpy-discussion] Fortran was dead ... [was Re: rewriting NumPy code in C or C++ or similar]

2011-03-15 Thread Sturla Molden
Den 14.03.2011 23:10, skrev Matthieu Brucher: Intel Fortran is an excellent Fortran compiler. Why is Fortran still better than C and C++? - some rules are different, like arrays passed to functions are ALWAYS supposed to be independent in Fortran, whereas in C, you have to add a restrict

Re: [Numpy-discussion] Fortran was dead ... [was Re: rewriting NumPy code in C or C++ or similar]

2011-03-15 Thread Matthieu Brucher
C++ templates maks binaries almost impossible to debug. Never had an issue with this and all my number crunching code is done through metaprogramming (with vectorization, cache blocking...) So I have a lot of complex template structures, and debugging them is easy. Then, if someone doesn't

Re: [Numpy-discussion] Fortran was dead ... [was Re: rewriting NumPy code in C or C++ or similar]

2011-03-15 Thread Sturla Molden
Den 14.03.2011 23:10, skrev Matthieu Brucher: - Fortran 95 has an excellent array support, which is not currently available in C/C++ (perhaps with ArBB?) In C++ you can actually make array libraries that behave almost like a Fortran compiler (cf. Blitz++, Intel Array Building Blocks), but

Re: [Numpy-discussion] Fortran was dead ... [was Re: rewriting NumPy code in C or C++ or similar]

2011-03-15 Thread Charles R Harris
On Tue, Mar 15, 2011 at 9:12 AM, Sturla Molden stu...@molden.no wrote: Den 14.03.2011 23:10, skrev Matthieu Brucher: - Fortran 95 has an excellent array support, which is not currently available in C/C++ (perhaps with ArBB?) In C++ you can actually make array libraries that behave almost

Re: [Numpy-discussion] Fortran was dead ... [was Re: rewriting NumPy code in C or C++ or similar]

2011-03-15 Thread Christopher Barker
On 3/15/11 8:33 AM, Charles R Harris wrote: There really isn't a satisfactory array library for C++. The fact that every couple of years there is another project to produce one testifies to that fact. And I think not just the fact that there is not one, but that perhaps C++ the language, or

Re: [Numpy-discussion] Fortran was dead ... [was Re: rewriting NumPy code in C or C++ or similar]

2011-03-15 Thread Sturla Molden
Den 15.03.2011 16:33, skrev Charles R Harris: There really isn't a satisfactory array library for C++. The fact that every couple of years there is another project to produce one testifies to that fact. In order to be competitive against Fortran 95, an array library for C++ must do all the

Re: [Numpy-discussion] Fortran was dead ... [was Re: rewriting NumPy code in C or C++ or similar]

2011-03-15 Thread Sturla Molden
Den 15.03.2011 17:10, skrev Christopher Barker: I've been slowly arriving to the conclusion that that is no place for C++ in programming. If you really need to twiddle bits, use C. If you need high performance numerics, use Fortran. If you need high level complex data structures, use Python

Re: [Numpy-discussion] Fortran was dead ... [was Re: rewriting NumPy code in C or C++ or similar]

2011-03-15 Thread Sturla Molden
Den 15.03.2011 16:05, skrev Matthieu Brucher: BTW, instead of Blitzz++, you have vigra and Eigen that are the new equivalent libraries, and you may want to keep an eye on Intel's ArBB. Intel's ArBB is interesting. But in order for this to work, there must be an idustry standard that other

Re: [Numpy-discussion] Fortran was dead ... [was Re: rewriting NumPy code in C or C++ or similar]

2011-03-15 Thread Yung-Yu Chen
On Tue, Mar 15, 2011 at 12:10, Christopher Barker chris.bar...@noaa.govwrote: On 3/15/11 8:33 AM, Charles R Harris wrote: There really isn't a satisfactory array library for C++. The fact that every couple of years there is another project to produce one testifies to that fact. And I

Re: [Numpy-discussion] Fortran was dead ... [was Re: rewriting NumPy code in C or C++ or similar]

2011-03-15 Thread Sturla Molden
Den 15.03.2011 18:01, skrev Yung-Yu Chen: I really love the capabilities Fortran provides for quick array operations, especially floating-points. What I think Fortran is still lacking is better support of C pointers and structures. Fortran 90 has user defined types, but they are not ABI

Re: [Numpy-discussion] Fortran was dead ... [was Re: rewriting NumPy code in C or C++ or similar]

2011-03-15 Thread Neal Becker
Charles R Harris wrote: On Tue, Mar 15, 2011 at 9:12 AM, Sturla Molden stu...@molden.no wrote: Den 14.03.2011 23:10, skrev Matthieu Brucher: - Fortran 95 has an excellent array support, which is not currently available in C/C++ (perhaps with ArBB?) In C++ you can actually make array

[Numpy-discussion] Fortran was dead ... [was Re: rewriting NumPy code in C or C++ or similar]

2011-03-14 Thread Sebastian Haase
On Mon, Mar 14, 2011 at 9:24 PM, Ondrej Certik ond...@certik.cz wrote: Hi Sturla, On Tue, Mar 8, 2011 at 6:25 AM, Sturla Molden stu...@molden.no wrote: Den 08.03.2011 05:05, skrev Dan Halbert: Thanks, that's a good suggestion. I have not written Fortran since 1971, but it's come a long way.

Re: [Numpy-discussion] Fortran was dead ... [was Re: rewriting NumPy code in C or C++ or similar]

2011-03-14 Thread Matthieu Brucher
Hi, Intel Fortran is an excellent Fortran compiler. Why is Fortran still better than C and C++? - some rules are different, like arrays passed to functions are ALWAYS supposed to be independent in Fortran, whereas in C, you have to add a restrict keyword - due to the last fact, Fortran is a