/*! main */
int main(int argc, char** argv)
{
srand(time(NULL));
int N(5);
CPPL::dssmatrix A(N);
A(0,0)=1;
A(2,3)=2;
A(1,2)=3;
A(4,1)=4;
//A(3,2)+=0.1;
std::cout << "A =\n" << A << std::endl;
std::cout << "A(1,1)=" << A(1,1) << std::endl;
for(std::vector<CPPL::dcomponent>::const_iterator it=A.data.begin(); it!=A.data.end(); it++){
std::cout << "A(" << it->i << "," << it->j << ") =" << it->v << std::endl;
}
A.checkup();
A.put(3,4, 3.4);
A.del(1,4);
A.del(0);
std::cout << "A =\n" << A << std::endl;
const CPPL::dssmatrix B(A);
//// write/read ////
B.write( "tmp.txt" );
CPPL::dssmatrix C;
C.read( "tmp.txt" );
std::cout << "C-B =\n" << C-B << "<-Should be zero." << std::endl;
return 0;
}
/*! main */
int main(int argc, char** argv)
{
srand(time(NULL));
int M(5), N(3), CAP(4);
CPPL::dssmatrix A(M,N,CAP);
A.fput(0,0, 1.);
A.fput(3,2, 2.);
A.fput(1,2, 3.);
A.fput(4,1, 4.);
cout << "A =\n" << A << endl;
for(int c=0; c<A.vol; c++){
cout << "A(" << A.indx[c] << "," << A.jndx[c] << ") ="
<< A.array[c] << endl;
}
//A.put(1,2, 4.5);
//A.add(1,2, 0.1);
//A.sub(1,2, 0.1);
//A.mult(1,2, 10.);
//A.div(1,2, 10.);
//A.del(1,2);
A.fdel(2);
cout << "A =\n" << A << endl;
const CPPL::dssmatrix B(A);
//// write/read ////
B.write( "tmp.txt" );
CPPL::dssmatrix C;
C.read( "tmp.txt" );
cout << "C-B =\n" << C-B << "<-Should be zero." << endl;
return 0;
}