Thread View: sci.math.num-analysis
2 messages
2 total messages
Started by Han de Bruijn
Mon, 02 Nov 1992 08:57
RE: Complex Matrices: C source
Author: Han de Bruijn
Date: Mon, 02 Nov 1992 08:57
Date: Mon, 02 Nov 1992 08:57
26 lines
1794 bytes
1794 bytes
In article <mjerosc-311092221815@esemac.erenj.com> Michael Jerosch-Herold: > I know that lapack has routines which accomplish this, but > I prefer C code and before I try my luck with the Fortran to > C translator "f2c" I thought I first inquire if such code > can be found at an ftp site. Also "f2c" is not too good at > dealing with code related to matrices, I heard. Another possibility is to _interface_ Fortran with C. It requires some rudumentary knowledge of _both_ languages, but usually pays the effort. (I'm an experienced Fortran programmer but I learned a few bits of C.) For example NAG brings out the so called "C header files". If you include these in a C-program, then you can invoke routines from the NAG Fortran library (object codes), as if they were written in C. Why should you re-invent the wheel, huh? > Any ftp site which has this sort of code? In what language is > BLAS written? It's all in Fortran, of course. -- * Han de Bruijn; Applications&Graphics | "A little bit of Physics * No * TUD Computing Centre; P.O. Box 354 | would be NO idleness in * Oil * 2600 AJ Delft; The Netherlands. | Mathematics" (HdB). * for * E-mail: Han.deBruijn@RC.TUDelft.NL --| Fax: +31 15 78 37 87 ----* Blood
Re: Complex Matrices: C source
Author: M. J. Saltzman
Date: Mon, 02 Nov 1992 18:31
Date: Mon, 02 Nov 1992 18:31
36 lines
1966 bytes
1966 bytes
In article <1992Nov2.085721.5553@dutrun2.tudelft.nl> rcpshdb@dutrun2.tudelft.nl (Han de Bruijn) writes: >In article <mjerosc-311092221815@esemac.erenj.com> Michael Jerosch-Herold: >[...] > >> Any ftp site which has this sort of code? In what language is >> BLAS written? > >It's all in Fortran, of course. Well, they were certainly designed with Fortran in mind, so the interfaces assume things like "column-major" storage for 2-dimensional arrays, and the original, portable BLAS were written in Fortran. On the other hand, there are C implementations available as well, and there's no inherent reason why BLAS routines couldn't be written in any language, for the most natural vector and matrix data structures in that language. The BLAS is mainly a description of certain vector and matrix operations, not necessarily a particular implementation of that description. Also, I'm sure that most vendor-supplied BLAS are written in carefully hand-optimized assembler, taking account of cache, vector registers, pipelining, etc., available on the *particular* machine for which they are written. Most vendor-supplied BLAS *are* written with Fortran calling conventions. >-- >* Han de Bruijn; Applications&Graphics | "A little bit of Physics * No >* TUD Computing Centre; P.O. Box 354 | would be NO idleness in * Oil >* 2600 AJ Delft; The Netherlands. | Mathematics" (HdB). * for >* E-mail: Han.deBruijn@RC.TUDelft.NL --| Fax: +31 15 78 37 87 ----* Blood -- Matthew Saltzman Clemson University Math Sciences mjs@clemson.edu
Thread Navigation
This is a paginated view of messages in the thread with full content displayed inline.
Messages are displayed in chronological order, with the original post highlighted in green.
Use pagination controls to navigate through all messages in large threads.
Back to All Threads