/* //////////////////////////////////////////////////////////////////////////// -- testing any solver */ int main( int argc, char** argv ) { magma_int_t info = 0; TESTING_INIT(); magma_sopts zopts; magma_queue_t queue=NULL; magma_queue_create( /*devices[ opts->device ],*/ &queue ); float one = MAGMA_S_MAKE(1.0, 0.0); float zero = MAGMA_S_MAKE(0.0, 0.0); magma_s_matrix A={Magma_CSR}, B={Magma_CSR}, B_d={Magma_CSR}; magma_s_matrix x={Magma_CSR}, b={Magma_CSR}; int i=1; CHECK( magma_sparse_opts( argc, argv, &zopts, &i, queue )); B.blocksize = zopts.blocksize; B.alignment = zopts.alignment; if ( zopts.solver_par.solver != Magma_PCG && zopts.solver_par.solver != Magma_PGMRES && zopts.solver_par.solver != Magma_PBICGSTAB && zopts.solver_par.solver != Magma_ITERREF && zopts.solver_par.solver != Magma_LOBPCG ) zopts.precond_par.solver = Magma_NONE; CHECK( magma_ssolverinfo_init( &zopts.solver_par, &zopts.precond_par, queue )); while( i < argc ) { if ( strcmp("LAPLACE2D", argv[i]) == 0 && i+1 < argc ) { // Laplace test i++; magma_int_t laplace_size = atoi( argv[i] ); CHECK( magma_sm_5stencil( laplace_size, &A, queue )); } else { // file-matrix test CHECK( magma_s_csr_mtx( &A, argv[i], queue )); } printf( "\n# matrix info: %d-by-%d with %d nonzeros\n\n", (int) A.num_rows,(int) A.num_cols,(int) A.nnz ); // for the eigensolver case zopts.solver_par.ev_length = A.num_rows; CHECK( magma_seigensolverinfo_init( &zopts.solver_par, queue )); // scale matrix CHECK( magma_smscale( &A, zopts.scaling, queue )); CHECK( magma_smconvert( A, &B, Magma_CSR, zopts.output_format, queue )); CHECK( magma_smtransfer( B, &B_d, Magma_CPU, Magma_DEV, queue )); // vectors and initial guess CHECK( magma_svinit( &b, Magma_DEV, A.num_cols, 1, one, queue )); //magma_svinit( &x, Magma_DEV, A.num_cols, 1, one, queue ); //magma_s_spmv( one, B_d, x, zero, b, queue ); // b = A x //magma_smfree(&x, queue ); CHECK( magma_svinit( &x, Magma_DEV, A.num_cols, 1, zero, queue )); info = magma_s_solver( B_d, b, &x, &zopts, queue ); if( info != 0 ){ printf("error: solver returned: %s (%d).\n", magma_strerror( info ), info ); } magma_ssolverinfo( &zopts.solver_par, &zopts.precond_par, queue ); magma_smfree(&B_d, queue ); magma_smfree(&B, queue ); magma_smfree(&A, queue ); magma_smfree(&x, queue ); magma_smfree(&b, queue ); i++; } cleanup: magma_smfree(&B_d, queue ); magma_smfree(&B, queue ); magma_smfree(&A, queue ); magma_smfree(&x, queue ); magma_smfree(&b, queue ); magma_ssolverinfo_free( &zopts.solver_par, &zopts.precond_par, queue ); magma_queue_destroy( queue ); TESTING_FINALIZE(); return info; }
/* //////////////////////////////////////////////////////////////////////////// -- testing any solver */ int main( int argc, char** argv ) { magma_int_t info = 0; TESTING_INIT(); magma_sopts zopts; magma_queue_t queue=NULL; magma_queue_create( &queue ); float one = MAGMA_S_MAKE(1.0, 0.0); float zero = MAGMA_S_MAKE(0.0, 0.0); magma_s_matrix A={Magma_CSR}, B={Magma_CSR}, B_d={Magma_CSR}; magma_s_matrix x={Magma_CSR}, b={Magma_CSR}; int i=1; CHECK( magma_sparse_opts( argc, argv, &zopts, &i, queue )); B.blocksize = zopts.blocksize; B.alignment = zopts.alignment; CHECK( magma_ssolverinfo_init( &zopts.solver_par, &zopts.precond_par, queue )); while( i < argc ) { if ( strcmp("LAPLACE2D", argv[i]) == 0 && i+1 < argc ) { // Laplace test i++; magma_int_t laplace_size = atoi( argv[i] ); CHECK( magma_sm_5stencil( laplace_size, &A, queue )); } else { // file-matrix test CHECK( magma_s_csr_mtx( &A, argv[i], queue )); } // for the eigensolver case zopts.solver_par.ev_length = A.num_cols; CHECK( magma_seigensolverinfo_init( &zopts.solver_par, queue )); // scale matrix CHECK( magma_smscale( &A, zopts.scaling, queue )); // preconditioner if ( zopts.solver_par.solver != Magma_ITERREF ) { CHECK( magma_s_precondsetup( A, b, &zopts.solver_par, &zopts.precond_par, queue ) ); } CHECK( magma_smconvert( A, &B, Magma_CSR, zopts.output_format, queue )); printf( "\n%% matrix info: %d-by-%d with %d nonzeros\n\n", int(A.num_rows), int(A.num_cols), int(A.nnz) ); printf("matrixinfo = [ \n"); printf("%% size (m x n) || nonzeros (nnz) || nnz/m || stored nnz\n"); printf("%%======================================================================" "======%%\n"); printf(" %8d %8d %10d %4d %10d\n", int(B.num_rows), int(B.num_cols), int(B.true_nnz), int(B.true_nnz/B.num_rows), int(B.nnz) ); printf("%%======================================================================" "======%%\n"); printf("];\n"); CHECK( magma_smtransfer( B, &B_d, Magma_CPU, Magma_DEV, queue )); // vectors and initial guess CHECK( magma_svinit( &b, Magma_DEV, A.num_rows, 1, one, queue )); //magma_svinit( &x, Magma_DEV, A.num_cols, 1, one, queue ); //magma_s_spmv( one, B_d, x, zero, b, queue ); // b = A x //magma_smfree(&x, queue ); CHECK( magma_svinit( &x, Magma_DEV, A.num_cols, 1, zero, queue )); info = magma_s_solver( B_d, b, &x, &zopts, queue ); if( info != 0 ) { printf("%%error: solver returned: %s (%d).\n", magma_strerror( info ), int(info) ); } printf("data = [\n"); magma_ssolverinfo( &zopts.solver_par, &zopts.precond_par, queue ); printf("];\n\n"); printf("precond_info = [\n"); printf("%% setup runtime\n"); printf(" %.6f %.6f\n", zopts.precond_par.setuptime, zopts.precond_par.runtime ); printf("];\n\n"); magma_smfree(&B_d, queue ); magma_smfree(&B, queue ); magma_smfree(&A, queue ); magma_smfree(&x, queue ); magma_smfree(&b, queue ); i++; } cleanup: magma_smfree(&B_d, queue ); magma_smfree(&B, queue ); magma_smfree(&A, queue ); magma_smfree(&x, queue ); magma_smfree(&b, queue ); magma_ssolverinfo_free( &zopts.solver_par, &zopts.precond_par, queue ); magma_queue_destroy( queue ); TESTING_FINALIZE(); return info; }