From: tiller@solace.me.uiuc.edu (Mike Tiller)
Subject: Bug
Organization: GNUs Not Usenet
Distribution: gnu
Date: Mon, 11 Jan 1993 07:51:44 GMT

You compiler tells me the following should be reported.  I was
running this with gcc-2.3.3 on a sparc-sunos4.1 and this code
used to work with gcc-2.2.2 (I think thats the correct version).

# 1 "Example1.C"
# 1 "Scalar.h" 1



 
 

# 1 "/usr/local/gnu/lib/g++-include/math.h" 1 3
 
 



















#pragma interface











# 64 "/usr/local/gnu/lib/g++-include/math.h" 3

extern "C" {

double  acos(double);
double  acosh(double);
double  asin(double);
double  asinh(double);
double  atan(double);
double  atan2(double, double);
double  atanh(double);
double  cbrt(double);
double  ceil(double);
double  copysign(double,double);
double  cos(double);
double  cosh(double);
double  drem(double,double);
double  erf(double);
double  erfc(double);
double  exp(double);
double  expm1(double);
double  fabs(double);
double  finite(double);
double  floor(double);
double	fmod(double, double);
double  frexp(double, int*);
double  gamma(double);
double  hypot(double,double);
double  infnan(int);

 
int     isinf(double);
int     isnan(double);

double  j0(double);
double  j1(double);
double  jn(int, double);
double  ldexp(double, int);
double  lgamma(double);
double  log(double);
double  log10(double);
double  log1p(double);
double  logb(double);
double  modf(double, double*);
double  pow(double, double);
double  rint(double);
double  scalb(double, int);
double  sin(double);
double  sinh(double);
double  sqrt(double);
double  tan(double);
double  tanh(double);
double  y0(double);
double  y1(double);
double  yn(int, double);
double	exp10(double);
void	sincos(double,double*,double*);

double aint(double);
double anint(double);
int irint(double);
int nint(double);
}



 

 


struct libm_exception
{
  int type;
  char* name;
  double arg1, arg2, retval;
};








extern "C" int matherr(libm_exception*);



# 1 "/usr/local/gnu/lib/gcc-lib/sparc-sun-sunos4.1/2.3.3/include/float.h" 1 3
 


 

    


    


    


    


    


    


    


    


    


    


    



    


    


    


    


    


    


    


    


    



    


    


    



union __convert_long_double {
  int __convert_long_double_i[4];
  long double __convert_long_double_d;
};


    


    


    


    


    


    




# 152 "/usr/local/gnu/lib/g++-include/math.h" 2 3


 









 







 

















































# 7 "Scalar.h" 2

# 1 "/wren7/tiller/Include/Gen/Array++.h" 1
 
 
 
 
 





# 1 "/usr/local/gnu/lib/g++-include/stdio.h" 1 3
 
 


















#pragma interface







 
 
 

# 1 "/usr/local/gnu/lib/g++-include/_G_config.h" 1 3
  






typedef long _G_clock_t;
typedef short _G_dev_t;
typedef long   _G_fpos_t;
typedef unsigned short _G_gid_t;
typedef unsigned long _G_ino_t;
typedef unsigned short _G_mode_t;
typedef short _G_nlink_t;
typedef long _G_off_t;
typedef int _G_pid_t;
typedef int _G_ptrdiff_t;
typedef int _G_sigset_t;
typedef int _G_size_t;
typedef long _G_time_t;
typedef unsigned short _G_uid_t;
typedef unsigned short _G_wchar_t;
typedef int   _G_ssize_t;
typedef char*   _G_va_list;











 






 
 

# 33 "/usr/local/gnu/lib/g++-include/stdio.h" 2 3


extern "C" {

































# 1 "/usr/local/gnu/lib/gcc-lib/sparc-sun-sunos4.1/2.3.3/include/stdio.h" 1 3
 




extern	struct	_iobuf {
	int	_cnt;
	unsigned char *_ptr;
	unsigned char *_base;
	int	_bufsiz;
	short	_flag;
	char	_file;		 
} _iob[];








































extern struct _iobuf 	*__hide_fopen ();
extern struct _iobuf 	*__hide_fdopen ();
extern struct _iobuf 	*__hide_freopen ();
extern struct _iobuf 	*__hide_popen ();
extern struct _iobuf 	*tmpfile();
extern long	__hide_ftell ();
extern char	*fgets();
extern char	*gets();
extern char	*__hide_sprintf ();
extern char	*ctermid();
extern char	*cuserid();
extern char	*__hide_tempnam ();
extern char	*tmpnam();

 
extern	int	__hide_printf (char *, ...);
extern	int	__hide_fprintf (struct _iobuf  *, char *, ...);






# 69 "/usr/local/gnu/lib/g++-include/stdio.h" 2 3







































}

