extern "C" magma_int_t magma_z_precond( magma_z_sparse_matrix A, magma_z_vector b, magma_z_vector *x, magma_z_preconditioner *precond, magma_queue_t queue ) { // set up precond parameters as solver parameters magma_z_solver_par psolver_par; psolver_par.epsilon = precond->epsilon; psolver_par.maxiter = precond->maxiter; psolver_par.restart = precond->restart; psolver_par.verbose = 0; switch( precond->solver ) { case Magma_CG: magma_zcg_res( A, b, x, &psolver_par, queue );break; case Magma_BICGSTAB: magma_zbicgstab( A, b, x, &psolver_par, queue );break; case Magma_GMRES: magma_zgmres( A, b, x, &psolver_par, queue );break; case Magma_JACOBI: magma_zjacobi( A, b, x, &psolver_par, queue );break; case Magma_BAITER: magma_zbaiter( A, b, x, &psolver_par, queue );break; } return MAGMA_SUCCESS; }
/* //////////////////////////////////////////////////////////////////////////// -- running magma_zgmres */ int main( int argc, char** argv) { TESTING_INIT(); magma_z_solver_par solver_par; magma_z_preconditioner precond_par; solver_par.epsilon = 10e-16; solver_par.maxiter = 1000; solver_par.restart = 30; solver_par.num_eigenvalues = 0; solver_par.ortho = Magma_CGS; solver_par.verbose = 0; int format = 0; int ortho = 0; int scale = 0; magma_scale_t scaling = Magma_NOSCALE; magma_z_sparse_matrix A, B, B_d; magma_z_vector x, b; B.blocksize = 8; B.alignment = 8; magmaDoubleComplex one = MAGMA_Z_MAKE(1.0, 0.0); magmaDoubleComplex zero = MAGMA_Z_MAKE(0.0, 0.0); B.storage_type = Magma_CSR; int i; for( i = 1; i < argc; ++i ) { if ( strcmp("--format", argv[i]) == 0 ) { format = atoi( argv[++i] ); switch( format ) { case 0: B.storage_type = Magma_CSR; break; case 1: B.storage_type = Magma_ELL; break; case 2: B.storage_type = Magma_ELLRT; break; case 3: B.storage_type = Magma_SELLP; break; } }else if ( strcmp("--mscale", argv[i]) == 0 ) { scale = atoi( argv[++i] ); switch( scale ) { case 0: scaling = Magma_NOSCALE; break; case 1: scaling = Magma_UNITDIAG; break; case 2: scaling = Magma_UNITROW; break; } }else if ( strcmp("--blocksize", argv[i]) == 0 ) { B.blocksize = atoi( argv[++i] ); }else if ( strcmp("--alignment", argv[i]) == 0 ) { B.alignment = atoi( argv[++i] ); }else if ( strcmp("--verbose", argv[i]) == 0 ) { solver_par.verbose = atoi( argv[++i] ); } else if ( strcmp("--ortho", argv[i]) == 0 ) { ortho = atoi( argv[++i] ); switch( ortho ) { case 0: solver_par.ortho = Magma_CGS; break; case 1: solver_par.ortho = Magma_MGS; break; case 2: solver_par.ortho = Magma_FUSED_CGS; break; } } else if ( strcmp("--restart", argv[i]) == 0 ) { solver_par.restart = atoi( argv[++i] ); } else if ( strcmp("--maxiter", argv[i]) == 0 ) { solver_par.maxiter = atoi( argv[++i] ); } else if ( strcmp("--tol", argv[i]) == 0 ) { sscanf( argv[++i], "%lf", &solver_par.epsilon ); } else break; } printf( "\n# usage: ./run_zgmres" " [ --format %d (0=CSR, 1=ELL 2=ELLRT, 3=SELLP)" " [ --blocksize %d --alignment %d ]" " --mscale %d (0=no, 1=unitdiag, 2=unitrownrm)" " --verbose %d (0=summary, k=details every k iterations)" " --restart %d --maxiter %d --tol %.2e" " --ortho %d (0=CGS, 1=MGS, 2=FUSED_CGS) ]" " matrices \n\n", format, (int) B.blocksize, (int) B.alignment, (int) scale, (int) solver_par.verbose, (int) solver_par.restart, (int) solver_par.maxiter, solver_par.epsilon, ortho ); magma_zsolverinfo_init( &solver_par, &precond_par ); while( i < argc ){ magma_z_csr_mtx( &A, argv[i] ); printf( "\n# matrix info: %d-by-%d with %d nonzeros\n\n", (int) A.num_rows,(int) A.num_cols,(int) A.nnz ); // scale matrix magma_zmscale( &A, scaling ); magma_z_mconvert( A, &B, Magma_CSR, B.storage_type ); magma_z_mtransfer( B, &B_d, Magma_CPU, Magma_DEV ); // vectors and initial guess magma_z_vinit( &b, Magma_DEV, A.num_cols, one ); magma_z_vinit( &x, Magma_DEV, A.num_cols, one ); magma_z_spmv( one, B_d, x, zero, b ); // b = A x magma_z_vfree(&x); magma_z_vinit( &x, Magma_DEV, A.num_cols, zero ); magma_zgmres( B_d, b, &x, &solver_par ); magma_zsolverinfo( &solver_par, &precond_par ); magma_z_mfree(&B_d); magma_z_mfree(&B); magma_z_mfree(&A); magma_z_vfree(&x); magma_z_vfree(&b); i++; } magma_zsolverinfo_free( &solver_par, &precond_par ); TESTING_FINALIZE(); return 0; }