Added series expansion for functions (classical) Li and S around x==0.

[ginac.git] / check / check_matrices.cpp

diff --git a/check/check_matrices.cpp b/check/check_matrices.cpp
index 465c244ca1a00f6a0733473d5bd2e872c109e636..29f75bcad124b914567d2f0050dc6e3a7ceada99 100644 (file)
--- a/check/check_matrices.cpp
+++ b/check/check_matrices.cpp
@@ -1,9 +1,11 @@
 /** @file check_matrices.cpp
  *
 /** @file check_matrices.cpp
  *

- *  Here we test manipulations on GiNaC's symbolic matrices. */
+ *  Here we test manipulations on GiNaC's symbolic matrices.  They are a
+ *  well-tried resource for cross-checking the underlying symbolic
+ *  manipulations. */

 
 /*
 
 /*

- *  GiNaC Copyright (C) 1999-2002 Johannes Gutenberg University Mainz, Germany
+ *  GiNaC Copyright (C) 1999-2005 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
  *
  *  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


@@ -17,19 +19,19 @@
  *
  *  You should have received a copy of the GNU General Public License
  *  along with this program; if not, write to the Free Software
  *
  *  You should have received a copy of the GNU General Public License
  *  along with this program; if not, write to the Free Software

- *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
+ *  Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA

  */
 
 #include "checks.h"
 
 /* determinants of some sparse symbolic matrices with coefficients in
  * an integral domain. */
  */
 
 #include "checks.h"
 
 /* determinants of some sparse symbolic matrices with coefficients in
  * an integral domain. */

-static unsigned integdom_matrix_determinants(void)
+static unsigned integdom_matrix_determinants()

 {
        unsigned result = 0;
        symbol a("a");
        
 {
        unsigned result = 0;
        symbol a("a");
        

-       for (unsigned size=3; size<20; ++size) {
+       for (unsigned size=3; size<22; ++size) {

                matrix A(size,size);
                // populate one element in each row:
                for (unsigned r=0; r<size-1; ++r)
                matrix A(size,size);
                // populate one element in each row:
                for (unsigned r=0; r<size-1; ++r)


@@ -51,12 +53,12 @@ static unsigned integdom_matrix_determinants(void)
 
 /* determinants of some symbolic matrices with multivariate rational function
  * coefficients. */
 
 /* determinants of some symbolic matrices with multivariate rational function
  * coefficients. */

-static unsigned rational_matrix_determinants(void)
+static unsigned rational_matrix_determinants()

 {
        unsigned result = 0;
        symbol a("a"), b("b"), c("c");
        
 {
        unsigned result = 0;
        symbol a("a"), b("b"), c("c");
        

-       for (unsigned size=3; size<8; ++size) {
+       for (unsigned size=3; size<9; ++size) {

                matrix A(size,size);
                for (unsigned r=0; r<size-1; ++r) {
                        // populate one or two elements in each row:
                matrix A(size,size);
                for (unsigned r=0; r<size-1; ++r) {
                        // populate one or two elements in each row:


@@ -85,12 +87,12 @@ static unsigned rational_matrix_determinants(void)
 }
 
 /* Some quite funny determinants with functions and stuff like that inside. */
 }
 
 /* Some quite funny determinants with functions and stuff like that inside. */

-static unsigned funny_matrix_determinants(void)
+static unsigned funny_matrix_determinants()

 {
        unsigned result = 0;
        symbol a("a"), b("b"), c("c");
        
 {
        unsigned result = 0;
        symbol a("a"), b("b"), c("c");
        

-       for (unsigned size=3; size<7; ++size) {
+       for (unsigned size=3; size<8; ++size) {

                matrix A(size,size);
                for (unsigned co=0; co<size-1; ++co) {
                        // populate one or two elements in each row:
                matrix A(size,size);
                for (unsigned co=0; co<size-1; ++co) {
                        // populate one or two elements in each row:


@@ -119,12 +121,12 @@ static unsigned funny_matrix_determinants(void)
 }
 
 /* compare results from different determinant algorithms.*/
 }
 
 /* compare results from different determinant algorithms.*/

-static unsigned compare_matrix_determinants(void)
+static unsigned compare_matrix_determinants()

 {
        unsigned result = 0;
        symbol a("a");
        
 {
        unsigned result = 0;
        symbol a("a");
        

-       for (unsigned size=2; size<7; ++size) {
+       for (unsigned size=2; size<8; ++size) {

                matrix A(size,size);
                for (unsigned co=0; co<size; ++co) {
                        for (unsigned ro=0; ro<size; ++ro) {
                matrix A(size,size);
                for (unsigned co=0; co<size; ++co) {
                        for (unsigned ro=0; ro<size; ++ro) {


@@ -156,12 +158,12 @@ static unsigned compare_matrix_determinants(void)
        return result;
 }
 
        return result;
 }
 

-static unsigned symbolic_matrix_inverse(void)
+static unsigned symbolic_matrix_inverse()

 {
        unsigned result = 0;
        symbol a("a"), b("b"), c("c");
        
 {
        unsigned result = 0;
        symbol a("a"), b("b"), c("c");
        

-       for (unsigned size=2; size<5; ++size) {
+       for (unsigned size=2; size<6; ++size) {

                matrix A(size,size);
                do {
                        for (unsigned co=0; co<size; ++co) {
                matrix A(size,size);
                do {
                        for (unsigned co=0; co<size; ++co) {


@@ -193,7 +195,7 @@ static unsigned symbolic_matrix_inverse(void)
        return result;
 }
 
        return result;
 }
 

-unsigned check_matrices(void)
+unsigned check_matrices()

 {
        unsigned result = 0;
        
 {
        unsigned result = 0;