DEAL::3 5 (3.0,15.)
DEAL::3 6 (4.0,18.)
DEAL::3 7 (5.0,21.)
-DEAL::3 8 (6.0,24.)
+DEAL::3 8 (0.0,0.0)
DEAL::3 9 (0.0,0.0)
DEAL::4 0 (0.0,0.0)
DEAL::4 1 (0.0,0.0)
DEAL::4 5 (6.0,20.)
DEAL::4 6 (8.0,24.)
DEAL::4 7 (10.,28.)
-DEAL::4 8 (12.,32.)
+DEAL::4 8 (0.0,0.0)
DEAL::4 9 (0.0,0.0)
DEAL::5 0 (0.0,0.0)
DEAL::5 1 (0.0,0.0)
DEAL::5 5 (9.0,25.)
DEAL::5 6 (12.,30.)
DEAL::5 7 (15.,35.)
-DEAL::5 8 (18.,40.)
+DEAL::5 8 (0.0,0.0)
DEAL::5 9 (0.0,0.0)
DEAL::6 0 (0.0,0.0)
DEAL::6 1 (0.0,0.0)
DEAL::6 5 (12.,30.)
DEAL::6 6 (16.,36.)
DEAL::6 7 (20.,42.)
-DEAL::6 8 (24.,48.)
+DEAL::6 8 (0.0,0.0)
DEAL::6 9 (0.0,0.0)
DEAL::7 0 (0.0,0.0)
DEAL::7 1 (0.0,0.0)
DEAL::7 5 (15.,35.)
DEAL::7 6 (20.,42.)
DEAL::7 7 (25.,49.)
-DEAL::7 8 (30.,56.)
+DEAL::7 8 (0.0,0.0)
DEAL::7 9 (0.0,0.0)
DEAL::8 0 (0.0,0.0)
DEAL::8 1 (0.0,0.0)
DEAL::3 5 (3.0,15.)
DEAL::3 6 (4.0,18.)
DEAL::3 7 (5.0,21.)
-DEAL::3 8 (6.0,24.)
+DEAL::3 8 (0.0,0.0)
DEAL::3 9 (0.0,0.0)
DEAL::4 0 (0.0,0.0)
DEAL::4 1 (0.0,0.0)
DEAL::4 5 (6.0,20.)
DEAL::4 6 (8.0,24.)
DEAL::4 7 (10.,28.)
-DEAL::4 8 (12.,32.)
+DEAL::4 8 (0.0,0.0)
DEAL::4 9 (0.0,0.0)
DEAL::5 0 (0.0,0.0)
DEAL::5 1 (0.0,0.0)
DEAL::5 5 (9.0,25.)
DEAL::5 6 (12.,30.)
DEAL::5 7 (15.,35.)
-DEAL::5 8 (18.,40.)
+DEAL::5 8 (0.0,0.0)
DEAL::5 9 (0.0,0.0)
DEAL::6 0 (0.0,0.0)
DEAL::6 1 (0.0,0.0)
DEAL::6 5 (12.,30.)
DEAL::6 6 (16.,36.)
DEAL::6 7 (20.,42.)
-DEAL::6 8 (24.,48.)
+DEAL::6 8 (0.0,0.0)
DEAL::6 9 (0.0,0.0)
DEAL::7 0 (0.0,0.0)
DEAL::7 1 (0.0,0.0)
DEAL::7 5 (15.,35.)
DEAL::7 6 (20.,42.)
DEAL::7 7 (25.,49.)
-DEAL::7 8 (30.,56.)
+DEAL::7 8 (0.0,0.0)
DEAL::7 9 (0.0,0.0)
DEAL::8 0 (0.0,0.0)
DEAL::8 1 (0.0,0.0)
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);
for (unsigned int i=0; i<w.size(); ++i)
w(i) = w(i)+std::complex<number>(1.);
deallog << m.matrix_scalar_product (v,w) << std::endl;
}
-
-DEAL::(-3.6e+04,3.2e+03)
-DEAL::(-3.6e+04,3.2e+03)
+DEAL::(-2.3e+04,1.3e+03)
+DEAL::(-2.3e+04,1.3e+03)
-DEAL::2.2e+02
-DEAL::2.2e+02
+DEAL::1.9e+02
+DEAL::1.9e+02
-DEAL::2.3e+02
-DEAL::2.3e+02
+DEAL::1.7e+02
+DEAL::1.7e+02
-DEAL::1.7e+02
-DEAL::1.7e+02
+DEAL::1.4e+02
+DEAL::1.4e+02
-(2.0,12.)(3.0,15.)(4.0,18.)(5.0,21.)(6.0,24.)
-(4.0,16.)(6.0,20.)(8.0,24.)(10.,28.)(12.,32.)
-(6.0,20.)(9.0,25.)(12.,30.)(15.,35.)(18.,40.)
-(8.0,24.)(12.,30.)(16.,36.)(20.,42.)(24.,48.)
-(10.,28.)(15.,35.)(20.,42.)(25.,49.)(30.,56.)
-(2.0,12.)(3.0,15.)(4.0,18.)(5.0,21.)(6.0,24.)
-(4.0,16.)(6.0,20.)(8.0,24.)(10.,28.)(12.,32.)
-(6.0,20.)(9.0,25.)(12.,30.)(15.,35.)(18.,40.)
-(8.0,24.)(12.,30.)(16.,36.)(20.,42.)(24.,48.)
-(10.,28.)(15.,35.)(20.,42.)(25.,49.)(30.,56.)
+(2.0,12.)(3.0,15.)(4.0,18.)(5.0,21.)
+(4.0,16.)(6.0,20.)(8.0,24.)(10.,28.)
+(6.0,20.)(9.0,25.)(12.,30.)(15.,35.)
+(8.0,24.)(12.,30.)(16.,36.)(20.,42.)
+(10.,28.)(15.,35.)(20.,42.)(25.,49.)
+(2.0,12.)(3.0,15.)(4.0,18.)(5.0,21.)
+(4.0,16.)(6.0,20.)(8.0,24.)(10.,28.)
+(6.0,20.)(9.0,25.)(12.,30.)(15.,35.)
+(8.0,24.)(12.,30.)(16.,36.)(20.,42.)
+(10.,28.)(15.,35.)(20.,42.)(25.,49.)
DEAL::3 2 (0.0,0.0)
DEAL::3 3 (0.0,0.0)
DEAL::3 4 (30.,70.)
-DEAL::3 5 (36.,80.)
+DEAL::3 5 (0.0,0.0)
DEAL::3 6 (0.0,0.0)
DEAL::3 7 (0.0,0.0)
DEAL::3 8 (0.0,0.0)
DEAL::4 2 (0.0,0.0)
DEAL::4 3 (0.0,0.0)
DEAL::4 4 (40.,84.)
-DEAL::4 5 (48.,96.)
+DEAL::4 5 (0.0,0.0)
DEAL::4 6 (0.0,0.0)
DEAL::4 7 (0.0,0.0)
DEAL::4 8 (0.0,0.0)
DEAL::5 2 (0.0,0.0)
DEAL::5 3 (0.0,0.0)
DEAL::5 4 (50.,98.)
-DEAL::5 5 (60.,1.1e+02)
+DEAL::5 5 (0.0,0.0)
DEAL::5 6 (0.0,0.0)
DEAL::5 7 (0.0,0.0)
DEAL::5 8 (0.0,0.0)
DEAL::3 2 (0.0,0.0)
DEAL::3 3 (0.0,0.0)
DEAL::3 4 (30.,70.)
-DEAL::3 5 (36.,80.)
+DEAL::3 5 (0.0,0.0)
DEAL::3 6 (0.0,0.0)
DEAL::3 7 (0.0,0.0)
DEAL::3 8 (0.0,0.0)
DEAL::4 2 (0.0,0.0)
DEAL::4 3 (0.0,0.0)
DEAL::4 4 (40.,84.)
-DEAL::4 5 (48.,96.)
+DEAL::4 5 (0.0,0.0)
DEAL::4 6 (0.0,0.0)
DEAL::4 7 (0.0,0.0)
DEAL::4 8 (0.0,0.0)
DEAL::5 2 (0.0,0.0)
DEAL::5 3 (0.0,0.0)
DEAL::5 4 (50.,98.)
-DEAL::5 5 (60.,1.1e+02)
+DEAL::5 5 (0.0,0.0)
DEAL::5 6 (0.0,0.0)
DEAL::5 7 (0.0,0.0)
DEAL::5 8 (0.0,0.0)
DEAL::4 1 (0.0,0.0)
DEAL::4 2 (0.0,0.0)
DEAL::4 3 (0.0,0.0)
-DEAL::4 4 (36.,80.)
-DEAL::4 5 (48.,96.)
-DEAL::4 6 (60.,1.1e+02)
+DEAL::4 4 (0.0,0.0)
+DEAL::4 5 (0.0,0.0)
+DEAL::4 6 (0.0,0.0)
DEAL::4 7 (0.0,0.0)
DEAL::4 8 (0.0,0.0)
DEAL::4 9 (0.0,0.0)
DEAL::4 1 (0.0,0.0)
DEAL::4 2 (0.0,0.0)
DEAL::4 3 (0.0,0.0)
-DEAL::4 4 (36.,80.)
-DEAL::4 5 (48.,96.)
-DEAL::4 6 (60.,1.1e+02)
+DEAL::4 4 (0.0,0.0)
+DEAL::4 5 (0.0,0.0)
+DEAL::4 6 (0.0,0.0)
DEAL::4 7 (0.0,0.0)
DEAL::4 8 (0.0,0.0)
DEAL::4 9 (0.0,0.0)
{
FullMatrix<std::complex<number> > m;
make_complex_matrix (m);
- m.add_col (2, 3.141, 4);
+ m.add_col (2, 3.141, 3);
print_matrix (m);
}
-
DEAL::0 0 (2.0,12.)
DEAL::0 1 (3.0,15.)
-DEAL::0 2 (23.,93.)
+DEAL::0 2 (20.,84.)
DEAL::0 3 (5.0,21.)
-DEAL::0 4 (6.0,24.)
DEAL::1 0 (4.0,16.)
DEAL::1 1 (6.0,20.)
-DEAL::1 2 (46.,1.2e+02)
+DEAL::1 2 (39.,1.1e+02)
DEAL::1 3 (10.,28.)
-DEAL::1 4 (12.,32.)
DEAL::2 0 (6.0,20.)
DEAL::2 1 (9.0,25.)
-DEAL::2 2 (69.,1.6e+02)
+DEAL::2 2 (59.,1.4e+02)
DEAL::2 3 (15.,35.)
-DEAL::2 4 (18.,40.)
DEAL::3 0 (8.0,24.)
DEAL::3 1 (12.,30.)
-DEAL::3 2 (91.,1.9e+02)
+DEAL::3 2 (79.,1.7e+02)
DEAL::3 3 (20.,42.)
-DEAL::3 4 (24.,48.)
DEAL::4 0 (10.,28.)
DEAL::4 1 (15.,35.)
-DEAL::4 2 (1.1e+02,2.2e+02)
+DEAL::4 2 (99.,2.0e+02)
DEAL::4 3 (25.,49.)
-DEAL::4 4 (30.,56.)
DEAL::0 0 (2.0,12.)
DEAL::0 1 (3.0,15.)
-DEAL::0 2 (23.,93.)
+DEAL::0 2 (20.,84.)
DEAL::0 3 (5.0,21.)
-DEAL::0 4 (6.0,24.)
DEAL::1 0 (4.0,16.)
DEAL::1 1 (6.0,20.)
-DEAL::1 2 (46.,1.2e+02)
+DEAL::1 2 (39.,1.1e+02)
DEAL::1 3 (10.,28.)
-DEAL::1 4 (12.,32.)
DEAL::2 0 (6.0,20.)
DEAL::2 1 (9.0,25.)
-DEAL::2 2 (69.,1.6e+02)
+DEAL::2 2 (59.,1.4e+02)
DEAL::2 3 (15.,35.)
-DEAL::2 4 (18.,40.)
DEAL::3 0 (8.0,24.)
DEAL::3 1 (12.,30.)
-DEAL::3 2 (91.,1.9e+02)
+DEAL::3 2 (79.,1.7e+02)
DEAL::3 3 (20.,42.)
-DEAL::3 4 (24.,48.)
DEAL::4 0 (10.,28.)
DEAL::4 1 (15.,35.)
-DEAL::4 2 (1.1e+02,2.2e+02)
+DEAL::4 2 (99.,2.0e+02)
DEAL::4 3 (25.,49.)
-DEAL::4 4 (30.,56.)
{
FullMatrix<std::complex<number> > m;
make_complex_matrix (m);
- m.add_col (2, 3.141, 4, 2.718, 3);
+ m.add_col (2, 3.141, 1, 2.718, 3);
print_matrix (m);
}
-
DEAL::0 0 (2.0,12.)
DEAL::0 1 (3.0,15.)
-DEAL::0 2 (36.,1.5e+02)
+DEAL::0 2 (27.,1.2e+02)
DEAL::0 3 (5.0,21.)
-DEAL::0 4 (6.0,24.)
DEAL::1 0 (4.0,16.)
DEAL::1 1 (6.0,20.)
-DEAL::1 2 (73.,2.0e+02)
+DEAL::1 2 (54.,1.6e+02)
DEAL::1 3 (10.,28.)
-DEAL::1 4 (12.,32.)
DEAL::2 0 (6.0,20.)
DEAL::2 1 (9.0,25.)
-DEAL::2 2 (1.1e+02,2.5e+02)
+DEAL::2 2 (81.,2.0e+02)
DEAL::2 3 (15.,35.)
-DEAL::2 4 (18.,40.)
DEAL::3 0 (8.0,24.)
DEAL::3 1 (12.,30.)
-DEAL::3 2 (1.5e+02,3.0e+02)
+DEAL::3 2 (1.1e+02,2.4e+02)
DEAL::3 3 (20.,42.)
-DEAL::3 4 (24.,48.)
DEAL::4 0 (10.,28.)
DEAL::4 1 (15.,35.)
-DEAL::4 2 (1.8e+02,3.5e+02)
+DEAL::4 2 (1.4e+02,2.9e+02)
DEAL::4 3 (25.,49.)
-DEAL::4 4 (30.,56.)
DEAL::0 0 (2.0,12.)
DEAL::0 1 (3.0,15.)
-DEAL::0 2 (36.,1.5e+02)
+DEAL::0 2 (27.,1.2e+02)
DEAL::0 3 (5.0,21.)
-DEAL::0 4 (6.0,24.)
DEAL::1 0 (4.0,16.)
DEAL::1 1 (6.0,20.)
-DEAL::1 2 (73.,2.0e+02)
+DEAL::1 2 (54.,1.6e+02)
DEAL::1 3 (10.,28.)
-DEAL::1 4 (12.,32.)
DEAL::2 0 (6.0,20.)
DEAL::2 1 (9.0,25.)
-DEAL::2 2 (1.1e+02,2.5e+02)
+DEAL::2 2 (81.,2.0e+02)
DEAL::2 3 (15.,35.)
-DEAL::2 4 (18.,40.)
DEAL::3 0 (8.0,24.)
DEAL::3 1 (12.,30.)
-DEAL::3 2 (1.5e+02,3.0e+02)
+DEAL::3 2 (1.1e+02,2.4e+02)
DEAL::3 3 (20.,42.)
-DEAL::3 4 (24.,48.)
DEAL::4 0 (10.,28.)
DEAL::4 1 (15.,35.)
-DEAL::4 2 (1.8e+02,3.5e+02)
+DEAL::4 2 (1.4e+02,2.9e+02)
DEAL::4 3 (25.,49.)
-DEAL::4 4 (30.,56.)
{
FullMatrix<std::complex<number> > m;
make_complex_matrix (m);
- m.swap_col (2, 4);
+ m.swap_col (2, 3);
print_matrix (m);
}
-
DEAL::0 0 (2.0,12.)
DEAL::0 1 (3.0,15.)
-DEAL::0 2 (6.0,24.)
-DEAL::0 3 (5.0,21.)
-DEAL::0 4 (4.0,18.)
+DEAL::0 2 (5.0,21.)
+DEAL::0 3 (4.0,18.)
DEAL::1 0 (4.0,16.)
DEAL::1 1 (6.0,20.)
-DEAL::1 2 (12.,32.)
-DEAL::1 3 (10.,28.)
-DEAL::1 4 (8.0,24.)
+DEAL::1 2 (10.,28.)
+DEAL::1 3 (8.0,24.)
DEAL::2 0 (6.0,20.)
DEAL::2 1 (9.0,25.)
-DEAL::2 2 (18.,40.)
-DEAL::2 3 (15.,35.)
-DEAL::2 4 (12.,30.)
+DEAL::2 2 (15.,35.)
+DEAL::2 3 (12.,30.)
DEAL::3 0 (8.0,24.)
DEAL::3 1 (12.,30.)
-DEAL::3 2 (24.,48.)
-DEAL::3 3 (20.,42.)
-DEAL::3 4 (16.,36.)
+DEAL::3 2 (20.,42.)
+DEAL::3 3 (16.,36.)
DEAL::4 0 (10.,28.)
DEAL::4 1 (15.,35.)
-DEAL::4 2 (30.,56.)
-DEAL::4 3 (25.,49.)
-DEAL::4 4 (20.,42.)
+DEAL::4 2 (25.,49.)
+DEAL::4 3 (20.,42.)
DEAL::0 0 (2.0,12.)
DEAL::0 1 (3.0,15.)
-DEAL::0 2 (6.0,24.)
-DEAL::0 3 (5.0,21.)
-DEAL::0 4 (4.0,18.)
+DEAL::0 2 (5.0,21.)
+DEAL::0 3 (4.0,18.)
DEAL::1 0 (4.0,16.)
DEAL::1 1 (6.0,20.)
-DEAL::1 2 (12.,32.)
-DEAL::1 3 (10.,28.)
-DEAL::1 4 (8.0,24.)
+DEAL::1 2 (10.,28.)
+DEAL::1 3 (8.0,24.)
DEAL::2 0 (6.0,20.)
DEAL::2 1 (9.0,25.)
-DEAL::2 2 (18.,40.)
-DEAL::2 3 (15.,35.)
-DEAL::2 4 (12.,30.)
+DEAL::2 2 (15.,35.)
+DEAL::2 3 (12.,30.)
DEAL::3 0 (8.0,24.)
DEAL::3 1 (12.,30.)
-DEAL::3 2 (24.,48.)
-DEAL::3 3 (20.,42.)
-DEAL::3 4 (16.,36.)
+DEAL::3 2 (20.,42.)
+DEAL::3 3 (16.,36.)
DEAL::4 0 (10.,28.)
DEAL::4 1 (15.,35.)
-DEAL::4 2 (30.,56.)
-DEAL::4 3 (25.,49.)
-DEAL::4 4 (20.,42.)
+DEAL::4 2 (25.,49.)
+DEAL::4 3 (20.,42.)
-DEAL::0 (-4.1e+02,4.1e+02)
-DEAL::1 (-5.0e+02,6.2e+02)
-DEAL::2 (-5.9e+02,8.4e+02)
-DEAL::3 (-6.8e+02,1.0e+03)
-DEAL::4 (-7.6e+02,1.3e+03)
+DEAL::0 (-2.7e+02,2.5e+02)
+DEAL::1 (-3.3e+02,3.8e+02)
+DEAL::2 (-4.0e+02,5.1e+02)
+DEAL::3 (-4.6e+02,6.4e+02)
+DEAL::4 (-5.2e+02,7.7e+02)
DEAL::0 (1.0,3.0)
DEAL::1 (2.0,4.0)
DEAL::2 (3.0,5.0)
DEAL::3 (4.0,6.0)
-DEAL::4 (5.0,7.0)
-DEAL::0 (-4.1e+02,4.1e+02)
-DEAL::1 (-5.0e+02,6.2e+02)
-DEAL::2 (-5.9e+02,8.3e+02)
-DEAL::3 (-6.8e+02,1.0e+03)
-DEAL::4 (-7.7e+02,1.2e+03)
+DEAL::0 (-2.7e+02,2.5e+02)
+DEAL::1 (-3.4e+02,3.8e+02)
+DEAL::2 (-4.0e+02,5.0e+02)
+DEAL::3 (-4.6e+02,6.3e+02)
+DEAL::4 (-5.3e+02,7.6e+02)
DEAL::0 (1.0,3.0)
DEAL::1 (2.0,4.0)
DEAL::2 (3.0,5.0)
DEAL::3 (4.0,6.0)
-DEAL::4 (5.0,7.0)
-DEAL::0 (-4.1e+02,4.1e+02)
-DEAL::1 (-5.0e+02,6.2e+02)
-DEAL::2 (-5.9e+02,8.4e+02)
-DEAL::3 (-6.8e+02,1.0e+03)
-DEAL::4 (-7.6e+02,1.3e+03)
+DEAL::0 (-2.7e+02,2.5e+02)
+DEAL::1 (-3.3e+02,3.8e+02)
+DEAL::2 (-4.0e+02,5.1e+02)
+DEAL::3 (-4.6e+02,6.4e+02)
+DEAL::4 (-5.2e+02,7.7e+02)
DEAL::0 (1.0,3.0)
DEAL::1 (2.0,4.0)
DEAL::2 (3.0,5.0)
DEAL::3 (4.0,6.0)
-DEAL::4 (5.0,7.0)
-DEAL::0 (-4.1e+02,4.1e+02)
-DEAL::1 (-5.0e+02,6.2e+02)
-DEAL::2 (-5.9e+02,8.3e+02)
-DEAL::3 (-6.8e+02,1.0e+03)
-DEAL::4 (-7.7e+02,1.2e+03)
+DEAL::0 (-2.7e+02,2.5e+02)
+DEAL::1 (-3.4e+02,3.8e+02)
+DEAL::2 (-4.0e+02,5.0e+02)
+DEAL::3 (-4.6e+02,6.3e+02)
+DEAL::4 (-5.3e+02,7.6e+02)
DEAL::0 (1.0,3.0)
DEAL::1 (2.0,4.0)
DEAL::2 (3.0,5.0)
DEAL::3 (4.0,6.0)
-DEAL::4 (5.0,7.0)
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_add (v, w);
print_vector (v);
print_vector (w);
}
-
-DEAL::0 (-4.1e+02,4.1e+02)
-DEAL::1 (-5.0e+02,6.2e+02)
-DEAL::2 (-5.9e+02,8.4e+02)
-DEAL::3 (-6.8e+02,1.0e+03)
-DEAL::4 (-7.6e+02,1.3e+03)
+DEAL::0 (-2.7e+02,2.5e+02)
+DEAL::1 (-3.3e+02,3.8e+02)
+DEAL::2 (-4.0e+02,5.1e+02)
+DEAL::3 (-4.6e+02,6.4e+02)
+DEAL::4 (-5.2e+02,7.7e+02)
DEAL::0 (1.0,3.0)
DEAL::1 (2.0,4.0)
DEAL::2 (3.0,5.0)
DEAL::3 (4.0,6.0)
-DEAL::4 (5.0,7.0)
-DEAL::0 (-4.1e+02,4.1e+02)
-DEAL::1 (-5.0e+02,6.2e+02)
-DEAL::2 (-5.9e+02,8.4e+02)
-DEAL::3 (-6.8e+02,1.0e+03)
-DEAL::4 (-7.6e+02,1.3e+03)
+DEAL::0 (-2.7e+02,2.5e+02)
+DEAL::1 (-3.3e+02,3.8e+02)
+DEAL::2 (-4.0e+02,5.1e+02)
+DEAL::3 (-4.6e+02,6.4e+02)
+DEAL::4 (-5.2e+02,7.7e+02)
DEAL::0 (1.0,3.0)
DEAL::1 (2.0,4.0)
DEAL::2 (3.0,5.0)
DEAL::3 (4.0,6.0)
-DEAL::4 (5.0,7.0)
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_domain_vector (v);
+ make_complex_range_vector (w);
m.Tvmult (v, w, true);
print_vector (v);
print_vector (v);
print_vector (w);
}
-
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_domain_vector (v);
+ make_complex_range_vector (w);
m.Tvmult_add (v, w);
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> > v, w;
- make_complex_vector (v);
- make_complex_vector (w);
+ make_complex_range_vector (v);
+ make_complex_range_vector (w);
m.precondition_Jacobi (v, w, 3.141);
print_vector (v);
print_vector (w);
}
-
FullMatrix<std::complex<number> > m;
make_complex_matrix (m);
Vector<std::complex<number> > v, w, x;
- make_complex_vector (v);
- make_complex_vector (w);
- make_complex_vector (x);
+ make_complex_range_vector (v);
+ make_complex_domain_vector (w);
+ make_complex_range_vector (x);
m.residual (v, w, x);
print_vector (v);
print_vector (w);
print_vector (x);
}
-
-DEAL::0 (4.1e+02,-4.1e+02)
-DEAL::1 (5.0e+02,-6.2e+02)
-DEAL::2 (5.9e+02,-8.2e+02)
-DEAL::3 (6.8e+02,-1.0e+03)
-DEAL::4 (7.8e+02,-1.2e+03)
+DEAL::0 (2.7e+02,-2.4e+02)
+DEAL::1 (3.4e+02,-3.7e+02)
+DEAL::2 (4.0e+02,-5.0e+02)
+DEAL::3 (4.7e+02,-6.3e+02)
+DEAL::4 (5.3e+02,-7.5e+02)
DEAL::0 (1.0,3.0)
DEAL::1 (2.0,4.0)
DEAL::2 (3.0,5.0)
DEAL::3 (4.0,6.0)
-DEAL::4 (5.0,7.0)
DEAL::0 (1.0,3.0)
DEAL::1 (2.0,4.0)
DEAL::2 (3.0,5.0)
DEAL::3 (4.0,6.0)
DEAL::4 (5.0,7.0)
-DEAL::0 (4.1e+02,-4.1e+02)
-DEAL::1 (5.0e+02,-6.2e+02)
-DEAL::2 (5.9e+02,-8.2e+02)
-DEAL::3 (6.8e+02,-1.0e+03)
-DEAL::4 (7.8e+02,-1.2e+03)
+DEAL::0 (2.7e+02,-2.4e+02)
+DEAL::1 (3.4e+02,-3.7e+02)
+DEAL::2 (4.0e+02,-5.0e+02)
+DEAL::3 (4.7e+02,-6.3e+02)
+DEAL::4 (5.3e+02,-7.5e+02)
DEAL::0 (1.0,3.0)
DEAL::1 (2.0,4.0)
DEAL::2 (3.0,5.0)
DEAL::3 (4.0,6.0)
-DEAL::4 (5.0,7.0)
DEAL::0 (1.0,3.0)
DEAL::1 (2.0,4.0)
DEAL::2 (3.0,5.0)
DEAL::3 5 (3.0,0.0)
DEAL::3 6 (4.0,0.0)
DEAL::3 7 (5.0,0.0)
-DEAL::3 8 (6.0,0.0)
+DEAL::3 8 (0.0,0.0)
DEAL::3 9 (0.0,0.0)
DEAL::4 0 (0.0,0.0)
DEAL::4 1 (0.0,0.0)
DEAL::4 5 (6.0,0.0)
DEAL::4 6 (8.0,0.0)
DEAL::4 7 (10.,0.0)
-DEAL::4 8 (12.,0.0)
+DEAL::4 8 (0.0,0.0)
DEAL::4 9 (0.0,0.0)
DEAL::5 0 (0.0,0.0)
DEAL::5 1 (0.0,0.0)
DEAL::5 5 (9.0,0.0)
DEAL::5 6 (12.,0.0)
DEAL::5 7 (15.,0.0)
-DEAL::5 8 (18.,0.0)
+DEAL::5 8 (0.0,0.0)
DEAL::5 9 (0.0,0.0)
DEAL::6 0 (0.0,0.0)
DEAL::6 1 (0.0,0.0)
DEAL::6 5 (12.,0.0)
DEAL::6 6 (16.,0.0)
DEAL::6 7 (20.,0.0)
-DEAL::6 8 (24.,0.0)
+DEAL::6 8 (0.0,0.0)
DEAL::6 9 (0.0,0.0)
DEAL::7 0 (0.0,0.0)
DEAL::7 1 (0.0,0.0)
DEAL::7 5 (15.,0.0)
DEAL::7 6 (20.,0.0)
DEAL::7 7 (25.,0.0)
-DEAL::7 8 (30.,0.0)
+DEAL::7 8 (0.0,0.0)
DEAL::7 9 (0.0,0.0)
DEAL::8 0 (0.0,0.0)
DEAL::8 1 (0.0,0.0)
DEAL::3 5 (3.0,0.0)
DEAL::3 6 (4.0,0.0)
DEAL::3 7 (5.0,0.0)
-DEAL::3 8 (6.0,0.0)
+DEAL::3 8 (0.0,0.0)
DEAL::3 9 (0.0,0.0)
DEAL::4 0 (0.0,0.0)
DEAL::4 1 (0.0,0.0)
DEAL::4 5 (6.0,0.0)
DEAL::4 6 (8.0,0.0)
DEAL::4 7 (10.,0.0)
-DEAL::4 8 (12.,0.0)
+DEAL::4 8 (0.0,0.0)
DEAL::4 9 (0.0,0.0)
DEAL::5 0 (0.0,0.0)
DEAL::5 1 (0.0,0.0)
DEAL::5 5 (9.0,0.0)
DEAL::5 6 (12.,0.0)
DEAL::5 7 (15.,0.0)
-DEAL::5 8 (18.,0.0)
+DEAL::5 8 (0.0,0.0)
DEAL::5 9 (0.0,0.0)
DEAL::6 0 (0.0,0.0)
DEAL::6 1 (0.0,0.0)
DEAL::6 5 (12.,0.0)
DEAL::6 6 (16.,0.0)
DEAL::6 7 (20.,0.0)
-DEAL::6 8 (24.,0.0)
+DEAL::6 8 (0.0,0.0)
DEAL::6 9 (0.0,0.0)
DEAL::7 0 (0.0,0.0)
DEAL::7 1 (0.0,0.0)
DEAL::7 5 (15.,0.0)
DEAL::7 6 (20.,0.0)
DEAL::7 7 (25.,0.0)
-DEAL::7 8 (30.,0.0)
+DEAL::7 8 (0.0,0.0)
DEAL::7 9 (0.0,0.0)
DEAL::8 0 (0.0,0.0)
DEAL::8 1 (0.0,0.0)
FullMatrix<std::complex<number> > m;
make_matrix (m);
Vector<std::complex<number> > v,w;
- make_vector (v);
- make_vector (w);
+ make_range_vector (v);
+ make_domain_vector (w);
for (unsigned int i=0; i<w.size(); ++i)
w(i) = w(i)+std::complex<number>(1.);
deallog << m.matrix_scalar_product (v,w) << std::endl;
}
-
-DEAL::(5.0e+03,0.0)
-DEAL::(5.0e+03,0.0)
+DEAL::(3.0e+03,0.0)
+DEAL::(3.0e+03,0.0)
-DEAL::90.
-DEAL::90.
+DEAL::75.
+DEAL::75.
-DEAL::1.0e+02
-DEAL::1.0e+02
+DEAL::70.
+DEAL::70.
-DEAL::70.
-DEAL::70.
+DEAL::54.
+DEAL::54.
-(2.0,0.0)(3.0,0.0)(4.0,0.0)(5.0,0.0)(6.0,0.0)
-(4.0,0.0)(6.0,0.0)(8.0,0.0)(10.,0.0)(12.,0.0)
-(6.0,0.0)(9.0,0.0)(12.,0.0)(15.,0.0)(18.,0.0)
-(8.0,0.0)(12.,0.0)(16.,0.0)(20.,0.0)(24.,0.0)
-(10.,0.0)(15.,0.0)(20.,0.0)(25.,0.0)(30.,0.0)
-(2.0,0.0)(3.0,0.0)(4.0,0.0)(5.0,0.0)(6.0,0.0)
-(4.0,0.0)(6.0,0.0)(8.0,0.0)(10.,0.0)(12.,0.0)
-(6.0,0.0)(9.0,0.0)(12.,0.0)(15.,0.0)(18.,0.0)
-(8.0,0.0)(12.,0.0)(16.,0.0)(20.,0.0)(24.,0.0)
-(10.,0.0)(15.,0.0)(20.,0.0)(25.,0.0)(30.,0.0)
+(2.0,0.0)(3.0,0.0)(4.0,0.0)(5.0,0.0)
+(4.0,0.0)(6.0,0.0)(8.0,0.0)(10.,0.0)
+(6.0,0.0)(9.0,0.0)(12.,0.0)(15.,0.0)
+(8.0,0.0)(12.,0.0)(16.,0.0)(20.,0.0)
+(10.,0.0)(15.,0.0)(20.,0.0)(25.,0.0)
+(2.0,0.0)(3.0,0.0)(4.0,0.0)(5.0,0.0)
+(4.0,0.0)(6.0,0.0)(8.0,0.0)(10.,0.0)
+(6.0,0.0)(9.0,0.0)(12.,0.0)(15.,0.0)
+(8.0,0.0)(12.,0.0)(16.,0.0)(20.,0.0)
+(10.,0.0)(15.,0.0)(20.,0.0)(25.,0.0)
DEAL::3 2 (0.0,0.0)
DEAL::3 3 (0.0,0.0)
DEAL::3 4 (30.,0.0)
-DEAL::3 5 (36.,0.0)
+DEAL::3 5 (0.0,0.0)
DEAL::3 6 (0.0,0.0)
DEAL::3 7 (0.0,0.0)
DEAL::3 8 (0.0,0.0)
DEAL::4 2 (0.0,0.0)
DEAL::4 3 (0.0,0.0)
DEAL::4 4 (40.,0.0)
-DEAL::4 5 (48.,0.0)
+DEAL::4 5 (0.0,0.0)
DEAL::4 6 (0.0,0.0)
DEAL::4 7 (0.0,0.0)
DEAL::4 8 (0.0,0.0)
DEAL::5 2 (0.0,0.0)
DEAL::5 3 (0.0,0.0)
DEAL::5 4 (50.,0.0)
-DEAL::5 5 (60.,0.0)
+DEAL::5 5 (0.0,0.0)
DEAL::5 6 (0.0,0.0)
DEAL::5 7 (0.0,0.0)
DEAL::5 8 (0.0,0.0)
DEAL::3 2 (0.0,0.0)
DEAL::3 3 (0.0,0.0)
DEAL::3 4 (30.,0.0)
-DEAL::3 5 (36.,0.0)
+DEAL::3 5 (0.0,0.0)
DEAL::3 6 (0.0,0.0)
DEAL::3 7 (0.0,0.0)
DEAL::3 8 (0.0,0.0)
DEAL::4 2 (0.0,0.0)
DEAL::4 3 (0.0,0.0)
DEAL::4 4 (40.,0.0)
-DEAL::4 5 (48.,0.0)
+DEAL::4 5 (0.0,0.0)
DEAL::4 6 (0.0,0.0)
DEAL::4 7 (0.0,0.0)
DEAL::4 8 (0.0,0.0)
DEAL::5 2 (0.0,0.0)
DEAL::5 3 (0.0,0.0)
DEAL::5 4 (50.,0.0)
-DEAL::5 5 (60.,0.0)
+DEAL::5 5 (0.0,0.0)
DEAL::5 6 (0.0,0.0)
DEAL::5 7 (0.0,0.0)
DEAL::5 8 (0.0,0.0)
DEAL::4 1 (0.0,0.0)
DEAL::4 2 (0.0,0.0)
DEAL::4 3 (0.0,0.0)
-DEAL::4 4 (36.,0.0)
-DEAL::4 5 (48.,0.0)
-DEAL::4 6 (60.,0.0)
+DEAL::4 4 (0.0,0.0)
+DEAL::4 5 (0.0,0.0)
+DEAL::4 6 (0.0,0.0)
DEAL::4 7 (0.0,0.0)
DEAL::4 8 (0.0,0.0)
DEAL::4 9 (0.0,0.0)
DEAL::4 1 (0.0,0.0)
DEAL::4 2 (0.0,0.0)
DEAL::4 3 (0.0,0.0)
-DEAL::4 4 (36.,0.0)
-DEAL::4 5 (48.,0.0)
-DEAL::4 6 (60.,0.0)
+DEAL::4 4 (0.0,0.0)
+DEAL::4 5 (0.0,0.0)
+DEAL::4 6 (0.0,0.0)
DEAL::4 7 (0.0,0.0)
DEAL::4 8 (0.0,0.0)
DEAL::4 9 (0.0,0.0)
{
FullMatrix<std::complex<number> > m;
make_matrix (m);
- m.add_col (2, 3.141, 4);
+ m.add_col (2, 3.141, 3);
print_matrix (m);
}
-
DEAL::0 0 (2.0,0.0)
DEAL::0 1 (3.0,0.0)
-DEAL::0 2 (23.,0.0)
+DEAL::0 2 (20.,0.0)
DEAL::0 3 (5.0,0.0)
-DEAL::0 4 (6.0,0.0)
DEAL::1 0 (4.0,0.0)
DEAL::1 1 (6.0,0.0)
-DEAL::1 2 (46.,0.0)
+DEAL::1 2 (39.,0.0)
DEAL::1 3 (10.,0.0)
-DEAL::1 4 (12.,0.0)
DEAL::2 0 (6.0,0.0)
DEAL::2 1 (9.0,0.0)
-DEAL::2 2 (69.,0.0)
+DEAL::2 2 (59.,0.0)
DEAL::2 3 (15.,0.0)
-DEAL::2 4 (18.,0.0)
DEAL::3 0 (8.0,0.0)
DEAL::3 1 (12.,0.0)
-DEAL::3 2 (91.,0.0)
+DEAL::3 2 (79.,0.0)
DEAL::3 3 (20.,0.0)
-DEAL::3 4 (24.,0.0)
DEAL::4 0 (10.,0.0)
DEAL::4 1 (15.,0.0)
-DEAL::4 2 (1.1e+02,0.0)
+DEAL::4 2 (99.,0.0)
DEAL::4 3 (25.,0.0)
-DEAL::4 4 (30.,0.0)
DEAL::0 0 (2.0,0.0)
DEAL::0 1 (3.0,0.0)
-DEAL::0 2 (23.,0.0)
+DEAL::0 2 (20.,0.0)
DEAL::0 3 (5.0,0.0)
-DEAL::0 4 (6.0,0.0)
DEAL::1 0 (4.0,0.0)
DEAL::1 1 (6.0,0.0)
-DEAL::1 2 (46.,0.0)
+DEAL::1 2 (39.,0.0)
DEAL::1 3 (10.,0.0)
-DEAL::1 4 (12.,0.0)
DEAL::2 0 (6.0,0.0)
DEAL::2 1 (9.0,0.0)
-DEAL::2 2 (69.,0.0)
+DEAL::2 2 (59.,0.0)
DEAL::2 3 (15.,0.0)
-DEAL::2 4 (18.,0.0)
DEAL::3 0 (8.0,0.0)
DEAL::3 1 (12.,0.0)
-DEAL::3 2 (91.,0.0)
+DEAL::3 2 (79.,0.0)
DEAL::3 3 (20.,0.0)
-DEAL::3 4 (24.,0.0)
DEAL::4 0 (10.,0.0)
DEAL::4 1 (15.,0.0)
-DEAL::4 2 (1.1e+02,0.0)
+DEAL::4 2 (99.,0.0)
DEAL::4 3 (25.,0.0)
-DEAL::4 4 (30.,0.0)
{
FullMatrix<std::complex<number> > m;
make_matrix (m);
- m.add_col (2, 3.141, 4, 2.718, 3);
+ m.add_col (2, 3.141, 1, 2.718, 3);
print_matrix (m);
}
-
DEAL::0 0 (2.0,0.0)
DEAL::0 1 (3.0,0.0)
-DEAL::0 2 (36.,0.0)
+DEAL::0 2 (27.,0.0)
DEAL::0 3 (5.0,0.0)
-DEAL::0 4 (6.0,0.0)
DEAL::1 0 (4.0,0.0)
DEAL::1 1 (6.0,0.0)
-DEAL::1 2 (73.,0.0)
+DEAL::1 2 (54.,0.0)
DEAL::1 3 (10.,0.0)
-DEAL::1 4 (12.,0.0)
DEAL::2 0 (6.0,0.0)
DEAL::2 1 (9.0,0.0)
-DEAL::2 2 (1.1e+02,0.0)
+DEAL::2 2 (81.,0.0)
DEAL::2 3 (15.,0.0)
-DEAL::2 4 (18.,0.0)
DEAL::3 0 (8.0,0.0)
DEAL::3 1 (12.,0.0)
-DEAL::3 2 (1.5e+02,0.0)
+DEAL::3 2 (1.1e+02,0.0)
DEAL::3 3 (20.,0.0)
-DEAL::3 4 (24.,0.0)
DEAL::4 0 (10.,0.0)
DEAL::4 1 (15.,0.0)
-DEAL::4 2 (1.8e+02,0.0)
+DEAL::4 2 (1.4e+02,0.0)
DEAL::4 3 (25.,0.0)
-DEAL::4 4 (30.,0.0)
DEAL::0 0 (2.0,0.0)
DEAL::0 1 (3.0,0.0)
-DEAL::0 2 (36.,0.0)
+DEAL::0 2 (27.,0.0)
DEAL::0 3 (5.0,0.0)
-DEAL::0 4 (6.0,0.0)
DEAL::1 0 (4.0,0.0)
DEAL::1 1 (6.0,0.0)
-DEAL::1 2 (73.,0.0)
+DEAL::1 2 (54.,0.0)
DEAL::1 3 (10.,0.0)
-DEAL::1 4 (12.,0.0)
DEAL::2 0 (6.0,0.0)
DEAL::2 1 (9.0,0.0)
-DEAL::2 2 (1.1e+02,0.0)
+DEAL::2 2 (81.,0.0)
DEAL::2 3 (15.,0.0)
-DEAL::2 4 (18.,0.0)
DEAL::3 0 (8.0,0.0)
DEAL::3 1 (12.,0.0)
-DEAL::3 2 (1.5e+02,0.0)
+DEAL::3 2 (1.1e+02,0.0)
DEAL::3 3 (20.,0.0)
-DEAL::3 4 (24.,0.0)
DEAL::4 0 (10.,0.0)
DEAL::4 1 (15.,0.0)
-DEAL::4 2 (1.8e+02,0.0)
+DEAL::4 2 (1.4e+02,0.0)
DEAL::4 3 (25.,0.0)
-DEAL::4 4 (30.,0.0)
{
FullMatrix<std::complex<number> > m;
make_matrix (m);
- m.swap_col (2, 4);
+ m.swap_col (2, 3);
print_matrix (m);
}
-
DEAL::0 0 (2.0,0.0)
DEAL::0 1 (3.0,0.0)
-DEAL::0 2 (6.0,0.0)
-DEAL::0 3 (5.0,0.0)
-DEAL::0 4 (4.0,0.0)
+DEAL::0 2 (5.0,0.0)
+DEAL::0 3 (4.0,0.0)
DEAL::1 0 (4.0,0.0)
DEAL::1 1 (6.0,0.0)
-DEAL::1 2 (12.,0.0)
-DEAL::1 3 (10.,0.0)
-DEAL::1 4 (8.0,0.0)
+DEAL::1 2 (10.,0.0)
+DEAL::1 3 (8.0,0.0)
DEAL::2 0 (6.0,0.0)
DEAL::2 1 (9.0,0.0)
-DEAL::2 2 (18.,0.0)
-DEAL::2 3 (15.,0.0)
-DEAL::2 4 (12.,0.0)
+DEAL::2 2 (15.,0.0)
+DEAL::2 3 (12.,0.0)
DEAL::3 0 (8.0,0.0)
DEAL::3 1 (12.,0.0)
-DEAL::3 2 (24.,0.0)
-DEAL::3 3 (20.,0.0)
-DEAL::3 4 (16.,0.0)
+DEAL::3 2 (20.,0.0)
+DEAL::3 3 (16.,0.0)
DEAL::4 0 (10.,0.0)
DEAL::4 1 (15.,0.0)
-DEAL::4 2 (30.,0.0)
-DEAL::4 3 (25.,0.0)
-DEAL::4 4 (20.,0.0)
+DEAL::4 2 (25.,0.0)
+DEAL::4 3 (20.,0.0)
DEAL::0 0 (2.0,0.0)
DEAL::0 1 (3.0,0.0)
-DEAL::0 2 (6.0,0.0)
-DEAL::0 3 (5.0,0.0)
-DEAL::0 4 (4.0,0.0)
+DEAL::0 2 (5.0,0.0)
+DEAL::0 3 (4.0,0.0)
DEAL::1 0 (4.0,0.0)
DEAL::1 1 (6.0,0.0)
-DEAL::1 2 (12.,0.0)
-DEAL::1 3 (10.,0.0)
-DEAL::1 4 (8.0,0.0)
+DEAL::1 2 (10.,0.0)
+DEAL::1 3 (8.0,0.0)
DEAL::2 0 (6.0,0.0)
DEAL::2 1 (9.0,0.0)
-DEAL::2 2 (18.,0.0)
-DEAL::2 3 (15.,0.0)
-DEAL::2 4 (12.,0.0)
+DEAL::2 2 (15.,0.0)
+DEAL::2 3 (12.,0.0)
DEAL::3 0 (8.0,0.0)
DEAL::3 1 (12.,0.0)
-DEAL::3 2 (24.,0.0)
-DEAL::3 3 (20.,0.0)
-DEAL::3 4 (16.,0.0)
+DEAL::3 2 (20.,0.0)
+DEAL::3 3 (16.,0.0)
DEAL::4 0 (10.,0.0)
DEAL::4 1 (15.,0.0)
-DEAL::4 2 (30.,0.0)
-DEAL::4 3 (25.,0.0)
-DEAL::4 4 (20.,0.0)
+DEAL::4 2 (25.,0.0)
+DEAL::4 3 (20.,0.0)
FullMatrix<std::complex<number> > m;
make_matrix (m);
Vector<std::complex<number> > v, w;
- make_vector (v);
- make_vector (w);
+ make_range_vector (v);
+ make_domain_vector (w);
m.vmult (v, w, true);
print_vector (v);
print_vector (v);
print_vector (w);
}
-
-DEAL::0 (71.,0.0)
-DEAL::1 (1.4e+02,0.0)
-DEAL::2 (2.1e+02,0.0)
-DEAL::3 (2.8e+02,0.0)
-DEAL::4 (3.6e+02,0.0)
+DEAL::0 (41.,0.0)
+DEAL::1 (82.,0.0)
+DEAL::2 (1.2e+02,0.0)
+DEAL::3 (1.6e+02,0.0)
+DEAL::4 (2.0e+02,0.0)
DEAL::0 (1.0,0.0)
DEAL::1 (2.0,0.0)
DEAL::2 (3.0,0.0)
DEAL::3 (4.0,0.0)
-DEAL::4 (5.0,0.0)
-DEAL::0 (70.,0.0)
-DEAL::1 (1.4e+02,0.0)
-DEAL::2 (2.1e+02,0.0)
-DEAL::3 (2.8e+02,0.0)
-DEAL::4 (3.5e+02,0.0)
+DEAL::0 (40.,0.0)
+DEAL::1 (80.,0.0)
+DEAL::2 (1.2e+02,0.0)
+DEAL::3 (1.6e+02,0.0)
+DEAL::4 (2.0e+02,0.0)
DEAL::0 (1.0,0.0)
DEAL::1 (2.0,0.0)
DEAL::2 (3.0,0.0)
DEAL::3 (4.0,0.0)
-DEAL::4 (5.0,0.0)
-DEAL::0 (71.,0.0)
-DEAL::1 (1.4e+02,0.0)
-DEAL::2 (2.1e+02,0.0)
-DEAL::3 (2.8e+02,0.0)
-DEAL::4 (3.6e+02,0.0)
+DEAL::0 (41.,0.0)
+DEAL::1 (82.,0.0)
+DEAL::2 (1.2e+02,0.0)
+DEAL::3 (1.6e+02,0.0)
+DEAL::4 (2.0e+02,0.0)
DEAL::0 (1.0,0.0)
DEAL::1 (2.0,0.0)
DEAL::2 (3.0,0.0)
DEAL::3 (4.0,0.0)
-DEAL::4 (5.0,0.0)
-DEAL::0 (70.,0.0)
-DEAL::1 (1.4e+02,0.0)
-DEAL::2 (2.1e+02,0.0)
-DEAL::3 (2.8e+02,0.0)
-DEAL::4 (3.5e+02,0.0)
+DEAL::0 (40.,0.0)
+DEAL::1 (80.,0.0)
+DEAL::2 (1.2e+02,0.0)
+DEAL::3 (1.6e+02,0.0)
+DEAL::4 (2.0e+02,0.0)
DEAL::0 (1.0,0.0)
DEAL::1 (2.0,0.0)
DEAL::2 (3.0,0.0)
DEAL::3 (4.0,0.0)
-DEAL::4 (5.0,0.0)
FullMatrix<std::complex<number> > m;
make_matrix (m);
Vector<std::complex<number> > v, w;
- make_vector (v);
- make_vector (w);
+ make_range_vector (v);
+ make_domain_vector (w);
m.vmult_add (v, w);
print_vector (v);
print_vector (w);
}
-
-DEAL::0 (71.,0.0)
-DEAL::1 (1.4e+02,0.0)
-DEAL::2 (2.1e+02,0.0)
-DEAL::3 (2.8e+02,0.0)
-DEAL::4 (3.6e+02,0.0)
+DEAL::0 (41.,0.0)
+DEAL::1 (82.,0.0)
+DEAL::2 (1.2e+02,0.0)
+DEAL::3 (1.6e+02,0.0)
+DEAL::4 (2.0e+02,0.0)
DEAL::0 (1.0,0.0)
DEAL::1 (2.0,0.0)
DEAL::2 (3.0,0.0)
DEAL::3 (4.0,0.0)
-DEAL::4 (5.0,0.0)
-DEAL::0 (71.,0.0)
-DEAL::1 (1.4e+02,0.0)
-DEAL::2 (2.1e+02,0.0)
-DEAL::3 (2.8e+02,0.0)
-DEAL::4 (3.6e+02,0.0)
+DEAL::0 (41.,0.0)
+DEAL::1 (82.,0.0)
+DEAL::2 (1.2e+02,0.0)
+DEAL::3 (1.6e+02,0.0)
+DEAL::4 (2.0e+02,0.0)
DEAL::0 (1.0,0.0)
DEAL::1 (2.0,0.0)
DEAL::2 (3.0,0.0)
DEAL::3 (4.0,0.0)
-DEAL::4 (5.0,0.0)
FullMatrix<std::complex<number> > m;
make_matrix (m);
Vector<std::complex<number> > v, w, x;
- make_vector (v);
- make_vector (w);
- make_vector (x);
+ make_range_vector (v);
+ make_domain_vector (w);
+ make_range_vector (x);
m.residual (v, w, x);
print_vector (v);
print_vector (w);
print_vector (x);
}
-
-DEAL::0 (-69.,0.0)
-DEAL::1 (-1.4e+02,0.0)
-DEAL::2 (-2.1e+02,0.0)
-DEAL::3 (-2.8e+02,0.0)
-DEAL::4 (-3.4e+02,0.0)
+DEAL::0 (-39.,0.0)
+DEAL::1 (-78.,0.0)
+DEAL::2 (-1.2e+02,0.0)
+DEAL::3 (-1.6e+02,0.0)
+DEAL::4 (-2.0e+02,0.0)
DEAL::0 (1.0,0.0)
DEAL::1 (2.0,0.0)
DEAL::2 (3.0,0.0)
DEAL::3 (4.0,0.0)
-DEAL::4 (5.0,0.0)
DEAL::0 (1.0,0.0)
DEAL::1 (2.0,0.0)
DEAL::2 (3.0,0.0)
DEAL::3 (4.0,0.0)
DEAL::4 (5.0,0.0)
-DEAL::0 (-69.,0.0)
-DEAL::1 (-1.4e+02,0.0)
-DEAL::2 (-2.1e+02,0.0)
-DEAL::3 (-2.8e+02,0.0)
-DEAL::4 (-3.4e+02,0.0)
+DEAL::0 (-39.,0.0)
+DEAL::1 (-78.,0.0)
+DEAL::2 (-1.2e+02,0.0)
+DEAL::3 (-1.6e+02,0.0)
+DEAL::4 (-2.0e+02,0.0)
DEAL::0 (1.0,0.0)
DEAL::1 (2.0,0.0)
DEAL::2 (3.0,0.0)
DEAL::3 (4.0,0.0)
-DEAL::4 (5.0,0.0)
DEAL::0 (1.0,0.0)
DEAL::1 (2.0,0.0)
DEAL::2 (3.0,0.0)
FullMatrix<std::complex<number> > m;
make_matrix (m);
Vector<std::complex<number> > x, y;
- make_vector (x);
- make_vector (y);
+ make_range_vector (x);
+ make_domain_vector (y);
m.forward (x, y);
print_vector (x);
print_vector (y);
}
-
DEAL::1 (0.0,0.0)
DEAL::2 (0.0,0.0)
DEAL::3 (0.0,0.0)
-DEAL::4 (0.0,0.0)
+DEAL::4 (5.0,0.0)
DEAL::0 (1.0,0.0)
DEAL::1 (2.0,0.0)
DEAL::2 (3.0,0.0)
DEAL::3 (4.0,0.0)
-DEAL::4 (5.0,0.0)
DEAL::0 (0.50,0.0)
DEAL::1 (0.0,0.0)
DEAL::2 (0.0,0.0)
DEAL::3 (0.0,0.0)
-DEAL::4 (0.0,0.0)
+DEAL::4 (5.0,0.0)
DEAL::0 (1.0,0.0)
DEAL::1 (2.0,0.0)
DEAL::2 (3.0,0.0)
DEAL::3 (4.0,0.0)
-DEAL::4 (5.0,0.0)
FullMatrix<std::complex<number> > m;
make_matrix (m);
Vector<std::complex<number> > x, y;
- make_vector (x);
- make_vector (y);
+ make_range_vector (x);
+ make_domain_vector (y);
+ for (std::complex<number> &entry : y)
+ entry += std::complex<number>(1.0, 0.0);
m.backward (x, y);
print_vector (x);
print_vector (y);
}
-
-DEAL::0 (0.0,0.0)
-DEAL::1 (0.0,0.0)
-DEAL::2 (0.0,0.0)
-DEAL::3 (0.0,0.0)
-DEAL::4 (0.17,0.0)
-DEAL::0 (1.0,0.0)
-DEAL::1 (2.0,0.0)
-DEAL::2 (3.0,0.0)
-DEAL::3 (4.0,0.0)
+DEAL::0 (0.25,0.0)
+DEAL::1 (0.056,0.0)
+DEAL::2 (0.021,0.0)
+DEAL::3 (0.25,0.0)
DEAL::4 (5.0,0.0)
-DEAL::0 (0.0,0.0)
-DEAL::1 (0.0,0.0)
-DEAL::2 (0.0,0.0)
-DEAL::3 (0.0,0.0)
-DEAL::4 (0.17,0.0)
-DEAL::0 (1.0,0.0)
-DEAL::1 (2.0,0.0)
-DEAL::2 (3.0,0.0)
-DEAL::3 (4.0,0.0)
+DEAL::0 (2.0,0.0)
+DEAL::1 (3.0,0.0)
+DEAL::2 (4.0,0.0)
+DEAL::3 (5.0,0.0)
+DEAL::0 (0.25,0.0)
+DEAL::1 (0.056,0.0)
+DEAL::2 (0.021,0.0)
+DEAL::3 (0.25,0.0)
DEAL::4 (5.0,0.0)
+DEAL::0 (2.0,0.0)
+DEAL::1 (3.0,0.0)
+DEAL::2 (4.0,0.0)
+DEAL::3 (5.0,0.0)
-DEAL::90.
-DEAL::90.
+DEAL::75.
+DEAL::75.
-DEAL::1.0e+02
-DEAL::1.0e+02
+DEAL::70.
+DEAL::70.
-DEAL::70.
-DEAL::70.
+DEAL::54.
+DEAL::54.
-2.0 3.0 4.0 5.0 6.0
-4.0 6.0 8.0 10. 12.
-6.0 9.0 12. 15. 18.
-8.0 12. 16. 20. 24.
-10. 15. 20. 25. 30.
-2.0 3.0 4.0 5.0 6.0
-4.0 6.0 8.0 10. 12.
-6.0 9.0 12. 15. 18.
-8.0 12. 16. 20. 24.
-10. 15. 20. 25. 30.
+2.0 3.0 4.0 5.0
+4.0 6.0 8.0 10.
+6.0 9.0 12. 15.
+8.0 12. 16. 20.
+10. 15. 20. 25.
+2.0 3.0 4.0 5.0
+4.0 6.0 8.0 10.
+6.0 9.0 12. 15.
+8.0 12. 16. 20.
+10. 15. 20. 25.
-DEAL::0 0 0
-DEAL::0 1 0
-DEAL::0 2 0
-DEAL::0 3 0
-DEAL::0 4 0
-DEAL::0 5 0
-DEAL::0 6 0
-DEAL::0 7 0
-DEAL::0 8 0
-DEAL::0 9 0
-DEAL::1 0 0
-DEAL::1 1 0
-DEAL::1 2 0
-DEAL::1 3 0
-DEAL::1 4 0
-DEAL::1 5 0
-DEAL::1 6 0
-DEAL::1 7 0
-DEAL::1 8 0
-DEAL::1 9 0
-DEAL::2 0 0
-DEAL::2 1 0
-DEAL::2 2 0
-DEAL::2 3 0
-DEAL::2 4 0
-DEAL::2 5 0
-DEAL::2 6 0
-DEAL::2 7 0
-DEAL::2 8 0
-DEAL::2 9 0
-DEAL::3 0 0
-DEAL::3 1 0
-DEAL::3 2 0
-DEAL::3 3 0
+DEAL::0 0 0.0
+DEAL::0 1 0.0
+DEAL::0 2 0.0
+DEAL::0 3 0.0
+DEAL::0 4 0.0
+DEAL::0 5 0.0
+DEAL::0 6 0.0
+DEAL::0 7 0.0
+DEAL::0 8 0.0
+DEAL::0 9 0.0
+DEAL::1 0 0.0
+DEAL::1 1 0.0
+DEAL::1 2 0.0
+DEAL::1 3 0.0
+DEAL::1 4 0.0
+DEAL::1 5 0.0
+DEAL::1 6 0.0
+DEAL::1 7 0.0
+DEAL::1 8 0.0
+DEAL::1 9 0.0
+DEAL::2 0 0.0
+DEAL::2 1 0.0
+DEAL::2 2 0.0
+DEAL::2 3 0.0
+DEAL::2 4 0.0
+DEAL::2 5 0.0
+DEAL::2 6 0.0
+DEAL::2 7 0.0
+DEAL::2 8 0.0
+DEAL::2 9 0.0
+DEAL::3 0 0.0
+DEAL::3 1 0.0
+DEAL::3 2 0.0
+DEAL::3 3 0.0
DEAL::3 4 30.
-DEAL::3 5 36.
-DEAL::3 6 0
-DEAL::3 7 0
-DEAL::3 8 0
-DEAL::3 9 0
-DEAL::4 0 0
-DEAL::4 1 0
-DEAL::4 2 0
-DEAL::4 3 0
+DEAL::3 5 0.0
+DEAL::3 6 0.0
+DEAL::3 7 0.0
+DEAL::3 8 0.0
+DEAL::3 9 0.0
+DEAL::4 0 0.0
+DEAL::4 1 0.0
+DEAL::4 2 0.0
+DEAL::4 3 0.0
DEAL::4 4 40.
-DEAL::4 5 48.
-DEAL::4 6 0
-DEAL::4 7 0
-DEAL::4 8 0
-DEAL::4 9 0
-DEAL::5 0 0
-DEAL::5 1 0
-DEAL::5 2 0
-DEAL::5 3 0
+DEAL::4 5 0.0
+DEAL::4 6 0.0
+DEAL::4 7 0.0
+DEAL::4 8 0.0
+DEAL::4 9 0.0
+DEAL::5 0 0.0
+DEAL::5 1 0.0
+DEAL::5 2 0.0
+DEAL::5 3 0.0
DEAL::5 4 50.
-DEAL::5 5 60.
-DEAL::5 6 0
-DEAL::5 7 0
-DEAL::5 8 0
-DEAL::5 9 0
-DEAL::6 0 0
-DEAL::6 1 0
-DEAL::6 2 0
-DEAL::6 3 0
-DEAL::6 4 0
-DEAL::6 5 0
-DEAL::6 6 0
-DEAL::6 7 0
-DEAL::6 8 0
-DEAL::6 9 0
-DEAL::7 0 0
-DEAL::7 1 0
-DEAL::7 2 0
-DEAL::7 3 0
-DEAL::7 4 0
-DEAL::7 5 0
-DEAL::7 6 0
-DEAL::7 7 0
-DEAL::7 8 0
-DEAL::7 9 0
-DEAL::8 0 0
-DEAL::8 1 0
-DEAL::8 2 0
-DEAL::8 3 0
-DEAL::8 4 0
-DEAL::8 5 0
-DEAL::8 6 0
-DEAL::8 7 0
-DEAL::8 8 0
-DEAL::8 9 0
-DEAL::9 0 0
-DEAL::9 1 0
-DEAL::9 2 0
-DEAL::9 3 0
-DEAL::9 4 0
-DEAL::9 5 0
-DEAL::9 6 0
-DEAL::9 7 0
-DEAL::9 8 0
-DEAL::9 9 0
-DEAL::0 0 0
-DEAL::0 1 0
-DEAL::0 2 0
-DEAL::0 3 0
-DEAL::0 4 0
-DEAL::0 5 0
-DEAL::0 6 0
-DEAL::0 7 0
-DEAL::0 8 0
-DEAL::0 9 0
-DEAL::1 0 0
-DEAL::1 1 0
-DEAL::1 2 0
-DEAL::1 3 0
-DEAL::1 4 0
-DEAL::1 5 0
-DEAL::1 6 0
-DEAL::1 7 0
-DEAL::1 8 0
-DEAL::1 9 0
-DEAL::2 0 0
-DEAL::2 1 0
-DEAL::2 2 0
-DEAL::2 3 0
-DEAL::2 4 0
-DEAL::2 5 0
-DEAL::2 6 0
-DEAL::2 7 0
-DEAL::2 8 0
-DEAL::2 9 0
-DEAL::3 0 0
-DEAL::3 1 0
-DEAL::3 2 0
-DEAL::3 3 0
+DEAL::5 5 0.0
+DEAL::5 6 0.0
+DEAL::5 7 0.0
+DEAL::5 8 0.0
+DEAL::5 9 0.0
+DEAL::6 0 0.0
+DEAL::6 1 0.0
+DEAL::6 2 0.0
+DEAL::6 3 0.0
+DEAL::6 4 0.0
+DEAL::6 5 0.0
+DEAL::6 6 0.0
+DEAL::6 7 0.0
+DEAL::6 8 0.0
+DEAL::6 9 0.0
+DEAL::7 0 0.0
+DEAL::7 1 0.0
+DEAL::7 2 0.0
+DEAL::7 3 0.0
+DEAL::7 4 0.0
+DEAL::7 5 0.0
+DEAL::7 6 0.0
+DEAL::7 7 0.0
+DEAL::7 8 0.0
+DEAL::7 9 0.0
+DEAL::8 0 0.0
+DEAL::8 1 0.0
+DEAL::8 2 0.0
+DEAL::8 3 0.0
+DEAL::8 4 0.0
+DEAL::8 5 0.0
+DEAL::8 6 0.0
+DEAL::8 7 0.0
+DEAL::8 8 0.0
+DEAL::8 9 0.0
+DEAL::9 0 0.0
+DEAL::9 1 0.0
+DEAL::9 2 0.0
+DEAL::9 3 0.0
+DEAL::9 4 0.0
+DEAL::9 5 0.0
+DEAL::9 6 0.0
+DEAL::9 7 0.0
+DEAL::9 8 0.0
+DEAL::9 9 0.0
+DEAL::0 0 0.0
+DEAL::0 1 0.0
+DEAL::0 2 0.0
+DEAL::0 3 0.0
+DEAL::0 4 0.0
+DEAL::0 5 0.0
+DEAL::0 6 0.0
+DEAL::0 7 0.0
+DEAL::0 8 0.0
+DEAL::0 9 0.0
+DEAL::1 0 0.0
+DEAL::1 1 0.0
+DEAL::1 2 0.0
+DEAL::1 3 0.0
+DEAL::1 4 0.0
+DEAL::1 5 0.0
+DEAL::1 6 0.0
+DEAL::1 7 0.0
+DEAL::1 8 0.0
+DEAL::1 9 0.0
+DEAL::2 0 0.0
+DEAL::2 1 0.0
+DEAL::2 2 0.0
+DEAL::2 3 0.0
+DEAL::2 4 0.0
+DEAL::2 5 0.0
+DEAL::2 6 0.0
+DEAL::2 7 0.0
+DEAL::2 8 0.0
+DEAL::2 9 0.0
+DEAL::3 0 0.0
+DEAL::3 1 0.0
+DEAL::3 2 0.0
+DEAL::3 3 0.0
DEAL::3 4 30.
-DEAL::3 5 36.
-DEAL::3 6 0
-DEAL::3 7 0
-DEAL::3 8 0
-DEAL::3 9 0
-DEAL::4 0 0
-DEAL::4 1 0
-DEAL::4 2 0
-DEAL::4 3 0
+DEAL::3 5 0.0
+DEAL::3 6 0.0
+DEAL::3 7 0.0
+DEAL::3 8 0.0
+DEAL::3 9 0.0
+DEAL::4 0 0.0
+DEAL::4 1 0.0
+DEAL::4 2 0.0
+DEAL::4 3 0.0
DEAL::4 4 40.
-DEAL::4 5 48.
-DEAL::4 6 0
-DEAL::4 7 0
-DEAL::4 8 0
-DEAL::4 9 0
-DEAL::5 0 0
-DEAL::5 1 0
-DEAL::5 2 0
-DEAL::5 3 0
+DEAL::4 5 0.0
+DEAL::4 6 0.0
+DEAL::4 7 0.0
+DEAL::4 8 0.0
+DEAL::4 9 0.0
+DEAL::5 0 0.0
+DEAL::5 1 0.0
+DEAL::5 2 0.0
+DEAL::5 3 0.0
DEAL::5 4 50.
-DEAL::5 5 60.
-DEAL::5 6 0
-DEAL::5 7 0
-DEAL::5 8 0
-DEAL::5 9 0
-DEAL::6 0 0
-DEAL::6 1 0
-DEAL::6 2 0
-DEAL::6 3 0
-DEAL::6 4 0
-DEAL::6 5 0
-DEAL::6 6 0
-DEAL::6 7 0
-DEAL::6 8 0
-DEAL::6 9 0
-DEAL::7 0 0
-DEAL::7 1 0
-DEAL::7 2 0
-DEAL::7 3 0
-DEAL::7 4 0
-DEAL::7 5 0
-DEAL::7 6 0
-DEAL::7 7 0
-DEAL::7 8 0
-DEAL::7 9 0
-DEAL::8 0 0
-DEAL::8 1 0
-DEAL::8 2 0
-DEAL::8 3 0
-DEAL::8 4 0
-DEAL::8 5 0
-DEAL::8 6 0
-DEAL::8 7 0
-DEAL::8 8 0
-DEAL::8 9 0
-DEAL::9 0 0
-DEAL::9 1 0
-DEAL::9 2 0
-DEAL::9 3 0
-DEAL::9 4 0
-DEAL::9 5 0
-DEAL::9 6 0
-DEAL::9 7 0
-DEAL::9 8 0
-DEAL::9 9 0
+DEAL::5 5 0.0
+DEAL::5 6 0.0
+DEAL::5 7 0.0
+DEAL::5 8 0.0
+DEAL::5 9 0.0
+DEAL::6 0 0.0
+DEAL::6 1 0.0
+DEAL::6 2 0.0
+DEAL::6 3 0.0
+DEAL::6 4 0.0
+DEAL::6 5 0.0
+DEAL::6 6 0.0
+DEAL::6 7 0.0
+DEAL::6 8 0.0
+DEAL::6 9 0.0
+DEAL::7 0 0.0
+DEAL::7 1 0.0
+DEAL::7 2 0.0
+DEAL::7 3 0.0
+DEAL::7 4 0.0
+DEAL::7 5 0.0
+DEAL::7 6 0.0
+DEAL::7 7 0.0
+DEAL::7 8 0.0
+DEAL::7 9 0.0
+DEAL::8 0 0.0
+DEAL::8 1 0.0
+DEAL::8 2 0.0
+DEAL::8 3 0.0
+DEAL::8 4 0.0
+DEAL::8 5 0.0
+DEAL::8 6 0.0
+DEAL::8 7 0.0
+DEAL::8 8 0.0
+DEAL::8 9 0.0
+DEAL::9 0 0.0
+DEAL::9 1 0.0
+DEAL::9 2 0.0
+DEAL::9 3 0.0
+DEAL::9 4 0.0
+DEAL::9 5 0.0
+DEAL::9 6 0.0
+DEAL::9 7 0.0
+DEAL::9 8 0.0
+DEAL::9 9 0.0
-DEAL::0 0 0
-DEAL::0 1 0
-DEAL::0 2 0
-DEAL::0 3 0
-DEAL::0 4 0
-DEAL::0 5 0
-DEAL::0 6 0
-DEAL::0 7 0
-DEAL::0 8 0
-DEAL::0 9 0
-DEAL::1 0 0
-DEAL::1 1 0
-DEAL::1 2 0
-DEAL::1 3 0
-DEAL::1 4 0
-DEAL::1 5 0
-DEAL::1 6 0
-DEAL::1 7 0
-DEAL::1 8 0
-DEAL::1 9 0
-DEAL::2 0 0
-DEAL::2 1 0
-DEAL::2 2 0
-DEAL::2 3 0
-DEAL::2 4 0
-DEAL::2 5 0
-DEAL::2 6 0
-DEAL::2 7 0
-DEAL::2 8 0
-DEAL::2 9 0
-DEAL::3 0 0
-DEAL::3 1 0
-DEAL::3 2 0
-DEAL::3 3 0
+DEAL::0 0 0.0
+DEAL::0 1 0.0
+DEAL::0 2 0.0
+DEAL::0 3 0.0
+DEAL::0 4 0.0
+DEAL::0 5 0.0
+DEAL::0 6 0.0
+DEAL::0 7 0.0
+DEAL::0 8 0.0
+DEAL::0 9 0.0
+DEAL::1 0 0.0
+DEAL::1 1 0.0
+DEAL::1 2 0.0
+DEAL::1 3 0.0
+DEAL::1 4 0.0
+DEAL::1 5 0.0
+DEAL::1 6 0.0
+DEAL::1 7 0.0
+DEAL::1 8 0.0
+DEAL::1 9 0.0
+DEAL::2 0 0.0
+DEAL::2 1 0.0
+DEAL::2 2 0.0
+DEAL::2 3 0.0
+DEAL::2 4 0.0
+DEAL::2 5 0.0
+DEAL::2 6 0.0
+DEAL::2 7 0.0
+DEAL::2 8 0.0
+DEAL::2 9 0.0
+DEAL::3 0 0.0
+DEAL::3 1 0.0
+DEAL::3 2 0.0
+DEAL::3 3 0.0
DEAL::3 4 30.
DEAL::3 5 40.
DEAL::3 6 50.
-DEAL::3 7 0
-DEAL::3 8 0
-DEAL::3 9 0
-DEAL::4 0 0
-DEAL::4 1 0
-DEAL::4 2 0
-DEAL::4 3 0
-DEAL::4 4 36.
-DEAL::4 5 48.
-DEAL::4 6 60.
-DEAL::4 7 0
-DEAL::4 8 0
-DEAL::4 9 0
-DEAL::5 0 0
-DEAL::5 1 0
-DEAL::5 2 0
-DEAL::5 3 0
-DEAL::5 4 0
-DEAL::5 5 0
-DEAL::5 6 0
-DEAL::5 7 0
-DEAL::5 8 0
-DEAL::5 9 0
-DEAL::6 0 0
-DEAL::6 1 0
-DEAL::6 2 0
-DEAL::6 3 0
-DEAL::6 4 0
-DEAL::6 5 0
-DEAL::6 6 0
-DEAL::6 7 0
-DEAL::6 8 0
-DEAL::6 9 0
-DEAL::7 0 0
-DEAL::7 1 0
-DEAL::7 2 0
-DEAL::7 3 0
-DEAL::7 4 0
-DEAL::7 5 0
-DEAL::7 6 0
-DEAL::7 7 0
-DEAL::7 8 0
-DEAL::7 9 0
-DEAL::8 0 0
-DEAL::8 1 0
-DEAL::8 2 0
-DEAL::8 3 0
-DEAL::8 4 0
-DEAL::8 5 0
-DEAL::8 6 0
-DEAL::8 7 0
-DEAL::8 8 0
-DEAL::8 9 0
-DEAL::9 0 0
-DEAL::9 1 0
-DEAL::9 2 0
-DEAL::9 3 0
-DEAL::9 4 0
-DEAL::9 5 0
-DEAL::9 6 0
-DEAL::9 7 0
-DEAL::9 8 0
-DEAL::9 9 0
-DEAL::0 0 0
-DEAL::0 1 0
-DEAL::0 2 0
-DEAL::0 3 0
-DEAL::0 4 0
-DEAL::0 5 0
-DEAL::0 6 0
-DEAL::0 7 0
-DEAL::0 8 0
-DEAL::0 9 0
-DEAL::1 0 0
-DEAL::1 1 0
-DEAL::1 2 0
-DEAL::1 3 0
-DEAL::1 4 0
-DEAL::1 5 0
-DEAL::1 6 0
-DEAL::1 7 0
-DEAL::1 8 0
-DEAL::1 9 0
-DEAL::2 0 0
-DEAL::2 1 0
-DEAL::2 2 0
-DEAL::2 3 0
-DEAL::2 4 0
-DEAL::2 5 0
-DEAL::2 6 0
-DEAL::2 7 0
-DEAL::2 8 0
-DEAL::2 9 0
-DEAL::3 0 0
-DEAL::3 1 0
-DEAL::3 2 0
-DEAL::3 3 0
+DEAL::3 7 0.0
+DEAL::3 8 0.0
+DEAL::3 9 0.0
+DEAL::4 0 0.0
+DEAL::4 1 0.0
+DEAL::4 2 0.0
+DEAL::4 3 0.0
+DEAL::4 4 0.0
+DEAL::4 5 0.0
+DEAL::4 6 0.0
+DEAL::4 7 0.0
+DEAL::4 8 0.0
+DEAL::4 9 0.0
+DEAL::5 0 0.0
+DEAL::5 1 0.0
+DEAL::5 2 0.0
+DEAL::5 3 0.0
+DEAL::5 4 0.0
+DEAL::5 5 0.0
+DEAL::5 6 0.0
+DEAL::5 7 0.0
+DEAL::5 8 0.0
+DEAL::5 9 0.0
+DEAL::6 0 0.0
+DEAL::6 1 0.0
+DEAL::6 2 0.0
+DEAL::6 3 0.0
+DEAL::6 4 0.0
+DEAL::6 5 0.0
+DEAL::6 6 0.0
+DEAL::6 7 0.0
+DEAL::6 8 0.0
+DEAL::6 9 0.0
+DEAL::7 0 0.0
+DEAL::7 1 0.0
+DEAL::7 2 0.0
+DEAL::7 3 0.0
+DEAL::7 4 0.0
+DEAL::7 5 0.0
+DEAL::7 6 0.0
+DEAL::7 7 0.0
+DEAL::7 8 0.0
+DEAL::7 9 0.0
+DEAL::8 0 0.0
+DEAL::8 1 0.0
+DEAL::8 2 0.0
+DEAL::8 3 0.0
+DEAL::8 4 0.0
+DEAL::8 5 0.0
+DEAL::8 6 0.0
+DEAL::8 7 0.0
+DEAL::8 8 0.0
+DEAL::8 9 0.0
+DEAL::9 0 0.0
+DEAL::9 1 0.0
+DEAL::9 2 0.0
+DEAL::9 3 0.0
+DEAL::9 4 0.0
+DEAL::9 5 0.0
+DEAL::9 6 0.0
+DEAL::9 7 0.0
+DEAL::9 8 0.0
+DEAL::9 9 0.0
+DEAL::0 0 0.0
+DEAL::0 1 0.0
+DEAL::0 2 0.0
+DEAL::0 3 0.0
+DEAL::0 4 0.0
+DEAL::0 5 0.0
+DEAL::0 6 0.0
+DEAL::0 7 0.0
+DEAL::0 8 0.0
+DEAL::0 9 0.0
+DEAL::1 0 0.0
+DEAL::1 1 0.0
+DEAL::1 2 0.0
+DEAL::1 3 0.0
+DEAL::1 4 0.0
+DEAL::1 5 0.0
+DEAL::1 6 0.0
+DEAL::1 7 0.0
+DEAL::1 8 0.0
+DEAL::1 9 0.0
+DEAL::2 0 0.0
+DEAL::2 1 0.0
+DEAL::2 2 0.0
+DEAL::2 3 0.0
+DEAL::2 4 0.0
+DEAL::2 5 0.0
+DEAL::2 6 0.0
+DEAL::2 7 0.0
+DEAL::2 8 0.0
+DEAL::2 9 0.0
+DEAL::3 0 0.0
+DEAL::3 1 0.0
+DEAL::3 2 0.0
+DEAL::3 3 0.0
DEAL::3 4 30.
DEAL::3 5 40.
DEAL::3 6 50.
-DEAL::3 7 0
-DEAL::3 8 0
-DEAL::3 9 0
-DEAL::4 0 0
-DEAL::4 1 0
-DEAL::4 2 0
-DEAL::4 3 0
-DEAL::4 4 36.
-DEAL::4 5 48.
-DEAL::4 6 60.
-DEAL::4 7 0
-DEAL::4 8 0
-DEAL::4 9 0
-DEAL::5 0 0
-DEAL::5 1 0
-DEAL::5 2 0
-DEAL::5 3 0
-DEAL::5 4 0
-DEAL::5 5 0
-DEAL::5 6 0
-DEAL::5 7 0
-DEAL::5 8 0
-DEAL::5 9 0
-DEAL::6 0 0
-DEAL::6 1 0
-DEAL::6 2 0
-DEAL::6 3 0
-DEAL::6 4 0
-DEAL::6 5 0
-DEAL::6 6 0
-DEAL::6 7 0
-DEAL::6 8 0
-DEAL::6 9 0
-DEAL::7 0 0
-DEAL::7 1 0
-DEAL::7 2 0
-DEAL::7 3 0
-DEAL::7 4 0
-DEAL::7 5 0
-DEAL::7 6 0
-DEAL::7 7 0
-DEAL::7 8 0
-DEAL::7 9 0
-DEAL::8 0 0
-DEAL::8 1 0
-DEAL::8 2 0
-DEAL::8 3 0
-DEAL::8 4 0
-DEAL::8 5 0
-DEAL::8 6 0
-DEAL::8 7 0
-DEAL::8 8 0
-DEAL::8 9 0
-DEAL::9 0 0
-DEAL::9 1 0
-DEAL::9 2 0
-DEAL::9 3 0
-DEAL::9 4 0
-DEAL::9 5 0
-DEAL::9 6 0
-DEAL::9 7 0
-DEAL::9 8 0
-DEAL::9 9 0
+DEAL::3 7 0.0
+DEAL::3 8 0.0
+DEAL::3 9 0.0
+DEAL::4 0 0.0
+DEAL::4 1 0.0
+DEAL::4 2 0.0
+DEAL::4 3 0.0
+DEAL::4 4 0.0
+DEAL::4 5 0.0
+DEAL::4 6 0.0
+DEAL::4 7 0.0
+DEAL::4 8 0.0
+DEAL::4 9 0.0
+DEAL::5 0 0.0
+DEAL::5 1 0.0
+DEAL::5 2 0.0
+DEAL::5 3 0.0
+DEAL::5 4 0.0
+DEAL::5 5 0.0
+DEAL::5 6 0.0
+DEAL::5 7 0.0
+DEAL::5 8 0.0
+DEAL::5 9 0.0
+DEAL::6 0 0.0
+DEAL::6 1 0.0
+DEAL::6 2 0.0
+DEAL::6 3 0.0
+DEAL::6 4 0.0
+DEAL::6 5 0.0
+DEAL::6 6 0.0
+DEAL::6 7 0.0
+DEAL::6 8 0.0
+DEAL::6 9 0.0
+DEAL::7 0 0.0
+DEAL::7 1 0.0
+DEAL::7 2 0.0
+DEAL::7 3 0.0
+DEAL::7 4 0.0
+DEAL::7 5 0.0
+DEAL::7 6 0.0
+DEAL::7 7 0.0
+DEAL::7 8 0.0
+DEAL::7 9 0.0
+DEAL::8 0 0.0
+DEAL::8 1 0.0
+DEAL::8 2 0.0
+DEAL::8 3 0.0
+DEAL::8 4 0.0
+DEAL::8 5 0.0
+DEAL::8 6 0.0
+DEAL::8 7 0.0
+DEAL::8 8 0.0
+DEAL::8 9 0.0
+DEAL::9 0 0.0
+DEAL::9 1 0.0
+DEAL::9 2 0.0
+DEAL::9 3 0.0
+DEAL::9 4 0.0
+DEAL::9 5 0.0
+DEAL::9 6 0.0
+DEAL::9 7 0.0
+DEAL::9 8 0.0
+DEAL::9 9 0.0
{
FullMatrix<number> m;
make_matrix (m);
- m.add_col (2, 3.141, 4);
+ m.add_col (2, 3.141, 3);
print_matrix (m);
}
-
DEAL::0 0 2.0
DEAL::0 1 3.0
-DEAL::0 2 23.
+DEAL::0 2 20.
DEAL::0 3 5.0
-DEAL::0 4 6.0
DEAL::1 0 4.0
DEAL::1 1 6.0
-DEAL::1 2 46.
+DEAL::1 2 39.
DEAL::1 3 10.
-DEAL::1 4 12.
DEAL::2 0 6.0
DEAL::2 1 9.0
-DEAL::2 2 69.
+DEAL::2 2 59.
DEAL::2 3 15.
-DEAL::2 4 18.
DEAL::3 0 8.0
DEAL::3 1 12.
-DEAL::3 2 91.
+DEAL::3 2 79.
DEAL::3 3 20.
-DEAL::3 4 24.
DEAL::4 0 10.
DEAL::4 1 15.
-DEAL::4 2 1.1e+02
+DEAL::4 2 99.
DEAL::4 3 25.
-DEAL::4 4 30.
DEAL::0 0 2.0
DEAL::0 1 3.0
-DEAL::0 2 23.
+DEAL::0 2 20.
DEAL::0 3 5.0
-DEAL::0 4 6.0
DEAL::1 0 4.0
DEAL::1 1 6.0
-DEAL::1 2 46.
+DEAL::1 2 39.
DEAL::1 3 10.
-DEAL::1 4 12.
DEAL::2 0 6.0
DEAL::2 1 9.0
-DEAL::2 2 69.
+DEAL::2 2 59.
DEAL::2 3 15.
-DEAL::2 4 18.
DEAL::3 0 8.0
DEAL::3 1 12.
-DEAL::3 2 91.
+DEAL::3 2 79.
DEAL::3 3 20.
-DEAL::3 4 24.
DEAL::4 0 10.
DEAL::4 1 15.
-DEAL::4 2 1.1e+02
+DEAL::4 2 99.
DEAL::4 3 25.
-DEAL::4 4 30.
{
FullMatrix<number> m;
make_matrix (m);
- m.add_col (2, 3.141, 4, 2.718, 3);
+ m.add_col (2, 3.141, 1, 2.718, 3);
print_matrix (m);
}
-
DEAL::0 0 2.0
DEAL::0 1 3.0
-DEAL::0 2 36.
+DEAL::0 2 27.
DEAL::0 3 5.0
-DEAL::0 4 6.0
DEAL::1 0 4.0
DEAL::1 1 6.0
-DEAL::1 2 73.
+DEAL::1 2 54.
DEAL::1 3 10.
-DEAL::1 4 12.
DEAL::2 0 6.0
DEAL::2 1 9.0
-DEAL::2 2 1.1e+02
+DEAL::2 2 81.
DEAL::2 3 15.
-DEAL::2 4 18.
DEAL::3 0 8.0
DEAL::3 1 12.
-DEAL::3 2 1.5e+02
+DEAL::3 2 1.1e+02
DEAL::3 3 20.
-DEAL::3 4 24.
DEAL::4 0 10.
DEAL::4 1 15.
-DEAL::4 2 1.8e+02
+DEAL::4 2 1.4e+02
DEAL::4 3 25.
-DEAL::4 4 30.
DEAL::0 0 2.0
DEAL::0 1 3.0
-DEAL::0 2 36.
+DEAL::0 2 27.
DEAL::0 3 5.0
-DEAL::0 4 6.0
DEAL::1 0 4.0
DEAL::1 1 6.0
-DEAL::1 2 73.
+DEAL::1 2 54.
DEAL::1 3 10.
-DEAL::1 4 12.
DEAL::2 0 6.0
DEAL::2 1 9.0
-DEAL::2 2 1.1e+02
+DEAL::2 2 81.
DEAL::2 3 15.
-DEAL::2 4 18.
DEAL::3 0 8.0
DEAL::3 1 12.
-DEAL::3 2 1.5e+02
+DEAL::3 2 1.1e+02
DEAL::3 3 20.
-DEAL::3 4 24.
DEAL::4 0 10.
DEAL::4 1 15.
-DEAL::4 2 1.8e+02
+DEAL::4 2 1.4e+02
DEAL::4 3 25.
-DEAL::4 4 30.
{
FullMatrix<number> m;
make_matrix (m);
- m.swap_col (2, 4);
+ m.swap_col (2, 3);
print_matrix (m);
}
-
DEAL::0 0 2.0
DEAL::0 1 3.0
-DEAL::0 2 6.0
-DEAL::0 3 5.0
-DEAL::0 4 4.0
+DEAL::0 2 5.0
+DEAL::0 3 4.0
DEAL::1 0 4.0
DEAL::1 1 6.0
-DEAL::1 2 12.
-DEAL::1 3 10.
-DEAL::1 4 8.0
+DEAL::1 2 10.
+DEAL::1 3 8.0
DEAL::2 0 6.0
DEAL::2 1 9.0
-DEAL::2 2 18.
-DEAL::2 3 15.
-DEAL::2 4 12.
+DEAL::2 2 15.
+DEAL::2 3 12.
DEAL::3 0 8.0
DEAL::3 1 12.
-DEAL::3 2 24.
-DEAL::3 3 20.
-DEAL::3 4 16.
+DEAL::3 2 20.
+DEAL::3 3 16.
DEAL::4 0 10.
DEAL::4 1 15.
-DEAL::4 2 30.
-DEAL::4 3 25.
-DEAL::4 4 20.
+DEAL::4 2 25.
+DEAL::4 3 20.
DEAL::0 0 2.0
DEAL::0 1 3.0
-DEAL::0 2 6.0
-DEAL::0 3 5.0
-DEAL::0 4 4.0
+DEAL::0 2 5.0
+DEAL::0 3 4.0
DEAL::1 0 4.0
DEAL::1 1 6.0
-DEAL::1 2 12.
-DEAL::1 3 10.
-DEAL::1 4 8.0
+DEAL::1 2 10.
+DEAL::1 3 8.0
DEAL::2 0 6.0
DEAL::2 1 9.0
-DEAL::2 2 18.
-DEAL::2 3 15.
-DEAL::2 4 12.
+DEAL::2 2 15.
+DEAL::2 3 12.
DEAL::3 0 8.0
DEAL::3 1 12.
-DEAL::3 2 24.
-DEAL::3 3 20.
-DEAL::3 4 16.
+DEAL::3 2 20.
+DEAL::3 3 16.
DEAL::4 0 10.
DEAL::4 1 15.
-DEAL::4 2 30.
-DEAL::4 3 25.
-DEAL::4 4 20.
+DEAL::4 2 25.
+DEAL::4 3 20.
FullMatrix<number> m;
make_matrix (m);
Vector<number> v, w;
- make_vector (v);
- make_vector (w);
+ make_range_vector (v);
+ make_domain_vector (w);
m.vmult (v, w, true);
print_vector (v);
print_vector (v);
print_vector (w);
}
-
-DEAL::0 71.
-DEAL::1 1.4e+02
-DEAL::2 2.1e+02
-DEAL::3 2.8e+02
-DEAL::4 3.6e+02
+DEAL::0 41.
+DEAL::1 82.
+DEAL::2 1.2e+02
+DEAL::3 1.6e+02
+DEAL::4 2.0e+02
DEAL::0 1.0
DEAL::1 2.0
DEAL::2 3.0
DEAL::3 4.0
-DEAL::4 5.0
-DEAL::0 70.
-DEAL::1 1.4e+02
-DEAL::2 2.1e+02
-DEAL::3 2.8e+02
-DEAL::4 3.5e+02
+DEAL::0 40.
+DEAL::1 80.
+DEAL::2 1.2e+02
+DEAL::3 1.6e+02
+DEAL::4 2.0e+02
DEAL::0 1.0
DEAL::1 2.0
DEAL::2 3.0
DEAL::3 4.0
-DEAL::4 5.0
-DEAL::0 71.
-DEAL::1 1.4e+02
-DEAL::2 2.1e+02
-DEAL::3 2.8e+02
-DEAL::4 3.6e+02
+DEAL::0 41.
+DEAL::1 82.
+DEAL::2 1.2e+02
+DEAL::3 1.6e+02
+DEAL::4 2.0e+02
DEAL::0 1.0
DEAL::1 2.0
DEAL::2 3.0
DEAL::3 4.0
-DEAL::4 5.0
-DEAL::0 70.
-DEAL::1 1.4e+02
-DEAL::2 2.1e+02
-DEAL::3 2.8e+02
-DEAL::4 3.5e+02
+DEAL::0 40.
+DEAL::1 80.
+DEAL::2 1.2e+02
+DEAL::3 1.6e+02
+DEAL::4 2.0e+02
DEAL::0 1.0
DEAL::1 2.0
DEAL::2 3.0
DEAL::3 4.0
-DEAL::4 5.0
FullMatrix<number> m;
make_matrix (m);
Vector<number> v, w;
- make_vector (v);
- make_vector (w);
+ make_range_vector (v);
+ make_domain_vector (w);
m.vmult_add (v, w);
print_vector (v);
print_vector (w);
}
-
-DEAL::0 71.
-DEAL::1 1.4e+02
-DEAL::2 2.1e+02
-DEAL::3 2.8e+02
-DEAL::4 3.6e+02
+DEAL::0 41.
+DEAL::1 82.
+DEAL::2 1.2e+02
+DEAL::3 1.6e+02
+DEAL::4 2.0e+02
DEAL::0 1.0
DEAL::1 2.0
DEAL::2 3.0
DEAL::3 4.0
-DEAL::4 5.0
-DEAL::0 71.
-DEAL::1 1.4e+02
-DEAL::2 2.1e+02
-DEAL::3 2.8e+02
-DEAL::4 3.6e+02
+DEAL::0 41.
+DEAL::1 82.
+DEAL::2 1.2e+02
+DEAL::3 1.6e+02
+DEAL::4 2.0e+02
DEAL::0 1.0
DEAL::1 2.0
DEAL::2 3.0
DEAL::3 4.0
-DEAL::4 5.0
FullMatrix<number> m;
make_matrix (m);
Vector<number> v, w, x;
- make_vector (v);
- make_vector (w);
- make_vector (x);
+ make_range_vector (v);
+ make_domain_vector (w);
+ make_range_vector (x);
m.residual (v, w, x);
print_vector (v);
print_vector (w);
print_vector (x);
}
-
-DEAL::0 -69.
-DEAL::1 -1.4e+02
-DEAL::2 -2.1e+02
-DEAL::3 -2.8e+02
-DEAL::4 -3.4e+02
+DEAL::0 -39.
+DEAL::1 -78.
+DEAL::2 -1.2e+02
+DEAL::3 -1.6e+02
+DEAL::4 -2.0e+02
DEAL::0 1.0
DEAL::1 2.0
DEAL::2 3.0
DEAL::3 4.0
-DEAL::4 5.0
DEAL::0 1.0
DEAL::1 2.0
DEAL::2 3.0
DEAL::3 4.0
DEAL::4 5.0
-DEAL::0 -69.
-DEAL::1 -1.4e+02
-DEAL::2 -2.1e+02
-DEAL::3 -2.8e+02
-DEAL::4 -3.4e+02
+DEAL::0 -39.
+DEAL::1 -78.
+DEAL::2 -1.2e+02
+DEAL::3 -1.6e+02
+DEAL::4 -2.0e+02
DEAL::0 1.0
DEAL::1 2.0
DEAL::2 3.0
DEAL::3 4.0
-DEAL::4 5.0
DEAL::0 1.0
DEAL::1 2.0
DEAL::2 3.0
-DEAL::0 0 0
-DEAL::0 1 0
-DEAL::0 2 0
-DEAL::0 3 0
-DEAL::0 4 0
-DEAL::0 5 0
-DEAL::0 6 0
-DEAL::0 7 0
-DEAL::0 8 0
-DEAL::0 9 0
-DEAL::1 0 0
-DEAL::1 1 0
-DEAL::1 2 0
-DEAL::1 3 0
-DEAL::1 4 0
-DEAL::1 5 0
-DEAL::1 6 0
-DEAL::1 7 0
-DEAL::1 8 0
-DEAL::1 9 0
-DEAL::2 0 0
-DEAL::2 1 0
-DEAL::2 2 0
-DEAL::2 3 0
-DEAL::2 4 0
-DEAL::2 5 0
-DEAL::2 6 0
-DEAL::2 7 0
-DEAL::2 8 0
-DEAL::2 9 0
-DEAL::3 0 0
-DEAL::3 1 0
-DEAL::3 2 0
-DEAL::3 3 0
+DEAL::0 0 0.0
+DEAL::0 1 0.0
+DEAL::0 2 0.0
+DEAL::0 3 0.0
+DEAL::0 4 0.0
+DEAL::0 5 0.0
+DEAL::0 6 0.0
+DEAL::0 7 0.0
+DEAL::0 8 0.0
+DEAL::0 9 0.0
+DEAL::1 0 0.0
+DEAL::1 1 0.0
+DEAL::1 2 0.0
+DEAL::1 3 0.0
+DEAL::1 4 0.0
+DEAL::1 5 0.0
+DEAL::1 6 0.0
+DEAL::1 7 0.0
+DEAL::1 8 0.0
+DEAL::1 9 0.0
+DEAL::2 0 0.0
+DEAL::2 1 0.0
+DEAL::2 2 0.0
+DEAL::2 3 0.0
+DEAL::2 4 0.0
+DEAL::2 5 0.0
+DEAL::2 6 0.0
+DEAL::2 7 0.0
+DEAL::2 8 0.0
+DEAL::2 9 0.0
+DEAL::3 0 0.0
+DEAL::3 1 0.0
+DEAL::3 2 0.0
+DEAL::3 3 0.0
DEAL::3 4 2.0
DEAL::3 5 3.0
DEAL::3 6 4.0
DEAL::3 7 5.0
-DEAL::3 8 6.0
-DEAL::3 9 0
-DEAL::4 0 0
-DEAL::4 1 0
-DEAL::4 2 0
-DEAL::4 3 0
+DEAL::3 8 0.0
+DEAL::3 9 0.0
+DEAL::4 0 0.0
+DEAL::4 1 0.0
+DEAL::4 2 0.0
+DEAL::4 3 0.0
DEAL::4 4 4.0
DEAL::4 5 6.0
DEAL::4 6 8.0
DEAL::4 7 10.
-DEAL::4 8 12.
-DEAL::4 9 0
-DEAL::5 0 0
-DEAL::5 1 0
-DEAL::5 2 0
-DEAL::5 3 0
+DEAL::4 8 0.0
+DEAL::4 9 0.0
+DEAL::5 0 0.0
+DEAL::5 1 0.0
+DEAL::5 2 0.0
+DEAL::5 3 0.0
DEAL::5 4 6.0
DEAL::5 5 9.0
DEAL::5 6 12.
DEAL::5 7 15.
-DEAL::5 8 18.
-DEAL::5 9 0
-DEAL::6 0 0
-DEAL::6 1 0
-DEAL::6 2 0
-DEAL::6 3 0
+DEAL::5 8 0.0
+DEAL::5 9 0.0
+DEAL::6 0 0.0
+DEAL::6 1 0.0
+DEAL::6 2 0.0
+DEAL::6 3 0.0
DEAL::6 4 8.0
DEAL::6 5 12.
DEAL::6 6 16.
DEAL::6 7 20.
-DEAL::6 8 24.
-DEAL::6 9 0
-DEAL::7 0 0
-DEAL::7 1 0
-DEAL::7 2 0
-DEAL::7 3 0
+DEAL::6 8 0.0
+DEAL::6 9 0.0
+DEAL::7 0 0.0
+DEAL::7 1 0.0
+DEAL::7 2 0.0
+DEAL::7 3 0.0
DEAL::7 4 10.
DEAL::7 5 15.
DEAL::7 6 20.
DEAL::7 7 25.
-DEAL::7 8 30.
-DEAL::7 9 0
-DEAL::8 0 0
-DEAL::8 1 0
-DEAL::8 2 0
-DEAL::8 3 0
-DEAL::8 4 0
-DEAL::8 5 0
-DEAL::8 6 0
-DEAL::8 7 0
-DEAL::8 8 0
-DEAL::8 9 0
-DEAL::9 0 0
-DEAL::9 1 0
-DEAL::9 2 0
-DEAL::9 3 0
-DEAL::9 4 0
-DEAL::9 5 0
-DEAL::9 6 0
-DEAL::9 7 0
-DEAL::9 8 0
-DEAL::9 9 0
-DEAL::0 0 0
-DEAL::0 1 0
-DEAL::0 2 0
-DEAL::0 3 0
-DEAL::0 4 0
-DEAL::0 5 0
-DEAL::0 6 0
-DEAL::0 7 0
-DEAL::0 8 0
-DEAL::0 9 0
-DEAL::1 0 0
-DEAL::1 1 0
-DEAL::1 2 0
-DEAL::1 3 0
-DEAL::1 4 0
-DEAL::1 5 0
-DEAL::1 6 0
-DEAL::1 7 0
-DEAL::1 8 0
-DEAL::1 9 0
-DEAL::2 0 0
-DEAL::2 1 0
-DEAL::2 2 0
-DEAL::2 3 0
-DEAL::2 4 0
-DEAL::2 5 0
-DEAL::2 6 0
-DEAL::2 7 0
-DEAL::2 8 0
-DEAL::2 9 0
-DEAL::3 0 0
-DEAL::3 1 0
-DEAL::3 2 0
-DEAL::3 3 0
+DEAL::7 8 0.0
+DEAL::7 9 0.0
+DEAL::8 0 0.0
+DEAL::8 1 0.0
+DEAL::8 2 0.0
+DEAL::8 3 0.0
+DEAL::8 4 0.0
+DEAL::8 5 0.0
+DEAL::8 6 0.0
+DEAL::8 7 0.0
+DEAL::8 8 0.0
+DEAL::8 9 0.0
+DEAL::9 0 0.0
+DEAL::9 1 0.0
+DEAL::9 2 0.0
+DEAL::9 3 0.0
+DEAL::9 4 0.0
+DEAL::9 5 0.0
+DEAL::9 6 0.0
+DEAL::9 7 0.0
+DEAL::9 8 0.0
+DEAL::9 9 0.0
+DEAL::0 0 0.0
+DEAL::0 1 0.0
+DEAL::0 2 0.0
+DEAL::0 3 0.0
+DEAL::0 4 0.0
+DEAL::0 5 0.0
+DEAL::0 6 0.0
+DEAL::0 7 0.0
+DEAL::0 8 0.0
+DEAL::0 9 0.0
+DEAL::1 0 0.0
+DEAL::1 1 0.0
+DEAL::1 2 0.0
+DEAL::1 3 0.0
+DEAL::1 4 0.0
+DEAL::1 5 0.0
+DEAL::1 6 0.0
+DEAL::1 7 0.0
+DEAL::1 8 0.0
+DEAL::1 9 0.0
+DEAL::2 0 0.0
+DEAL::2 1 0.0
+DEAL::2 2 0.0
+DEAL::2 3 0.0
+DEAL::2 4 0.0
+DEAL::2 5 0.0
+DEAL::2 6 0.0
+DEAL::2 7 0.0
+DEAL::2 8 0.0
+DEAL::2 9 0.0
+DEAL::3 0 0.0
+DEAL::3 1 0.0
+DEAL::3 2 0.0
+DEAL::3 3 0.0
DEAL::3 4 2.0
DEAL::3 5 3.0
DEAL::3 6 4.0
DEAL::3 7 5.0
-DEAL::3 8 6.0
-DEAL::3 9 0
-DEAL::4 0 0
-DEAL::4 1 0
-DEAL::4 2 0
-DEAL::4 3 0
+DEAL::3 8 0.0
+DEAL::3 9 0.0
+DEAL::4 0 0.0
+DEAL::4 1 0.0
+DEAL::4 2 0.0
+DEAL::4 3 0.0
DEAL::4 4 4.0
DEAL::4 5 6.0
DEAL::4 6 8.0
DEAL::4 7 10.
-DEAL::4 8 12.
-DEAL::4 9 0
-DEAL::5 0 0
-DEAL::5 1 0
-DEAL::5 2 0
-DEAL::5 3 0
+DEAL::4 8 0.0
+DEAL::4 9 0.0
+DEAL::5 0 0.0
+DEAL::5 1 0.0
+DEAL::5 2 0.0
+DEAL::5 3 0.0
DEAL::5 4 6.0
DEAL::5 5 9.0
DEAL::5 6 12.
DEAL::5 7 15.
-DEAL::5 8 18.
-DEAL::5 9 0
-DEAL::6 0 0
-DEAL::6 1 0
-DEAL::6 2 0
-DEAL::6 3 0
+DEAL::5 8 0.0
+DEAL::5 9 0.0
+DEAL::6 0 0.0
+DEAL::6 1 0.0
+DEAL::6 2 0.0
+DEAL::6 3 0.0
DEAL::6 4 8.0
DEAL::6 5 12.
DEAL::6 6 16.
DEAL::6 7 20.
-DEAL::6 8 24.
-DEAL::6 9 0
-DEAL::7 0 0
-DEAL::7 1 0
-DEAL::7 2 0
-DEAL::7 3 0
+DEAL::6 8 0.0
+DEAL::6 9 0.0
+DEAL::7 0 0.0
+DEAL::7 1 0.0
+DEAL::7 2 0.0
+DEAL::7 3 0.0
DEAL::7 4 10.
DEAL::7 5 15.
DEAL::7 6 20.
DEAL::7 7 25.
-DEAL::7 8 30.
-DEAL::7 9 0
-DEAL::8 0 0
-DEAL::8 1 0
-DEAL::8 2 0
-DEAL::8 3 0
-DEAL::8 4 0
-DEAL::8 5 0
-DEAL::8 6 0
-DEAL::8 7 0
-DEAL::8 8 0
-DEAL::8 9 0
-DEAL::9 0 0
-DEAL::9 1 0
-DEAL::9 2 0
-DEAL::9 3 0
-DEAL::9 4 0
-DEAL::9 5 0
-DEAL::9 6 0
-DEAL::9 7 0
-DEAL::9 8 0
-DEAL::9 9 0
+DEAL::7 8 0.0
+DEAL::7 9 0.0
+DEAL::8 0 0.0
+DEAL::8 1 0.0
+DEAL::8 2 0.0
+DEAL::8 3 0.0
+DEAL::8 4 0.0
+DEAL::8 5 0.0
+DEAL::8 6 0.0
+DEAL::8 7 0.0
+DEAL::8 8 0.0
+DEAL::8 9 0.0
+DEAL::9 0 0.0
+DEAL::9 1 0.0
+DEAL::9 2 0.0
+DEAL::9 3 0.0
+DEAL::9 4 0.0
+DEAL::9 5 0.0
+DEAL::9 6 0.0
+DEAL::9 7 0.0
+DEAL::9 8 0.0
+DEAL::9 9 0.0