extern "C" {

int    fclose(struct _iobuf *);
struct _iobuf *  fdopen(int, const char*);
int    fflush(struct _iobuf *);
int    fgetc(struct _iobuf *);
char*  fgets(char*, int, struct _iobuf  *);
struct _iobuf *  fopen(const char*, const char*);
int    fprintf(struct _iobuf *, const char* ...);
int    fputc(int, struct _iobuf *);
int    fputs(const char*, struct _iobuf *);
_G_size_t  fread(void*, _G_size_t , _G_size_t , struct _iobuf *);



struct _iobuf *  freopen(const char*, const char*, struct _iobuf *);

int    fscanf(struct _iobuf *, const char* ...);
int    fseek(struct _iobuf *, long, int);
long   ftell(struct _iobuf  *);
_G_size_t  fwrite(const void*, _G_size_t , _G_size_t , struct _iobuf *);
char*  gets(char*);
int    getw(struct _iobuf *);
int    pclose(struct _iobuf *);
void   perror(const char*);
struct _iobuf *  popen(const char*, const char*);
int    printf(const char* ...);
int    puts(const char*);
int    putw(int, struct _iobuf *);
int    rewind(struct _iobuf *);
int    scanf(const char* ...);
void   setbuf(struct _iobuf *, char*);
void   setbuffer(struct _iobuf *, char*, int);
int    setlinebuf(struct _iobuf *);
int    setvbuf(struct _iobuf *, char*, int, _G_size_t );
int    sscanf(char*, const char* ...);
struct _iobuf *  tmpfile();
int    ungetc(int, struct _iobuf *);
int    vfprintf (...) ;
int    vprintf (...) ;
char*  sprintf (...) ;
char*  vsprintf (...) ;

extern int _filbuf (...) ;
extern int _flsbuf (...) ;

}


















# 11 "/wren7/tiller/Include/Gen/Array++.h" 2



template <class T>
class	Array
{
	protected:
		T	*array;
		int	size;

	public:
		void	SetSize(int s)
		{
			T	*temp;
			int	i;

			if (array == 0 )
			{
				array = new T[s];
				for(i=0;i<s;i++)
					array[i] = T();
			}
			else
			{
				temp = array;
				array = new T[s];
				for(i=0;i<s;i++)
					if (i<size)
						array[i] = temp[i];
					else
						array[i] = T();

				delete [] temp;
			}
			size = s;
		}
		Array()
		{
			size = 0;
			array = 0 ;
		}
		Array(const T& a)
		{
			size = 0;
			array = 0 ;
			operator+=(a);
		}
		Array(const T& a, const T& b)
		{
			size = 0;
			array = 0 ;
			operator+=(a);
			operator+=(b);
		}
		Array(const T& a, const T& b, const T& c)
		{
			size = 0;
			array = 0 ;
			operator+=(a);
			operator+=(b);
			operator+=(c);
		}
		Array(const Array &a)
		{
			int	i;

			size = a.size;
			array = new T[size];
			for(i=0;i<size;i++)
				array[i] = a.array[i];
		}
		Array operator=(const Array &a)
		{
			int	i;

			size = a.size;
			array = new T[size];
			for(i=0;i<size;i++)
				array[i] = a.array[i];
			return(a);
		}
		int	OutOfRange(int i)
		{
			if (i-1 >= size || i-1 < 0)
				return(1);
			else
				return(0);
		}
		inline int Size() const { return(size); }
		inline T& Elem(int i) {
			if (i-1 >= size)
				SetSize(i);
			return(array[i-1]);
		}
		inline T GetCopy(int i) const
		{
			if (i-1 >= size)
				return(T());
			return(array[i-1]);
		}
		inline void Append(const T& a)
		{
			SetSize(size+1);
			array[size-1] = a;
		}
		inline T Remove(int n)
		{
			T	ret = T();
			if (n>size || n<1)
				return(ret);
			else
			{
				T	*temp;
				int	i, j;

				temp = array;
				array = new T[size-1];
				for(i=0,j=0;i<size;i++)
				{
					if (i==n-1)
					{
						ret = temp[i];
					}
					else
					{
						array[j] = temp[i];
						j++;
					}
				}
				delete [] temp;
				size--;
			}

			return(ret);
		}
		inline void Append(T *a)
		{
			SetSize(size+1);
			array[size-1] = *a;
		}
		inline T operator()(int i) const { return(GetCopy(i)); }
		inline T& operator[](int i) { return(Elem(i)); }
		inline T operator+=(const T& a)
		{
			SetSize(size+1);
			array[size-1] = a;
			return(a);
		}
		inline T *operator+=(T *a)
		{
			SetSize(size+1);
			array[size-1] = *a;
			return(a);
		}
		virtual ~Array()
		{
			delete [] array;
		}
};


