-funroll-loops -felide-constructors -fnonnull-objects $(INCLUDE)
ifeq ($(shell uname),SunOS)
-CXXFLAGS := -V2.7.2.3 $(CXXFLAGS)
-CXXFLAGS.g:= -V2.7.2.3 $(CXXFLAGS.g)
+#CXXFLAGS := -V2.7.2.3 $(CXXFLAGS)
+#CXXFLAGS.g:= -V2.7.2.3 $(CXXFLAGS.g)
+INCLUDE2.8 = -I/usr/local/source/libstdc++-2.8.0/libstdc++ \
+ -I/usr/local/source/libstdc++-2.8.0/libio
+CXXFLAGS := $(CXXFLAGS) -fno-rtti -fno-exceptions $(INCLUDE2.8)
+CXXFLAGS.g:= $(CXXFLAGS.g) -W -fno-rtti -fno-exceptions $(INCLUDE2.8)
endif
lib: lib.g.a lib.a
+
lib.a: ../lib/liblac.a($(o-files))
lib.g.a: ../lib/liblac.g.a($(go-files))
-dFMatrix::dFMatrix (const dFMatrix &m)
+dFMatrix::dFMatrix (const dFMatrix &m)
{
init (m.dim_image, m.dim_range);
for (int i=0; i!=dim_image*dim_range; ++i)
t01 = M.el(0,1)*M.el(2,0),
t04 = M.el(0,2)*M.el(2,0),
t07 = 1.0/(t4*M.el(2,2)-t6*M.el(2,1)-t8*M.el(2,2)+
- t00*M.el(2,1)+t01*M.el(1,2)-t04*M.el(1,1));
+ t00*M.el(2,1)+t01*M.el(1,2)-t04*M.el(1,1));
el(0,0) = (M.el(1,1)*M.el(2,2)-M.el(1,2)*M.el(2,1))*t07;
el(0,1) = -(M.el(0,1)*M.el(2,2)-M.el(0,2)*M.el(2,1))*t07;
el(0,2) = -(-M.el(0,1)*M.el(1,2)+M.el(0,2)*M.el(1,1))*t07;
{}
-dVector::dVector(int n)
+dVector::dVector(int n) :
+ dim(n),
+ maxdim(n),
+ val(0)
{
Assert (n>0, ExcInvalidNumber(n));
- dim = n;
- maxdim = n;
if (n)
{
val = new double[maxdim];
Assert (val != 0, ExcOutOfMemory());
clear ();
}
- else
- val = 0;
}
-dVector::dVector(const dVector& v)
+dVector::dVector(const dVector& v) :
+ VectorBase(v),
+ dim(v.n()),
+ maxdim(v.n()),
+ val(0)
{
- dim = v.n();
- maxdim = v.n();
-
if (dim)
{
val = new double[maxdim];
for (int i=0; i<dim; i++)
val[i] = v.val[i];
}
- else
- val = 0;
}
-iVector::iVector()
+iVector::iVector() :
+ dim(1),
+ maxdim(1),
+ val(new int[1])
{
- dim=maxdim=1;
- val = new int[1];
Assert (val != 0, ExcOutOfMemory());
clear ();
}
-iVector::iVector(int n)
+
+iVector::iVector(int n) :
+ dim(n),
+ maxdim(n),
+ val(new int[n])
{
Assert (n>0, ExcInvalidNumber(n));
- dim = n;
- maxdim = n;
- val = new int[maxdim];
Assert (val != 0, ExcOutOfMemory());
clear ();
}
-iVector::iVector(const iVector& v)
+iVector::iVector(const iVector& v) :
+ dim(v.dim),
+ maxdim(v.maxdim),
+ val(0)
{
reinit(v.dim,1);
int i;
Assert (n>0, ExcInvalidNumber(n));
if (n>maxdim)
{
- delete[] val;
+ if (val != 0) delete[] val;
val = new int[n];
Assert (val != 0, ExcOutOfMemory());
maxdim = n;