paweld/fpc
1
0
Fork 0
mirror of https://gitlab.com/freepascal.org/fpc/source.git synced 2025-08-13 12:19:18 +02:00
Code Issues Packages Projects Releases Wiki Activity

* small changes

Browse Source
This commit is contained in:
marco 2000-01-25 20:22:31 +00:00
parent 6a7e4fcb89
commit 5f4ae3e48d
1 changed files with 115 additions and 109 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

224
packages/numlib/numlib.txt
View File

@ -1,111 +1,117 @@
{ {
$Id$ $Id$
This file is part of the Numlib package. This file is part of the Numlib package.
Copyright (c) 1986-2000 by Copyright (c) 1986-2000 by
Kees van Ginneken, Wil Kortsmit and Loek van Reij of the Kees van Ginneken, Wil Kortsmit and Loek van Reij of the
Computational centre of the Eindhoven University of Technology Computational centre of the Eindhoven University of Technology
FPC port Code by Marco van de Voort (marco@freepascal.org) FPC port Code by Marco van de Voort (marco@freepascal.org)
Documentation by Michael van Canneyt (Michael@freepascal.org) Documentation by Michael van Canneyt (Michael@freepascal.org)
This is an internal document with information collected during porting This is an internal document with information collected during porting
numlib to FPC numlib to FPC
See the file COPYING.FPC, included in this distribution, See the file COPYING.FPC, included in this distribution,
for details about the copyright. for details about the copyright.
This program is distributed in the hope that it will be useful, This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
**********************************************************************} **********************************************************************}
NumLib ported.txt, internals or "Developpers docs". NumLib ported.txt, internals or "Developpers docs".
ARBFLOAT, basic Floating point type. ARBFLOAT, basic Floating point type.
----------------------------------- -----------------------------------
In the FPC revision instead of picking a certain floating point type, In the FPC revision instead of picking a certain floating point type,
a new type "ArbFloat" is defined which is used as floating point type a new type "ArbFloat" is defined which is used as floating point type
throughout the entire library. If the floating point type is changed, throughout the entire library. If the floating point type is changed,
define or undefine ArbExtended and add the machineconstants change to define or undefine ArbExtended and add the machineconstants change to
the type selected. the type selected.
This allows IEEE Double (64bit) and Extended(80bit). This allows IEEE Double (64bit) and Extended(80bit).
ARBINT, basic INTEGER type. ARBINT, basic INTEGER type.
----------------------------------- -----------------------------------
Because in plain FPC mode Integer =16-bits (for TP compatibility), and in Because in plain FPC mode Integer =16-bits (for TP compatibility), and in
Delphi 32-bits, I changed all integers to ArbInt. Delphi 32-bits, I changed all integers to ArbInt.
The basic idea is the same as ArbFloat, but it is less consequently used, The basic idea is the same as ArbFloat, but it is less consequently used,
mainly because some typecastings of pointers to words existed. These mainly because some typecastings of pointers to words existed. These
typecastings should never be 16-bits in FPC, so all local variables are typecastings should never be 16-bits in FPC, so all local variables are
longint. (which is currently always 32bits) longint. (which is currently always 32bits)
VECTOR or MATRIX as ArbFloat. VECTOR or MATRIX as ArbFloat.
----------------------------------- -----------------------------------
NumLib often passes Matrices and Vectors as one ArbFloat + some integer NumLib often passes Matrices and Vectors as one ArbFloat + some integer
values, then maps the following pmatrix type over it, and accesses it as an values, then maps the following pmatrix type over it, and accesses it as an
array or vector: array or vector:
procedure dosomething(var invalue:ArbFloat); procedure dosomething(var invalue:ArbFloat);
type Row=ARRAY[0..maxelements] OF ArbFloat; type Row=ARRAY[0..maxelements] OF ArbFloat;
Matrix=Array[0..maxelements] OF ^ROW; Matrix=Array[0..maxelements] OF ^ROW;
pmatrix=^matrix; pmatrix=^matrix;
Var pa : pmatrix; Var pa : pmatrix;
begin begin
pa=@invalue; pa=@invalue;
pa[x]^[y]:=valuexy pa[x]^[y]:=valuexy
END; END;
The calling side looks like this: The calling side looks like this:
VAR L : ARRAY[0..1999] OF ArbFloat; VAR L : ARRAY[0..1999] OF ArbFloat;
DoSomething(L[0]); DoSomething(L[0]);
----------- -----------
Questions that remain open/incompleteneses in the package we got: Questions that remain open/incompleteneses in the package we got:
- Typ, mdt and Dsl,spl are undocumented. Typ is quite understandable though. - Typ, mdt and Dsl,spl are undocumented. Typ is quite understandable though.
Mdt and dsl contain probably procedures that, in earlier version were Mdt and dsl contain probably procedures that, in earlier version were
used as locals in some unit. When the procedures were also used in other units, used as locals in some unit. When the procedures were also used in other units,
they were moved to a different unit, but the documentation wasn't extended. they were moved to a different unit, but the documentation wasn't extended.
SPL is different in many ways. Contains comments (including some english SPL is different in many ways. Contains comments (including some english
ones) ones)
- All procedures with an extra l appended to the name in unit SLE are - All procedures with an extra l appended to the name in unit SLE are
undocumented. undocumented.
- The archive we got seems to be a copy of the working directory of the author, - The archive we got seems to be a copy of the working directory of the author,
a snapshot during a never finished restructure. (probably matlab and a snapshot during a never finished restructure. (probably matlab and
similar programs took over) similar programs took over)
The sources/finished directory was the "new" archive, which generated a The sources/finished directory was the "new" archive, which generated a
.dll. All graphics using routines, and new units were never finished. .dll. All graphics using routines, and new units were never finished.
The problem is that also the documentation was never finished. The problem is that also the documentation was never finished.
- How to implement a less ugly calling convention, without loosing speed. - How to implement a less ugly calling convention, without loosing speed.
(which still can be important in nummerics)? (which still can be important in nummerics)?
Other remarks: Other remarks:
- Spe needs some constants recalculated to get full precision for extended, - Spe needs some constants recalculated to get full precision for extended,
the files to calc the constants aren't included. (copied from some reference the files to calc the constants aren't included. (copied from some reference
book?) Some units have literature references in the documentation. Spe book?) Some units have literature references in the documentation. Spe
(Murphy's law) doesn't. (Murphy's law) doesn't. Some other units (INT) also have this problem.
- Found out what MDT does. MDT is the core of DET. Det just reads a vector,
------------- reformats it to form a matrix, and then passes it to MDT.
Some translation problems: Such a vector only contains the required fields to build a certain kind of
matrix. (e.g. for a bandmatrix, the diagonals). MDT=Matrix determinant.
I left the term "Bandmatrix" untranslated, and give the mathematical definition
here, and hope you know what it is in english: -------------
Some translation problems:
If A is a "n x n" bandmatrix with leftbound l, and rightbound r then
I left the term "Bandmatrix" untranslated, and give the mathematical definition
Aij=0 if j<i-l or j>i+r here, and hope you know what it is in english:
{
If A is a "n x n" bandmatrix with leftbound l, and rightbound r then
Aij=0 if j<i-l or j>i+r
{
$Log$ $Log$
Revision 1.1 2000-01-24 22:11:28 marco Revision 1.2 2000-01-25 20:22:31 marco
* small changes
Revision 1.1 2000/01/24 22:11:28 marco
* initial version * initial version
}
}