# 8 "Scalar.h" 2

# 1 "adExpr.h" 1



 
 

# 1 "/wren7/tiller/Include/GC.h" 1


extern	"C" {
# 1 "/usr/local/gnu/lib/g++-include/stdio.h" 1 3
 
 
















# 174 "/usr/local/gnu/lib/g++-include/stdio.h" 3

# 4 "/wren7/tiller/Include/GC.h" 2

};

 
 
 
 

class	GC
{
	private:
		class	GC	*GC_next;
		class	GC	*GC_prev;
		static	class	GC	*head;
		static	class	GC	*tail;

	protected:
		int	ref_count;

	public:
		static	int		GC_size;
		GC();
		inline	void	Ref() { ref_count++; }
		inline	void	UnRef() { ref_count--; }
		inline	int		AmountOfGarbage() { return(GC_size); }
		static	void	Clean();
		static	void	Print();	 
		GC		*RemoveMe();
		~GC();
};


# 7 "adExpr.h" 2


class	adExpression : public GC
{
	private:
	public:
		adExpression() { }
		virtual	~adExpression() { }
		virtual double	Value() = 0;
		virtual	double	D(const class Scalar&) = 0;
		virtual void	Print() = 0;
};

class	adConstExpr : public adExpression
{
	private:
		double	v;
	public:
		adConstExpr(double);
		double	Value();
		double	D(const class Scalar&);
		void Print() { printf(" <([%d] Const %g)> ", ref_count, v); }
};

class	adBinExpr : public adExpression
{
	protected:
		adExpression	*e1, *e2;
	public:
		adBinExpr(Scalar&, Scalar&);
		void Print()
		{
			printf(" <([%d] Bin ", ref_count);
			e1->Print();
			e2->Print();
			printf(")> ");
		}
		~adBinExpr();
};

class	adAddExpr : public adBinExpr
{
	public:
		adAddExpr(Scalar&, Scalar&);
		double	Value();
		double	D(const class Scalar&);
};

class	adMultExpr : public adBinExpr
{
	public:
		adMultExpr(Scalar&, Scalar&);
		double	Value();
		double	D(const class Scalar&);
};

class	adSubExpr : public adBinExpr
{
	public:
		adSubExpr(Scalar&, Scalar&);
		double	Value();
		double	D(const class Scalar&);
};

class	adDivExpr : public adBinExpr
{
	public:
		adDivExpr(Scalar&, Scalar&);
		double	Value();
		double	D(const class Scalar&);
};

class	adPowExpr : public adExpression
{
	protected:
		adExpression	*b;
		double	exp;
	public:
		adPowExpr(Scalar&, double);
		double	Value();
		double	D(const class Scalar&);
		void Print()
		{
			printf(" <([%d] Pow ", ref_count);
			b->Print();
			printf(" %g)> ", exp);
		}
		~adPowExpr();
};

class	adExtExpr : public adExpression
{
	typedef	double	(*DF)(const Scalar &, const Scalar &);
	protected:
		adExpression *b;
		DF		df;
	public:
		adExtExpr(Scalar&, DF);
		double	Value() { return(b->Value()); }
		double	D(const Scalar &e);
		void Print()
		{
			printf(" <([%d] External )> ", ref_count);
		}
};

 
Scalar operator+(double, Scalar&);
Scalar operator+(Scalar&, double);
Scalar operator+(Scalar&, Scalar&);

 
Scalar operator-(double, Scalar&);
Scalar operator-(Scalar&, double);
Scalar operator-(Scalar&, Scalar&);

 
Scalar operator*(double, Scalar&);
Scalar operator*(Scalar&, double);
Scalar operator*(Scalar&, Scalar&);

 
Scalar operator/(double, Scalar&);
Scalar operator/(Scalar&, double);
Scalar operator/(Scalar&, Scalar&);

 
Scalar	pow(Scalar&, double);


# 9 "Scalar.h" 2


class	Scalar
{
	private:
		adExpression	*p;

	public:
		Scalar();
		Scalar(double);
		Scalar(adExpression *);
		Scalar(const Scalar &);
		inline	double Value() const
		{
			return(p->Value());
		}
		inline	double	D(const Scalar& e) const
		{
			return(p->D(e));
		}
		void	Print()
		{
			p->Print();
			printf("\n");
		}
		void operator=(double);
		~Scalar();
		inline	operator adExpression*() const { return(p); }
		inline	operator double() const { return(p->Value()); }

