Fortran - ? open SVD f90 code

Author

? open SVD f90 code

Cheng Cosine

2004-10-19, 3:58 pm

Though Fortran Recipes offers SVD code but that only works for real var.

I tried to d/l codes from netlib but my compiler, LF90v45, does not

compile it, but only gives me a series of errors. The package from

netlib contains many files so I dare not to trace them.

What are other site(s) I can d/l SVD alg that deals both real and complex

variables?

Thanks,
by Cheng Cosine
Oct/19/2k4 UT

Jugoslav Dujic

2004-10-19, 3:58 pm

Cheng Cosine wrote:
| Though Fortran Recipes offers SVD code but that only works for real var.
|
| I tried to d/l codes from netlib but my compiler, LF90v45, does not
|
| compile it, but only gives me a series of errors. The package from
|
| netlib contains many files so I dare not to trace them.
|
| What are other site(s) I can d/l SVD alg that deals both real and complex
|
| variables?

Which "netlib package"?

I'm not too familiar with SVD, but I did download some related LAPACK routines
(www.netlib.org/lapack) which compiled & worked just fine. Go to

http://www.netlib.org/lapack/individualroutines.html

choose your datatype, then search for "svd" there.

--
Jugoslav
___________
www.geocities.com/jdujic

Please reply to the newsgroup.
You can find my real e-mail on my home page above.

Cheng Cosine

2004-10-21, 3:58 am

"Jugoslav Dujic" <jdujic@yahoo.com> wrote in message
news:2tkoaaF2291mmU1@uni-berlin.de...
> Cheng Cosine wrote:
> | Though Fortran Recipes offers SVD code but that only works for real var.
> |
> | I tried to d/l codes from netlib but my compiler, LF90v45, does not
> |
> | compile it, but only gives me a series of errors. The package from
> |
> | netlib contains many files so I dare not to trace them.
> |
> | What are other site(s) I can d/l SVD alg that deals both real and
complex
> |
> | variables?
>
> Which "netlib package"?
>
> I'm not too familiar with SVD, but I did download some related LAPACK
routines
> (www.netlib.org/lapack) which compiled & worked just fine. Go to
>
> http://www.netlib.org/lapack/individualroutines.html
>
> choose your datatype, then search for "svd" there.
>

I guess we have the same files. One of the file is called

zgesvd.f90

I guess it's also from LAPACK you mentioned. But LF90v45 does not like it.
Neither does the free compiler in our school's computer center.

Jugoslav Dujic

2004-10-21, 8:57 am

Cheng Cosine wrote:
| "Jugoslav Dujic" <jdujic@yahoo.com> wrote in message
| news:2tkoaaF2291mmU1@uni-berlin.de...
|| Cheng Cosine wrote:
||| Though Fortran Recipes offers SVD code but that only works for real var.
|||
||| I tried to d/l codes from netlib but my compiler, LF90v45, does not
|||
||| compile it, but only gives me a series of errors. The package from
|||
||| netlib contains many files so I dare not to trace them.
|||
||| What are other site(s) I can d/l SVD alg that deals both real and complex
|||
||| variables?
||
|| Which "netlib package"?
||
|| I'm not too familiar with SVD, but I did download some related LAPACK
|| routines (www.netlib.org/lapack) which compiled & worked just fine. Go to
||
|| http://www.netlib.org/lapack/individualroutines.html
||
|| choose your datatype, then search for "svd" there.
||
|
| I guess we have the same files. One of the file is called
|
| zgesvd.f90
|
| I guess it's also from LAPACK you mentioned. But LF90v45 does not like it.
| Neither does the free compiler in our school's computer center.

Hmmm. LAPACK is one of the most nicely & most standardly written codes I've
seen.
Could it possibly be that Lahey doesn't support complex*16 (or at least
that particular syntax)?

--
Jugoslav
___________
www.geocities.com/jdujic

Please reply to the newsgroup.
You can find my real e-mail on my home page above.

Janne Blomqvist

2004-10-21, 8:57 am

On 2004-10-21, Jugoslav Dujic <jdujic@yahoo.com> wrote:
> Cheng Cosine wrote:
>| I guess we have the same files. One of the file is called
>|
>| zgesvd.f90

zgesvd.f, that is. LAPACK is fixed form F77 after all, not F90 free
form. Perhaps, if the file suffix has been changed the compiler assumes
it's free form.

