check/lsolve_onedim.cpp

   1 /** @file lsolve_onedim.cpp
   2  *
   3  * This test routine does some simple checks on solving a polynomial for a
   4  * variable. */
   5
   6 /*
   7  *  GiNaC Copyright (C) 1999 Johannes Gutenberg University Mainz, Germany
   8  *
   9  *  This program is free software; you can redistribute it and/or modify
  10  *  it under the terms of the GNU General Public License as published by
  11  *  the Free Software Foundation; either version 2 of the License, or
  12  *  (at your option) any later version.
  13  *
  14  *  This program is distributed in the hope that it will be useful,
  15  *  but WITHOUT ANY WARRANTY; without even the implied warranty of
  16  *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  17  *  GNU General Public License for more details.
  18  *
  19  *  You should have received a copy of the GNU General Public License
  20  *  along with this program; if not, write to the Free Software
  21  *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
  22  */
  23
  24 #include <ginac/ginac.h>
  25 using namespace GiNaC;
  26
  27 unsigned lsolve_onedim(void)
  28 {
  29     unsigned result = 0;
  30     symbol x("x");
  31     ex eq, aux;
  32
  33     cout << "checking linear solve..." << flush;
  34     clog << "---------linear solve:" << endl;
  35
  36     eq = (3*x+5 == numeric(8));
  37     aux = lsolve(eq,x);
  38     if (aux != 1) {
  39         result++;
  40         clog << "solution of 3*x+5==8 erroneously returned "
  41              << aux << endl;
  42     }
  43
  44     if (! result) {
  45         cout << " passed ";
  46         clog << "(no output)" << endl;
  47     } else {
  48         cout << " failed ";
  49     }
  50
  51     return result;
  52 }