		int	operator==(Scalar &e) { return(e.p==p); }
		void operator+=(Scalar&);
		void operator-=(Scalar&);
		void operator*=(Scalar&);
		void operator/=(Scalar&);

		static	Scalar Default();
		static	Scalar Invalid();

		friend	void	operator=(Scalar &, const Scalar &);
		friend	adExpression;
		friend	adConstExpr;
		friend	adAddExpr;
		friend	adMultExpr;
};

inline	void	operator=(Scalar &e, const Scalar &s)
{
	e.p->UnRef();
	e.p = s.p;
	e.p->Ref();
}

class	ScVector : public Array<Scalar>
{
	public:
		ScVector() : Array<Scalar>() { }
		ScVector(double *v, int n)
		{
			int	i;
			
			for(i=0;i<n;i++)
				Append(Scalar(v[i]));
		}
		int	Offset(Scalar& e)
		{
			int	i;

			for(i=0;i<size;i++)
				if (array[i]==e)
					return(i+1);
			return(-1);
		}
};

# 98 "Scalar.h"



# 1 "Example1.C" 2

# 1 "/usr/local/gnu/lib/g++-include/iostream.h" 1 3
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 



#pragma interface



# 1 "/usr/local/gnu/lib/g++-include/streambuf.h" 1 3
 
 
 
 
 
 
 
 
 
 
 
 
 
 





#pragma interface


   
















class ostream; class streambuf; class backupbuf;

 







extern "C" int __underflow(streambuf*);
extern "C" int __overflow(streambuf*, int);

typedef _G_off_t streamoff;
typedef _G_off_t streampos;  

typedef unsigned long __fmtflags;
typedef unsigned char __iostate;

struct _ios_fields {  
    streambuf *_strbuf;
    ostream* _tie;
    int _width;
    __fmtflags _flags;
    _G_wchar_t _fill;
    __iostate _state;
    __iostate _exceptions;
    int _precision;
};















# 95 "/usr/local/gnu/lib/g++-include/streambuf.h" 3


class ios : public _ios_fields {
  public:
    typedef __fmtflags fmtflags;
    typedef int iostate;
    typedef int openmode;
    enum io_state {
	goodbit = 0 ,
	eofbit = 1 ,
	failbit = 2 ,
	badbit = 4  };
    enum open_mode {
	in = 1 ,
	out = 2 ,
	ate = 4 ,
	app = 8 ,
	trunc = 16 ,
	nocreate = 32 ,
	noreplace = 64 ,
	bin = 128  };
    enum seek_dir { beg, cur, end};
     
    enum { skipws=01, left=02, right=04, internal=010,
	   dec=020, oct=040, hex=0100,
	   showbase=0200, showpoint=0400, uppercase=01000, showpos=02000,
	   scientific=04000, fixed=010000, unitbuf=020000, stdio=040000,
	   dont_close=0x80000000  
	   };
    enum {  
	basefield=dec+oct+hex,
	floatfield = scientific+fixed,
	adjustfield = left+right+internal
    };

# 138 "/usr/local/gnu/lib/g++-include/streambuf.h" 3


    ostream* tie() const { return _tie; }
    ostream* tie(ostream* val) { ostream* save=_tie; _tie=val; return save; }

     
    _G_wchar_t fill() const { return (_G_wchar_t)_fill; }
    _G_wchar_t fill(_G_wchar_t newf)
	{_G_wchar_t oldf = (_G_wchar_t)_fill; _fill = (char)newf; return oldf;}
    fmtflags flags() const { return _flags; }
    fmtflags flags(fmtflags new_val) {
	fmtflags old_val = _flags; _flags = new_val; return old_val; }
    int precision() const { return _precision; }
    int precision(int newp) {
	unsigned short oldp = _precision; _precision = (unsigned short)newp;
	return oldp; }
    fmtflags setf(fmtflags val) {
	fmtflags oldbits = _flags;
	_flags |= val; return oldbits; }
    fmtflags setf(fmtflags val, fmtflags mask) {
	fmtflags oldbits = _flags;
	_flags = (_flags & ~mask) | (val & mask); return oldbits; }
    fmtflags unsetf(fmtflags mask) {
	fmtflags oldbits = _flags & mask;
	_flags &= ~mask; return oldbits; }
    int width() const { return _width; }
    int width(int val) { int save = _width; _width = val; return save; }




    void _throw_failure() { }


