| -∞ | NAN | yes
+ * )
+ *
+ * References:
+ * cephes, $(LINK http://en.wikipedia.org/wiki/Gamma_function)
*/
+/* Documentation prepared by Don Clugston */
real tgamma(real x)
{
- version (linux)
- return std.c.math.tgammal(x);
- else
- throw new NotImplemented("tgamma");
+ return std.c.math.tgammal(x);
+
+ // Use etc.gamma.tgamma for those C systems that are missing it
}
/**************************************
@@ -1561,3 +1595,70 @@
assert(feqrel(-real.max,real.infinity)==0);
assert(feqrel(real.max,-real.max)==0);
}
+
+
+/***********************************
+ * Evaluate polynomial A(x) = a0 + a1x + a2x² + a3x³ ...
+ *
+ * Uses Horner's rule A(x) = a0 + x(a1 + x(a2 + x(a3 + ...)))
+ * Params:
+ * A = array of coefficients a0, a1, etc.
+ */
+real poly(real x, real[] A)
+in
+{
+ assert(A.length > 0);
+}
+body
+{
+ version (D_InlineAsm)
+ {
+ asm // assembler by W. Bright
+ {
+ // EDX = (A.length - 1) * real.sizeof
+ mov ECX,A[EBP] ; // ECX = A.length
+ dec ECX ;
+ lea EDX,[ECX][ECX*8] ;
+ add EDX,ECX ;
+ add EDX,A+4[EBP] ;
+ fld real ptr [EDX] ; // ST0 = coeff[ECX]
+ jecxz return_ST ;
+ fld x[EBP] ; // ST0 = x
+ fxch ST(1) ; // ST1 = x, ST0 = r
+ align 4 ;
+ L2: fmul ST,ST(1) ; // r *= x
+ fld real ptr -10[EDX] ;
+ sub EDX,10 ; // deg--
+ faddp ST(1),ST ;
+ dec ECX ;
+ jne L2 ;
+ fxch ST(1) ; // ST1 = r, ST0 = x
+ fstp ST(0) ; // dump x
+ align 4 ;
+ return_ST: ;
+ ;
+ }
+ }
+ else
+ {
+ i = A.length - 1;
+ real r = A[i];
+ while (--i >= 0)
+ {
+ r *= x;
+ r += A[i];
+ }
+ return r;
+ }
+}
+
+unittest
+{
+ debug (math) printf("math.poly.unittest\n");
+ real x = 3.1;
+ static real pp[] = [56.1, 32.7, 6];
+
+ assert( poly(x, pp) == (56.1L + (32.7L + 6L * x) * x) );
+}
+
+
diff -uNr dmd-0.135/dmd/src/phobos/std/md5.d dmd-0.136/dmd/src/phobos/std/md5.d
--- dmd-0.135/dmd/src/phobos/std/md5.d 2005-10-03 13:02:50.000000000 +0200
+++ dmd-0.136/dmd/src/phobos/std/md5.d 2005-10-17 00:40:30.000000000 +0200
@@ -2,6 +2,69 @@
* Derived from the RSA Data Security, Inc. MD5 Message-Digest Algorithm.
*/
+/**
+ * Computes MD5 digests of arbitrary data. MD5 digests are 16 byte quantities that are like a checksum or crc, but are more robust.
+ *
+ * There are two ways to do this. The first does it all in one function call to
+ * sum(). The second is for when the data is buffered.
+ *
+ * Bugs:
+ * MD5 digests have been demonstrated to not be unique.
+ *
+ * Author:
+ * The routines and algorithms are derived from the
+ * $(I RSA Data Security, Inc. MD5 Message-Digest Algorithm).
+ *
+ * References:
+ * $(LINK2 http://en.wikipedia.org/wiki/Md5, Wikipedia on MD5)
+ *
+ * Macros:
+ * WIKI = StdMd5
+ */
+
+/++++++++++++++++++++++++++++++++
+ Example:
+
+--------------------
+// This code is derived from the
+// RSA Data Security, Inc. MD5 Message-Digest Algorithm.
+
+import std.md5;
+import std.string;
+import std.c.stdio;
+import std.stdio;
+
+int main(char[][] args)
+{
+ foreach (char[] arg; args)
+ MDFile(arg);
+ return 0;
+}
+
+/* Digests a file and prints the result. */
+void MDFile(char[] filename)
+{
+ FILE* file;
+ MD5_CTX context;
+ int len;
+ ubyte[4 * 1024] buffer;
+ ubyte digest[16];
+
+ if ((file = fopen(std.string.toStringz(filename), "rb")) == null)
+ writefln("%s can't be opened", filename);
+ else
+ {
+ context.start();
+ while ((len = fread(buffer, 1, buffer.size, file)) != 0)
+ context.update(buffer[0 .. len]);
+ context.finish(digest);
+ fclose(file);
+
+ writefln("MD5 (%s) = %s", filename, digestToString(digest));
+ }
+}
+--------------------
+ +/
/* Copyright (C) 1991-2, RSA Data Security, Inc. Created 1991. All
rights reserved.
@@ -30,7 +93,51 @@
import std.string;
-/* MD5 context. */
+/***************************************
+ * Computes MD5 digest of array of data.
+ */
+
+void sum(ubyte[16] digest, void[] data)
+{
+ MD5_CTX context;
+
+ context.start();
+ context.update(data);
+ context.finish(digest);
+}
+
+/******************
+ * Prints a message digest in hexadecimal to stdout.
+ */
+void printDigest(ubyte digest[16])
+{
+ foreach (ubyte u; digest)
+ printf("%02x", u);
+}
+
+/****************************************
+ * Converts MD5 digest to a string.
+ */
+
+char[] digestToString(ubyte[16] digest)
+{
+ char[] result = new char[32];
+ int i;
+
+ foreach (ubyte u; digest)
+ {
+ result[i] = std.string.hexdigits[u >> 4];
+ result[i + 1] = std.string.hexdigits[u & 15];
+ i += 2;
+ }
+ return result;
+}
+
+/**
+ * Holds context of MD5 computation.
+ *
+ * Used when data to be digested is buffered.
+ */
struct MD5_CTX
{
uint state[4] = /* state (ABCD) */
@@ -108,14 +215,15 @@
a += b;
}
- /* MD5 initialization. Begins an MD5 operation, writing a new context.
+ /**
+ * MD5 initialization. Begins an MD5 operation, writing a new context.
*/
void start()
{
*this = MD5_CTX.init;
}
- /* MD5 block update operation. Continues an MD5 message-digest
+ /** MD5 block update operation. Continues an MD5 message-digest
operation, processing another message block, and updating the
context.
*/
@@ -151,8 +259,8 @@
memcpy(&buffer[index], &input[i], inputLen-i);
}
- /* MD5 finalization. Ends an MD5 message-digest operation, writing the
- * the message digest and zeroizing the context.
+ /** MD5 finalization. Ends an MD5 message-digest operation, writing the
+ * the message to digest and zeroing the context.
*/
void finish(ubyte[16] digest) /* message digest */
{
@@ -334,27 +442,6 @@
}
}
-/*************************************
- * Computes MD5 digest of array of data
- */
-
-void sum(ubyte[16] digest, void[] data)
-{
- MD5_CTX context;
-
- context.start();
- context.update(data);
- context.finish(digest);
-}
-
-/* Prints a message digest in hexadecimal.
- */
-void printDigest(ubyte digest[16])
-{
- foreach (ubyte u; digest)
- printf("%02x", u);
-}
-
unittest
{
debug(md5) printf("std.md5.unittest\n");
@@ -384,5 +471,7 @@
"1234567890123456789012345678901234567890");
assert(digest == cast(ubyte[])x"57edf4a22be3c955ac49da2e2107b67a");
+ assert(digestToString(cast(ubyte[16])x"c3fcd3d76192e4007dfb496cca67e13b")
+ == "C3FCD3D76192E4007DFB496CCA67E13B");
}
diff -uNr dmd-0.135/dmd/src/phobos/std/zip.d dmd-0.136/dmd/src/phobos/std/zip.d
--- dmd-0.135/dmd/src/phobos/std/zip.d 2005-10-03 13:02:50.000000000 +0200
+++ dmd-0.136/dmd/src/phobos/std/zip.d 2005-10-17 00:40:30.000000000 +0200
@@ -1,4 +1,20 @@
+/**
+ * Read/write data in the $(LINK2 http://www.info-_zip.org, zip archive) format.
+ * Makes use of the etc.c.zlib compression library.
+ *
+ * Bugs:
+ * $(UL
+ * $(LI Multi-disk zips not supported.)
+ * $(LI Only Zip version 20 formats are supported.)
+ * $(LI Only supports compression modes 0 (no compression) and 8 (deflate).)
+ * $(LI Does not support encryption.)
+ * )
+ *
+ * Macros:
+ * WIKI = StdZip
+ */
+
module std.zip;
private import std.zlib;
@@ -7,9 +23,8 @@
//debug=print;
-/* Throw this on errors.
+/** Thrown on error.
*/
-
class ZipException : Exception
{
this(char[] msg)
@@ -18,27 +33,36 @@
}
}
+/**
+ * A member of the ZipArchive.
+ */
class ArchiveMember
{
- ushort madeVersion = 20;
- ushort extractVersion = 20;
- ushort flags;
- ushort compressionMethod;
- std.date.DosFileTime time;
- uint crc32;
- uint compressedSize;
- uint expandedSize;
- ushort diskNumber;
- ushort internalAttributes;
- uint externalAttributes;
+ ushort madeVersion = 20; /// Read Only
+ ushort extractVersion = 20; /// Read Only
+ ushort flags; /// Read/Write: normally set to 0
+ ushort compressionMethod; /// Read/Write: 0 for compression, 8 for deflate
+ std.date.DosFileTime time; /// Read/Write: Last modified time of the member. It's in the DOS date/time format.
+ uint crc32; /// Read Only: cyclic redundancy check (CRC) value
+ uint compressedSize; /// Read Only: size of data of member in compressed form.
+ uint expandedSize; /// Read Only: size of data of member in expanded form.
+ ushort diskNumber; /// Read Only: should be 0.
+ ushort internalAttributes; /// Read/Write
+ uint externalAttributes; /// Read/Write
private uint offset;
+ /**
+ * Read/Write: Usually the file name of the archive member; it is used to
+ * index the archive directory for the member. Each member must have a unique
+ * name[]. Do not change without removing member from the directory first.
+ */
char[] name;
- ubyte[] extra;
- char[] comment;
- ubyte[] compressedData;
- ubyte[] expandedData;
+
+ ubyte[] extra; /// Read/Write: extra data for this member.
+ char[] comment; /// Read/Write: comment associated with this member.
+ ubyte[] compressedData; /// Read Only: data of member in compressed form.
+ ubyte[] expandedData; /// Read/Write: data of member in uncompressed form.
debug(print)
{
@@ -60,16 +84,33 @@
}
}
+/**
+ * Object representing the entire archive.
+ * ZipArchives are collections of ArchiveMembers.
+ */
class ZipArchive
{
- ubyte[] data; // representing the entire Zip contents
+ ubyte[] data; /// Read Only: array representing the entire contents of the archive.
uint endrecOffset;
- uint diskNumber;
- uint diskStartDir;
- uint numEntries;
- uint totalEntries;
- char[] comment;
+ uint diskNumber; /// Read Only: 0 since multi-disk zip archives are not supported.
+ uint diskStartDir; /// Read Only: 0 since multi-disk zip archives are not supported.
+ uint numEntries; /// Read Only: number of ArchiveMembers in the directory.
+ uint totalEntries; /// Read Only: same as totalEntries.
+ char[] comment; /// Read/Write: the archive comment. Must be less than 65536 bytes in length.
+
+ /**
+ * Read Only: array indexed by the name of each member of the archive.
+ * Example:
+ * All the members of the archive can be accessed with a foreach loop:
+ * --------------------
+ * ZipArchive archive = new ZipArchive(data);
+ * foreach (ArchiveMember am; archive.directory)
+ * {
+ * writefln("member name is '%s'", am.name);
+ * }
+ * --------------------
+ */
ArchiveMember[char[]] directory;
debug (print)
@@ -86,23 +127,36 @@
/* ============ Creating a new archive =================== */
- /* Constructor to use when creating a new archive.
+ /** Constructor to use when creating a new archive.
*/
-
this()
{
}
+ /** Add de to the archive.
+ */
void addMember(ArchiveMember de)
{
directory[de.name] = de;
}
+ /** Delete de from the archive.
+ */
void deleteMember(ArchiveMember de)
{
directory.remove(de.name);
}
+ /**
+ * Construct an archive out of the current members of the archive.
+ *
+ * Fills in the properties data[], diskNumber, diskStartDir, numEntries,
+ * totalEntries, and directory[].
+ * For each ArchiveMember, fills in properties crc32, compressedSize,
+ * compressedData[].
+ *
+ * Returns: array representing the entire archive.
+ */
void[] build()
{ uint i;
uint directoryOffset;
@@ -224,7 +278,19 @@
/* ============ Reading an existing archive =================== */
- /* Constructor to use when reading an existing archive.
+ /**
+ * Constructor to use when reading an existing archive.
+ *
+ * Fills in the properties data[], diskNumber, diskStartDir, numEntries,
+ * totalEntries, comment[], and directory[].
+ * For each ArchiveMember, fills in
+ * properties madeVersion, extractVersion, flags, compressionMethod, time,
+ * crc32, compressedSize, expandedSize, compressedData[], diskNumber,
+ * internalAttributes, externalAttributes, name[], extra[], comment[].
+ * Use expand() to get the expanded data for each ArchiveMember.
+ *
+ * Params:
+ * buffer = the entire contents of the archive.
*/
this(void[] buffer)
@@ -324,6 +390,13 @@
throw new ZipException("invalid directory entry 3");
}
+ /*****
+ * Decompress the contents of archive member de and return the expanded
+ * data.
+ *
+ * Fills in properties extractVersion, flags, compressionMethod, time,
+ * crc32, compressedSize, expandedSize, expandedData[], name[], extra[].
+ */
ubyte[] expand(ArchiveMember de)
{ uint namelen;
uint extralen;
diff -uNr dmd-0.135/dmd/src/phobos/std/zlib.d dmd-0.136/dmd/src/phobos/std/zlib.d
--- dmd-0.135/dmd/src/phobos/std/zlib.d 2005-10-03 13:02:50.000000000 +0200
+++ dmd-0.136/dmd/src/phobos/std/zlib.d 2005-10-17 00:40:30.000000000 +0200
@@ -1,5 +1,11 @@
-/* zlib.d
- * A D interface to the zlib library, www.gzip.org/zlib
+/**
+ * Compress/decompress data using the $(LINK2 http://www._zlib.net, zlib library).
+ *
+ * References:
+ * $(LINK2 http://en.wikipedia.org/wiki/Zlib, Wikipedia)
+ *
+ * Macros:
+ * WIKI = StdZlib
*/
@@ -45,6 +51,8 @@
}
/**************************************************
+ * Compute the Adler32 checksum of the data in buf[]. adler is the starting
+ * value when computing a cumulative checksum.
*/
uint adler32(uint adler, void[] buf)
@@ -65,6 +73,8 @@
}
/*********************************
+ * Compute the CRC32 checksum of the data in buf[]. crc is the starting value
+ * when computing a cumulative checksum.
*/
uint crc32(uint crc, void[] buf)
@@ -85,6 +95,11 @@
}
/*********************************************
+ * Compresses the data in srcbuf[] using compression _level level.
+ * The default value
+ * for level is 6, legal values are 1..9, with 1 being the least compression
+ * and 9 being the most.
+ * Returns the compressed data.
*/
void[] compress(void[] srcbuf, int level)
@@ -111,6 +126,7 @@
}
/*********************************************
+ * ditto
*/
void[] compress(void[] buf)
@@ -119,19 +135,14 @@
}
/*********************************************
+ * Decompresses the data in srcbuf[].
+ * Params: destlen = size of the uncompressed data.
+ * It need not be accurate, but the decompression will be faster if the exact
+ * size is supplied.
+ * Returns: the decompressed data.
*/
-void[] uncompress(void[] srcbuf)
-{
- return uncompress(srcbuf, 0, 15);
-}
-
-void[] uncompress(void[] srcbuf, uint destlen)
-{
- return uncompress(srcbuf, destlen, 15);
-}
-
-void[] uncompress(void[] srcbuf, uint destlen, int winbits)
+void[] uncompress(void[] srcbuf, uint destlen = 0u, int winbits = 15)
{
int err;
void[] destbuf;
@@ -212,6 +223,7 @@
+/
/*********************************************
+ * Used when the data to be compressed is not all in one buffer.
*/
class Compress
@@ -231,10 +243,10 @@
}
public:
- this()
- {
- }
+ /**
+ * Construct. level is the same as for D.zlib.compress().
+ */
this(int level)
in
{
@@ -245,6 +257,11 @@
this.level = level;
}
+ /// ditto
+ this()
+ {
+ }
+
~this()
{ int err;
@@ -257,6 +274,11 @@
}
}
+ /**
+ * Compress the data in buf and return the compressed data.
+ * The buffers
+ * returned from successive calls to this should be concatenated together.
+ */
void[] compress(void[] buf)
{ int err;
void[] destbuf;
@@ -291,12 +313,25 @@
return destbuf;
}
- void[] flush()
- {
- return flush(Z_FINISH);
- }
-
- void[] flush(int mode)
+ /***
+ * Compress and return any remaining data.
+ * The returned data should be appended to that returned by compress().
+ * Params:
+ * mode = one of the following:
+ * $(DL
+ $(DT Z_SYNC_FLUSH )
+ $(DD Syncs up flushing to the next byte boundary.
+ Used when more data is to be compressed later on.)
+ $(DT Z_FULL_FLUSH )
+ $(DD Syncs up flushing to the next byte boundary.
+ Used when more data is to be compressed later on,
+ and the decompressor needs to be restartable at this
+ point.)
+ $(DT Z_FINISH)
+ $(DD (default) Used when finished compressing the data. )
+ )
+ */
+ void[] flush(int mode = Z_FINISH)
in
{
assert(mode == Z_FINISH || mode == Z_SYNC_FLUSH || mode == Z_FULL_FLUSH);
@@ -333,6 +368,10 @@
}
}
+/******
+ * Used when the data to be decompressed is not all in one buffer.
+ */
+
class UnCompress
{
private:
@@ -351,15 +390,20 @@
}
public:
- this()
- {
- }
+ /**
+ * Construct. destbufsize is the same as for D.zlib.uncompress().
+ */
this(uint destbufsize)
{
this.destbufsize = destbufsize;
}
+ /** ditto */
+ this()
+ {
+ }
+
~this()
{ int err;
@@ -373,6 +417,11 @@
done = 1;
}
+ /**
+ * Decompress the data in buf and return the decompressed data.
+ * The buffers returned from successive calls to this should be concatenated
+ * together.
+ */
void[] uncompress(void[] buf)
in
{
@@ -414,6 +463,11 @@
return destbuf;
}
+ /**
+ * Decompress and return any remaining data.
+ * The returned data should be appended to that returned by uncompress().
+ * The UnCompress object cannot be used further.
+ */
void[] flush()
in
{
diff -uNr dmd-0.135/dmd/src/phobos/std.ddoc dmd-0.136/dmd/src/phobos/std.ddoc
--- dmd-0.135/dmd/src/phobos/std.ddoc 2005-10-03 13:02:50.000000000 +0200
+++ dmd-0.136/dmd/src/phobos/std.ddoc 2005-10-17 00:40:30.000000000 +0200
@@ -16,16 +16,17 @@
-
+
-
-Last update $(DATETIME)
-
|
@@ -97,21 +98,18 @@
|
-$(TITLE)
+$(TITLE)
$(BODY)
|
-Feedback and Comments
- Add feedback and comments regarding this
- page.
+
+Copyright © 1999-$(YEAR) by Digital Mars, All Rights Reserved |
+Page generated by $(LINK2 http://www.digitalmars.com/d/ddoc.html, Ddoc).
+
-
-$(SMALL Copyright © 1999-$(YEAR) by Digital Mars, All Rights Reserved
-Page generated by $(LINK2 http://www.digitalmars.com/d/ddoc.html, Ddoc).)
-
|