home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #31 / NN_1992_31.iso / spool / sci / math / numanal / 3678 < prev    next >
Encoding:
Text File  |  1992-12-25  |  2.0 KB  |  52 lines
  1. Newsgroups: sci.math.num-analysis
  2. Path: sparky!uunet!cs.utexas.edu!torn!skule.ecf!drill.me!steinman
  3. From: steinman@me.utoronto.ca (David Steinman)
  4. Subject: Re: LU decomposition on very large matrix.
  5. Message-ID: <BztzMq.59z@me.utoronto.ca>
  6. Sender: news@me.utoronto.ca (News Reader)
  7. Organization: UofT Mechanical Engineering
  8. References: <HCOBB.92Dec23143856@fly2.berkeley.edu>
  9. Distribution: inet
  10. Date: Fri, 25 Dec 1992 20:01:38 GMT
  11. Lines: 39
  12.  
  13. hcobb@fly2.berkeley.edu (Henry J. Cobb) writes:
  14.  
  15. >I am interested in C code for the manipulation of a VERY large matrix.
  16. >Basically, what I need to do is perform an LU decomposition on a matrix
  17. >too large to operate on the whole thing at once.  Therefore, I need an
  18. >algorithm that allows me to operate on chunks of the matrix at a time,
  19. >and which allows an in-place LU decomposition, so that I can minimize
  20. >the amount of disk storage necessary for the computations.  Any
  21. >help is GREATLY appreciated.
  22.  
  23. Assuming your matrix results from a finite analysis technique, you
  24. should look into an out-of-core skyline matrix solver.  The only
  25. caveat is that i/o time can be of the same order as solve time.  I have a
  26. reference to a published Fortran code, which could be messily ported to 
  27. C via f2c (PD fortran to C converter) or by hand (with lotsa goto's).
  28.  
  29.     AUTHOR = "G. Dhatt and G. Touzot",
  30.     TITLE = "The finite element method displayed",
  31.     ADDRESS = "Chichester",
  32.     PUBLISHER = "John Wiley \& Sons",
  33.     YEAR = "1984"}
  34.  
  35. If you go with a skyline solver, you should also consider Sloan's
  36. algorithm, which reorders nodes to minimize the skyline.
  37.  
  38.     AUTHOR = "S. W. Sloan",
  39.     TITLE = "A fortran program for profile and wavefront reduction",
  40.     JOURNAL = "Int J Num Meth Eng",
  41.     YEAR = "1989",
  42.     VOLUME = "28",
  43.     PAGES = "2651--2679"}
  44.  
  45. I believe the code provided is in Fortran, but you should be able to 
  46. C it clearly (or, again, use f2c).
  47.  
  48. If you hear about any C-codes for the above, please let me know.  I'm
  49. not thrilled about being Fortran-dependent for such things.
  50. --
  51. Dave!
  52.