    streambuf* rdbuf() const { return _strbuf; }
    void clear(iostate state = 0) {
	_state = _strbuf ? state : state|badbit;
	if (_state & _exceptions) _throw_failure(); }
    void set(iostate flag) { _state |= flag;
	if (_state & _exceptions) _throw_failure(); }
    int good() const { return _state == 0; }
    int eof() const { return _state & ios::eofbit; }
    int fail() const { return _state & (ios::badbit|ios::failbit); }
    int bad() const { return _state & ios::badbit; }
    iostate rdstate() const { return _state; }
    operator void*() const { return fail() ? (void*)0 : (void*)(-1); }
    int operator!() const { return fail(); }
    iostate exception(iostate enable) {
	iostate old = _exceptions; _exceptions = enable;
	if (_state & _exceptions) _throw_failure();
	return old; }

    static int sync_with_stdio(int on);
    static void sync_with_stdio() { sync_with_stdio(1); }









  protected:
    ios(streambuf* sb = 0, ostream* tie = 0);
    virtual ~ios();
    void init(streambuf* sb) { _state=0; _strbuf=sb; }
};




typedef ios::seek_dir _seek_dir;


 
 
 
 
 




















 
 
class streammarker {
    friend class streambuf;



    friend int __underflow(streambuf*);

    struct streammarker *_next;   
    streambuf *_sbuf;  
    streampos _spos;  
    void set_streampos(streampos sp) { _spos = sp; }
    void set_offset(int offset) { _pos = offset; _spos = (streampos)(-2); }
     
     
    int _pos;
  public:
    streammarker(streambuf *sb);
    ~streammarker();
    int saving() { return  _spos == -2; }
    int delta(streammarker&);
    int delta();
};

struct __streambuf {
     
    int _flags;		 
    char* _gptr;	 
    char* _egptr;	 
    char* _eback;	 
    char* _pbase;	 
    char* _pptr;	 
    char* _epptr;	 
    char* _base;	 
    char* _ebuf;	 
    struct streambuf *_chain;

     
    friend class streammarker;
    char *_other_gbase;  
    char *_aux_limit;   
    char *_other_egptr;  
    streammarker *_markers;


     
    unsigned short _cur_column;
    char _unused;
    char _shortbuf[1];
};

extern unsigned __adjust_column(unsigned start, const char *line, int count);

struct streambuf : private __streambuf {
    friend class ios;
    friend class istream;
    friend class ostream;
    friend class streammarker;



    friend int __underflow(streambuf*);

  protected:
    static streambuf* _list_all;  
    streambuf*& xchain() { return _chain; }
    void _un_link();
    void _link_in();
    char* gptr() const { return _gptr; }
    char* pptr() const { return _pptr; }
    char* egptr() const { return _egptr; }
    char* epptr() const { return _epptr; }
    char* pbase() const { return _pbase; }
    char* eback() const { return _eback; }
    char* base() const { return _base; }
    char* ebuf() const { return _ebuf; }
    int blen() const { return _ebuf - _base; }
    void xput_char(char c) { *_pptr++ = c; }
    int xflags() { return _flags; }
    int xflags(int f) { int fl = _flags; _flags = f; return fl; }
    void xsetflags(int f) { _flags |= f; }
    void xsetflags(int f, int mask) { _flags = (_flags & ~mask) | (f & mask); }
    void gbump(int n) { _gptr += n; }
    void pbump(int n) { _pptr += n; }
    void setb(char* b, char* eb, int a=0);
    void setp(char* p, char* ep) { _pbase=_pptr=p; _epptr=ep; }
    void setg(char* eb, char* g, char *eg) { _eback=eb; _gptr=g; _egptr=eg; }
    char *shortbuf() { return _shortbuf; }

    int in_backup() { return _flags & 0x100 ; }
     
    char *Gbase() { return in_backup() ? _other_gbase : _eback; }
     
    char *eGptr() { return in_backup() ? _other_egptr : _egptr; }
     
    char *Bbase() { return in_backup() ? _eback : _other_gbase; }
    char *Bptr() { return _aux_limit; }
     
    char *eBptr() { return in_backup() ? _egptr : _other_egptr; }
    char *Nbase() { return _other_gbase; }
    char *eNptr() { return _other_egptr; }
    int have_backup() { return _other_gbase != 0 ; }
    int have_markers() { return _markers != 0 ; }
    int _least_marker();
    void switch_to_main_get_area();
    void switch_to_backup_area();
    void free_backup_area();
    void unsave_markers();  
    int put_mode() { return _flags & 0x800 ; }
    int switch_to_get_mode();
    
