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 a

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-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 ? > Comp

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 Ma

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 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

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: * c

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 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 ___ NumPy-D

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-15 Thread Neal Becker
Charles R Harris wrote: > On Tue, Mar 15, 2011 at 9:12 AM, Sturla Molden 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 librar

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 Yung-Yu Chen
On Tue, Mar 15, 2011 at 12:10, Christopher Barker wrote: > 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 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 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 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: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

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 Charles R Harris
On Tue, Mar 15, 2011 at 9:12 AM, Sturla Molden 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 like a > Fo

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 th

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 w

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 > rest

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

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 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 Fortran, whereas in C, you

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 langua

[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 wrote: > Hi Sturla, > > On Tue, Mar 8, 2011 at 6:25 AM, Sturla Molden 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. I was a little worried