- Changes to */Makefile.in triggered by stupid automake-version-thingie.

[ginac.git] / check / matrix_checks.cpp

diff --git a/check/matrix_checks.cpp b/check/matrix_checks.cpp
index 4b973ef314cb8297f9e3a1d997fad9351057f7b1..9e907ec365127aa63a70091d15a44e9664813079 100644 (file)
--- a/check/matrix_checks.cpp
+++ b/check/matrix_checks.cpp
@@ -1,8 +1,9 @@
 /** @file matrix_checks.cpp
  *
- *  Here we test manipulations on GiNaC's symbolic matrices.
- *
- *  GiNaC Copyright (C) 1999 Johannes Gutenberg University Mainz, Germany
+ *  Here we test manipulations on GiNaC's symbolic matrices. */
+
+/*
+ *  GiNaC Copyright (C) 1999-2000 Johannes Gutenberg University Mainz, Germany
  *
  *  This program is free software; you can redistribute it and/or modify
  *  it under the terms of the GNU General Public License as published by
@@ -20,7 +21,11 @@
  */
 
 #include <stdexcept>
-#include <ginac/ginac.h>
+#include "ginac.h"
+
+#ifndef NO_NAMESPACE_GINAC
+using namespace GiNaC;
+#endif // ndef NO_NAMESPACE_GINAC
 
 static unsigned matrix_determinants(void)
 {
@@ -121,10 +126,10 @@ static unsigned matrix_invert2(void)
     matrix m_i = m.inverse();
     ex det = m.determinant().expand();
     
-    if ( (normal(m_i(0,0)*det) != d) ||
-         (normal(m_i(0,1)*det) != -b) ||
-         (normal(m_i(1,0)*det) != -c) ||
-         (normal(m_i(1,1)*det) != a) ) {
+    if ((normal(m_i(0,0)*det) != d) ||
+        (normal(m_i(0,1)*det) != -b) ||
+        (normal(m_i(1,0)*det) != -c) ||
+        (normal(m_i(1,1)*det) != a)) {
         clog << "inversion of 2x2 matrix " << m
              << " erroneously returned " << m_i << endl;
         return 1;
@@ -144,15 +149,15 @@ static unsigned matrix_invert3(void)
     matrix m_i = m.inverse();
     ex det = m.determinant().normal().expand();
     
-    if ( (normal(m_i(0,0)*det) != (e*i-f*h)) ||
-         (normal(m_i(0,1)*det) != (c*h-b*i)) ||
-         (normal(m_i(0,2)*det) != (b*f-c*e)) ||
-         (normal(m_i(1,0)*det) != (f*g-d*i)) ||
-         (normal(m_i(1,1)*det) != (a*i-c*g)) ||
-         (normal(m_i(1,2)*det) != (c*d-a*f)) ||
-         (normal(m_i(2,0)*det) != (d*h-e*g)) ||
-         (normal(m_i(2,1)*det) != (b*g-a*h)) ||
-         (normal(m_i(2,2)*det) != (a*e-b*d)) ) {
+    if ((normal(m_i(0,0)*det) != (e*i-f*h)) ||
+        (normal(m_i(0,1)*det) != (c*h-b*i)) ||
+        (normal(m_i(0,2)*det) != (b*f-c*e)) ||
+        (normal(m_i(1,0)*det) != (f*g-d*i)) ||
+        (normal(m_i(1,1)*det) != (a*i-c*g)) ||
+        (normal(m_i(1,2)*det) != (c*d-a*f)) ||
+        (normal(m_i(2,0)*det) != (d*h-e*g)) ||
+        (normal(m_i(2,1)*det) != (b*g-a*h)) ||
+        (normal(m_i(2,2)*det) != (a*e-b*d))) {
         clog << "inversion of 3x3 matrix " << m
              << " erroneously returned " << m_i << endl;
         return 1;
@@ -225,7 +230,7 @@ unsigned matrix_checks(void)
     result += matrix_invert3();
     result += matrix_misc();
     
-    if (! result) {
+    if (!result) {
         cout << " passed ";
         clog << "(no output)" << endl;
     } else {