Re: Matrix exponentiation in FORTRANThis is a discussion on Re: Matrix exponentiation in FORTRAN within the pl1 forums in Programming Languages category; "Arjan" <arjan.van.dijk @ rivm.nl> wrote in message news:6997dd3e-520c-4ebc-bbed-9e81c7e0a325 @ h25g2000hsf.googlegroups.com... > > I have used the FORTRAN-version of Expokit to estimate the > simultaneous decay > of 300 types of decaying atoms, where daughter-nuclides arise from > their mothers, > with good results. The code halts when you offer a zero-matrix, but > the author > of expokit has explained that this is a "feature". Modification of the > code was easy. > > Arjan Beware using this hotchpotch of Fortran 90 and F77. The codes use the CMPLX function on double precision values, and do not specify the kind. ... | ||||||||
| | ||||||||
| ||||||||
| Register | FAQ | Calendar | Search | Today's Posts | Mark Forums Read |
|
#1
03-31-2008, 09:23 AM
| |||
| |||
 Re: Matrix exponentiation in FORTRAN "Arjan" <arjan.van.dijk@rivm.nl> wrote in message news:6997dd3e-520c-4ebc-bbed-9e81c7e0a325@h25g2000hsf.googlegroups.com... > > I have used the FORTRAN-version of Expokit to estimate the > simultaneous decay > of 300 types of decaying atoms, where daughter-nuclides arise from > their mothers, > with good results. The code halts when you offer a zero-matrix, but > the author > of expokit has explained that this is a "feature". Modification of the > code was easy. > > Arjan Beware using this hotchpotch of Fortran 90 and F77. The codes use the CMPLX function on double precision values, and do not specify the kind. This converts to a single-precision result. As well, there are other dubious practices of using DBLE(A*B) where apparently DPROD should have been used, or DBLE(A)*DBLE(B).  |
 |
| Thread Tools | |
| Display Modes | |
Linear Mode
