Initial commit.

1 files changed, 99 insertions, 0 deletions

diff --git a/ic-reals-6.3/base/Vector.c b/ic-reals-6.3/base/Vector.c
new file mode 100644
index 0000000..51cef26
--- /dev/null
+++ b/ic-reals-6.3/base/Vector.c
@@ -0,0 +1,99 @@
+/*
+ * Copyright (C) 2000, Imperial College
+ *
+ * This file is part of the Imperial College Exact Real Arithmetic Library.
+ * See the copyright notice included in the distribution for conditions
+ * of use.
+ */
+
+#include <stdio.h>
+#include "real.h"
+#include "real-impl.h"
+
+/*
+ * Functions for allocating and managing vector LFTs in the heap.
+ */
+
+Vec *
+allocVec()
+{
+	Vec *vec;
+
+	if ((vec = (Vec *) malloc (sizeof(Vec))) == NULL) 
+		Error(FATAL, E_INT, "allocVec", "malloc failed");
+
+
+#ifdef DAVINCI
+	newNodeId(vec);
+#else
+#ifdef TRACE
+	newNodeId(vec);
+#endif
+#endif
+
+
+	vec->tag.type = VECTOR;
+	vec->tag.dumped = FALSE;
+	vec->strm = (Real) NULL;
+
+#ifdef DAVINCI
+	beginGraphUpdate();
+	newNode(vec, VECTOR);
+	endGraphUpdate();
+#endif
+
+	return vec;
+}
+
+/*
+ * The next family of functions are one level up. They allocate structures
+ * in the heap for lfts and also initialize both the entries in the lft
+ * as well as the pointers to the heap objects to which the lft is being
+ * applied.
+ */
+
+/*
+ * Create a vector with small entries. A vector is a rational. We might
+ * wish to disallow 0/0.  
+ */
+Real
+vector_Int(int a, int b)
+{
+	Vec *vec;
+
+	if ((a == 0) && (b == 0))
+		Error(FATAL, E_INT, "vector_Int", "zero column vector");
+
+	vec = allocVec();
+	mpz_init_set_si(vec->vec[0], a);
+	mpz_init_set_si(vec->vec[1], b);
+	canonVector(vec->vec);
+
+	if (vectorSign(vec->vec) == 0)
+		vec->tag.isSigned = TRUE;
+	else
+		vec->tag.isSigned = FALSE;
+
+	return (Real) vec;
+}
+
+Real
+vector_Z(mpz_t a, mpz_t b)
+{
+	Vec *vec;
+
+	if ((mpz_sgn(a) == 0) && (mpz_sgn(b) == 0))
+		Error(FATAL, E_INT, "vector_Z", "zero column vector");
+
+	vec = allocVec();
+	mpz_init_set(vec->vec[0], a);
+	mpz_init_set(vec->vec[1], b);
+	canonVector(vec->vec);
+
+	if (vectorSign(vec->vec) == 0)
+		vec->tag.isSigned = TRUE;
+	else
+		vec->tag.isSigned = FALSE;
+
+	return (Real) vec;
+}