【刘庆源码共享】稀疏线性系统求解算法MGMRES(m) 之 基础类（Double类，封装double）

/* 
 * Copyright (c) 2009 湖南师范大学数计院 一心飞翔项目组
 * All Right Reserved
 *
 * 文件名：Double.h
 * 摘  要：定义一个与double对应的类，有助于捕捉处理 除数为0 的异常
 *
 * 作  者：刘 庆
 * 完成日期：2009年8月24日
 *
*/

#include <math.h>
#include <iostream>
#include <fstream>
#include "MyException.h"
#include "define.h"

#ifndef _DOUBLE__
#define _DOUBLE__

 

class Double
{
 double data;
public:
 Double()
 {
  data = 0.0;
 }
 Double(double d)
 {
  data = d;
 }
 ~Double(){ }

 double GetData()const
 {
  return data;
 }

 template <typename T> Double& operator = (const T t)
 {
  data = t;
  return *this;
 }
 Double& operator = (const Double& d)
 {
  data = d.data;
  return *this;
 }

 template <typename T> double operator / (const T t)
 {
  if (t == 0)
  {
   cout << "Divide By Zero Exception" << endl;
   return 0;
   throw DivideByZeroException("Divide By Zero Exception");
  }
  else
  {
   return data/t;
  }
 }
 double operator / (const Double& t)
 {
  if ( t.data==0 )
  {
   cout << "Divide By Zero Exception" << endl;
   return 0;
   throw DivideByZeroException("Divide By Zero Exception");
  }
  else
  {
   return data/t.data;
  }
 }
 template <typename T> double operator * (const T t)
 {
  return data*t;
 }
 double operator * (const Double& d)
 {
  return data*d.data;
 }
 template <typename T> double operator - (const T t)
 {
  return data - t;
 }
 double operator - (const Double& d)
 {
  return data - d.data;
 }
 template <typename T> double operator + (const T t)
 {
  return data + t;
 }
 double operator + (const Double& t)
 {
  return data + t.data;
 }
 template <typename T> Double& operator += (const T t)
 {
  data = data + t;
  return *this;
 }
 Double& operator += (const Double& t)
 {
  data = data + t.data;
  return *this;
 }
 template <typename T> Double& operator *= (const T t)
 {
  data = data*t;
  return *this;
 }
 Double& operator *= (const Double& t)
 {
  data = data * t.data;
  return *this;
 }
 template <typename T> Double& operator -= (const T t)
 {
  data = data - t;
  return *this;
 }
 Double& operator -= (const Double& t)
 {
  data = data - t.data;
  return *this;
 }
 template <typename T> Double& operator /= (const T t)
 {
  if (t == 0)
   throw DivideByZeroException("Divide By Zero Exception");
  else
  {
   data = data/t;
   return *this;
  }
 }
 Double& operator /= (const Double& t)
 {
  if(t.data==0)
   throw DivideByZeroException("Divide By Zero Exception");
  else
  {
   data = data/t.data;
   return *this;
  }
  return *this;
 }
 template <typename T> bool operator > (const T t) const
 {
  if(data > t)
   return true;
  else
   return false;
 }
 bool operator > (const Double& d)
 {
  if (data > d.data)
   return true;
  else
   return false;
 }
 template <typename T> bool operator < (const T t) const
 {
  if (data < t)
   return true;
  else
   return false;
 }
 bool operator < (const Double& d)
 {
  if (data < d.data)
   return true;
  else
   return false;
 } 
 template <typename T> bool operator >= (const T t) const
 {
  if (data >= t)
   return true;
  else
   return false;
 }
 bool operator >= (const Double& d)
 {
  if (data >= d.data)
   return true;
  else
   return false;
 }
 template <typename T> bool operator <= (const T t) const
 {
  if (data <= t)
   return true;
  else
   return false;
 }
 bool operator <= (const Double& d)
 {
  if (data<=d.data)
   return true;
  else
   return false;
 }
 template <typename T> bool operator == (const T t) const
 {
  if (t == 0)
  {
   if(fabs(data)<=_PRECIOUS)
    return true;
   else
    return false;
  }
  if (data == t)
   return true;
  else
   return false;
 }

 template <typename T> bool operator == (const T t)
 {
  if (t == 0)
  {
   if (fabs(data) <= _PRECIOUS)
    return true;
   else
    return false;
  }
  if (data == t)
   return true;
  else
   return false;
 }
 bool operator == (const Double& d)
 {
  if (data == d.data)
   return true;
  else
   return false;
 }
 template <typename T> bool operator != (const T t)const
 {
  if (t == 0)
  {
   if (fabs(data) > _PRECIOUS)
    return true;
   else
    return false;
  }
  if (data != t)
   return true;
  else
   return false;
 }
 template <typename T> bool operator != (const T t)
 {
  if (t == 0)
  {
   if (fabs(data) > _PRECIOUS)
    return true;
   else
    return false;
  }
  if (data != t)
   return true;
  else
   return false;
 }
 bool operator != (const Double& d)
 {
  if (data != d.data)
   return true;
  else
   return false;
 }
 friend ostream& operator << (ostream& os, Double& d)
 {
  os << d.data;
  return os;
 }
 friend istream& operator >> (istream& is, Double& d)
 {
  is >> d.data;
  return is;
 }
 friend ofstream& operator << (ofstream& ofs, Double& d)
 {
  ofs << d.data;
  return ofs;
 }
 friend ifstream& operator >> (ifstream& ifs, Double& d)
 {
  ifs >> d.data;
  return ifs;
 }
 bool EqualZero()
 {
  if (fabs(data) < _PRECIOUS)
   return true;
  else
   return false;
 }
 bool NotEqualZero()
 {
  return !EqualZero();
 }
};

 

#endif