home *** CD-ROM | disk | FTP | other *** search
- lable by anonymous ftp from some ftp sites (try archie to discover
- them). The files are usually called memmove.c or memset.c.
- Some ftp sites which currently (Jan '94) have a version of these files are
- munnari.oz.au (in Australia), ftp.uu.net, gatekeeper.dec.com (USA), and
- unix.hensa.ac.uk (in the UK). The directory in which you will find
- memmove.c and memset.c typically looks like .../bsd-sources/lib/libc/...
-
- There are two further machine-dependent quantities that should be set.
- These are machine epsilon or the unit roundoff for double precision
- arithmetic, and the maximum value produced by the rand() routine, which is
- used in rand_vec() and rand_mat().
-
-
- The current definitions of these are
-
- #define MACHEPS 2.2e-16
- #define MAX_RAND 2.147483648e9
-
- The value of MACHEPS should be correct for all IEEE standard double
- precision arithmetic.
-
- However, ANSI C's <float.h> contains #define'd quantities DBL_EPSILON
- and RAND_MAX, so if you have an ANSI C compiler and headers, replace the
- above two lines of machine.h with
-
- #include <float.h>
- /* for Real == float */
- #define MACHEPS DBL_EPSILON
- #define MAX_RAND RAND_MAX
-
- The default value given for MAX_RAND is 2^31 , as the Pyramid 9810 and
- the SPARC 2's both have 32 bit words. There is a program macheps.c which
- is included in your source files which computes and prints out the value of
- MACHEPS for your machine.
-
- Some other macros control some aspects of Meschach. One of these is
- SEGMENTED which should be #define'd if you are working with a machine or
- compiler that does not allow large arrays to be allocated. For example,
- the most common memory models for MS-DOS compilers do not allow more than
- 64Kbyte to be allocated in one block. This limits square matrices to be no
- more than 9090 . Inserting #define SEGMENTED 1 into machine.h will mean
- that matrices are allocated a row at a time.
-
-
-
- 4. SAMPLE TESTS
-
- There are several programs for checking Meschach called torture
- (source: torture.c) for the dense routines, sptort (source: sptort.c) for
- the sparse routines, ztorture (source ztorture.c) for a complex version of
- torture, memtort (source memtort.c) for memory allocation/deallocation,
- itertort (source itertort.c) for iterative methods, mfuntort (source
- mfuntort.c) for computing powers of dense matrices, iotort (source
- iotort.c) for I/O routines. These can be compiled using make by "make
- torture", "make sptort", etc. The programs are part of meschach0.shar.
-
-
- 5. OTHER PROBLEMS
-
- Meschach is not a commercial package, so we do not guarantee that
- everything will be perfect or will install smoothly. Inevitably there will
- be unforeseen problems. If you come across any bugs or inconsistencies, please
- let us know. If you need to modify the results of the configure script, or
- need to construct your own machine.h and makefile's, please send them to
- us. A number of people sent us the machine dependent files for Meschach 1.1,
- but with the use of configure, and the new information needed for version
- 1.2, these machine dependent files don't have quite the right information.
- Hopefully, though, they are redundant. Non-Unix platforms at present
- require ``manual'' installation. Because of the variety of platforms
- (MS-DOS, Macintosh, VAX/VMS, Prime, Amiga, Atari, ....) this is left up to
- the users of these platforms. We hope that you can use the dâ2tiê
- 0 á ™â -é „å€àŴ‘íŴ-í ‘è ½è ëЙâ0 á -éŴ‡íŴ-íŴ‡íŴ-í â‘ á á á ½è ëЙâ Uã
- 0 á -éŴ‡íŴ-íŴ‡íŴ-í á á á ½è ëЙâ0 á -éŴ‡íŴ-íŴ‡íŴ-í á á á ½è ëЙâpGâ WãÊÿÿª „å Ŵ‘í ã á ê“å ’å€Ŵà(0å0›å€1◰à ’›í€!î ’’íŴî ’Ŵí €â“å Qáðÿÿ⇩p á* ê
- 0 á ™â -é “å`â‘ á‘çŴà Ŵ’íŴ-í‘ç€à ‘è ½è ëЙâ0 á -éŴ‡íŴ-íŴ‡íŴ-í á á á ½è rror messages etc., will use the new names. Also note that the new
- iterative routines have a very different calling sequence. If you need the
- old iterative routines, they are in oldmeschach.shar.
-
- If you wish to let us know what you have done, etc., our email
- addresses are
-
- david.stewart@anu.edu.au
- zbigniew.leyk@anu.edu.au
-
- Good luck!
-
-
- ACKNOWLEDGMENTS
-
-
- Many people have helped in various ways with ideas and suggestions.
- Needless to say, the bugs are all ours! But these people should be thanked
- for their encouragement etc. These include a number of people at
- University of Queensland: Graeme Chandler, David De Wit, Martin Sharry,
- Michael Forbes, Phil Kilby, John Holt, Phil Pollett and Tony Watts. At the
- Australian National University: Mike Osborne, Steve Roberts, Margaret Kahn
- and Teresa Leyk. Karen George of the University of Canberra has been a
- source of both ideas and encouragement. Email has become significant part
- of work, and many people have pointed out bugs, inconsistencies and
- improvements to Meschach by email. These people include Ajay Shah of the
- University of Southern California, Dov Grobgeld of the Weizmann Institute,
- John Edstrom of the University of Calgary, Eric Grosse, one of the netlib
- organisers, Ole Saether of Oslo, Norway, Alfred Thiele and Pierre
- Asselin of Carnegie-Mellon Univeristy, Daniel Polani of the University of
- Mainz, Marian Slodicka of Slovakia, Kaifu Wu of Pomona, Hidetoshi
- Shimodaira of the University of Tokyo, Eng Siong of Edinburgh, Hirokawa Rui
- of the University of Tokyo, Marko Slyz of the University of Michigan, and
- Brook Milligan of the University of Texas. This list is only partial, and
- there are many others who have corresponded with us on details about
- Meschach and the like. Finally our thanks go to all those that have had to
- struggle with compilers and other things to get Meschach to work.
-
-
-
-
-
- FileDataŵcts ³
-
- u Eýÿÿ2õ ®L÷◰‹À; è TeX output 1993.08.10:1622⇧ ÿÿÿÿ ©R³™ ýbM ☓2™ ýflQ2™‹ ºÎ⇦ó#ÖÆ áH G® cmss17ÎMeschach‰fl+°Ç™‹ ‡È¢ó)ŵÎR6 ff cmss12ÔMatrix„³/Computations›in›C‰ ®C™™™™™‹ ¥jhfi²×óKñ`y