    streambuf(int flags=0);
  public:
    static int flush_all();
    static void flush_all_linebuffered();  
    virtual int underflow() = 0;  
    virtual int overflow(int c = (-1) ) = 0;  
    virtual int doallocate();
    virtual streampos seekoff(streamoff, _seek_dir, int mode=ios::in|ios::out);
    virtual streampos seekpos(streampos pos, int mode = ios::in|ios::out);
    int seekmark(streammarker& mark, int delta = 0);
    int sputbackc(char c);
    int sungetc();
    virtual ~streambuf();
    int unbuffered() { return _flags & 2  ? 1 : 0; }
    int linebuffered() { return _flags & 0x200  ? 1 : 0; }
    void unbuffered(int i)
	{ if (i) _flags |= 2 ; else _flags &= ~2 ; }
    void linebuffered(int i)
	{ if (i) _flags |= 0x200 ; else _flags &= ~0x200 ; }
    int allocate() {  
	if (base() || unbuffered()) return 0;
	else return doallocate(); }
     
    void allocbuf() { if (base() == 0 ) doallocbuf(); }
    void doallocbuf();
    virtual int sync();
    virtual int pbackfail(int c);
    virtual streambuf* setbuf(char* p, int len);
    int in_avail() { return _egptr - _gptr; }
    int out_waiting() { return _pptr - _pbase; }
    virtual int xsputn(const char* s, int n);
    int sputn(const char* s, int n) { return xsputn(s, n); }
    int padn(char pad, int n);  
    virtual int xsgetn(char* s, int n);
    int sgetn(char* s, int n) { return xsgetn(s, n); }
    int ignore(int);
    virtual int get_column();
    virtual int set_column(int);
    long sgetline(char* buf, _G_size_t n, char delim, int putback_delim);
    int sbumpc() {
	if (_gptr >= _egptr && __underflow(this) == (-1) ) return (-1) ;
	else return *(unsigned char*)_gptr++; }
    int sgetc() {
	if (_gptr >= _egptr && __underflow(this) == (-1) ) return (-1) ;
	else return *(unsigned char*)_gptr; }
    int snextc() {
	if (_gptr >= _egptr && __underflow(this) == (-1) ) return (-1) ;
	return _gptr++, sgetc(); }
    int sputc(int c) {
	if (_pptr >= _epptr) return __overflow(this, (unsigned char)c);
	else return *_pptr++ = c, (unsigned char)c; }
    void stossc() { if (_gptr < _egptr) _gptr++; }
    int vscan(char const *fmt0, _G_va_list ap, ios* stream = 0 );
    int scan(char const *fmt0 ...);
    int vform(char const *fmt0, _G_va_list ap);
    int form(char const *fmt0 ...);




};

 
 

 
 
 
 

class backupbuf : public streambuf {
    friend class streammarker;
  protected:
    backupbuf(int flags=0) : streambuf(flags|0x4000 ) { }
  public:
    virtual int pbackfail(int c);
    virtual int underflow();
    virtual int overflow(int c = (-1) );
};

struct __file_fields {
    short _fileno;
    int _blksize;
    _G_fpos_t  _offset;
 
};

class filebuf : public backupbuf {
  protected:
    struct __file_fields _fb;
    void init();
  public:
    static const int openprot;  
    filebuf();
    filebuf(int fd);
    filebuf(int fd, char* p, int len);
    ~filebuf();
    filebuf* attach(int fd);
    filebuf* open(const char *filename, const char *mode);
    filebuf* open(const char *filename, ios::openmode mode, int prot = 0664);
    virtual int underflow();
    virtual int overflow(int c = (-1) );
    int is_open() const { return _fb._fileno >= 0; }
    int fd() const { return is_open() ? _fb._fileno : (-1) ; }
    filebuf* close();
    virtual int doallocate();
    virtual streampos seekoff(streamoff, _seek_dir, int mode=ios::in|ios::out);
    virtual streambuf* setbuf(char* p, int len);
    int xsputn(const char* s, int n);
    int xsgetn(char* s, int n);
    virtual int sync();
  protected:  
 
    int is_reading() { return eback() != egptr(); }
    char* cur_ptr() { return is_reading() ?  gptr() : pptr(); }
     
    char* file_ptr() { return eGptr(); }
    int do_write(const char *data, int to_do);
    int do_flush() { return do_write(_pbase, _pptr-_pbase); }
     
    virtual _G_ssize_t sys_read(char* buf, _G_size_t size);
    virtual _G_fpos_t  sys_seek(_G_fpos_t , _seek_dir);
    virtual _G_ssize_t sys_write(const void*, long);
    virtual int sys_stat(void*);  
    virtual int sys_close();
};

inline ios::ios(streambuf* sb  , ostream* tie  ) {
		_state = sb ? ios::goodbit : ios::badbit; _exceptions=0;
		_strbuf=sb; _tie = tie; _width=0; _fill=' ';
		_flags=ios::skipws|ios::dec; _precision=6; }
