mirror of
https://github.com/QB64Official/qb64.git
synced 2024-09-28 11:17:47 +00:00
930 lines
24 KiB
Markdown
930 lines
24 KiB
Markdown
|
**QB64** uses C++ to compile a BAS program into an executable program. The following is a list of the **C [FUNCTION](FUNCTION)s** that can be used.
|
||
|
|
||
|
**Note: C++ Header files should be placed in the QB64 folder and are not required after a program is compiled.**
|
||
|
|
||
|
|
||
|
## **C++ Variable Types**
|
||
|
|
||
|
|
||
|
| C Name | Description | Size | QB64 Type |
|
||
|
| ----------- | ----------------------- | ----------------- | ------------ |
|
||
|
|[_BYTE](_BYTE) | Character or small integer. | 1 byte | [_BYTE](_BYTE) |
|
||
|
|[INTEGER](INTEGER) | Short Integer(Word) | 2 byte | [INTEGER](INTEGER)|
|
||
|
|[LONG](LONG) | Integer(Dword) | 4 byte | [LONG](LONG) |
|
||
|
|[LONG](LONG) | Int32, Long integer or Long | 4 byte | [LONG](LONG) |
|
||
|
|[_INTEGER64](_INTEGER64) | Long long (Qword) | 8 byte | [_INTEGER64](_INTEGER64) |
|
||
|
|[Boolean](Boolean) | Boolean value true or false. | 1 byte | [_BYTE](_BYTE) |
|
||
|
|[SINGLE](SINGLE) | Floating point number | 4 byte | [SINGLE](SINGLE) |
|
||
|
|[DOUBLE](DOUBLE) | Double precision floating. | 8 byte | [DOUBLE](DOUBLE) |
|
||
|
|[_FLOAT](_FLOAT) | Long double precision float | 10 byte | [_FLOAT](_FLOAT) |
|
||
|
|[Unicode](Unicode) | Wide character([Unicode](Unicode)) | 2 or 4 |
|
||
|
|[_OFFSET](_OFFSET) | void pointer(void *) | ANY | [_OFFSET](_OFFSET)|
|
||
|
|
||
|
|
||
|
## C Functions and Subs
|
||
|
|
||
|
```vb
|
||
|
|
||
|
DECLARE LIBRARY
|
||
|
**'ctime.h**
|
||
|
FUNCTION clock () 'arithmetic type elapsed processor representing time.
|
||
|
FUNCTION difftime# (BYVAL time2 AS _UNSIGNED LONG, BYVAL time1 AS _UNSIGNED LONG)
|
||
|
'seconds between time2 and time1
|
||
|
**'ctype.h**
|
||
|
FUNCTION isalnum% (BYVAL c AS INTEGER) 'is an alphabet letter(isalpha(c) or isdigit(c))
|
||
|
FUNCTION isalpha% (BYVAL c AS INTEGER) 'is letter (isupper(c) or islower(c))
|
||
|
FUNCTION isdigit% (BYVAL c AS INTEGER) 'is a decimal digit
|
||
|
FUNCTION isgraph% (BYVAL c AS INTEGER) 'is a printing character other than space
|
||
|
FUNCTION islower% (BYVAL c AS INTEGER) 'is a lower-case letter
|
||
|
FUNCTION isprint% (BYVAL c AS INTEGER) 'is printing character. ASCII: &H20 (" ") to &H7E (~)
|
||
|
FUNCTION ispunct% (BYVAL c AS INTEGER) 'is printing character other than space, letter, digit
|
||
|
FUNCTION isspace% (BYVAL c AS INTEGER) 'is space, formfeed, newline, return, tab, vertical tab
|
||
|
FUNCTION isupper% (BYVAL c AS INTEGER) 'is upper-case letter
|
||
|
FUNCTION isxdigit% (BYVAL c AS INTEGER)'is a hexdecimal digit character(0 thru 9 or A thru F)
|
||
|
FUNCTION tolower% (BYVAL c AS INTEGER) 'return lower-case equivalent
|
||
|
FUNCTION toupper% (BYVAL c AS INTEGER) 'return upper-case equivalent
|
||
|
|
||
|
**'math.h**
|
||
|
FUNCTION acos# (BYVAL x AS DOUBLE)
|
||
|
FUNCTION asin# (BYVAL x AS DOUBLE)
|
||
|
FUNCTION atan# (BYVAL x AS DOUBLE) 'arc tangent of x does not designate the quadrant
|
||
|
FUNCTION atan2# (BYVAL y AS DOUBLE, BYVAL x AS DOUBLE) 'arc-tangent of y/x designates quadrant
|
||
|
FUNCTION cosh# (BYVAL x AS DOUBLE)
|
||
|
FUNCTION ldexp# (BYVAL base AS DOUBLE, BYVAL exponent AS INTEGER) 'base times 2 to exponent
|
||
|
FUNCTION pow# (BYVAL base AS DOUBLE, BYVAL exponent AS DOUBLE) 'base number raised to exponent
|
||
|
FUNCTION sinh# (BYVAL x AS DOUBLE)
|
||
|
FUNCTION tanh# (BYVAL x AS DOUBLE)
|
||
|
|
||
|
**'stdio.h** return non-zero INTEGER on failure.
|
||
|
FUNCTION remove% (filename AS STRING) 'removes specified file
|
||
|
FUNCTION rename% (oldname AS STRING, newname AS STRING) 'renames file oldname to newname.
|
||
|
|
||
|
**'stdlib.h**
|
||
|
FUNCTION atol& (str AS STRING) 'convert string to Long (errno not necessarily set)
|
||
|
FUNCTION atoi% (str AS STRING) 'convert string to Integer (errno not necessarily set)
|
||
|
FUNCTION atof# (str AS STRING) 'convert string to Double (errno not necessarily set)
|
||
|
FUNCTION rand& () 'random number
|
||
|
SUB srand (BYVAL seed AS _UNSIGNED LONG) 'random seeded number
|
||
|
|
||
|
**'string.h**
|
||
|
FUNCTION memchr& (BasePtr, value AS STRING, BYVAL Bytes AS LONG) ' returns pointer to first
|
||
|
'occurance of string in a set number of bytes at the memory block pointer designated.
|
||
|
|
||
|
FUNCTION memcmp (pointer1, pointer2, BYVAL Bytes AS _UNSIGNED LONG )
|
||
|
'compares pointer bytes. Returns 0 if match, positive if ptr1>ptr2, negative if ptr1<ptr2
|
||
|
|
||
|
SUB memcpy (DestPtr, SourcePtr, BYVAL Bytes AS _UNSIGNED LONG)
|
||
|
'copies the number of source bytes from source pointer to destination pointer
|
||
|
|
||
|
SUB memmove (DestPtr AS var_TYPE, SourcePtr AS var_TYPE, BYVAL bytesize AS LONG)
|
||
|
'moves a variable pointer value to destination from a source as a type byte size.
|
||
|
|
||
|
SUB memset (pointer AS **var_TYPE**, BYVAL value AS _UNSIGNED LONG, BYVAL nbytes AS _UNSIGNED LONG)
|
||
|
'The **var_TYPE** MUST match the type of the argument passed! Fills a block of memory. Sets the
|
||
|
'first bytes of the block of memory at pointer to the specified value as a character.
|
||
|
|
||
|
FUNCTION strcspn% (str1 AS STRING, str2 AS STRING)
|
||
|
'length of prefix of str1 consisting of characters not in str2.
|
||
|
|
||
|
FUNCTION strcmp% (str1 AS STRING, str2 AS STRING) 'compares str1 with str2, negative value if
|
||
|
'str1<str2, zero if str1=str2, positive if str1>str2
|
||
|
|
||
|
FUNCTION strncmp% (str1 AS STRING, str2 AS STRING, BYVAL Bytes AS INTEGER)'compares first byte
|
||
|
'characters of str1 and str2, negative if str1<str2, 0 if str1=str2, positive if str1>str2
|
||
|
|
||
|
SUB strncpy (dest AS STRING, source AS STRING, BYVAL Bytes AS INTEGER) 'Copies first
|
||
|
'bytes of source to destination. If source ends before number of bytes, dest padded with 0's
|
||
|
END DECLARE
|
||
|
|
||
|
```
|
||
|
|
||
|
How to use the [SUB](SUB) **memmove** to transfer [TYPE](TYPE) data when using Libraries.
|
||
|
|
||
|
```vb
|
||
|
|
||
|
TYPE a 'Note: the TYPE must be placed before the DECLARE LIBRARY if used in it!
|
||
|
b AS DOUBLE
|
||
|
c AS LONG
|
||
|
END TYPE
|
||
|
|
||
|
DECLARE LIBRARY
|
||
|
SUB memmove (Dest AS a, Source AS a, BYVAL bytesize AS LONG)
|
||
|
END DECLARE
|
||
|
|
||
|
|
||
|
DIM d(10) AS a
|
||
|
d(0).b = 1.5
|
||
|
d(0).c = 99
|
||
|
|
||
|
memmove d(10), d(0), LEN(d()) 'LEN gives the total byte size of the TYPE
|
||
|
|
||
|
PRINT d(10).b
|
||
|
PRINT d(10).c
|
||
|
|
||
|
```
|
||
|
|
||
|
```text
|
||
|
|
||
|
1.5
|
||
|
99
|
||
|
```
|
||
|
|
||
|
> *Explanation:* When a [TYPE](TYPE) variable is moved to another variable or array index, all [TYPE](TYPE) dot values are moved with it.
|
||
|
|
||
|
Creating different **memset** functions for each variable type to be used.
|
||
|
|
||
|
```vb
|
||
|
|
||
|
DECLARE LIBRARY
|
||
|
FUNCTION memsetB& ALIAS **memset** (p AS _BYTE, BYVAL c AS _UNSIGNED LONG, BYVAL n AS _UNSIGNED LONG)
|
||
|
FUNCTION memsetI& ALIAS **memset** (p AS INTEGER, BYVAL c AS _UNSIGNED LONG, BYVAL n AS _UNSIGNED LONG)
|
||
|
FUNCTION memsetL& ALIAS **memset** (p AS LONG, BYVAL c AS _UNSIGNED LONG, BYVAL n AS _UNSIGNED LONG)
|
||
|
END DECLARE
|
||
|
|
||
|
DIM Barray(1 TO 10) AS _BYTE
|
||
|
res& = memsetB(Barray(1), 65, 5)
|
||
|
|
||
|
FOR i = 1 TO 10
|
||
|
PRINT Barray(i); 'display BYTE array decimal values
|
||
|
NEXT
|
||
|
PRINT: PRINT
|
||
|
|
||
|
DIM Larray(1 TO 10) AS LONG
|
||
|
res& = memsetL(Larray(1), 65, 5)
|
||
|
|
||
|
FOR i = 1 TO 10
|
||
|
PRINT Larray(i); 'displays LONG array decimal values
|
||
|
NEXT i
|
||
|
PRINT
|
||
|
FOR i = 1 TO 10
|
||
|
PRINT " " + HEX$(Larray(i)); 'displays each byte value &H41 = 65
|
||
|
NEXT i
|
||
|
|
||
|
```
|
||
|
<sub>Code example by stylin</sub>
|
||
|
|
||
|
> *Explanation:* When 5 bytes are put into a 4 byte LONG array value, the fifth byte goes into the next array element.
|
||
|
|
||
|
## Bit Casting
|
||
|
|
||
|
Header file: *Cast.h*
|
||
|
|
||
|
```cpp
|
||
|
|
||
|
|
||
|
float bitcast(int t)
|
||
|
{
|
||
|
return *(float*)&t;
|
||
|
}
|
||
|
int swap_endian(unsigned int k)
|
||
|
{
|
||
|
return (k>>24)|((k<<8) & 0x00FF0000)|((k>>8) & 0x0000FF00)|(k<<24);
|
||
|
}
|
||
|
|
||
|
```
|
||
|
|
||
|
<sub>Courtesy of Darth Who</sub>
|
||
|
|
||
|
```vb
|
||
|
|
||
|
DECLARE LIBRARY "Cast"
|
||
|
FUNCTION bitcast## (BYVAL t AS LONG)
|
||
|
FUNCTION swap_endian (BYVAL k AS _UNSIGNED LONG)
|
||
|
END DECLARE
|
||
|
Value& = &HE7750340
|
||
|
PRINT HEX$(Value&)
|
||
|
PRINT bitcast##(Value&)
|
||
|
PRINT bitcast##(swap_endian(Value&))
|
||
|
|
||
|
```
|
||
|
|
||
|
## Fast Math
|
||
|
|
||
|
**Fastmath.h** header file. Library to speed up program calculations. Use with [DECLARE LIBRARY](DECLARE-LIBRARY) "Fastmath"
|
||
|
|
||
|
```cpp
|
||
|
|
||
|
using namespace std;
|
||
|
unsigned long xrander=123456789, yrander=362436069, zrander = 521288629;
|
||
|
unsigned long trander;
|
||
|
double compbase = 2.3025850929940456840179914546844d;
|
||
|
double basecomp = 1.0d;
|
||
|
int shift1 = 1;
|
||
|
int shift2 = 5;
|
||
|
int shift3 = 16;
|
||
|
|
||
|
float Fast_Sqrt(float val) //log2(n) approximation //this is not used in the LN approximation
|
||
|
{
|
||
|
//float chk
|
||
|
union
|
||
|
{
|
||
|
int tmp;
|
||
|
float val;
|
||
|
} vals;
|
||
|
vals.val = val;
|
||
|
vals.tmp -= 1<<23; // Remove last bit so 1.0 gives 1.0
|
||
|
// tmp is now an approximation to logbase2(val)
|
||
|
vals.tmp >>= 1; // divide by 2
|
||
|
vals.tmp += 1<<29; // add 64 to exponent: (e+127)/2 =(e/2)+63;
|
||
|
//vals.tmp = (1<<29) + (vals.tmp >> 1) - (1<<22) + 0x4C000; // I am working on a way to improve this value
|
||
|
// that represents (e/2)-64 but want e/2
|
||
|
return vals.val;
|
||
|
}
|
||
|
|
||
|
float Fast_InvSqrt( float number ) // originally by Silicon Graphics slightly higher accuracy variant below.
|
||
|
{
|
||
|
long i;
|
||
|
float x2, y;
|
||
|
const float threehalfs = 1.5F;
|
||
|
x2 = number * 0.5F;
|
||
|
y = number;
|
||
|
i = * ( long * ) &y; // evil floating point bit level hacking [sic]
|
||
|
i = 0x5f375a86 - ( i >> 1 ); // what the fuck? [sic] original silicon graphics constant: 0x5f3759df
|
||
|
y = * ( float * ) &i;
|
||
|
y = y * ( threehalfs - ( x2 * y * y ) ); // 1st iteration
|
||
|
//y = y * ( threehalfs - ( x2 * y * y ) ); // 2nd iteration, this can be removed
|
||
|
return y;
|
||
|
}
|
||
|
|
||
|
double Fast_Pow(double a, double b) //fastpower originally developed by Martin Ankerl
|
||
|
{
|
||
|
int tmp = (*(1 + (int *)&a));
|
||
|
int tmp2 = (int)(b * (tmp - 1072632447) + 1072632447);
|
||
|
double p = 0.0;
|
||
|
*(1 + (int * )&p) = tmp2;
|
||
|
//p = p * a / 2.71828F ; failed attempt to auto correct the accuracy
|
||
|
return tmp;
|
||
|
}
|
||
|
|
||
|
double Fast_Exp(double y) //2.87921
|
||
|
{
|
||
|
double d;
|
||
|
//*((int*)(&d) + 0) = 0;
|
||
|
*((int*)(&d) + 1) = (int)(1512775 * y + 1072632447);
|
||
|
return d;
|
||
|
}
|
||
|
|
||
|
|
||
|
double Fast_XLnX(double x) //Borchardt's algorithm only accurate close to the origin
|
||
|
{
|
||
|
union
|
||
|
{
|
||
|
int tmp;
|
||
|
float val;
|
||
|
} vals;
|
||
|
vals.val = x;
|
||
|
vals.tmp -= 1<<23; // Remove last bit so 1.0 gives 1.0
|
||
|
vals.tmp >>= 1; // divide by 2
|
||
|
vals.tmp += 1<<29; // add 64 to exponent: (e+127)/2 =(e/2)+63;
|
||
|
return 6.0F * x * (x - 1.0F) / (x + 4.0F * vals.val ++);
|
||
|
}
|
||
|
|
||
|
double Fast_LnX(double x)
|
||
|
{
|
||
|
union
|
||
|
{
|
||
|
int tmp;
|
||
|
float val;
|
||
|
} vals;
|
||
|
vals.val = x;
|
||
|
vals.tmp -= 1<<23; // Remove last bit so 1.0 gives 1.0
|
||
|
vals.tmp >>= 1; // divide by 2
|
||
|
vals.tmp += 1<<29; // add 64 to exponent: (e+127)/2 =(e/2)+63;
|
||
|
return 6.0F * x * (x - 1.0F) / (x * (x + 4.0F * vals.val ++));
|
||
|
}
|
||
|
|
||
|
double Fast_Log10 (double x)
|
||
|
{
|
||
|
union
|
||
|
{
|
||
|
int tmp;
|
||
|
float val;
|
||
|
} vals;
|
||
|
vals.val = x;
|
||
|
vals.tmp -= 1<<23; // Remove last bit so 1.0 gives 1.0
|
||
|
vals.tmp >>= 1; // divide by 2
|
||
|
vals.tmp += 1<<29; // add 64 to exponent: (e+127)/2 =(e/2)+63;
|
||
|
return 6.0F * x * (x - 1.0F) / (x * 2.302585092994F * (x + 4.0F * vals.val ++));
|
||
|
}
|
||
|
|
||
|
|
||
|
double Fast_LogPi (double x)
|
||
|
{
|
||
|
union
|
||
|
{
|
||
|
int tmp;
|
||
|
float val;
|
||
|
} vals;
|
||
|
vals.val = x;
|
||
|
vals.tmp -= 1<<23; // Remove last bit so 1.0 gives 1.0
|
||
|
vals.tmp >>= 1; // divide by 2
|
||
|
vals.tmp += 1<<29; // add 64 to exponent: (e+127)/2 =(e/2)+63;
|
||
|
return 6.0F * x * (x - 1.0F) / (x * 1.1447298858494F * (x + 4.0F * vals.val ++));
|
||
|
}
|
||
|
|
||
|
double Fast_Sin(double x) //currently only supports range between -3 * pi and 3 * pi
|
||
|
{
|
||
|
double sine;
|
||
|
//always wrap input angle to -PI..PI I know it is abit of a compicated algorithm but hey
|
||
|
int Piintdiv = x * 0.15915494309189F;
|
||
|
double Pix = x - 6.2831853071796F * Piintdiv;
|
||
|
if (Pix < -3.14159265F)
|
||
|
Pix += 6.28318531F;
|
||
|
else
|
||
|
if (Pix > 3.14159265F)
|
||
|
Pix -= 6.28318531F;
|
||
|
if (Pix < 0)
|
||
|
sine = (0.405284735F * Pix + 1.27323954F) * Pix;
|
||
|
//sine = 1.27323954 * x + .405284735 * x * x; // I do so love horner form
|
||
|
else
|
||
|
sine = (1.27323954 - 0.405284735 * Pix) * Pix;
|
||
|
return sine;
|
||
|
}
|
||
|
|
||
|
|
||
|
double Fast_Cos(double x)
|
||
|
{
|
||
|
double Cosine;
|
||
|
int Piintdiv = x * 0.15915494309189F;
|
||
|
double Pix = x - 6.2831853071796F * Piintdiv;
|
||
|
Pix += 1.57079632;
|
||
|
if (Pix < -3.14159265F)
|
||
|
Pix += 6.28318531F;
|
||
|
else
|
||
|
if (Pix > 3.14159265F)
|
||
|
Pix -= 6.28318531F;
|
||
|
if (Pix < 0)
|
||
|
Cosine = (1.27323954 + 0.405284735 * Pix) * Pix;
|
||
|
else
|
||
|
Cosine = (1.27323954 - 0.405284735 * Pix) * Pix;
|
||
|
return Cosine;
|
||
|
}
|
||
|
|
||
|
double Fast_Tan(double x)
|
||
|
{
|
||
|
double Cosine;
|
||
|
double sine;
|
||
|
//always wrap input angle to -PI..PI
|
||
|
if (x < -3.14159265F)
|
||
|
x += 6.28318531F;
|
||
|
else
|
||
|
if (x > 3.14159265F)
|
||
|
x -= 6.28318531F;
|
||
|
//compute sine
|
||
|
if (x < 0)
|
||
|
sine = (0.405284735F * x + 1.27323954F) * x; // I do so love horner form the alernative on the nex line is slower
|
||
|
//sine = 1.27323954 * x + .405284735 * x * x;
|
||
|
else
|
||
|
sine = (1.27323954 - 0.405284735 * x) * x;
|
||
|
x += 1.57079632;
|
||
|
if (x < 0)
|
||
|
Cosine = (1.27323954 + 0.405284735 * x) * x;
|
||
|
else
|
||
|
Cosine = (1.27323954 - 0.405284735 * x) * x;
|
||
|
return sine / Cosine;
|
||
|
}
|
||
|
|
||
|
double Fast_CoTan(double x)
|
||
|
{
|
||
|
double Cosine;
|
||
|
double sine;
|
||
|
//always wrap input angle to -PI..PI
|
||
|
if (x < -3.14159265F)
|
||
|
x += 6.28318531F;
|
||
|
else
|
||
|
if (x > 3.14159265F)
|
||
|
x -= 6.28318531F;
|
||
|
//compute sine
|
||
|
if (x < 0)
|
||
|
sine = (0.405284735F * x + 1.27323954F) * x; // I do so love horner form the alernative on the nex line is slower
|
||
|
//sine = 1.27323954 * x + .405284735 * x * x;
|
||
|
else
|
||
|
sine = (1.27323954 - 0.405284735 * x) * x;
|
||
|
x += 1.57079632;
|
||
|
if (x < 0)
|
||
|
Cosine = (1.27323954 + 0.405284735 * x) * x;
|
||
|
else
|
||
|
Cosine = (1.27323954 - 0.405284735 * x) * x;
|
||
|
return Cosine / sine;
|
||
|
}
|
||
|
|
||
|
double Fast_Sec(double x)
|
||
|
{
|
||
|
double Cosine;
|
||
|
x += 1.57079632;
|
||
|
if (x > 3.14159265)
|
||
|
x -= 6.28318531;
|
||
|
if (x < 0)
|
||
|
Cosine = (1.27323954 + 0.405284735 * x) * x;
|
||
|
else
|
||
|
Cosine = (1.27323954 - 0.405284735 * x) * x;
|
||
|
return 1 / Cosine;
|
||
|
}
|
||
|
|
||
|
double Fast_Csc(double x)
|
||
|
{
|
||
|
double sine;
|
||
|
//always wrap input angle to -PI..PI
|
||
|
if (x < -3.14159265F)
|
||
|
x += 6.28318531F;
|
||
|
else
|
||
|
if (x > 3.14159265F)
|
||
|
x -= 6.28318531F;
|
||
|
//compute sine
|
||
|
if (x < 0)
|
||
|
sine = (0.405284735F * x + 1.27323954F) * x;
|
||
|
//sine = 1.27323954 * x + .405284735 * x * x; // I do so love horner form
|
||
|
else
|
||
|
sine = (1.27323954 - 0.405284735 * x) * x;
|
||
|
return 1 / sine;
|
||
|
}
|
||
|
|
||
|
|
||
|
// |error| < 0.005
|
||
|
float Fast_Atan2(float y, float x)
|
||
|
{
|
||
|
float PIBY2_FLOAT = 1.5707963F;
|
||
|
float PI_FLOAT = 3.14159265F;
|
||
|
if (x ## 0.0f)
|
||
|
{
|
||
|
if (y > 0.0f) return PIBY2_FLOAT;
|
||
|
if (y ## 0.0f) return 0.0f;
|
||
|
return 0 - PIBY2_FLOAT;
|
||
|
}
|
||
|
float atan;
|
||
|
float z = y/x;
|
||
|
if ( fabsf( z ) < 1.0f )
|
||
|
{
|
||
|
atan = z/(1.0f + 0.28f*z*z);
|
||
|
if ( x < 0.0f )
|
||
|
{
|
||
|
if ( y < 0.0f ) return atan - PI_FLOAT;
|
||
|
return atan + PI_FLOAT;
|
||
|
}
|
||
|
}
|
||
|
else
|
||
|
{
|
||
|
atan = PIBY2_FLOAT - z/(z*z + 0.28f);
|
||
|
if ( y < 0.0f )
|
||
|
return atan - PI_FLOAT;
|
||
|
}
|
||
|
return atan;
|
||
|
}
|
||
|
|
||
|
float Fast_Atan(float y)
|
||
|
{
|
||
|
float x = 1.0F;
|
||
|
float PIBY2_FLOAT = 1.5707963F;
|
||
|
float PI_FLOAT = 3.14159265F;
|
||
|
float atan;
|
||
|
float z = y/x;
|
||
|
if (fabsf( z ) < 1.0f )
|
||
|
{
|
||
|
atan = z/(1.0f + 0.28f*z*z);
|
||
|
if ( x < 0.0f )
|
||
|
{
|
||
|
if ( y < 0.0f ) return atan - PI_FLOAT;
|
||
|
return atan + PI_FLOAT;
|
||
|
}
|
||
|
}
|
||
|
else
|
||
|
{
|
||
|
atan = PIBY2_FLOAT - z/(z*z + 0.28f);
|
||
|
if ( y < 0.0f )
|
||
|
return atan - PI_FLOAT;
|
||
|
}
|
||
|
return atan;
|
||
|
}
|
||
|
|
||
|
double Fast_ACos(double x)
|
||
|
{
|
||
|
float retval = (-0.69813170079773212F * x * x - 0.87266462599716477F) * x + 1.5707963267948966F;
|
||
|
return retval;
|
||
|
}
|
||
|
|
||
|
double Fast_ASin(double x)
|
||
|
{
|
||
|
float retval = (0.69813170079773212F * x * x + 0.87266462599716477F) * x;
|
||
|
return retval;
|
||
|
}
|
||
|
|
||
|
double Fast_SinH(double x)
|
||
|
{
|
||
|
double y = -x;
|
||
|
double k;
|
||
|
double d;
|
||
|
//*((int*)(&d) + 0) = 0;
|
||
|
*((int*)(&d) + 1) = (int)(1512775 * y + 1072632447);
|
||
|
//*((int*)(&k) + 0) = 0;
|
||
|
*((int*)(&k) + 1) = (int)(1512775 * x + 1072632447);
|
||
|
return (k - d) / 2;
|
||
|
}
|
||
|
|
||
|
double Fast_CosH(double x)
|
||
|
{
|
||
|
double y = -x;
|
||
|
double k;
|
||
|
double d;
|
||
|
//*((int*)(&d) + 0) = 0;
|
||
|
*((int*)(&d) + 1) = (int)(1512775 * y + 1072632447);
|
||
|
//*((int*)(&k) + 0) = 0;
|
||
|
*((int*)(&k) + 1) = (int)(1512775 * x + 1072632447);
|
||
|
return (k + d) / 2;
|
||
|
}
|
||
|
|
||
|
double Fast_TanH(double x) // from http://www.musicdsp.org/showone.php?id=238
|
||
|
{
|
||
|
//double xt;
|
||
|
if (x < -3)
|
||
|
{
|
||
|
return -1;
|
||
|
}
|
||
|
else if (x > 3)
|
||
|
{
|
||
|
return 1;
|
||
|
}
|
||
|
else
|
||
|
{
|
||
|
double xz = x * x;
|
||
|
return x * (27 + xz) / (27 + 9 * xz);
|
||
|
}
|
||
|
//return xt
|
||
|
}
|
||
|
|
||
|
double Fast_ATanH(double x)
|
||
|
{
|
||
|
double mlnx = 1-x;
|
||
|
double plnx = 1+x;
|
||
|
union
|
||
|
{
|
||
|
int tmp;
|
||
|
float val;
|
||
|
} vals;
|
||
|
union
|
||
|
{
|
||
|
int tmp;
|
||
|
float val;
|
||
|
} vals2;
|
||
|
vals.val = plnx;
|
||
|
vals.tmp -= 1<<23; // Remove last bit so 1.0 gives 1.0
|
||
|
vals.tmp >>= 1; // divide by 2
|
||
|
vals.tmp += 1<<29; // add 64 to exponent: (e+127)/2 =(e/2)+63;
|
||
|
vals2.val = mlnx;
|
||
|
vals2.tmp -= 1<<23; // Remove last bit so 1.0 gives 1.0
|
||
|
vals2.tmp >>= 1; // divide by 2
|
||
|
vals2.tmp += 1<<29; // add 64 to exponent: (e+127)/2 =(e/2)+63;
|
||
|
return (6.0F * mlnx * (mlnx - 1.0F) / (mlnx * (mlnx + 4.0F * vals.val ++)) - 6.0F * plnx * (plnx - 1.0F) / (plnx * (plnx + 4.0F * vals2.val ++))) * 0.5f;
|
||
|
}
|
||
|
|
||
|
double Fast_ACosH(double val)
|
||
|
{
|
||
|
double zsqrz;
|
||
|
//float chk
|
||
|
union
|
||
|
{
|
||
|
int tmp;
|
||
|
float val;
|
||
|
} vals;
|
||
|
union
|
||
|
{
|
||
|
int tmp;
|
||
|
float val;
|
||
|
} val2s;
|
||
|
vals.val = val + 1;
|
||
|
vals.tmp -= 1<<23; // Remove last bit so 1.0 gives 1.0
|
||
|
// tmp is now an approximation to logbase2(val)
|
||
|
vals.tmp >>= 1; // divide by 2
|
||
|
vals.tmp += 1<<29; // add 64 to exponent: (e+127)/2 =(e/2)+63;
|
||
|
val2s.val = val - 1;
|
||
|
val2s.tmp -= 1<<23; // Remove last bit so 1.0 gives 1.0
|
||
|
// tmp is now an approximation to logbase2(val)
|
||
|
val2s.tmp >>= 1; // divide by 2
|
||
|
val2s.tmp += 1<<29; // add 64 to exponent: (e+127)/2 =(e/2)+63;
|
||
|
zsqrz = val + vals.val * val2s.val;
|
||
|
return log(zsqrz);
|
||
|
}
|
||
|
|
||
|
double Fast_ASinH(double x)
|
||
|
{
|
||
|
union
|
||
|
{
|
||
|
int tmp;
|
||
|
float val;
|
||
|
} vals;
|
||
|
vals.val = x * x + 1;
|
||
|
vals.tmp -= 1<<23; // Remove last bit so 1.0 gives 1.0
|
||
|
// tmp is now an approximation to logbase2(val)
|
||
|
vals.tmp >>= 1; // divide by 2
|
||
|
vals.tmp += 1<<29; // add 64 to exponent: (e+127)/2 =(e/2)+63;
|
||
|
return log(x + vals.val);
|
||
|
}
|
||
|
|
||
|
double Fast_ASecH(double x)
|
||
|
{
|
||
|
double zsqrz;
|
||
|
//float chk
|
||
|
union
|
||
|
{
|
||
|
int tmp;
|
||
|
float val;
|
||
|
} vals;
|
||
|
union
|
||
|
{
|
||
|
int tmp;
|
||
|
float val;
|
||
|
} val2s;
|
||
|
vals.val = 1 / x + 1;
|
||
|
vals.tmp -= 1<<23; // Remove last bit so 1.0 gives 1.0
|
||
|
// tmp is now an approximation to logbase2(val)
|
||
|
vals.tmp >>= 1; // divide by 2
|
||
|
vals.tmp += 1<<29; // add 64 to exponent: (e+127)/2 =(e/2)+63;
|
||
|
val2s.val = 1 / x - 1;
|
||
|
val2s.tmp -= 1<<23; // Remove last bit so 1.0 gives 1.0
|
||
|
// tmp is now an approximation to logbase2(val)
|
||
|
val2s.tmp >>= 1; // divide by 2
|
||
|
val2s.tmp += 1<<29; // add 64 to exponent: (e+127)/2 =(e/2)+63;
|
||
|
zsqrz = 1 / x + vals.val * val2s.val;
|
||
|
return log(zsqrz);
|
||
|
}
|
||
|
|
||
|
double Fast_ACscH(double val)
|
||
|
{
|
||
|
double x = 1 / val;
|
||
|
union
|
||
|
{
|
||
|
int tmp;
|
||
|
float val;
|
||
|
} vals;
|
||
|
vals.val = x * x + 1;
|
||
|
vals.tmp -= 1<<23; // Remove last bit so 1.0 gives 1.0
|
||
|
// tmp is now an approximation to logbase2(val)
|
||
|
vals.tmp >>= 1; // divide by 2
|
||
|
vals.tmp += 1<<29; // add 64 to exponent: (e+127)/2 =(e/2)+63;
|
||
|
return log(x + vals.val);
|
||
|
}
|
||
|
|
||
|
double Fast_ACotH(double valts)
|
||
|
{
|
||
|
double x = 1 / valts;
|
||
|
double mlnx = 1-x;
|
||
|
double plnx = 1+x;
|
||
|
union
|
||
|
{
|
||
|
int tmp;
|
||
|
float val;
|
||
|
} vals;
|
||
|
union
|
||
|
{
|
||
|
int tmp;
|
||
|
float val;
|
||
|
} vals2;
|
||
|
vals.val = plnx;
|
||
|
vals.tmp -= 1<<23; // Remove last bit so 1.0 gives 1.0
|
||
|
vals.tmp >>= 1; // divide by 2
|
||
|
vals.tmp += 1<<29; // add 64 to exponent: (e+127)/2 =(e/2)+63;
|
||
|
vals2.val = mlnx;
|
||
|
vals2.tmp -= 1<<23; // Remove last bit so 1.0 gives 1.0
|
||
|
vals2.tmp >>= 1; // divide by 2
|
||
|
vals2.tmp += 1<<29; // add 64 to exponent: (e+127)/2 =(e/2)+63;
|
||
|
return (6.0F * mlnx * (mlnx - 1.0F) / (mlnx * (mlnx + 4.0F * vals.val ++)) - 6.0F * plnx * (plnx - 1.0F) / (plnx * (plnx + 4.0F * vals2.val ++))) / 2;
|
||
|
}
|
||
|
|
||
|
void Fast_RandInit(float x,int i,int j, int k)
|
||
|
{
|
||
|
float ftl;
|
||
|
float ftl2;
|
||
|
float ftl3;
|
||
|
ftl = 10000 * sin(x*7);
|
||
|
xrander = * ( long * ) &ftl;
|
||
|
ftl2 = 10000 * cos(x*3);
|
||
|
yrander = * ( long * ) &ftl2;
|
||
|
ftl3 = 10000 * sin(x*2);
|
||
|
zrander = * ( long * ) &ftl3;
|
||
|
if ((i ## 0)||(j ## 0)||( k ## 0 ))
|
||
|
{
|
||
|
shift1 = 1;
|
||
|
shift2 = 5;
|
||
|
shift3 = 16;
|
||
|
}
|
||
|
else
|
||
|
{
|
||
|
shift1 = i; //8 * sin(x) + 9;
|
||
|
shift3 = k; //8 * cos(x) + 9;
|
||
|
shift2 = j; //(shift1 + shift3) / 2;
|
||
|
}
|
||
|
}
|
||
|
|
||
|
inline float Fast_Rand(void) //based off of George Marsaglia's XORSHIFT algorithms
|
||
|
{ //
|
||
|
xrander ^= xrander << shift1;
|
||
|
xrander ^= xrander >> shift2;
|
||
|
xrander ^= xrander << shift3;
|
||
|
trander = xrander;
|
||
|
xrander = yrander;
|
||
|
yrander = zrander;
|
||
|
zrander = trander ^ xrander ^ yrander;
|
||
|
return float (zrander)/4294967295;
|
||
|
}
|
||
|
|
||
|
inline float Fast_Sign(float f)
|
||
|
{
|
||
|
float r = 1.0f;
|
||
|
(int&)r |= ((int&)f & 0x80000000);
|
||
|
return r;
|
||
|
}
|
||
|
|
||
|
double High_ATanH(double x)
|
||
|
{
|
||
|
return (log(1 + x) - log(1 - x)) / 2;
|
||
|
}
|
||
|
|
||
|
double High_ACosH(double x)
|
||
|
{
|
||
|
return log(x + sqrt(x + 1) * sqrt(x - 1));
|
||
|
}
|
||
|
|
||
|
double High_ASinH(double x)
|
||
|
{
|
||
|
return log(x + sqrt(x * x + 1));
|
||
|
}
|
||
|
|
||
|
double High_ASecH(double x)
|
||
|
{
|
||
|
double xz = 1 / x;
|
||
|
return log(sqrt(xz - 1) * sqrt(xz + 1) + xz);
|
||
|
}
|
||
|
|
||
|
double High_ACscH(double x)
|
||
|
{
|
||
|
return log(sqrt(1 + 1 / (x * x)) + 1 / x);
|
||
|
}
|
||
|
|
||
|
double High_ACotH(double x)
|
||
|
{
|
||
|
double xz = 1 / x;
|
||
|
return (log(1 + xz) - log(1 - xz)) / 2;
|
||
|
}
|
||
|
|
||
|
inline float High_Sign(float f)
|
||
|
{
|
||
|
//this may or may not be used in QB64 in the future
|
||
|
if (((int&)f & 0x7FFFFFFF)## 0)
|
||
|
{
|
||
|
return 0.0f;
|
||
|
}
|
||
|
else
|
||
|
{
|
||
|
float r = 1.0f;
|
||
|
(int&)r |= ((int&)f & 0x80000000);
|
||
|
return r;
|
||
|
}
|
||
|
}
|
||
|
int Misc_TrailZCount(unsigned int v)
|
||
|
{
|
||
|
//unsigned int v; // 32-bit word input to count zero bits on right
|
||
|
unsigned int c; // c will be the number of zero bits on the right,
|
||
|
// so if v is 1101000 (base 2), then c will be 3
|
||
|
// NOTE: if 0 ## v, then c = 31.
|
||
|
if (v & 0x1)
|
||
|
{
|
||
|
// special case for odd v (assumed to happen half of the time)
|
||
|
c = 0;
|
||
|
}
|
||
|
else
|
||
|
{
|
||
|
c = 1;
|
||
|
if ((v & 0xffff) ## 0)
|
||
|
{
|
||
|
v >>= 16;
|
||
|
c += 16;
|
||
|
}
|
||
|
if ((v & 0xff) ## 0)
|
||
|
{
|
||
|
v >>= 8;
|
||
|
c += 8;
|
||
|
}
|
||
|
if ((v & 0xf) ## 0)
|
||
|
{
|
||
|
v >>= 4;
|
||
|
c += 4;
|
||
|
}
|
||
|
if ((v & 0x3) ## 0)
|
||
|
{
|
||
|
v >>= 2;
|
||
|
c += 2;
|
||
|
}
|
||
|
c -= v & 0x1;
|
||
|
}
|
||
|
return c;
|
||
|
}
|
||
|
|
||
|
int Misc_Parity(int v)// is the number of 1 bits odd? true if so
|
||
|
{
|
||
|
//unsigned int v; // 32-bit word
|
||
|
v ^= v >> 1;
|
||
|
v ^= v >> 2;
|
||
|
v = (v & 0x11111111U) * 0x11111111U;
|
||
|
return (v >> 28) & 1;
|
||
|
}
|
||
|
|
||
|
int Misc_BitSet(int v)
|
||
|
//u//nsigned int v; // count the number of bits set in v
|
||
|
{
|
||
|
unsigned int c; // c accumulates the total bits set in v
|
||
|
for (c = 0; v; c++)
|
||
|
{
|
||
|
v &= v - 1; // clear the least significant bit set
|
||
|
}
|
||
|
return c;
|
||
|
}
|
||
|
|
||
|
long Misc_FloatToLong(float y) //1132462080
|
||
|
{
|
||
|
return * ( long * ) &y;
|
||
|
}
|
||
|
|
||
|
float Misc_LongToFloat(long y)
|
||
|
{
|
||
|
return * ( float * ) &y;
|
||
|
}
|
||
|
|
||
|
long Misc_UnSetRMBit(long x)
|
||
|
{
|
||
|
return x & (x - 1);
|
||
|
}
|
||
|
|
||
|
```
|
||
|
|
||
|
<sub>Library created by DarthWho</sub>
|
||
|
|
||
|
**[Latest Fast Math Update](http://dl.dropbox.com/u/12359848/fastmath.h)**
|
||
|
|
||
|
Using **trailzcount** to speed up the process of finding the prime factors of a number.
|
||
|
|
||
|
```vb
|
||
|
|
||
|
DECLARE LIBRARY "fastmath"
|
||
|
FUNCTION TrailZCount% ALIAS Misc_TrailZCount (BYVAL vals AS LONG)
|
||
|
END DECLARE
|
||
|
INPUT "enter a number: ", xt&
|
||
|
tabs = TrailZCount%(xt&)
|
||
|
PRINT "the prime factors of"; xt&; " are:";
|
||
|
IF tabs THEN
|
||
|
FOR i = 1 TO tabs
|
||
|
PRINT 2;
|
||
|
NEXT i
|
||
|
END IF
|
||
|
x& = xt& / (2 ^ tabs)
|
||
|
c& = 3
|
||
|
WHILE c& <= x&
|
||
|
WHILE x& MOD c& = 0
|
||
|
x& = x& / c&
|
||
|
PRINT c&;
|
||
|
WEND
|
||
|
c& = c& + 1
|
||
|
WEND
|
||
|
PRINT
|
||
|
|
||
|
```
|
||
|
|
||
|
Using **parity** to detect single bit errors simulated using **unsetrmbit**.
|
||
|
|
||
|
```vb
|
||
|
|
||
|
'parity checking code which may be used in order to detect a download error:
|
||
|
'will also use the bitset function to simulate a download error .bitset unsets the lowest set bit
|
||
|
'parity checking code which may be used in order to detect a download error:
|
||
|
'will also use the bitset function to simulate a download error .bitset unsets the lowest set bit
|
||
|
DECLARE LIBRARY "fastmath"
|
||
|
FUNCTION parity& ALIAS Misc_Parity (BYVAL val AS LONG)
|
||
|
FUNCTION unsetrmbit& ALIAS Misc_UnSetRMBit (BYVAL vals AS LONG)
|
||
|
END DECLARE
|
||
|
TYPE DLDATA
|
||
|
par AS _BIT
|
||
|
dataa AS LONG
|
||
|
END TYPE
|
||
|
DIM values(1) 'here is where the data is created
|
||
|
g& = 8
|
||
|
h& = 9
|
||
|
values.par(0) = parity&(g&)
|
||
|
values.dataa(0) = g&
|
||
|
values.par(1) = parity&(h&)
|
||
|
values.dataa(1) = h& 'transfer and error in values.dataa(1)
|
||
|
values.dataa(1) = unsetrmbit&(values.dataa(1)) 'checking for single bit errors in download
|
||
|
FOR i = 0 TO 1
|
||
|
IF values.par(i) = parity&(values.dataa(i)) THEN
|
||
|
PRINT "no single bit errors detected in in data packet"; i
|
||
|
ELSE
|
||
|
PRINT "error found in data packet"; i
|
||
|
END IF
|
||
|
NEXT
|
||
|
|
||
|
```
|
||
|
|
||
|
## See Also
|
||
|
|
||
|
* [DECLARE LIBRARY](DECLARE-LIBRARY), [BYVAL](BYVAL)
|
||
|
* [_OFFSET](_OFFSET), [_OFFSET (function)](_OFFSET-(function)) (lp, ptr and p names)
|
||
|
* [DLL Libraries](DLL-Libraries), [Windows Libraries](Windows-Libraries)
|
||
|
* [Libraries](Libraries)
|