int main()
{
double ar[SIZE_ar] = { 0, 1, 2, 3, 4, 5 };
double ar2[SIZE_ar2] = { 8,7,6,5,4,3,2,1,0};
DoubleVector dv1(SIZE_ar, ar);
DoubleVector dv2(SIZE_ar2, ar2);
DoubleVector dv3;
dv3 = dv1;
double p = dv1[1];
DoubleVector v;
v = dv2 + dv1 + dv1 - dv2 - dv1;
v.print();
dv1 += dv3 += dv2 -= dv2;
dv1.print();
dv1 = dv3 * 2 * 4 ;
dv1 *= 2;
dv3 = dv1 / 2;
dv3 /= 2;
double scalar = dv3% dv1;
BoolVector a(2);
cout<<"number in binary ...00000010 --> in array pos[31] bit is "<<a[31]
<<endl<<"in array pos[3] - bit is "<<a[30]<<endl;
a.set(30);
cout << "Nothing changes because in pos 30 it has 1. The number is the same (2) :" << a.get_number()<<endl;
a.set(31);
cout << "Change on pos 31 from 0 to 1. The new number is 3:" << a.get_number()<<endl;
BoolVector b(55);
b.set(28);
cout << "We are working with this number 00..110111 in binary (decimal 55)."<<endl
<<"We want to make it 00.. 111111 exacly 63 in decimal : " << b.get_number() << endl;
b.get(28);
cout << "We are working with 00.. 111111 decimal 63 changing into 00..110111 decimal 55:" << b.get_number() << endl;
}
int main()
{
DoubleVector dv(5, 0), dv1(10, 1);
cout << dv.size() << endl;
dv.print();
dv1[5] = 5.5;
dv1.print();
dv = dv1;
cout << dv.min() << "\t" << dv.max() << "\t" << dv.norm_2() << endl;
dv.print();
cout << "Hello World!" << endl;
DoubleMatrix a(2, 0.0);
a(0, 0) = 1.0;
a(1, 1) = 2.5;
a[0][1] = 0.3;
cout << "a:" << endl;
a.print();
cout << "|a| = " << a.det2x2() << endl;
DoubleMatrix e = a * a.inverted2x2();
cout << "e:" << endl;
e.print();
DoubleMatrix b(2, 5, 1.0);
cout << "b:" << endl;
b.print();
DoubleMatrix m(a * b);
cout << "a * b" << endl;
m.print();
b = a;
cout << "b = a :" << endl;
b.print();
DoubleMatrix c(b.transpose());
cout << "c = b.traspose():" << endl;
c.print();
DoubleMatrix ccc(3, 0.0);
ccc(0, 1) = 0.8;
ccc(1, 0) = 0.2; ccc(1, 2) = 0.4;
ccc(2, 1) = 0.5;
ccc(0, 2) = 1.0;
cout << "ccc:" << endl;
ccc.print();
cout << "|ccc| = " << ccc.det3x3() << endl;
DoubleMatrix iccc = ccc.inverted();
cout << "ccc^(-1):" << endl;
iccc.print();
DoubleMatrix eee = ccc * iccc;
cout << "eee = ccc * ccc^(-1) :" << endl;
eee.print();
cout << "sparce" << endl;
MappedDoubleMatrix mdm(3);
// mdm(0, 1) = 0.5;
// mdm(1, 0) = 0.1; mdm(1, 2) = 0.3;
// mdm(2, 1) = 0.7;
mdm(0, 0) = 1.0; mdm(0, 1) = 1.0;
mdm(1, 0) = 1.0; mdm(1, 1) = 1.0; mdm(1, 2) = 2.0;
mdm(2, 1) = 2.0; mdm(2, 2) = 1.0;
mdm.print();
DoubleVector v3(3, 1.0);
v3[0] += 1.0;
cout << "v3: ";
v3.print();
DoubleVector ccc_v3(mdm * v3);
cout << "mdm * v3: ";
ccc_v3.print();
cout << "v * ccc_v3 = " << v3 * ccc_v3 << endl;
DoubleVector X = mdm.conjugateGradient(v3);
X.print();
cout << "RowDoubleMatrix" << endl;
RowDoubleMatrix rdm(mdm);
X = rdm.conjugateGradient(v3);
X.print();
X = rdm.preconditionedConjugateGradient(v3);
X.print();
cout << " !!! " << endl;
// MappedDoubleMatrix tst(2);
// tst(0,0) = 4.0; tst(0, 1) = 1.0;
// tst(1,0) = 1.0; tst(1, 1) = 3.0;
// cout << "test system: " << endl;
// tst.print();
// DoubleVector tsv(2);
// tsv[0] = 1.0; tsv[1] = 2.0;
// cout << endl;
// tsv.print();
// DoubleVector X = tst.conjugateGradient(tsv);
// X.print();
MappedDoubleMatrix chol(3);
chol(0, 0) = 4.0; chol(0, 1) = 12.0; chol(0, 2) = -16.0;
chol(1, 0)= 12.0; chol(1, 1) = 37.0; chol(1, 2) = -43.0;
chol(2,0) = -16.0; chol(2, 1) = -43.0; chol(2, 2) = 98.0;
chol.print();
DoubleVector bc(3, 0.0);
bc[0] = 0.0; bc[1] = 1.0; bc[2] = 2.0;
DoubleVector cxx = chol.cholesky(bc);
cxx.print();
return 0;
}