inline ios::~ios() {
    if (!(_flags & (unsigned int)ios::dont_close)) delete _strbuf; }


# 24 "/usr/local/gnu/lib/g++-include/iostream.h" 2 3


class istream; class ostream;
typedef ios& (*__manip)(ios&);
typedef istream& (*__imanip)(istream&);
typedef ostream& (*__omanip)(ostream&);

extern istream& ws(istream& ins);
extern ostream& flush(ostream& outs);
extern ostream& endl(ostream& outs);
extern ostream& ends(ostream& outs);

class ostream : virtual public ios
{
     
    void do_osfx();
  public:
    ostream() { }
    ostream(streambuf* sb, ostream* tied=0 );
    int opfx() {
	if (!good()) return 0; else { if (_tie) _tie->flush(); return 1;} }
    void osfx() { if (flags() & (ios::unitbuf|ios::stdio))
		      do_osfx(); }
    streambuf* ostreambuf() const { return _strbuf; }
    ostream& flush();
    ostream& put(char c) { _strbuf->sputc(c); return *this; }
    ostream& put(unsigned char c) { return put((char)c); }

    ostream& write(const char *s, int n);
    ostream& write(const unsigned char *s, int n) { return write((const char*)s, n);}

    ostream& put(signed char c) { return put((char)c); }
    ostream& write(const signed char *s, int n) { return write((const char*)s, n);}

    ostream& write(const void *s, int n) { return write((const char*)s, n);}
    ostream& seekp(streampos);
    ostream& seekp(streamoff, _seek_dir);
    streampos tellp();
    ostream& form(const char *format ...);
    ostream& vform(const char *format, _G_va_list args);

    ostream& operator<<(char c);
    ostream& operator<<(unsigned char c) { return (*this) << (char)c; }

    ostream& operator<<(signed char c) { return (*this) << (char)c; }

    ostream& operator<<(const char *s);
    ostream& operator<<(const unsigned char *s)
	{ return (*this) << (const char*)s; }

    ostream& operator<<(const signed char *s)
	{ return (*this) << (const char*)s; }

    ostream& operator<<(void *p);
    ostream& operator<<(int n);
    ostream& operator<<(unsigned int n);
    ostream& operator<<(long n);
    ostream& operator<<(unsigned long n);

    ostream& operator<<(long long n);
    ostream& operator<<(unsigned long long n);

    ostream& operator<<(short n) {return operator<<((int)n);}
    ostream& operator<<(unsigned short n) {return operator<<((unsigned int)n);}
    ostream& operator<<(double n);
    ostream& operator<<(float n) { return operator<<((double)n); }
    ostream& operator<<(__omanip func) { return (*func)(*this); }
    ostream& operator<<(__manip func) {(*func)(*this); return *this;}
    ostream& operator<<(streambuf*);
};

class istream : virtual public ios
{
     
    _G_ssize_t _gcount;

    int _skip_ws();
  public:
    istream() { _gcount = 0; }
    istream(streambuf* sb, ostream*tied=0 );
    streambuf* istreambuf() const { return _strbuf; }
    istream& get(char* ptr, int len, char delim = '\n');
    istream& get(unsigned char* ptr, int len, char delim = '\n')
	{ return get((char*)ptr, len, delim); }
    istream& get(char& c);
    istream& get(unsigned char& c) { return get((char&)c); }
    istream& getline(char* ptr, int len, char delim = '\n');
    istream& getline(unsigned char* ptr, int len, char delim = '\n')
	{ return getline((char*)ptr, len, delim); }

    istream& get(signed char& c)  { return get((char&)c); }
    istream& get(signed char* ptr, int len, char delim = '\n')
	{ return get((char*)ptr, len, delim); }
    istream& getline(signed char* ptr, int len, char delim = '\n')
	{ return getline((char*)ptr, len, delim); }

    istream& read(char *ptr, int n);
    istream& read(unsigned char *ptr, int n) { return read((char*)ptr, n); }

    istream& read(signed char *ptr, int n) { return read((char*)ptr, n); }

    istream& read(void *ptr, int n) { return read((char*)ptr, n); }
     
