Compare commits
2 Commits
35f051b5da
...
c3629d7150
| Author | SHA1 | Date | |
|---|---|---|---|
| c3629d7150 | |||
| 3c3d2716a5 |
@ -5,7 +5,8 @@ set(CMAKE_C_STANDARD 99)
|
|||||||
|
|
||||||
add_compile_options(-g -Wall -pedantic)
|
add_compile_options(-g -Wall -pedantic)
|
||||||
|
|
||||||
add_executable(njvm njvm.c)
|
add_executable(njvm njvm.c
|
||||||
|
njvm.h)
|
||||||
|
|
||||||
# Include directories for njvm executable
|
# Include directories for njvm executable
|
||||||
target_include_directories(njvm PRIVATE ${CMAKE_CURRENT_SOURCE_DIR}/bigint/build/include)
|
target_include_directories(njvm PRIVATE ${CMAKE_CURRENT_SOURCE_DIR}/bigint/build/include)
|
||||||
|
|||||||
97
GC.c
97
GC.c
@ -1,105 +1,14 @@
|
|||||||
#ifndef GC
|
#ifndef GC
|
||||||
#define GC
|
#define GC
|
||||||
|
|
||||||
#include <stdlib.h>
|
#include <stdlib.h>
|
||||||
#include <stdio.h>
|
#include <stdio.h>
|
||||||
#include <string.h>
|
#include <string.h>
|
||||||
#include <stdbool.h>
|
#include <stdbool.h>
|
||||||
#include "objref.c"
|
#include "objref.c"
|
||||||
#include "instruktion.c"
|
#include "instruktion.c"
|
||||||
|
#include "support.c"
|
||||||
|
#include "bigint.h"
|
||||||
|
|
||||||
bool halbspeicher1IsActiv = true;
|
|
||||||
unsigned char* halbspeicher1;
|
|
||||||
unsigned char* halbspeicher2;
|
|
||||||
unsigned char* freizeiger;
|
|
||||||
unsigned int groesse;
|
|
||||||
|
|
||||||
void initHeap(unsigned int size){
|
|
||||||
groesse = size*1024/2;
|
|
||||||
if((halbspeicher1 = malloc(groesse)) == NULL) perror("malloc");
|
|
||||||
if((halbspeicher2 = malloc(groesse)) == NULL) perror("malloc");
|
|
||||||
freizeiger = halbspeicher1;
|
|
||||||
}
|
|
||||||
ObjRef copyObjToFreeMem(ObjRef obj){
|
|
||||||
ObjRef copy = (ObjRef) halbspeicher2[freizeiger];
|
|
||||||
unsigned int size = (IS_PRIMITIVE(obj))? obj->size : GET_ELEMENT_COUNT(obj)* sizeof(void *);
|
|
||||||
memcpy((void *) halbspeicher2[freizeiger], obj, size);
|
|
||||||
freizeiger += size;
|
|
||||||
return copy;
|
|
||||||
}
|
|
||||||
void setBrokenHeart(ObjRef obj, bool val){
|
|
||||||
//TODO
|
|
||||||
}
|
|
||||||
bool getBrokenHeart(ObjRef obj){
|
|
||||||
//TODO
|
|
||||||
return true;
|
|
||||||
}
|
|
||||||
void setForwardPointer(ObjRef obj, void* ref){
|
|
||||||
//TODO
|
|
||||||
}
|
|
||||||
void *getForwardPointer(ObjRef obj){
|
|
||||||
//TODO
|
|
||||||
return malloc(0);
|
|
||||||
}
|
|
||||||
|
|
||||||
ObjRef reallocate(ObjRef orig){
|
|
||||||
ObjRef copy;
|
|
||||||
if(orig == NULL){
|
|
||||||
copy == NULL;
|
|
||||||
}else{
|
|
||||||
if(getBrokenHeart(orig)){
|
|
||||||
setForwardPointer(copy,orig);
|
|
||||||
}else{
|
|
||||||
copy = copyObjToFreeMem(orig);
|
|
||||||
setBrokenHeart(orig,true);
|
|
||||||
setForwardPointer(orig,copy);
|
|
||||||
}
|
|
||||||
}
|
|
||||||
return copy;
|
|
||||||
}
|
|
||||||
unsigned int getSize(ObjRef obj){
|
|
||||||
return (IS_PRIMITIVE(obj))? obj->size : GET_ELEMENT_COUNT(obj)* sizeof(void *);
|
|
||||||
}
|
|
||||||
void scan(void){
|
|
||||||
void *scan = halbspeicher1;
|
|
||||||
while (scan < freizeiger){
|
|
||||||
if(!IS_PRIMITIVE((ObjRef)scan)){
|
|
||||||
for (unsigned int i = 0; i < GET_ELEMENT_COUNT((ObjRef)scan); i += getSize(GET_REFS_PTR((ObjRef)scan)[i])) {
|
|
||||||
GET_REFS_PTR((ObjRef)scan)[i] = reallocate(GET_REFS_PTR((ObjRef)scan)[i]);
|
|
||||||
}
|
|
||||||
}
|
|
||||||
scan += getSize((ObjRef)scan);
|
|
||||||
}
|
|
||||||
}
|
|
||||||
void *alloc(unsigned int size){
|
|
||||||
if((&freizeiger - &halbspeicher1)/8 + size > groesse) {
|
|
||||||
printf("Heap Overflow\n");
|
|
||||||
exit(EXIT_FAILURE);
|
|
||||||
}
|
|
||||||
// Set broken heart to false
|
|
||||||
setBrokenHeart((ObjRef)freizeiger,false);
|
|
||||||
freizeiger += sizeof(bool) + sizeof(void *);
|
|
||||||
void *ptr = &freizeiger;
|
|
||||||
freizeiger += size;
|
|
||||||
return ptr;
|
|
||||||
}
|
|
||||||
void swap(void){
|
|
||||||
unsigned char *tmp = halbspeicher1;
|
|
||||||
halbspeicher1 = halbspeicher2;
|
|
||||||
halbspeicher2 = tmp;
|
|
||||||
freizeiger = halbspeicher1;
|
|
||||||
halbspeicher1IsActiv = !halbspeicher1IsActiv;
|
|
||||||
}
|
|
||||||
void triggerGC(void){
|
|
||||||
printf("GC triggered\n");
|
|
||||||
//TODO
|
|
||||||
}
|
|
||||||
void copy(){
|
|
||||||
ObjRef *test;
|
|
||||||
unsigned int objResSize = 2;
|
|
||||||
if((test = malloc(2*sizeof(test)))==NULL) perror("malloc");
|
|
||||||
for (int i = 0; i < objResSize; ++i) {
|
|
||||||
|
|
||||||
}
|
|
||||||
}
|
|
||||||
|
|
||||||
#endif
|
#endif
|
||||||
|
|||||||
22
bigint/Makefile
Normal file
22
bigint/Makefile
Normal file
@ -0,0 +1,22 @@
|
|||||||
|
#
|
||||||
|
# Makefile for big integer library and test
|
||||||
|
#
|
||||||
|
|
||||||
|
DIRS = src tst
|
||||||
|
|
||||||
|
all: clean
|
||||||
|
for i in $(DIRS) ; do \
|
||||||
|
$(MAKE) -C $$i install ; \
|
||||||
|
done
|
||||||
|
|
||||||
|
clean:
|
||||||
|
for i in $(DIRS) ; do \
|
||||||
|
$(MAKE) -C $$i clean ; \
|
||||||
|
done
|
||||||
|
rm -rf ./build
|
||||||
|
rm -f *~
|
||||||
|
|
||||||
|
dist: clean
|
||||||
|
(cd .. ; \
|
||||||
|
tar -cvf bigint.tar bigint ; \
|
||||||
|
gzip -f bigint.tar)
|
||||||
92
bigint/README
Normal file
92
bigint/README
Normal file
@ -0,0 +1,92 @@
|
|||||||
|
The "bigint" Package
|
||||||
|
====================
|
||||||
|
|
||||||
|
1. What is it?
|
||||||
|
--------------
|
||||||
|
|
||||||
|
This package implements a multiple-precision integer arithmetic package,
|
||||||
|
i.e., a collection of functions which can calculate with integers having
|
||||||
|
arbitrarily many digits. The algorithms are taken from [D. Knuth: The
|
||||||
|
Art of Computer Programming, Vol. 2, Seminumerical Algorithms], the
|
||||||
|
implementation language is C.
|
||||||
|
|
||||||
|
|
||||||
|
2. "Multiple Precision" - how does it work?
|
||||||
|
-------------------------------------------
|
||||||
|
|
||||||
|
Each integer number is represented as an array of digits. The array
|
||||||
|
is large enough to hold the number of digits necessary to represent
|
||||||
|
the number. Each digit occupies a single byte, so the number base of
|
||||||
|
this representation is 256. Addition, subtraction, and multiplication
|
||||||
|
work as we all have learned it: perform the desired operation digit
|
||||||
|
by digit, starting from the least significant digit, and observing
|
||||||
|
any "carries" from one place to the next higher one. Division is a
|
||||||
|
little bit more complicated because there is a certain amount of
|
||||||
|
guesswork involved. Knuth gives a formal treatment of this guesswork.
|
||||||
|
|
||||||
|
|
||||||
|
3. How do I use it?
|
||||||
|
-------------------
|
||||||
|
|
||||||
|
Because every big integer may have a differently sized array to hold
|
||||||
|
its digits, these structures are dynamically allocated on the heap of
|
||||||
|
the C runtime system, and accessed by pointers. If you want to perform
|
||||||
|
an arithmetic operation on one or two big integers, you have to load
|
||||||
|
the corresponding pointers into a structure called BIP ("Big Integer
|
||||||
|
Processor"), and call the arithmetic function. When the function has
|
||||||
|
returned, the pointer to the result of the operation can be found in
|
||||||
|
another component of the BIP. The following functions are available:
|
||||||
|
|
||||||
|
int bigSgn(void); /* sign */
|
||||||
|
int bigCmp(void); /* comparison */
|
||||||
|
void bigNeg(void); /* negation */
|
||||||
|
void bigAdd(void); /* addition */
|
||||||
|
void bigSub(void); /* subtraction */
|
||||||
|
void bigMul(void); /* multiplication */
|
||||||
|
void bigDiv(void); /* division */
|
||||||
|
void bigFromInt(int n); /* conversion int --> big */
|
||||||
|
int bigToInt(void); /* conversion big --> int */
|
||||||
|
void bigRead(FILE *in); /* read a big integer */
|
||||||
|
void bigPrint(FILE *out); /* print a big integer */
|
||||||
|
void bigDump(FILE *out, BigObjRef bigObjRef); /* dump a big integer */
|
||||||
|
|
||||||
|
Some of these functions accept or return ordinary integers. For the
|
||||||
|
exact definition of each function's interface, please see the comments
|
||||||
|
in the function's source.
|
||||||
|
|
||||||
|
|
||||||
|
4. What else is needed?
|
||||||
|
-----------------------
|
||||||
|
|
||||||
|
The library tries to detect fatal errors in using its functions (e.g.,
|
||||||
|
null pointers to operands) as well as internal errors (which "cannot
|
||||||
|
happen"). In either case a user-supplied error routine is called, which
|
||||||
|
is supposed to print an error message and then to terminate the program.
|
||||||
|
|
||||||
|
The library does not attempt to manage memory. For this purpose, it
|
||||||
|
relies on a user-supplied function "void* newPrimObject(int dataSize)",
|
||||||
|
which should allocate sufficiently many bytes and return a pointer to
|
||||||
|
a created object. In addition, the library needs to access a data field
|
||||||
|
- that implements the size-as-needed cocept - within the created object,
|
||||||
|
to actualy manage the data structures that represent the values of the
|
||||||
|
integer number. For this, the following function can be used:
|
||||||
|
|
||||||
|
void * getPrimObjectDataPointer(void * obj){
|
||||||
|
ObjRef oo = ((ObjRef) (obj));
|
||||||
|
return oo->data;
|
||||||
|
}
|
||||||
|
|
||||||
|
This function must also be a user-supplied function that must be
|
||||||
|
implemented by the program that wants to use this library.
|
||||||
|
|
||||||
|
For details about these three user-supplied functions take a look in
|
||||||
|
the file "support.c" in the directory "tst".
|
||||||
|
|
||||||
|
|
||||||
|
5. What is in the directory "tst"?
|
||||||
|
----------------------------------
|
||||||
|
|
||||||
|
Well, you may have guessed it already: these are test cases for the
|
||||||
|
library. You can learn how to link against the library by inspecting
|
||||||
|
the "Makefile" for the tests, and you can find a simple implementation
|
||||||
|
of the support library.
|
||||||
BIN
bigint/build/bin/testbip
Executable file
BIN
bigint/build/bin/testbip
Executable file
Binary file not shown.
60
bigint/build/include/bigint.h
Normal file
60
bigint/build/include/bigint.h
Normal file
@ -0,0 +1,60 @@
|
|||||||
|
/*
|
||||||
|
* bigint.h -- big integer library
|
||||||
|
*/
|
||||||
|
|
||||||
|
|
||||||
|
#ifndef _BIGINT_H_
|
||||||
|
#define _BIGINT_H_
|
||||||
|
|
||||||
|
|
||||||
|
/* object representation */
|
||||||
|
typedef void* BigObjRef;
|
||||||
|
|
||||||
|
|
||||||
|
#include <stdio.h>
|
||||||
|
|
||||||
|
|
||||||
|
typedef struct {
|
||||||
|
int nd; /* number of digits; array may be bigger */
|
||||||
|
/* nd = 0 exactly when number = 0 */
|
||||||
|
unsigned char sign; /* one of BIG_NEGATIVE or BIG_POSITIVE */
|
||||||
|
/* zero always has BIG_POSITIVE here */
|
||||||
|
unsigned char digits[1]; /* the digits proper; number base is 256 */
|
||||||
|
/* LS digit first; MS digit is not zero */
|
||||||
|
} Big;
|
||||||
|
|
||||||
|
#include "support.h"
|
||||||
|
|
||||||
|
|
||||||
|
/* big integer processor registers */
|
||||||
|
|
||||||
|
typedef struct {
|
||||||
|
BigObjRef op1; /* first (or single) operand */
|
||||||
|
BigObjRef op2; /* second operand (if present) */
|
||||||
|
BigObjRef res; /* result of operation */
|
||||||
|
BigObjRef rem; /* remainder in case of division */
|
||||||
|
} BIP;
|
||||||
|
|
||||||
|
extern BIP bip; /* registers of the processor */
|
||||||
|
|
||||||
|
|
||||||
|
/* big integer processor functions */
|
||||||
|
|
||||||
|
int bigSgn(void); /* sign */
|
||||||
|
int bigCmp(void); /* comparison */
|
||||||
|
void bigNeg(void); /* negation */
|
||||||
|
void bigAdd(void); /* addition */
|
||||||
|
void bigSub(void); /* subtraction */
|
||||||
|
void bigMul(void); /* multiplication */
|
||||||
|
void bigDiv(void); /* division */
|
||||||
|
|
||||||
|
void bigFromInt(int n); /* conversion int --> big */
|
||||||
|
int bigToInt(void); /* conversion big --> int */
|
||||||
|
|
||||||
|
void bigRead(FILE *in); /* read a big integer */
|
||||||
|
void bigPrint(FILE *out); /* print a big integer */
|
||||||
|
|
||||||
|
void bigDump(FILE *out, BigObjRef bigObjRef); /* dump a big integer object */
|
||||||
|
|
||||||
|
|
||||||
|
#endif /* _BIGINT_H_ */
|
||||||
15
bigint/build/include/support.h
Normal file
15
bigint/build/include/support.h
Normal file
@ -0,0 +1,15 @@
|
|||||||
|
/*
|
||||||
|
* support.h -- object representation and support functions
|
||||||
|
*/
|
||||||
|
|
||||||
|
|
||||||
|
#ifndef _SUPPORT_H_
|
||||||
|
#define _SUPPORT_H_
|
||||||
|
|
||||||
|
/* support functions */
|
||||||
|
|
||||||
|
void fatalError(char *msg); /* print a message and exit */
|
||||||
|
void * newPrimObject(int dataSize); /* create a new primitive object */
|
||||||
|
void * getPrimObjectDataPointer(void * primObject);
|
||||||
|
|
||||||
|
#endif /* _SUPPORT_H_ */
|
||||||
BIN
bigint/build/lib/libbigint.a
Normal file
BIN
bigint/build/lib/libbigint.a
Normal file
Binary file not shown.
26
bigint/src/Makefile
Normal file
26
bigint/src/Makefile
Normal file
@ -0,0 +1,26 @@
|
|||||||
|
#
|
||||||
|
# Makefile for big integer library
|
||||||
|
#
|
||||||
|
|
||||||
|
BUILD = ../build
|
||||||
|
|
||||||
|
CC = gcc
|
||||||
|
CFLAGS = -g -Wall
|
||||||
|
|
||||||
|
all: support.h bigint.h libbigint.a
|
||||||
|
|
||||||
|
install: all
|
||||||
|
mkdir -p $(BUILD)/include
|
||||||
|
cp support.h $(BUILD)/include
|
||||||
|
cp bigint.h $(BUILD)/include
|
||||||
|
mkdir -p $(BUILD)/lib
|
||||||
|
cp libbigint.a $(BUILD)/lib
|
||||||
|
|
||||||
|
libbigint.a: bigint.o
|
||||||
|
ar -crs libbigint.a bigint.o
|
||||||
|
|
||||||
|
bigint.o: bigint.c bigint.h support.h
|
||||||
|
$(CC) $(CFLAGS) -o bigint.o -c bigint.c
|
||||||
|
|
||||||
|
clean:
|
||||||
|
rm -f *~ bigint.o libbigint.a
|
||||||
987
bigint/src/bigint.c
Normal file
987
bigint/src/bigint.c
Normal file
@ -0,0 +1,987 @@
|
|||||||
|
/*
|
||||||
|
* bigint.c -- big integer library
|
||||||
|
*/
|
||||||
|
|
||||||
|
|
||||||
|
#include <stdio.h>
|
||||||
|
#include <ctype.h>
|
||||||
|
|
||||||
|
#include "bigint.h"
|
||||||
|
|
||||||
|
/**************************************************************/
|
||||||
|
|
||||||
|
/* debugging */
|
||||||
|
|
||||||
|
|
||||||
|
#define DIV_CHK_01 0
|
||||||
|
#define DIV_CHK_02 0
|
||||||
|
#define DIV_CHK_03 0
|
||||||
|
#define DIV_CHK_04 0
|
||||||
|
#define DIV_CHK_05 0
|
||||||
|
#define DIV_CHK_06 0
|
||||||
|
#define DIV_CHK_07 0
|
||||||
|
#define DIV_CHK_08 0
|
||||||
|
#define DIV_CHK_09 0
|
||||||
|
#define DIV_CHK_10 0
|
||||||
|
#define DIV_CHK_11 0
|
||||||
|
|
||||||
|
|
||||||
|
/**************************************************************/
|
||||||
|
|
||||||
|
/* big integer representation */
|
||||||
|
|
||||||
|
|
||||||
|
#define BIG_NEGATIVE ((unsigned char) 0)
|
||||||
|
#define BIG_POSITIVE ((unsigned char) 1)
|
||||||
|
|
||||||
|
|
||||||
|
#define BIG_PTR(bigObjRef) ((Big *) (getPrimObjectDataPointer(bigObjRef)))
|
||||||
|
#define GET_ND(bigObjRef) (BIG_PTR(bigObjRef)->nd)
|
||||||
|
#define SET_ND(bigObjRef, val) (BIG_PTR(bigObjRef)->nd = (val))
|
||||||
|
#define GET_SIGN(bigObjRef) (BIG_PTR(bigObjRef)->sign)
|
||||||
|
#define SET_SIGN(bigObjRef, val) (BIG_PTR(bigObjRef)->sign = (val))
|
||||||
|
#define GET_DIGIT(bigObjRef, i) (BIG_PTR(bigObjRef)->digits[i])
|
||||||
|
#define SET_DIGIT(bigObjRef, i, val) (BIG_PTR(bigObjRef)->digits[i] = (val))
|
||||||
|
|
||||||
|
|
||||||
|
/**************************************************************/
|
||||||
|
|
||||||
|
/* global data */
|
||||||
|
|
||||||
|
|
||||||
|
/*
|
||||||
|
* registers of the big integer processor
|
||||||
|
*/
|
||||||
|
BIP bip = {
|
||||||
|
NULL, /* op1 */
|
||||||
|
NULL, /* op2 */
|
||||||
|
NULL, /* res */
|
||||||
|
NULL, /* rem */
|
||||||
|
};
|
||||||
|
|
||||||
|
|
||||||
|
/**************************************************************/
|
||||||
|
|
||||||
|
|
||||||
|
/*
|
||||||
|
* construct a new big integer object
|
||||||
|
*
|
||||||
|
* number of digits is given by parameter
|
||||||
|
* a reference to a proper object is returned
|
||||||
|
* but no component of the big integer is set
|
||||||
|
*
|
||||||
|
* ATTENTION: All object references stored in
|
||||||
|
* places other than the bip registers may become
|
||||||
|
* invalid as soon as this function is called!
|
||||||
|
*/
|
||||||
|
static BigObjRef newBig(int nd) {
|
||||||
|
int dataSize;
|
||||||
|
BigObjRef bigObjRef;
|
||||||
|
|
||||||
|
dataSize = sizeof(int) + 1 + nd;
|
||||||
|
bigObjRef = newPrimObject(dataSize);
|
||||||
|
return bigObjRef;
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
/**************************************************************/
|
||||||
|
|
||||||
|
/* big integer unsigned arithmetic */
|
||||||
|
|
||||||
|
|
||||||
|
/*
|
||||||
|
* exchange bip.op1 and bip.op2
|
||||||
|
*/
|
||||||
|
static void bigXchg(void) {
|
||||||
|
BigObjRef tmp;
|
||||||
|
|
||||||
|
tmp = bip.op1;
|
||||||
|
bip.op1 = bip.op2;
|
||||||
|
bip.op2 = tmp;
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
/*
|
||||||
|
* big integer unsigned comparison
|
||||||
|
*
|
||||||
|
* operands in bip.op1 and bip.op2
|
||||||
|
* result is < 0, = 0, or > 0 if and only if the
|
||||||
|
* same relation holds for bip.op1 and bip.op2
|
||||||
|
*/
|
||||||
|
static int bigUcmp(void) {
|
||||||
|
int nd1;
|
||||||
|
int nd2;
|
||||||
|
int diff;
|
||||||
|
|
||||||
|
/* compare sizes */
|
||||||
|
nd1 = GET_ND(bip.op1);
|
||||||
|
nd2 = GET_ND(bip.op2);
|
||||||
|
if (nd1 != nd2) {
|
||||||
|
/* sizes are different: we know the bigger number */
|
||||||
|
return nd1 - nd2;
|
||||||
|
}
|
||||||
|
/* sizes are equal: we must look at the digits */
|
||||||
|
while (nd1--) {
|
||||||
|
diff = (int) GET_DIGIT(bip.op1, nd1) -
|
||||||
|
(int) GET_DIGIT(bip.op2, nd1);
|
||||||
|
if (diff != 0) {
|
||||||
|
return diff;
|
||||||
|
}
|
||||||
|
}
|
||||||
|
/* the numbers are equal */
|
||||||
|
return 0;
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
/*
|
||||||
|
* big integer unsigned addition
|
||||||
|
*
|
||||||
|
* operands in bip.op1 and bip.op2
|
||||||
|
* result in bip.res
|
||||||
|
*/
|
||||||
|
static void bigUadd(void) {
|
||||||
|
int nd1;
|
||||||
|
int nd2;
|
||||||
|
int i;
|
||||||
|
unsigned short carry;
|
||||||
|
unsigned short aux;
|
||||||
|
int xchgFlag;
|
||||||
|
|
||||||
|
/* make sure op1 has at least as many digits as op2 */
|
||||||
|
nd1 = GET_ND(bip.op1);
|
||||||
|
nd2 = GET_ND(bip.op2);
|
||||||
|
if (nd1 < nd2) {
|
||||||
|
/* exchange operands */
|
||||||
|
bigXchg();
|
||||||
|
i = nd1;
|
||||||
|
nd1 = nd2;
|
||||||
|
nd2 = i;
|
||||||
|
xchgFlag = 1;
|
||||||
|
} else {
|
||||||
|
/* don't exchange operands */
|
||||||
|
xchgFlag = 0;
|
||||||
|
}
|
||||||
|
/* allocate result */
|
||||||
|
bip.res = newBig(nd1 + 1);
|
||||||
|
/* copy op2 to result */
|
||||||
|
for (i = 0; i < nd2; i++) {
|
||||||
|
SET_DIGIT(bip.res, i, GET_DIGIT(bip.op2, i));
|
||||||
|
}
|
||||||
|
/* fill result with 0 up to size of op1 */
|
||||||
|
for (; i < nd1; i++) {
|
||||||
|
SET_DIGIT(bip.res, i, 0);
|
||||||
|
}
|
||||||
|
/* res = op1 + res */
|
||||||
|
carry = 0x00;
|
||||||
|
for (i = 0; i < nd1; i++) {
|
||||||
|
aux = (unsigned short) GET_DIGIT(bip.op1, i) +
|
||||||
|
(unsigned short) GET_DIGIT(bip.res, i) +
|
||||||
|
carry;
|
||||||
|
SET_DIGIT(bip.res, i, aux & 0xFF);
|
||||||
|
carry = aux >> 8;
|
||||||
|
}
|
||||||
|
SET_DIGIT(bip.res, i, carry);
|
||||||
|
/* determine actual size of result */
|
||||||
|
i = nd1 + 1;
|
||||||
|
while (--i >= 0 && GET_DIGIT(bip.res, i) == 0) ;
|
||||||
|
SET_ND(bip.res, i + 1);
|
||||||
|
/* restore operands */
|
||||||
|
if (xchgFlag) {
|
||||||
|
bigXchg();
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
/*
|
||||||
|
* big integer unsigned subtraction
|
||||||
|
*
|
||||||
|
* operands in bip.op1 and bip.op2
|
||||||
|
* result in bip.res, must not be negative
|
||||||
|
*/
|
||||||
|
static void bigUsub(void) {
|
||||||
|
int nd1;
|
||||||
|
int nd2;
|
||||||
|
int i;
|
||||||
|
unsigned short carry;
|
||||||
|
unsigned short aux;
|
||||||
|
|
||||||
|
/* op1 must have at least as many digits as op2 */
|
||||||
|
nd1 = GET_ND(bip.op1);
|
||||||
|
nd2 = GET_ND(bip.op2);
|
||||||
|
if (nd1 < nd2) {
|
||||||
|
/* unsigned subtraction would yield negative result */
|
||||||
|
fatalError("internal library error #1 - THIS SHOULD NEVER HAPPEN!");
|
||||||
|
}
|
||||||
|
/* allocate result */
|
||||||
|
bip.res = newBig(nd1);
|
||||||
|
/* copy op2 to result */
|
||||||
|
for (i = 0; i < nd2; i++) {
|
||||||
|
SET_DIGIT(bip.res, i, GET_DIGIT(bip.op2, i));
|
||||||
|
}
|
||||||
|
/* fill result with 0 up to size of op1 */
|
||||||
|
for (; i < nd1; i++) {
|
||||||
|
SET_DIGIT(bip.res, i, 0);
|
||||||
|
}
|
||||||
|
/* res = op1 - res */
|
||||||
|
carry = 0x01;
|
||||||
|
for (i = 0; i < nd1; i++) {
|
||||||
|
aux = (unsigned short) GET_DIGIT(bip.op1, i) -
|
||||||
|
(unsigned short) GET_DIGIT(bip.res, i) +
|
||||||
|
carry + 0xFF;
|
||||||
|
SET_DIGIT(bip.res, i, aux & 0xFF);
|
||||||
|
carry = aux >> 8;
|
||||||
|
}
|
||||||
|
if (carry != 0x01) {
|
||||||
|
/* unsigned subtraction would yield negative result */
|
||||||
|
fatalError("internal library error #2 - THIS SHOULD NEVER HAPPEN!");
|
||||||
|
}
|
||||||
|
/* determine actual size of result */
|
||||||
|
i = nd1;
|
||||||
|
while (--i >= 0 && GET_DIGIT(bip.res, i) == 0) ;
|
||||||
|
SET_ND(bip.res, i + 1);
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
/*
|
||||||
|
* big integer unsigned multiplication
|
||||||
|
*
|
||||||
|
* operands in bip.op1 and bip.op2
|
||||||
|
* result in bip.res
|
||||||
|
*/
|
||||||
|
static void bigUmul(void) {
|
||||||
|
int nd1;
|
||||||
|
int nd2;
|
||||||
|
int i, j, k;
|
||||||
|
unsigned short carry;
|
||||||
|
unsigned short aux;
|
||||||
|
|
||||||
|
/* get sizes of operands */
|
||||||
|
nd1 = GET_ND(bip.op1);
|
||||||
|
nd2 = GET_ND(bip.op2);
|
||||||
|
/* allocate result */
|
||||||
|
bip.res = newBig(nd1 + nd2);
|
||||||
|
/* reset lower nd1 digits of result */
|
||||||
|
for (i = 0; i < nd1; i++) {
|
||||||
|
SET_DIGIT(bip.res, i, 0);
|
||||||
|
}
|
||||||
|
/* res = op1 * op2 */
|
||||||
|
for (j = 0; j < nd2; j++) {
|
||||||
|
carry = 0x00;
|
||||||
|
for (k = j, i = 0; i < nd1; k++, i++) {
|
||||||
|
aux = (unsigned short) GET_DIGIT(bip.op1, i) *
|
||||||
|
(unsigned short) GET_DIGIT(bip.op2, j) +
|
||||||
|
(unsigned short) GET_DIGIT(bip.res, k) +
|
||||||
|
carry;
|
||||||
|
SET_DIGIT(bip.res, k, aux & 0xFF);
|
||||||
|
carry = aux >> 8;
|
||||||
|
}
|
||||||
|
SET_DIGIT(bip.res, k, carry);
|
||||||
|
}
|
||||||
|
/* determine actual size of result */
|
||||||
|
i = nd1 + nd2;
|
||||||
|
while (--i >= 0 && GET_DIGIT(bip.res, i) == 0) ;
|
||||||
|
SET_ND(bip.res, i + 1);
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
/*
|
||||||
|
* big integer unsigned division by single digit divisor
|
||||||
|
*
|
||||||
|
* dividend in bip.rem, divisor in parameter
|
||||||
|
* quotient in bip.rem, remainder is returned
|
||||||
|
*/
|
||||||
|
static unsigned char bigUdiv1(unsigned char divisor) {
|
||||||
|
BigObjRef tmp;
|
||||||
|
int nd;
|
||||||
|
int i;
|
||||||
|
unsigned short d, r;
|
||||||
|
unsigned short aux;
|
||||||
|
|
||||||
|
/* get size of dividend */
|
||||||
|
nd = GET_ND(bip.rem);
|
||||||
|
/* check for division by zero */
|
||||||
|
d = (unsigned short) divisor;
|
||||||
|
if (d == 0) {
|
||||||
|
fatalError("internal library error #3 - THIS SHOULD NEVER HAPPEN!");
|
||||||
|
}
|
||||||
|
/* allocate result */
|
||||||
|
tmp = newBig(nd);
|
||||||
|
/* tmp = dividend / divisor, r = dividend % divisor */
|
||||||
|
r = 0;
|
||||||
|
for (i = nd - 1; i >= 0; i--) {
|
||||||
|
aux = (r << 8) | (unsigned short) GET_DIGIT(bip.rem, i);
|
||||||
|
SET_DIGIT(tmp, i, aux / d);
|
||||||
|
r = aux % d;
|
||||||
|
}
|
||||||
|
/* determine actual size of quotient */
|
||||||
|
i = nd;
|
||||||
|
while (--i >= 0 && GET_DIGIT(tmp, i) == 0) ;
|
||||||
|
SET_ND(tmp, i + 1);
|
||||||
|
/* store quotient */
|
||||||
|
bip.rem = tmp;
|
||||||
|
/* return remainder */
|
||||||
|
return (unsigned char) r;
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
/*
|
||||||
|
* big integer unsigned division
|
||||||
|
*
|
||||||
|
* dividend in bip.op1, divisor in bip.op2
|
||||||
|
* quotient in bip.res, remainder in bip.rem
|
||||||
|
*/
|
||||||
|
static void bigUdiv(void) {
|
||||||
|
BigObjRef tmp;
|
||||||
|
int nd1;
|
||||||
|
int nd2;
|
||||||
|
int nd3;
|
||||||
|
int i, j, k, l;
|
||||||
|
unsigned char r;
|
||||||
|
unsigned short scale;
|
||||||
|
unsigned short carry;
|
||||||
|
unsigned short aux;
|
||||||
|
unsigned short qhat;
|
||||||
|
unsigned short v1, v2;
|
||||||
|
unsigned short uj0, uj1, uj2, two;
|
||||||
|
|
||||||
|
/* get sizes of operands */
|
||||||
|
nd1 = GET_ND(bip.op1);
|
||||||
|
nd2 = GET_ND(bip.op2);
|
||||||
|
/* check for division by zero */
|
||||||
|
if (nd2 == 0) {
|
||||||
|
fatalError("division by zero");
|
||||||
|
}
|
||||||
|
/* check for small dividend */
|
||||||
|
if (bigUcmp() < 0) {
|
||||||
|
/* res = 0 */
|
||||||
|
bip.res = newBig(0);
|
||||||
|
SET_ND(bip.res, 0);
|
||||||
|
/* rem = op1; BUT THIS HAS TO BE A COPY! */
|
||||||
|
bip.rem = newBig(nd1);
|
||||||
|
for (i = 0; i < nd1; i++) {
|
||||||
|
SET_DIGIT(bip.rem, i, GET_DIGIT(bip.op1, i));
|
||||||
|
}
|
||||||
|
SET_ND(bip.rem, nd1);
|
||||||
|
return;
|
||||||
|
}
|
||||||
|
/* check for single digit divisor */
|
||||||
|
if (nd2 == 1) {
|
||||||
|
/* yes - use simple division by single digit divisor */
|
||||||
|
bip.rem = bip.op1;
|
||||||
|
r = bigUdiv1(GET_DIGIT(bip.op2, 0));
|
||||||
|
bip.res = bip.rem;
|
||||||
|
if (r == 0) {
|
||||||
|
bip.rem = newBig(0);
|
||||||
|
SET_ND(bip.rem, 0);
|
||||||
|
} else {
|
||||||
|
bip.rem = newBig(1);
|
||||||
|
SET_ND(bip.rem, 1);
|
||||||
|
SET_DIGIT(bip.rem, 0, r);
|
||||||
|
}
|
||||||
|
return;
|
||||||
|
}
|
||||||
|
/*
|
||||||
|
* now for the general case
|
||||||
|
*/
|
||||||
|
#if DIV_CHK_01
|
||||||
|
printf("div_chk #01: division, general case\n");
|
||||||
|
printf(" dividend = ");
|
||||||
|
bigDump(stdout, bip.op1);
|
||||||
|
printf("\n");
|
||||||
|
printf(" divisor = ");
|
||||||
|
bigDump(stdout, bip.op2);
|
||||||
|
printf("\n");
|
||||||
|
#endif
|
||||||
|
/* determine scale factor for normalization */
|
||||||
|
scale = (unsigned short) 256 /
|
||||||
|
((unsigned short) GET_DIGIT(bip.op2, nd2 - 1) + 1);
|
||||||
|
#if DIV_CHK_02
|
||||||
|
printf("div_chk #02: scale factor = %02X\n", scale);
|
||||||
|
#endif
|
||||||
|
/* normalize dividend, result is in bip.rem */
|
||||||
|
bip.rem = newBig(nd1 + 1);
|
||||||
|
carry = 0x00;
|
||||||
|
for (i = 0; i < nd1; i++) {
|
||||||
|
aux = (unsigned short) GET_DIGIT(bip.op1, i) * scale +
|
||||||
|
carry;
|
||||||
|
SET_DIGIT(bip.rem, i, aux & 0xFF);
|
||||||
|
carry = aux >> 8;
|
||||||
|
}
|
||||||
|
SET_DIGIT(bip.rem, i, carry);
|
||||||
|
SET_ND(bip.rem, nd1 + 1);
|
||||||
|
#if DIV_CHK_03
|
||||||
|
printf("div_chk #03: normalized dividend = ");
|
||||||
|
bigDump(stdout, bip.rem);
|
||||||
|
printf("\n");
|
||||||
|
#endif
|
||||||
|
/* normalize divisor, result is in bip.res */
|
||||||
|
bip.res = newBig(nd2);
|
||||||
|
carry = 0x00;
|
||||||
|
for (i = 0; i < nd2; i++) {
|
||||||
|
aux = (unsigned short) GET_DIGIT(bip.op2, i) * scale +
|
||||||
|
carry;
|
||||||
|
SET_DIGIT(bip.res, i, aux & 0xFF);
|
||||||
|
carry = aux >> 8;
|
||||||
|
}
|
||||||
|
if (carry != 0x00) {
|
||||||
|
/* overflow in divisor normalization */
|
||||||
|
fatalError("internal library error #4 - THIS SHOULD NEVER HAPPEN!");
|
||||||
|
}
|
||||||
|
SET_ND(bip.res, nd2);
|
||||||
|
#if DIV_CHK_04
|
||||||
|
printf("div_chk #04: normalized divisor = ");
|
||||||
|
bigDump(stdout, bip.res);
|
||||||
|
printf("\n");
|
||||||
|
#endif
|
||||||
|
/* allocate quotient */
|
||||||
|
nd3 = nd1 - nd2 + 1;
|
||||||
|
tmp = newBig(nd3);
|
||||||
|
/* extract the two most significand digits of divisor */
|
||||||
|
v1 = (unsigned short) GET_DIGIT(bip.res, nd2 - 1);
|
||||||
|
v2 = (unsigned short) GET_DIGIT(bip.res, nd2 - 2);
|
||||||
|
/* loop on digits of dividend and compute digits of quotient */
|
||||||
|
/* j is index into dividend, k is index into quotient */
|
||||||
|
for (j = nd1, k = nd3 - 1; k >= 0; j--, k--) {
|
||||||
|
#if DIV_CHK_05
|
||||||
|
printf("div_chk #05: j = %d, k = %d\n", j, k);
|
||||||
|
#endif
|
||||||
|
/* calculate qhat */
|
||||||
|
uj0 = (unsigned short) GET_DIGIT(bip.rem, j);
|
||||||
|
uj1 = (unsigned short) GET_DIGIT(bip.rem, j - 1);
|
||||||
|
uj2 = (unsigned short) GET_DIGIT(bip.rem, j - 2);
|
||||||
|
two = (uj0 << 8) | uj1;
|
||||||
|
if (uj0 == v1) {
|
||||||
|
qhat = (unsigned short) 255;
|
||||||
|
#if DIV_CHK_06
|
||||||
|
printf("div_chk #06a: qhat = %02X\n", qhat);
|
||||||
|
#endif
|
||||||
|
} else {
|
||||||
|
qhat = two / v1;
|
||||||
|
#if DIV_CHK_06
|
||||||
|
printf("div_chk #06b: qhat = %02X\n", qhat);
|
||||||
|
#endif
|
||||||
|
}
|
||||||
|
while (qhat * v2 > (((two - qhat * v1) << 8) | uj2)) {
|
||||||
|
qhat--;
|
||||||
|
#if DIV_CHK_07
|
||||||
|
printf("div_chk #07: qhat decremented, is now %02X\n", qhat);
|
||||||
|
#endif
|
||||||
|
}
|
||||||
|
/* multiply and subtract */
|
||||||
|
/* l is index into dividend, i is index into divisor */
|
||||||
|
carry = 0xFF;
|
||||||
|
for (l = j - nd2, i = 0; i < nd2; l++, i++) {
|
||||||
|
aux = (unsigned short) GET_DIGIT(bip.rem, l) -
|
||||||
|
(unsigned short) GET_DIGIT(bip.res, i) * qhat +
|
||||||
|
carry + 0xFE01;
|
||||||
|
SET_DIGIT(bip.rem, l, aux & 0xFF);
|
||||||
|
carry = aux >> 8;
|
||||||
|
}
|
||||||
|
aux = (unsigned short) GET_DIGIT(bip.rem, l) +
|
||||||
|
carry + 0xFE01;
|
||||||
|
SET_DIGIT(bip.rem, l, aux & 0xFF);
|
||||||
|
carry = aux >> 8;
|
||||||
|
#if DIV_CHK_08
|
||||||
|
printf("div_chk #08: remainder = ");
|
||||||
|
bigDump(stdout, bip.rem);
|
||||||
|
printf("\n");
|
||||||
|
#endif
|
||||||
|
/* test remainder and possibly add back */
|
||||||
|
if (carry != 0xFF) {
|
||||||
|
/* qhat is one too large */
|
||||||
|
qhat--;
|
||||||
|
#if DIV_CHK_09
|
||||||
|
printf("div_chk #09: qhat final correction, is now %02X\n", qhat);
|
||||||
|
#endif
|
||||||
|
/* add back */
|
||||||
|
/* l is index into dividend, i is index into divisor */
|
||||||
|
carry = 0x00;
|
||||||
|
for (l = j - nd2, i = 0; i < nd2; l++, i++) {
|
||||||
|
aux = (unsigned short) GET_DIGIT(bip.rem, l) +
|
||||||
|
(unsigned short) GET_DIGIT(bip.res, i) +
|
||||||
|
carry;
|
||||||
|
SET_DIGIT(bip.rem, l, aux & 0xFF);
|
||||||
|
carry = aux >> 8;
|
||||||
|
}
|
||||||
|
aux = (unsigned short) GET_DIGIT(bip.rem, l) +
|
||||||
|
carry;
|
||||||
|
SET_DIGIT(bip.rem, l, aux & 0xFF);
|
||||||
|
carry = aux >> 8;
|
||||||
|
if (carry != 0x01) {
|
||||||
|
/* missing carry in add-back sum */
|
||||||
|
fatalError("internal library error #5 - THIS SHOULD NEVER HAPPEN!");
|
||||||
|
}
|
||||||
|
#if DIV_CHK_10
|
||||||
|
printf("div_chk #10: remainder = ");
|
||||||
|
bigDump(stdout, bip.rem);
|
||||||
|
printf("\n");
|
||||||
|
#endif
|
||||||
|
}
|
||||||
|
/* store quotient digit */
|
||||||
|
SET_DIGIT(tmp, k, qhat);
|
||||||
|
#if DIV_CHK_11
|
||||||
|
printf("div_chk #11: quotient digit = %02X\n", qhat);
|
||||||
|
#endif
|
||||||
|
}
|
||||||
|
/* finish quotient */
|
||||||
|
i = nd3;
|
||||||
|
while (--i >= 0 && GET_DIGIT(tmp, i) == 0) ;
|
||||||
|
SET_ND(tmp, i + 1);
|
||||||
|
bip.res = tmp;
|
||||||
|
/* finish and unnormalize remainder */
|
||||||
|
i = nd1 + 1;
|
||||||
|
while (--i >= 0 && GET_DIGIT(bip.rem, i) == 0) ;
|
||||||
|
SET_ND(bip.rem, i + 1);
|
||||||
|
r = bigUdiv1(scale);
|
||||||
|
if (r != 0) {
|
||||||
|
/* non-zero remainder in unnormalization */
|
||||||
|
fatalError("internal library error #6 - THIS SHOULD NEVER HAPPEN!");
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
/**************************************************************/
|
||||||
|
|
||||||
|
/* nil reference exception */
|
||||||
|
|
||||||
|
|
||||||
|
static void nilRefException(void) {
|
||||||
|
fatalError("big integer library detected illegal nil reference");
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
/**************************************************************/
|
||||||
|
|
||||||
|
/* big integer arithmetic */
|
||||||
|
|
||||||
|
|
||||||
|
/*
|
||||||
|
* big integer sign
|
||||||
|
*
|
||||||
|
* operand in bip.op1
|
||||||
|
* result is < 0, = 0, or > 0 if and only if
|
||||||
|
* the same relation holds for bip.op1
|
||||||
|
*/
|
||||||
|
int bigSgn(void) {
|
||||||
|
if (bip.op1 == NULL) {
|
||||||
|
nilRefException();
|
||||||
|
}
|
||||||
|
if (GET_ND(bip.op1) == 0) {
|
||||||
|
return 0;
|
||||||
|
}
|
||||||
|
if (GET_SIGN(bip.op1) == BIG_POSITIVE) {
|
||||||
|
return 1;
|
||||||
|
} else {
|
||||||
|
return -1;
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
/*
|
||||||
|
* big integer comparison
|
||||||
|
*
|
||||||
|
* operands in bip.op1 and bip.op2
|
||||||
|
* result is < 0, = 0, or > 0 if and only if the
|
||||||
|
* same relation holds for bip.op1 and bip.op2
|
||||||
|
*/
|
||||||
|
int bigCmp(void) {
|
||||||
|
if (bip.op1 == NULL ||
|
||||||
|
bip.op2 == NULL) {
|
||||||
|
nilRefException();
|
||||||
|
}
|
||||||
|
if (GET_SIGN(bip.op1) == BIG_POSITIVE) {
|
||||||
|
if (GET_SIGN(bip.op2) == BIG_POSITIVE) {
|
||||||
|
/* op1 >= 0 and op2 >= 0 */
|
||||||
|
return bigUcmp();
|
||||||
|
} else {
|
||||||
|
/* op1 >= 0 and op2 < 0 */
|
||||||
|
return 1;
|
||||||
|
}
|
||||||
|
} else {
|
||||||
|
if (GET_SIGN(bip.op2) == BIG_POSITIVE) {
|
||||||
|
/* op1 < 0 and op2 >= 0 */
|
||||||
|
return -1;
|
||||||
|
} else {
|
||||||
|
/* op1 < 0 and op2 < 0 */
|
||||||
|
return -bigUcmp();
|
||||||
|
}
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
/*
|
||||||
|
* big integer negation
|
||||||
|
*
|
||||||
|
* operand in bip.op1
|
||||||
|
* result in bip.res
|
||||||
|
*/
|
||||||
|
void bigNeg(void) {
|
||||||
|
int nd;
|
||||||
|
int i;
|
||||||
|
|
||||||
|
if (bip.op1 == NULL) {
|
||||||
|
nilRefException();
|
||||||
|
}
|
||||||
|
/* make copy of operand */
|
||||||
|
nd = GET_ND(bip.op1);
|
||||||
|
bip.res = newBig(nd);
|
||||||
|
for (i = 0; i < nd; i++) {
|
||||||
|
SET_DIGIT(bip.res, i, GET_DIGIT(bip.op1, i));
|
||||||
|
}
|
||||||
|
SET_ND(bip.res, nd);
|
||||||
|
/* store inverted sign */
|
||||||
|
if (GET_SIGN(bip.op1) == BIG_NEGATIVE || nd == 0) {
|
||||||
|
SET_SIGN(bip.res, BIG_POSITIVE);
|
||||||
|
} else {
|
||||||
|
SET_SIGN(bip.res, BIG_NEGATIVE);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
/*
|
||||||
|
* big integer addition
|
||||||
|
*
|
||||||
|
* operands in bip.op1 and bip.op2
|
||||||
|
* result in bip.res
|
||||||
|
*/
|
||||||
|
void bigAdd(void) {
|
||||||
|
if (bip.op1 == NULL ||
|
||||||
|
bip.op2 == NULL) {
|
||||||
|
nilRefException();
|
||||||
|
}
|
||||||
|
if (GET_SIGN(bip.op1) == BIG_POSITIVE) {
|
||||||
|
if (GET_SIGN(bip.op2) == BIG_POSITIVE) {
|
||||||
|
/* op1 >= 0 and op2 >= 0 */
|
||||||
|
bigUadd();
|
||||||
|
SET_SIGN(bip.res, BIG_POSITIVE);
|
||||||
|
} else {
|
||||||
|
/* op1 >= 0 and op2 < 0 */
|
||||||
|
if (bigUcmp() >= 0) {
|
||||||
|
/* |op1| >= |op2| */
|
||||||
|
bigUsub();
|
||||||
|
SET_SIGN(bip.res, BIG_POSITIVE);
|
||||||
|
} else {
|
||||||
|
/* |op1| < |op2| */
|
||||||
|
bigXchg();
|
||||||
|
bigUsub();
|
||||||
|
SET_SIGN(bip.res, BIG_NEGATIVE);
|
||||||
|
bigXchg();
|
||||||
|
}
|
||||||
|
}
|
||||||
|
} else {
|
||||||
|
if (GET_SIGN(bip.op2) == BIG_POSITIVE) {
|
||||||
|
/* op1 < 0 and op2 >= 0 */
|
||||||
|
if (bigUcmp() <= 0) {
|
||||||
|
/* |op1| <= |op2| */
|
||||||
|
bigXchg();
|
||||||
|
bigUsub();
|
||||||
|
SET_SIGN(bip.res, BIG_POSITIVE);
|
||||||
|
bigXchg();
|
||||||
|
} else {
|
||||||
|
/* |op1| > |op2| */
|
||||||
|
bigUsub();
|
||||||
|
SET_SIGN(bip.res, BIG_NEGATIVE);
|
||||||
|
}
|
||||||
|
} else {
|
||||||
|
/* op1 < 0 and op2 < 0 */
|
||||||
|
bigUadd();
|
||||||
|
SET_SIGN(bip.res, BIG_NEGATIVE);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
/*
|
||||||
|
* big integer subtraction
|
||||||
|
*
|
||||||
|
* operands in bip.op1 and bip.op2
|
||||||
|
* result in bip.res
|
||||||
|
*/
|
||||||
|
void bigSub(void) {
|
||||||
|
if (bip.op1 == NULL ||
|
||||||
|
bip.op2 == NULL) {
|
||||||
|
nilRefException();
|
||||||
|
}
|
||||||
|
if (GET_SIGN(bip.op1) == BIG_POSITIVE) {
|
||||||
|
if (GET_SIGN(bip.op2) == BIG_POSITIVE) {
|
||||||
|
/* op1 >= 0 and op2 >= 0 */
|
||||||
|
if (bigUcmp() >= 0) {
|
||||||
|
/* |op1| >= |op2| */
|
||||||
|
bigUsub();
|
||||||
|
SET_SIGN(bip.res, BIG_POSITIVE);
|
||||||
|
} else {
|
||||||
|
/* |op1| < |op2| */
|
||||||
|
bigXchg();
|
||||||
|
bigUsub();
|
||||||
|
SET_SIGN(bip.res, BIG_NEGATIVE);
|
||||||
|
bigXchg();
|
||||||
|
}
|
||||||
|
} else {
|
||||||
|
/* op1 >= 0 and op2 < 0 */
|
||||||
|
bigUadd();
|
||||||
|
SET_SIGN(bip.res, BIG_POSITIVE);
|
||||||
|
}
|
||||||
|
} else {
|
||||||
|
if (GET_SIGN(bip.op2) == BIG_POSITIVE) {
|
||||||
|
/* op1 < 0 and op2 >= 0 */
|
||||||
|
bigUadd();
|
||||||
|
SET_SIGN(bip.res, BIG_NEGATIVE);
|
||||||
|
} else {
|
||||||
|
/* op1 < 0 and op2 < 0 */
|
||||||
|
if (bigUcmp() <= 0) {
|
||||||
|
/* |op1| <= |op2| */
|
||||||
|
bigXchg();
|
||||||
|
bigUsub();
|
||||||
|
SET_SIGN(bip.res, BIG_POSITIVE);
|
||||||
|
bigXchg();
|
||||||
|
} else {
|
||||||
|
/* |op1| > |op2| */
|
||||||
|
bigUsub();
|
||||||
|
SET_SIGN(bip.res, BIG_NEGATIVE);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
/*
|
||||||
|
* big integer multiplication
|
||||||
|
*
|
||||||
|
* operands in bip.op1 and bip.op2
|
||||||
|
* result in bip.res
|
||||||
|
*/
|
||||||
|
void bigMul(void) {
|
||||||
|
if (bip.op1 == NULL ||
|
||||||
|
bip.op2 == NULL) {
|
||||||
|
nilRefException();
|
||||||
|
}
|
||||||
|
bigUmul();
|
||||||
|
if (GET_SIGN(bip.op1) == GET_SIGN(bip.op2) || GET_ND(bip.res) == 0) {
|
||||||
|
SET_SIGN(bip.res, BIG_POSITIVE);
|
||||||
|
} else {
|
||||||
|
SET_SIGN(bip.res, BIG_NEGATIVE);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
/*
|
||||||
|
* big integer division, truncating towards zero
|
||||||
|
*
|
||||||
|
* dividend in bip.op1, divisor in bip.op2
|
||||||
|
* quotient in bip.res, remainder in bip.rem
|
||||||
|
*/
|
||||||
|
void bigDiv(void) {
|
||||||
|
if (bip.op1 == NULL ||
|
||||||
|
bip.op2 == NULL) {
|
||||||
|
nilRefException();
|
||||||
|
}
|
||||||
|
bigUdiv();
|
||||||
|
if (GET_SIGN(bip.op1) == GET_SIGN(bip.op2) || GET_ND(bip.res) == 0) {
|
||||||
|
SET_SIGN(bip.res, BIG_POSITIVE);
|
||||||
|
} else {
|
||||||
|
SET_SIGN(bip.res, BIG_NEGATIVE);
|
||||||
|
}
|
||||||
|
if (GET_SIGN(bip.op1) == BIG_POSITIVE || GET_ND(bip.rem) == 0) {
|
||||||
|
SET_SIGN(bip.rem, BIG_POSITIVE);
|
||||||
|
} else {
|
||||||
|
SET_SIGN(bip.rem, BIG_NEGATIVE);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
/**************************************************************/
|
||||||
|
|
||||||
|
/* big integer conversions */
|
||||||
|
|
||||||
|
|
||||||
|
/*
|
||||||
|
* conversion int --> big
|
||||||
|
*
|
||||||
|
* operand in parameter
|
||||||
|
* result in bip.res
|
||||||
|
*/
|
||||||
|
void bigFromInt(int n) {
|
||||||
|
int i;
|
||||||
|
|
||||||
|
bip.res = newBig(sizeof(int));
|
||||||
|
if (n < 0) {
|
||||||
|
n = -n;
|
||||||
|
SET_SIGN(bip.res, BIG_NEGATIVE);
|
||||||
|
} else {
|
||||||
|
SET_SIGN(bip.res, BIG_POSITIVE);
|
||||||
|
}
|
||||||
|
for (i = 0; i < sizeof(int); i++) {
|
||||||
|
SET_DIGIT(bip.res, i, n & 0xFF);
|
||||||
|
n >>= 8;
|
||||||
|
}
|
||||||
|
while (--i >= 0 && GET_DIGIT(bip.res, i) == 0) ;
|
||||||
|
SET_ND(bip.res, i + 1);
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
/*
|
||||||
|
* conversion big --> int
|
||||||
|
*
|
||||||
|
* operand in bip.op1
|
||||||
|
* result is returned
|
||||||
|
*/
|
||||||
|
int bigToInt(void) {
|
||||||
|
int nd;
|
||||||
|
int i;
|
||||||
|
int res;
|
||||||
|
|
||||||
|
if (bip.op1 == NULL) {
|
||||||
|
nilRefException();
|
||||||
|
}
|
||||||
|
nd = GET_ND(bip.op1);
|
||||||
|
if (nd > 4 ||
|
||||||
|
(nd == 4 && GET_DIGIT(bip.op1, 3) >= 0x80)) {
|
||||||
|
fatalError("big integer too big for conversion to int");
|
||||||
|
}
|
||||||
|
res = 0;
|
||||||
|
for (i = nd - 1; i >= 0; i--) {
|
||||||
|
res <<= 8;
|
||||||
|
res |= (unsigned int) GET_DIGIT(bip.op1, i);
|
||||||
|
}
|
||||||
|
if (GET_SIGN(bip.op1) == BIG_NEGATIVE) {
|
||||||
|
res = -res;
|
||||||
|
}
|
||||||
|
return res;
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
/**************************************************************/
|
||||||
|
|
||||||
|
/* big integer I/O */
|
||||||
|
|
||||||
|
|
||||||
|
/*
|
||||||
|
* read a big integer
|
||||||
|
*
|
||||||
|
* stream to read from in parameter
|
||||||
|
* result in bip.res
|
||||||
|
*/
|
||||||
|
void bigRead(FILE *in) {
|
||||||
|
int c;
|
||||||
|
int positive;
|
||||||
|
|
||||||
|
c = fgetc(in);
|
||||||
|
while (isspace(c)) {
|
||||||
|
c = fgetc(in);
|
||||||
|
}
|
||||||
|
if (c == '-') {
|
||||||
|
positive = 0;
|
||||||
|
c = fgetc(in);
|
||||||
|
} else {
|
||||||
|
positive = 1;
|
||||||
|
if (c == '+') {
|
||||||
|
c = fgetc(in);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
if (!isdigit(c)) {
|
||||||
|
fatalError("no digits in input");
|
||||||
|
}
|
||||||
|
bigFromInt(10);
|
||||||
|
bip.rem = bip.res;
|
||||||
|
bigFromInt(0);
|
||||||
|
while (isdigit(c)) {
|
||||||
|
bip.op1 = bip.res;
|
||||||
|
bip.op2 = bip.rem;
|
||||||
|
bigUmul();
|
||||||
|
bip.op1 = bip.res;
|
||||||
|
bigFromInt(c - '0');
|
||||||
|
bip.op2 = bip.res;
|
||||||
|
bigUadd();
|
||||||
|
c = fgetc(in);
|
||||||
|
}
|
||||||
|
ungetc(c, in);
|
||||||
|
if (positive || GET_ND(bip.res) == 0) {
|
||||||
|
SET_SIGN(bip.res, BIG_POSITIVE);
|
||||||
|
} else {
|
||||||
|
SET_SIGN(bip.res, BIG_NEGATIVE);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
/*
|
||||||
|
* print a big integer
|
||||||
|
*
|
||||||
|
* stream to write to in parameter
|
||||||
|
* number to print in bip.op1
|
||||||
|
*/
|
||||||
|
void bigPrint(FILE *out) {
|
||||||
|
int nd;
|
||||||
|
unsigned char r;
|
||||||
|
int skipZero;
|
||||||
|
|
||||||
|
if (bip.op1 == NULL) {
|
||||||
|
nilRefException();
|
||||||
|
}
|
||||||
|
nd = GET_ND(bip.op1);
|
||||||
|
if (nd == 0) {
|
||||||
|
fprintf(out, "0");
|
||||||
|
return;
|
||||||
|
}
|
||||||
|
if (GET_SIGN(bip.op1) == BIG_NEGATIVE) {
|
||||||
|
fprintf(out, "-");
|
||||||
|
}
|
||||||
|
/* number of digits in base 10 = number of digits
|
||||||
|
in base 256 * log10(256), and log10(256) < 2.5 */
|
||||||
|
nd = 2 * nd + nd / 2;
|
||||||
|
bip.rem = bip.op1;
|
||||||
|
bigFromInt(10);
|
||||||
|
bip.op2 = bip.res;
|
||||||
|
bigFromInt(1);
|
||||||
|
while (nd != 0) {
|
||||||
|
bip.op1 = bip.res;
|
||||||
|
bigUmul();
|
||||||
|
nd--;
|
||||||
|
}
|
||||||
|
bip.op1 = bip.rem;
|
||||||
|
bip.op2 = bip.res;
|
||||||
|
skipZero = 1;
|
||||||
|
do {
|
||||||
|
bigUdiv();
|
||||||
|
if (GET_ND(bip.res) == 0) {
|
||||||
|
if (!skipZero) {
|
||||||
|
fprintf(out, "0");
|
||||||
|
}
|
||||||
|
} else {
|
||||||
|
if (GET_ND(bip.res) != 1) {
|
||||||
|
fatalError("internal library error #7 - THIS SHOULD NEVER HAPPEN!");
|
||||||
|
}
|
||||||
|
fprintf(out, "%c", GET_DIGIT(bip.res, 0) + '0');
|
||||||
|
skipZero = 0;
|
||||||
|
}
|
||||||
|
bip.op1 = bip.rem;
|
||||||
|
bip.rem = bip.op2;
|
||||||
|
r = bigUdiv1(10);
|
||||||
|
bip.op2 = bip.rem;
|
||||||
|
} while (r == 0);
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
/**************************************************************/
|
||||||
|
|
||||||
|
/* debugging */
|
||||||
|
|
||||||
|
|
||||||
|
/*
|
||||||
|
* dump a big integer object
|
||||||
|
*/
|
||||||
|
void bigDump(FILE *out, BigObjRef bigObjRef) {
|
||||||
|
int nd;
|
||||||
|
unsigned char sign;
|
||||||
|
int i;
|
||||||
|
|
||||||
|
if (bigObjRef == NULL) {
|
||||||
|
nilRefException();
|
||||||
|
}
|
||||||
|
nd = GET_ND(bigObjRef);
|
||||||
|
sign = GET_SIGN(bigObjRef);
|
||||||
|
fprintf(out, "[%d %c", nd, sign == BIG_POSITIVE ? '+' : '-');
|
||||||
|
for (i = 0; i < nd; i++) {
|
||||||
|
fprintf(out, " %02X", GET_DIGIT(bigObjRef, i));
|
||||||
|
}
|
||||||
|
fprintf(out, "]");
|
||||||
|
}
|
||||||
60
bigint/src/bigint.h
Normal file
60
bigint/src/bigint.h
Normal file
@ -0,0 +1,60 @@
|
|||||||
|
/*
|
||||||
|
* bigint.h -- big integer library
|
||||||
|
*/
|
||||||
|
|
||||||
|
|
||||||
|
#ifndef _BIGINT_H_
|
||||||
|
#define _BIGINT_H_
|
||||||
|
|
||||||
|
|
||||||
|
/* object representation */
|
||||||
|
typedef void* BigObjRef;
|
||||||
|
|
||||||
|
|
||||||
|
#include <stdio.h>
|
||||||
|
|
||||||
|
|
||||||
|
typedef struct {
|
||||||
|
int nd; /* number of digits; array may be bigger */
|
||||||
|
/* nd = 0 exactly when number = 0 */
|
||||||
|
unsigned char sign; /* one of BIG_NEGATIVE or BIG_POSITIVE */
|
||||||
|
/* zero always has BIG_POSITIVE here */
|
||||||
|
unsigned char digits[1]; /* the digits proper; number base is 256 */
|
||||||
|
/* LS digit first; MS digit is not zero */
|
||||||
|
} Big;
|
||||||
|
|
||||||
|
#include "support.h"
|
||||||
|
|
||||||
|
|
||||||
|
/* big integer processor registers */
|
||||||
|
|
||||||
|
typedef struct {
|
||||||
|
BigObjRef op1; /* first (or single) operand */
|
||||||
|
BigObjRef op2; /* second operand (if present) */
|
||||||
|
BigObjRef res; /* result of operation */
|
||||||
|
BigObjRef rem; /* remainder in case of division */
|
||||||
|
} BIP;
|
||||||
|
|
||||||
|
extern BIP bip; /* registers of the processor */
|
||||||
|
|
||||||
|
|
||||||
|
/* big integer processor functions */
|
||||||
|
|
||||||
|
int bigSgn(void); /* sign */
|
||||||
|
int bigCmp(void); /* comparison */
|
||||||
|
void bigNeg(void); /* negation */
|
||||||
|
void bigAdd(void); /* addition */
|
||||||
|
void bigSub(void); /* subtraction */
|
||||||
|
void bigMul(void); /* multiplication */
|
||||||
|
void bigDiv(void); /* division */
|
||||||
|
|
||||||
|
void bigFromInt(int n); /* conversion int --> big */
|
||||||
|
int bigToInt(void); /* conversion big --> int */
|
||||||
|
|
||||||
|
void bigRead(FILE *in); /* read a big integer */
|
||||||
|
void bigPrint(FILE *out); /* print a big integer */
|
||||||
|
|
||||||
|
void bigDump(FILE *out, BigObjRef bigObjRef); /* dump a big integer object */
|
||||||
|
|
||||||
|
|
||||||
|
#endif /* _BIGINT_H_ */
|
||||||
BIN
bigint/src/bigint.o
Normal file
BIN
bigint/src/bigint.o
Normal file
Binary file not shown.
BIN
bigint/src/libbigint.a
Normal file
BIN
bigint/src/libbigint.a
Normal file
Binary file not shown.
15
bigint/src/support.h
Normal file
15
bigint/src/support.h
Normal file
@ -0,0 +1,15 @@
|
|||||||
|
/*
|
||||||
|
* support.h -- object representation and support functions
|
||||||
|
*/
|
||||||
|
|
||||||
|
|
||||||
|
#ifndef _SUPPORT_H_
|
||||||
|
#define _SUPPORT_H_
|
||||||
|
|
||||||
|
/* support functions */
|
||||||
|
|
||||||
|
void fatalError(char *msg); /* print a message and exit */
|
||||||
|
void * newPrimObject(int dataSize); /* create a new primitive object */
|
||||||
|
void * getPrimObjectDataPointer(void * primObject);
|
||||||
|
|
||||||
|
#endif /* _SUPPORT_H_ */
|
||||||
28
bigint/tst/Makefile
Normal file
28
bigint/tst/Makefile
Normal file
@ -0,0 +1,28 @@
|
|||||||
|
#
|
||||||
|
# Makefile for big integer test
|
||||||
|
#
|
||||||
|
|
||||||
|
BUILD = ../build
|
||||||
|
|
||||||
|
CC = gcc
|
||||||
|
CFLAGS = -g -Wall -I$(BUILD)/include
|
||||||
|
LDFLAGS = -g -Wall -L$(BUILD)/lib
|
||||||
|
LDLIBS = -lbigint
|
||||||
|
|
||||||
|
all: testbip
|
||||||
|
|
||||||
|
install: all
|
||||||
|
mkdir -p $(BUILD)/bin
|
||||||
|
cp testbip $(BUILD)/bin
|
||||||
|
|
||||||
|
testbip: testbip.o support.o
|
||||||
|
$(CC) $(LDFLAGS) -o testbip testbip.o support.o $(LDLIBS)
|
||||||
|
|
||||||
|
testbip.o: testbip.c
|
||||||
|
$(CC) $(CFLAGS) -o testbip.o -c testbip.c
|
||||||
|
|
||||||
|
support.o: support.c
|
||||||
|
$(CC) $(CFLAGS) -o support.o -c support.c
|
||||||
|
|
||||||
|
clean:
|
||||||
|
rm -f *~ testbip.o support.o testbip
|
||||||
58
bigint/tst/support.c
Normal file
58
bigint/tst/support.c
Normal file
@ -0,0 +1,58 @@
|
|||||||
|
/*
|
||||||
|
* support.c -- support functions for big integer library
|
||||||
|
*/
|
||||||
|
|
||||||
|
|
||||||
|
#include <stdio.h>
|
||||||
|
#include <stdlib.h>
|
||||||
|
#include <stdbool.h>
|
||||||
|
|
||||||
|
#include "support.h"
|
||||||
|
|
||||||
|
typedef struct ObjRef{
|
||||||
|
bool brokenHeart;
|
||||||
|
struct ObjRef *forward_pointer;
|
||||||
|
unsigned int size;
|
||||||
|
unsigned char data[1];
|
||||||
|
} *ObjRef;
|
||||||
|
|
||||||
|
|
||||||
|
/*
|
||||||
|
* This routine is called in case a fatal error has occurred.
|
||||||
|
* It should print the error message and terminate the program.
|
||||||
|
*/
|
||||||
|
void fatalError(char *msg) {
|
||||||
|
printf("Fatal error: %s\n", msg);
|
||||||
|
exit(1);
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
/*
|
||||||
|
* This function is called whenever a new primitive object with
|
||||||
|
* a certain amount of internal memory is needed. It should return
|
||||||
|
* an object reference to a regular object, which contains a freely
|
||||||
|
* usable memory area of at least the requested size (measured in
|
||||||
|
* bytes). The memory area need not be initialized in any way.
|
||||||
|
*
|
||||||
|
* Note that this function may move all objects in memory at will
|
||||||
|
* (due to, e.g., garbage collection), as long as the pointers in
|
||||||
|
* the global "bip" structure point to the correct objects when
|
||||||
|
* the function returns.
|
||||||
|
*/
|
||||||
|
void * newPrimObject(int dataSize) {
|
||||||
|
ObjRef bigObjRef;
|
||||||
|
int size = sizeof(unsigned int) + dataSize * sizeof(unsigned char);
|
||||||
|
bigObjRef = malloc(size);
|
||||||
|
bigObjRef->size = size;
|
||||||
|
if (bigObjRef == NULL) {
|
||||||
|
fatalError("newPrimObject() got no memory");
|
||||||
|
}
|
||||||
|
bigObjRef->size = size;
|
||||||
|
return bigObjRef;
|
||||||
|
}
|
||||||
|
|
||||||
|
void * getPrimObjectDataPointer(void * obj){
|
||||||
|
ObjRef oo = ((ObjRef) (obj));
|
||||||
|
return oo->data;
|
||||||
|
}
|
||||||
|
|
||||||
BIN
bigint/tst/support.o
Normal file
BIN
bigint/tst/support.o
Normal file
Binary file not shown.
BIN
bigint/tst/testbip
Executable file
BIN
bigint/tst/testbip
Executable file
Binary file not shown.
1253
bigint/tst/testbip.c
Normal file
1253
bigint/tst/testbip.c
Normal file
File diff suppressed because it is too large
Load Diff
BIN
bigint/tst/testbip.o
Normal file
BIN
bigint/tst/testbip.o
Normal file
Binary file not shown.
189
njvm.c
189
njvm.c
@ -8,12 +8,17 @@
|
|||||||
#include "debugMenu.c"
|
#include "debugMenu.c"
|
||||||
#include "bigint.h"
|
#include "bigint.h"
|
||||||
#include "record.c"
|
#include "record.c"
|
||||||
#include "GC.c"
|
|
||||||
#include "SDA.c"
|
#include "SDA.c"
|
||||||
|
|
||||||
// Debug
|
// Debug
|
||||||
int debug = 0;
|
int debug = 0;
|
||||||
|
|
||||||
|
// Purge Flag
|
||||||
|
int purgeFlag = false;
|
||||||
|
|
||||||
|
// Heap Size
|
||||||
|
int argVheapSizeKiB = 8192;
|
||||||
|
|
||||||
// Program
|
// Program
|
||||||
struct program program;
|
struct program program;
|
||||||
|
|
||||||
@ -22,7 +27,7 @@ struct sda sda;
|
|||||||
|
|
||||||
// Stack
|
// Stack
|
||||||
struct stack stack;
|
struct stack stack;
|
||||||
#define SIZE 64
|
#define SIZE 10000
|
||||||
|
|
||||||
//Register
|
//Register
|
||||||
struct stack reg;
|
struct stack reg;
|
||||||
@ -46,7 +51,7 @@ void execute(struct program program) {
|
|||||||
ObjRef tempObj2;
|
ObjRef tempObj2;
|
||||||
int tempInt;
|
int tempInt;
|
||||||
for (i = 0; i < *program.size; ++i) {
|
for (i = 0; i < *program.size; ++i) {
|
||||||
if (debug == 1 || bp == i) debugMenu(fp, stack, &debug, i, &bp);
|
//if (debug == 1 || bp == i) debugMenu(fp, stack, &debug, i, &bp);
|
||||||
if(debug == 1) printf("(%i)",i);
|
if(debug == 1) printf("(%i)",i);
|
||||||
switch (program.program[i] >> 24) {
|
switch (program.program[i] >> 24) {
|
||||||
case HALT:
|
case HALT:
|
||||||
@ -315,9 +320,184 @@ void execute(struct program program) {
|
|||||||
void tests(void) {
|
void tests(void) {
|
||||||
}
|
}
|
||||||
|
|
||||||
|
/*
|
||||||
|
* Garbage Collection
|
||||||
|
*/
|
||||||
|
|
||||||
|
ObjRef *heap;
|
||||||
|
int heapSizeKiB;
|
||||||
|
char *memTargetPtr; //zielspeicherzeiger
|
||||||
|
char *freeHeapPtr; //freizeiger
|
||||||
|
char *halfHeapPtr; //halbspeicherende
|
||||||
|
int gcCount = 0;
|
||||||
|
|
||||||
|
#define MSB (1 << (8 * sizeof(unsigned int) - 1))
|
||||||
|
#define IS_PRIMITIVE(objRef) (((objRef)->size & MSB) == 0)
|
||||||
|
#define GET_SIZE(objRef)((objRef)->size& ~MSB) //Gets the size of an ObjRef
|
||||||
|
#define GET_REFS(objRef)((ObjRef *)((objRef)->data)) //Pointer on first ObjRef in Object
|
||||||
|
|
||||||
|
ObjRef copyObjectToFreeMem(ObjRef orig) {
|
||||||
|
//printf("%s%d\n", "Size in CpyObj: ", GET_SIZE(orig) );
|
||||||
|
int size;
|
||||||
|
if (IS_PRIMITIVE(orig)) {
|
||||||
|
size = GET_SIZE(orig);
|
||||||
|
} else {
|
||||||
|
size = sizeof(void *) + sizeof(unsigned int) + sizeof(bool) + GET_SIZE(orig) * sizeof(ObjRef *);
|
||||||
|
}
|
||||||
|
memcpy(freeHeapPtr, orig, size);
|
||||||
|
ObjRef oldPtr = (ObjRef) freeHeapPtr;
|
||||||
|
//printf("%s%d\n","size after memcpy ", GET_SIZE(oldPtr));
|
||||||
|
//printf("%s%d\n","actual Size: ", size);
|
||||||
|
//oldPtr->size = size;
|
||||||
|
freeHeapPtr += size;
|
||||||
|
return oldPtr;
|
||||||
|
}
|
||||||
|
|
||||||
|
ObjRef relocate(ObjRef orig) {
|
||||||
|
/* if (orig == NULL) {
|
||||||
|
printf("%s\n", "null");
|
||||||
|
} else {
|
||||||
|
printf("%s%d\n","Size in relocate ", GET_SIZE(orig));
|
||||||
|
}
|
||||||
|
*/
|
||||||
|
|
||||||
|
ObjRef copy;
|
||||||
|
if (orig == NULL) {
|
||||||
|
/* relocate(nil) = nil */
|
||||||
|
copy = NULL;
|
||||||
|
} else {
|
||||||
|
if (orig->brokenHeart == 1) {
|
||||||
|
/* Objekt ist bereits kopiert */
|
||||||
|
copy = orig->forwardPointer;
|
||||||
|
//copy->size =orig->size;
|
||||||
|
//printf("%s\n", "already copied");
|
||||||
|
} else {
|
||||||
|
/* Objekt muss noch kopiert werden*/
|
||||||
|
copy = copyObjectToFreeMem(orig);
|
||||||
|
orig->brokenHeart = true;
|
||||||
|
orig->forwardPointer = copy;
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
return copy;
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
void scan(void) {
|
||||||
|
char *scanPtr = memTargetPtr;
|
||||||
|
while (scanPtr < freeHeapPtr) {
|
||||||
|
/* es gibt noch Objekte, die
|
||||||
|
* gescannt werden müssen.
|
||||||
|
*/ ObjRef scanOr = (ObjRef) scanPtr;
|
||||||
|
// printf("%d\n", GET_SIZE(scanOr));
|
||||||
|
if (!IS_PRIMITIVE(scanOr)) {
|
||||||
|
/* Ist ein Verbundobjekt */
|
||||||
|
for (int i = 0; i < GET_SIZE(scanOr); i++) {
|
||||||
|
GET_REFS(scanOr)[i] = relocate(GET_REFS(scanOr)[i]);
|
||||||
|
}
|
||||||
|
scanPtr += sizeof(bool) + (GET_SIZE(scanOr) * sizeof(ObjRef *)) + sizeof(unsigned int) + sizeof(void *);
|
||||||
|
} else {
|
||||||
|
//printf("%s %d\n", "is primitive, size: ", GET_SIZE(scanOr));
|
||||||
|
scanPtr += GET_SIZE(scanOr);
|
||||||
|
|
||||||
|
}
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
void swap() {
|
||||||
|
if (memTargetPtr == (char *) heap) {
|
||||||
|
//printf("%s\n", "swapped to second half");
|
||||||
|
freeHeapPtr = halfHeapPtr;
|
||||||
|
memTargetPtr = halfHeapPtr;
|
||||||
|
halfHeapPtr = ((char *) heap) + (heapSizeKiB * 1024);
|
||||||
|
} else {
|
||||||
|
//printf("%s\n", "swapped to first half");
|
||||||
|
freeHeapPtr = ((char *) heap);
|
||||||
|
memTargetPtr = ((char *) heap);
|
||||||
|
halfHeapPtr = ((char *) heap) + ((heapSizeKiB * 1024) / 2);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
void garbageCollector() {
|
||||||
|
char *memToPurgePtr = halfHeapPtr - ((heapSizeKiB * 1024) / 2);
|
||||||
|
swap();
|
||||||
|
//printf("%s\n", "stack");
|
||||||
|
for (int i = 0; i < stack.current; i++) {
|
||||||
|
if (stack.stack[i].isObjRef) {
|
||||||
|
stack.stack[i].u.objRef = relocate(stack.stack[i].u.objRef);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
//printf("%s\n", "globals");
|
||||||
|
for (int i = 0; i < sda.size; i++) {
|
||||||
|
sda.sda[i] = relocate(sda.sda[i]);
|
||||||
|
}
|
||||||
|
//printf("%s\n", "retVal");
|
||||||
|
//returnValue = relocate(returnValue);
|
||||||
|
|
||||||
|
//printf("%s\n", "op1");
|
||||||
|
bip.op1 = relocate(bip.op1);
|
||||||
|
ObjRef b = bip.op1;
|
||||||
|
if(b!= NULL)
|
||||||
|
//printf("%d\n", GET_SIZE(b));
|
||||||
|
|
||||||
|
//printf("%s\n", "op2");
|
||||||
|
bip.op2 = relocate(bip.op2);
|
||||||
|
b = bip.op2;
|
||||||
|
if(b!= NULL)
|
||||||
|
//printf("%d\n", GET_SIZE(b));
|
||||||
|
|
||||||
|
//printf("%s\n", "res");
|
||||||
|
bip.res = relocate(bip.res);
|
||||||
|
b = bip.res;
|
||||||
|
if(b!= NULL)
|
||||||
|
//printf("%d\n", GET_SIZE(b));
|
||||||
|
|
||||||
|
//printf("%s\n", "rem");
|
||||||
|
bip.rem = relocate(bip.rem);
|
||||||
|
b = bip.rem;
|
||||||
|
if(b!= NULL)
|
||||||
|
//printf("%d\n", GET_SIZE(b));
|
||||||
|
|
||||||
|
|
||||||
|
scan();
|
||||||
|
if (purgeFlag) {
|
||||||
|
memset(memToPurgePtr, 0, heapSizeKiB * 1024 / 2);
|
||||||
|
}
|
||||||
|
gcCount++;
|
||||||
|
}
|
||||||
|
|
||||||
|
ObjRef alloc(unsigned int size) {
|
||||||
|
if (freeHeapPtr + size > halfHeapPtr) {
|
||||||
|
garbageCollector();
|
||||||
|
if (freeHeapPtr + size > halfHeapPtr) {
|
||||||
|
fprintf(stderr, "Heap is full\n");
|
||||||
|
exit(1);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
ObjRef or;
|
||||||
|
or = (ObjRef) freeHeapPtr;
|
||||||
|
freeHeapPtr = freeHeapPtr + size;
|
||||||
|
return or;
|
||||||
|
}
|
||||||
|
|
||||||
|
void heapAllocator(int n) {
|
||||||
|
heapSizeKiB = n;
|
||||||
|
if ((heap = malloc(n * 1024)) == NULL) {
|
||||||
|
perror("malloc");
|
||||||
|
exit(1);
|
||||||
|
}
|
||||||
|
memTargetPtr = (char *) heap;
|
||||||
|
freeHeapPtr = memTargetPtr;
|
||||||
|
halfHeapPtr = (memTargetPtr + (n * 1024) / 2);
|
||||||
|
}
|
||||||
|
|
||||||
|
/*
|
||||||
|
* Main
|
||||||
|
*/
|
||||||
|
|
||||||
int main(int argc, char *argv[]) {
|
int main(int argc, char *argv[]) {
|
||||||
// Init the Heap
|
// Init the Heap
|
||||||
initHeap(1000);
|
//initHeap(1000);
|
||||||
|
|
||||||
// Initialize the Stack
|
// Initialize the Stack
|
||||||
int size = SIZE;
|
int size = SIZE;
|
||||||
@ -372,6 +552,7 @@ int main(int argc, char *argv[]) {
|
|||||||
run = 1;
|
run = 1;
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
heapAllocator(argVheapSizeKiB);
|
||||||
// Init the sda
|
// Init the sda
|
||||||
ObjRef s[sizeSDA];
|
ObjRef s[sizeSDA];
|
||||||
sda.size = &sizeSDA;
|
sda.size = &sizeSDA;
|
||||||
|
|||||||
11
njvm.h
Normal file
11
njvm.h
Normal file
@ -0,0 +1,11 @@
|
|||||||
|
//
|
||||||
|
// Created by Elias Bennour on 28.01.24.
|
||||||
|
//
|
||||||
|
#include "objref.c"
|
||||||
|
|
||||||
|
#ifndef NINJA_NJVM_H
|
||||||
|
#define NINJA_NJVM_H
|
||||||
|
|
||||||
|
ObjRef alloc(unsigned int size);
|
||||||
|
|
||||||
|
#endif //NINJA_NJVM_H
|
||||||
4
objref.c
4
objref.c
@ -4,8 +4,8 @@
|
|||||||
|
|
||||||
typedef struct ObjRef{
|
typedef struct ObjRef{
|
||||||
// if brkokenHeart and forward_pointer added at this point bip throws errors
|
// if brkokenHeart and forward_pointer added at this point bip throws errors
|
||||||
// bool brokenHeart;
|
bool brokenHeart;
|
||||||
// struct ObjRef *forward_pointer;
|
struct ObjRef *forwardPointer;
|
||||||
unsigned int size;
|
unsigned int size;
|
||||||
unsigned char data[1];
|
unsigned char data[1];
|
||||||
} *ObjRef;
|
} *ObjRef;
|
||||||
|
|||||||
5
record.c
5
record.c
@ -5,12 +5,13 @@
|
|||||||
#define RECORD
|
#define RECORD
|
||||||
#include "stackslot.c"
|
#include "stackslot.c"
|
||||||
#include "instruktion.c"
|
#include "instruktion.c"
|
||||||
#include "GC.c"
|
#include "njvm.h"
|
||||||
|
|
||||||
ObjRef newRecord(int size){
|
ObjRef newRecord(int size){
|
||||||
ObjRef record;
|
ObjRef record;
|
||||||
unsigned int objSize;
|
unsigned int objSize;
|
||||||
objSize = sizeof(*record) + (size * sizeof(void *));
|
objSize = sizeof(*record) + (size * sizeof(void *));
|
||||||
if((record = malloc(objSize)) == NULL){
|
if((record = alloc(objSize)) == NULL){
|
||||||
perror("malloc");
|
perror("malloc");
|
||||||
}
|
}
|
||||||
record->size = MSB;
|
record->size = MSB;
|
||||||
|
|||||||
@ -1,10 +1,13 @@
|
|||||||
|
#ifndef SUPPORT_C
|
||||||
|
#define SUPPORT_C
|
||||||
|
|
||||||
#include "support.h"
|
#include "support.h"
|
||||||
#include <stdio.h>
|
#include <stdio.h>
|
||||||
#include <stdlib.h>
|
#include <stdlib.h>
|
||||||
|
|
||||||
#include "objref.c"
|
#include "objref.c"
|
||||||
#include "GC.c"
|
#include "njvm.h"
|
||||||
|
|
||||||
void fatalError(char *msg){
|
void fatalError(char *msg){
|
||||||
printf("Fatal error: %s\n", msg);
|
printf("Fatal error: %s\n", msg);
|
||||||
exit(1);
|
exit(1);
|
||||||
@ -13,7 +16,7 @@ void fatalError(char *msg){
|
|||||||
void * newPrimObject(int dataSize) {
|
void * newPrimObject(int dataSize) {
|
||||||
ObjRef bigObjRef;
|
ObjRef bigObjRef;
|
||||||
|
|
||||||
bigObjRef = malloc(sizeof(unsigned int) +
|
bigObjRef = alloc(sizeof(unsigned int) +
|
||||||
dataSize * sizeof(unsigned char));
|
dataSize * sizeof(unsigned char));
|
||||||
if (bigObjRef == NULL) {
|
if (bigObjRef == NULL) {
|
||||||
fatalError("newPrimObject() got no memory");
|
fatalError("newPrimObject() got no memory");
|
||||||
@ -26,3 +29,5 @@ void * getPrimObjectDataPointer(void * obj){
|
|||||||
ObjRef oo = ((ObjRef) (obj));
|
ObjRef oo = ((ObjRef) (obj));
|
||||||
return oo->data;
|
return oo->data;
|
||||||
}
|
}
|
||||||
|
|
||||||
|
#endif
|
||||||
Loading…
x
Reference in New Issue
Block a user