I am attaching the Fortran .f files rs.f & rg.f, if that would help. You can also get these routines (and more, such as the routines called by rg.f and rs.f) from "netlib" - http://www.netlib.org/eispack/index.html
The reason that you can find these from google search is that they are part of EISPACK which is a well-known linear algebra package, whereas wilder.f is a fortran routine written specifically for a particular package. You should contact the author of that package for the source code. For the source codes in "mypkg", look in the directory (in Windows, for example): C:\Program Files\R\R-2.6.1/src/library/mypkg/src For building R from source check the R-FAQ. Hope this helps, Ravi. ---------------------------------------------------------------------------- ------- Ravi Varadhan, Ph.D. Assistant Professor, The Center on Aging and Health Division of Geriatric Medicine and Gerontology Johns Hopkins University Ph: (410) 502-2619 Fax: (410) 614-9625 Email: [EMAIL PROTECTED] Webpage: http://www.jhsph.edu/agingandhealth/People/Faculty/Varadhan.html ---------------------------------------------------------------------------- -------- -----Original Message----- From: Shubha Vishwanath Karanth [mailto:[EMAIL PROTECTED] Sent: Friday, February 08, 2008 1:05 AM To: Ravi Varadhan; [EMAIL PROTECTED] Subject: RE: [R] .fortran code Hi Ravi, ...In the same fashion, I was trying to look at wilder.f of 'wilderSum' command of the TTR package. I even googled for the keyword, "wilder.f wilderSum" or " subroutine wilder", but did not find anything useful. Also I did not get when you said " So,you can look at these files if you had installed the source code". Where to look at these files(without going to google)...is there anything which I need to install more? BR, Shubha Shubha Karanth | Amba Research Ph +91 80 3980 8031 | Mob +91 94 4886 4510 Bangalore * Colombo * London * New York * San José * Singapore * www.ambaresearch.com -----Original Message----- From: Ravi Varadhan [mailto:[EMAIL PROTECTED] Sent: Thursday, February 07, 2008 8:56 PM To: Shubha Vishwanath Karanth; [EMAIL PROTECTED] Subject: RE: [R] .fortran code If you look at the source code for eigen() you will notice that the .Fortran calls a subroutine named either "rs.f" or "rg.f" depending on whether your matrix is real and symmetric or real and asymmetric, respectively. So, you can look at these files if you had installed the source code. Even if you had not installed the full source, you can google using keywrods "rs.f eigenvalue" or "rg.f eigenvalue" to see the Fortran code for these functions available from EISPACK. However, I am not sure how illuminating that would be since rs.f and rg.f are just calling routines that each calls a number of other subroutines. So, unless you are willing to spend a great deal of time, you may not be any wiser at the end of all this. Ravi. ---------------------------------------------------------------------------- ------- Ravi Varadhan, Ph.D. Assistant Professor, The Center on Aging and Health Division of Geriatric Medicine and Gerontology Johns Hopkins University Ph: (410) 502-2619 Fax: (410) 614-9625 Email: [EMAIL PROTECTED] Webpage: http://www.jhsph.edu/agingandhealth/People/Faculty/Varadhan.html ---------------------------------------------------------------------------- -------- -----Original Message----- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Shubha Vishwanath Karanth Sent: Thursday, February 07, 2008 9:42 AM To: [EMAIL PROTECTED] Subject: [R] .fortran code Hi, I need to look/understand what the ".fortran()" doing in say, the source code of the "eigen" command. How do I look into this? Thanks, Shubha This e-mail may contain confidential and/or privileged i...{{dropped:13}} ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. This e-mail may contain confidential and/or privileged information. If you are not the intended recipient (or have received this e-mail in error) please notify the sender immediately and destroy this e-mail. Any unauthorized copying, disclosure or distribution of the material in this e-mail is strictly forbidden. Any views or opinions presented are solely those of the author and do not necessarily represent those of Amba Holdings Inc., and/or its affiliates. Important additional terms relating to this email can be obtained at http://www.ambaresearch.com/disclaimer
______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