    istream& get(streambuf& sb, char delim = '\n');
    istream& gets(char **s, char delim = '\n');
    int ipfx(int need) {
	if (!good()) { set(ios::failbit); return 0; }
	if (_tie && (need == 0 || rdbuf()->in_avail() < need)) _tie->flush();
	if (!need && (flags() & ios::skipws)) return _skip_ws();
	return 1;
    }
    int ipfx0() {  
	if (!good()) { set(ios::failbit); return 0; }
	if (_tie) _tie->flush();
	if (flags() & ios::skipws) return _skip_ws();
	return 1;
    }
    int ipfx1() {  
	if (!good()) { set(ios::failbit); return 0; }
	if (_tie && rdbuf()->in_avail() == 0) _tie->flush();
	return 1;
    }
    int get() { if (!ipfx1()) return (-1) ;
		int ch = _strbuf->sbumpc();
		if (ch == (-1) ) set(ios::eofbit);
		return ch; }
    int peek() { if (!ipfx1()) return (-1) ;
		int ch = _strbuf->sgetc();
		if (ch == (-1) ) set(ios::eofbit);
		return ch; }
    _G_ssize_t gcount() { return _gcount; }
    istream& ignore(int n=1, int delim = (-1) );
    istream& seekg(streampos);
    istream& seekg(streamoff, _seek_dir);
    streampos tellg();
    istream& putback(char ch) {
	if (good() && _strbuf->sputbackc(ch) == (-1) ) clear(ios::badbit);
	return *this;}
    istream& unget() {
	if (good() && _strbuf->sungetc() == (-1) ) clear(ios::badbit);
	return *this;}
    istream& scan(const char *format ...);
    istream& vscan(const char *format, _G_va_list args);





    istream& operator>>(char*);
    istream& operator>>(unsigned char* p) { return operator>>((char*)p); }

    istream& operator>>(signed char*p) { return operator>>((char*)p); }

    istream& operator>>(char& c);
    istream& operator>>(unsigned char& c) {return operator>>((char&)c);}

    istream& operator>>(signed char& c) {return operator>>((char&)c);}

    istream& operator>>(int&);
    istream& operator>>(long&);

    istream& operator>>(long long&);

    istream& operator>>(short&);
    istream& operator>>(unsigned int&);
    istream& operator>>(unsigned long&);

    istream& operator>>(unsigned long long&);

    istream& operator>>(unsigned short&);
    istream& operator>>(float&);
    istream& operator>>(double&);
    istream& operator>>( __manip func) {(*func)(*this); return *this;}
    istream& operator>>(__imanip func) { return (*func)(*this); }
    istream& operator>>(streambuf*);
};


class iostream : public istream, public ostream
{
    _G_ssize_t _gcount;
  public:
    iostream() { _gcount = 0; }
    iostream(streambuf* sb, ostream*tied=0 );
};

extern istream cin;
extern ostream cout, cerr, clog;  

struct Iostream_init { } ;   

inline ios& dec(ios& i)
{ i.setf(ios::dec, ios::dec|ios::hex|ios::oct); return i; }
inline ios& hex(ios& i)
{ i.setf(ios::hex, ios::dec|ios::hex|ios::oct); return i; }
inline ios& oct(ios& i)
{ i.setf(ios::oct, ios::dec|ios::hex|ios::oct); return i; }


# 2 "Example1.C" 2


static	ScVector	Temperatures;
static	ScVector	DesignVars;

static	const	double	T0 = 0.0;
static	const	double	C1 = 1.0;
static	const	double	K = 1.0;

void	GetForwardSolution(ScVector& T)
{
 
	DesignVars[1] = -40.0*K/C1;

	T[1] = T0;
	T[2] = T0;
	T[3] = T0-(C1/K)*(DesignVars[1]/2.0);
	T[4] = T0-(C1/K)*(DesignVars[1]/2.0);
	T[5] = T0-(C1/K)*DesignVars[1];
	T[6] = T0-(C1/K)*DesignVars[1];
}

double	Sensitivity(Scalar &T, Scalar &b)
{
	return(T.D(b));
}

Scalar	func(ScVector&, ScVector&);

main()
{
	ScVector	T;
	ScVector	&b = DesignVars;
	Scalar	G;
	double	v, L;
	int	i, size;

	GetForwardSolution(Temperatures);
	size = Temperatures.Size();

	for(i=1;i<=size;i++)
		T[i] = new adExtExpr(Temperatures[i], Sensitivity);

	cout << "T[3].D(b[1]) = " << T[3].D(b[1]) << endl;
	G = func(T, b);

	cout << "i	dG/db" << endl;
	for(i=1;i<=DesignVars.Size();i++)
		cout << i << "\t" << G.D(b[i]) << "\n";

	L = b[1].Value();
	v = (C1/K)*(20+(C1*L)/(2*K));
	cout << v << endl;
}

Scalar
func(ScVector& T, ScVector& b)
{
	Scalar	G;

	G += pow(T[3]-T[2]-20.0, 2.0);
	
	return(G);
}

Newsgroups: gnu.g++.bug