>| I guess it's also from LAPACK you mentioned. But LF90v45 does not like it.
>| Neither does the free compiler in our school's computer center.
>
> Hmmm. LAPACK is one of the most nicely & most standardly written codes I've
> seen.
> Could it possibly be that Lahey doesn't support complex*16 (or at least
> that particular syntax)?

$ lf95 --version
Lahey/Fujitsu Fortran 95 Compiler Release L6.20b

$ lf95 --warn -c zgesvd.f
Encountered 0 errors, 0 warnings in file zgesvd.f.

--
Janne Blomqvist

Rich Townsend

2004-10-21, 4:03 pm

Jugoslav Dujic wrote:
> Cheng Cosine wrote:
> | "Jugoslav Dujic" <jdujic@yahoo.com> wrote in message
> | news:2tkoaaF2291mmU1@uni-berlin.de...
> || Cheng Cosine wrote:
> ||| Though Fortran Recipes offers SVD code but that only works for real var.
> |||
> ||| I tried to d/l codes from netlib but my compiler, LF90v45, does not
> |||
> ||| compile it, but only gives me a series of errors. The package from
> |||
> ||| netlib contains many files so I dare not to trace them.
> |||
> ||| What are other site(s) I can d/l SVD alg that deals both real and complex
> |||
> ||| variables?
> ||
> || Which "netlib package"?
> ||
> || I'm not too familiar with SVD, but I did download some related LAPACK
> || routines (www.netlib.org/lapack) which compiled & worked just fine. Go to
> ||
> || http://www.netlib.org/lapack/individualroutines.html
> ||
> || choose your datatype, then search for "svd" there.
> ||
> |
> | I guess we have the same files. One of the file is called
> |
> | zgesvd.f90
> |
> | I guess it's also from LAPACK you mentioned. But LF90v45 does not like it.
> | Neither does the free compiler in our school's computer center.
>
> Hmmm. LAPACK is one of the most nicely & most standardly written codes I've
> seen.
> Could it possibly be that Lahey doesn't support complex*16 (or at least
> that particular syntax)?
>

I've had no problems compiling LAPACK using the LF95 compiler (v 6.2c
for Linux). The LF90 product is pretty old.

Its worth mentioning that LAPACK is pure Fortran 77. There is a Fortran
95 wrapper to LAPACK, LAPACK95, but all of the routines contained in it
have the prefix la_. Therefore, I don't know where the 'zgesvd.f90'
routine is coming from -- but it doesn't belong to either library.

cheers,

Rich

--
Dr Richard H D Townsend
Bartol Research Institute
University of Delaware

[ Delete VOID for valid email address ]

Cheng Cosine

2004-10-21, 8:56 pm

I know LF90v45 is pretty old, but I am a student and I paid in full for
LF90v45 for couple hundred USD. It is not cheap for me to buy LF95. My
school does not teach Fortran anymore but only C.

That's why I am looking for something free or cheaper.

by Cheng Cosine
Oct/21/2k4 UT

Richard E Maine

2004-10-21, 8:56 pm

"Cheng Cosine" <acosine@ms13.url.com.tw> writes:

> I know LF90v45 is pretty old, but I am a student and I paid in full for
> LF90v45 for couple hundred USD. It is not cheap for me to buy LF95. My
> school does not teach Fortran anymore but only C.
>
> That's why I am looking for something free or cheaper.

You might get a little more constructive help if you were a little more
specific about your problem than that the compiler "does not like it"
and gave "a lot of error messages."

Frankly, I didn't bother to reply before because I figured that
someone who has posted as many questions here as you have ought to
know to post more information than that by now. I'd have probably
spoken up before if it was a newbie post.

I'd be pretty surprised if the compiler couldn't be convinced to compile
that code with no more than minor work, if that. But if you don't give
more information about the problem, you aren't going to get much more
in the way of useful help.

Though... I did notice the earlier post suggesting that you might be
using the wrong source form. That sure would produce a lot of error
messages... and be absolutely trivial to fix. Seemed to me like a
pretty plausible guess from very little data, but I heard noticed no
response at all to it. What are the odds that you already got told
the answer but didn't notice?

--
Richard Maine | Good judgment comes from experience;
email: my first.last at org.domain | experience comes from bad judgment.
org: nasa, domain: gov | -- Mark Twain

Sponsored Links