Exemple #1
0
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;
}
Exemple #2
0
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;
}