int main(int argc, char *argv[]) { SimpleMatrix m; assert (m.nrow() == 0); assert (m.ncol() == 0); m.resize(3, 2); assert (m.nrow() == 3); assert (m.ncol() == 2); m.zero(); assert (m[0][0] == 0.0); assert (m[2][1] == 0.0); std::vector<double> col(3, 1.0); m.appendCol(col); assert (m[0][0] == 0.0); assert (m[2][2] == 1.0); m.resize(2,3); std::vector<double> row(3, 2.0); m.appendRow(row); assert (m[0][0] == 0.0); assert (m[2][2] == 2.0); return 0; }
void NonlinearRelationWithSignInversed::computeJachlambda(double t, SiconosVector& x, SiconosVector& lambda, SimpleMatrix& D) { // double *h = &(*_jachlambda)(0,0); #ifdef SICONOS_DEBUG std::cout << "NonlinearRelationWithSignInversed::computeJachlambda " << " at " << " " << t << std::endl; #endif D.zero(); }
void NonlinearRelationReduced2::computeJachlambda(double t, SiconosVector& x, SiconosVector& lambda, SimpleMatrix& D) { D.zero(); }