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
  26 #ifndef NO_GINAC_NAMESPACE
  27 using namespace GiNaC;
  28 #endif // ndef NO_GINAC_NAMESPACE
  29
  30 unsigned lsolve_onedim(void)
  31 {
  32     unsigned result = 0;
  33     symbol x("x");
  34     ex eq, aux;
  35
  36     cout << "checking linear solve..." << flush;
  37     clog << "---------linear solve:" << endl;
  38
  39     eq = (3*x+5 == numeric(8));
  40     aux = lsolve(eq,x);
  41     if (aux != 1) {
  42         result++;
  43         clog << "solution of 3*x+5==8 erroneously returned "
  44              << aux << endl;
  45     }
  46
  47     if (!result) {
  48         cout << " passed ";
  49         clog << "(no output)" << endl;
  50     } else {
  51         cout << " failed ";
  52     }
  53
  54     return result;
  55 }