check ()
{
FullMatrix<std::complex<number> > m;
- make_complex_matrix (m);
+ make_complex_square_matrix (m);
Vector<std::complex<number> > v;
- make_complex_vector (v);
+ make_complex_range_vector (v);
deallog << m.matrix_norm_square (v) << std::endl;
}
-
check ()
{
FullMatrix<std::complex<number> > m;
- make_complex_matrix (m);
+ make_complex_square_matrix (m);
deallog << m.relative_symmetry_norm2() << std::endl;
}
-
check ()
{
FullMatrix<std::complex<number> > m;
- make_complex_matrix (m);
+ make_complex_square_matrix (m);
deallog << m.trace() << std::endl;
}
-
check ()
{
FullMatrix<std::complex<number> > m;
- make_complex_matrix (m);
+ make_complex_square_matrix (m);
m.diagadd (3.141);
print_matrix (m);
}
-
check ()
{
FullMatrix<std::complex<number> > m;
- make_complex_matrix (m);
+ make_complex_square_matrix (m);
m.symmetrize ();
print_matrix (m);
}
-
check ()
{
FullMatrix<std::complex<number> > m;
- make_complex_matrix (m);
+ make_complex_square_matrix (m);
m(0,0) = m(1,1) = m(2,2) = m(3,3) = m(4,4) = 50;
m.gauss_jordan ();
print_matrix (m);
}
-
check ()
{
FullMatrix<std::complex<number> > m, n, o;
- make_complex_matrix (m);
- make_complex_matrix (n);
- make_complex_matrix (o);
+ make_complex_square_matrix (m);
+ make_complex_square_matrix (n);
+ make_complex_square_matrix (o);
m.mmult (n, o, true);
print_matrix (m);
print_matrix (n);
print_matrix (o);
}
-
check ()
{
FullMatrix<std::complex<number> > m, n, o;
- make_complex_matrix (m);
- make_complex_matrix (n);
- make_complex_matrix (o);
+ make_complex_square_matrix (m);
+ make_complex_square_matrix (n);
+ make_complex_square_matrix (o);
m.Tmmult (n, o, true);
print_matrix (m);
print_matrix (n);
print_matrix (o);
}
-
FullMatrix<std::complex<number> > m;
make_complex_matrix (m);
Vector<std::complex<number> > v, w;
- make_complex_vector (v);
- make_complex_vector (w);
+ make_complex_range_vector (v);
+ make_complex_domain_vector (w);
m.vmult (v, w, true);
print_vector (v);
print_vector (v);
print_vector (w);
}
-
check ()
{
FullMatrix<std::complex<number> > m;
- make_complex_matrix (m);
+ make_complex_square_matrix (m);
Vector<std::complex<number> > x, y;
- make_complex_vector (x);
- make_complex_vector (y);
+ make_complex_range_vector (x);
+ make_complex_range_vector (y);
m.forward (x, y);
print_vector (x);
print_vector (y);
}
-
check ()
{
FullMatrix<std::complex<number> > m;
- make_complex_matrix (m);
+ make_complex_square_matrix (m);
Vector<std::complex<number> > x, y;
- make_complex_vector (x);
- make_complex_vector (y);
+ make_complex_range_vector (x);
+ make_complex_range_vector (y);
m.backward (x, y);
print_vector (x);
print_vector (y);
}
-
check ()
{
FullMatrix<std::complex<number> > m;
- make_matrix (m);
+ make_square_matrix (m);
Vector<std::complex<number> > v;
- make_vector (v);
+ make_range_vector (v);
deallog << m.matrix_norm_square (v) << std::endl;
}
-
check ()
{
FullMatrix<std::complex<number> > m;
- make_matrix (m);
+ make_square_matrix (m);
deallog << m.relative_symmetry_norm2() << std::endl;
}
-
check ()
{
FullMatrix<std::complex<number> > m;
- make_matrix (m);
+ make_square_matrix (m);
deallog << m.trace() << std::endl;
}
-
check ()
{
FullMatrix<std::complex<number> > m;
- make_matrix (m);
+ make_square_matrix (m);
m.diagadd (3.141);
print_matrix (m);
}
-
check ()
{
FullMatrix<std::complex<number> > m;
- make_matrix (m);
+ make_square_matrix (m);
m.symmetrize ();
print_matrix (m);
}
-
check ()
{
FullMatrix<std::complex<number> > m;
- make_matrix (m);
+ make_square_matrix (m);
m(0,0) = m(1,1) = m(2,2) = m(3,3) = m(4,4) = 50;
m.gauss_jordan ();
print_matrix (m);
}
-
check ()
{
FullMatrix<std::complex<number> > m, n, o;
- make_matrix (m);
- make_matrix (n);
- make_matrix (o);
+ make_square_matrix (m);
+ make_square_matrix (n);
+ make_square_matrix (o);
m.mmult (n, o, true);
print_matrix (m);
print_matrix (n);
print_matrix (o);
}
-
check ()
{
FullMatrix<std::complex<number> > m, n, o;
- make_matrix (m);
- make_matrix (n);
- make_matrix (o);
+ make_square_matrix (m);
+ make_square_matrix (n);
+ make_square_matrix (o);
m.Tmmult (n, o, true);
print_matrix (m);
print_matrix (n);
print_matrix (o);
}
-
check ()
{
FullMatrix<std::complex<number> > m;
- make_matrix (m);
+ make_square_matrix (m);
Vector<std::complex<number> > v, w;
- make_vector (v);
- make_vector (w);
+ make_range_vector (v);
+ make_range_vector (w);
m.precondition_Jacobi (v, w, 3.141);
print_vector (v);
print_vector (w);
}
-
check ()
{
FullMatrix<number> m;
- make_matrix (m);
+ make_square_matrix (m);
Vector<number> v;
- make_vector (v);
+ make_range_vector (v);
deallog << m.matrix_norm_square (v) << std::endl;
}
-
check ()
{
FullMatrix<number> m;
- make_matrix (m);
+ make_square_matrix (m);
Vector<number> v,w;
- make_vector (v);
- make_vector (w);
+ make_range_vector (v);
+ make_range_vector (w);
for (unsigned int i=0; i<w.size(); ++i)
w(i) = w(i)+1.;
deallog << m.matrix_scalar_product (v,w) << std::endl;
}
-
check ()
{
FullMatrix<number> m;
- make_matrix (m);
+ make_square_matrix (m);
deallog << m.relative_symmetry_norm2() << std::endl;
}
-
check ()
{
FullMatrix<number> m;
- make_matrix (m);
+ make_square_matrix (m);
deallog << m.trace() << std::endl;
}
-
check ()
{
FullMatrix<number> m;
- make_matrix (m);
+ make_square_matrix (m);
m.diagadd (3.141);
print_matrix (m);
}
-
check ()
{
FullMatrix<number> m;
- make_matrix (m);
+ make_square_matrix (m);
m.symmetrize ();
print_matrix (m);
}
-
check ()
{
FullMatrix<number> m;
- make_matrix (m);
+ make_square_matrix (m);
m(0,0) = m(1,1) = m(2,2) = m(3,3) = m(4,4) = 50;
m.gauss_jordan ();
print_matrix (m);
}
-
check ()
{
FullMatrix<number> m, n, o;
- make_matrix (m);
- make_matrix (n);
- make_matrix (o);
+ make_square_matrix (m);
+ make_square_matrix (n);
+ make_square_matrix (o);
m.mmult (n, o, true);
print_matrix (m);
print_matrix (n);
print_matrix (o);
}
-
check ()
{
FullMatrix<number> m, n, o;
- make_matrix (m);
- make_matrix (n);
- make_matrix (o);
+ make_square_matrix (m);
+ make_square_matrix (n);
+ make_square_matrix (o);
m.Tmmult (n, o, true);
print_matrix (m);
print_matrix (n);
print_matrix (o);
}
-
check ()
{
FullMatrix<number> m;
- make_matrix (m);
+ make_square_matrix (m);
Vector<number> v, w;
- make_vector (v);
- make_vector (w);
+ make_range_vector (v);
+ make_range_vector (w);
m.precondition_Jacobi (v, w, 3.141);
print_vector (v);
print_vector (w);
}
-
check ()
{
FullMatrix<number> m;
- make_matrix (m);
+ make_square_matrix (m);
Vector<number> x, y;
- make_vector (x);
- make_vector (y);
+ make_range_vector (x);
+ make_range_vector (y);
m.forward (x, y);
print_vector (x);
print_vector (y);
}
-
check ()
{
FullMatrix<number> m;
- make_matrix (m);
+ make_square_matrix (m);
Vector<number> x, y;
- make_vector (x);
- make_vector (y);
+ make_range_vector (x);
+ make_range_vector (y);
m.backward (x, y);
print_vector (x);
print_vector (y);
}
-
check ()
{
FullMatrix<number> m;
- make_matrix (m);
+ make_square_matrix (m);
for (typename FullMatrix<number>::const_iterator
p = m.begin(); p!=m.end(); ++p)
<< p->value()
<< std::endl;
}
-
check ()
{
FullMatrix<number> m;
- make_matrix (m);
+ make_square_matrix (m);
for (typename FullMatrix<number>::const_iterator
<< p->value()
<< std::endl;
}
-