int quaddobl_laursys_solver ( void ) { int fail,nbq,nbtasks,topdim,filter,factor,verbose; fail = syscon_read_quaddobl_Laurent_system(); fail = syscon_number_of_quaddobl_Laurentials(&nbq); printf("-> read %d polynomials\n", nbq); printf("\n"); read_solver_options(&nbtasks,&topdim,&filter,&factor,&verbose); printf("\nCalling the solver ...\n\n"); fail = quaddobl_laursys_solve(nbtasks,topdim,filter,factor,verbose); fail = quaddobl_laursys_write(topdim); return fail; }
void test_quaddobl_Newton_Laurent_step ( void ) { int fail,dim,len; printf("\nRunning Newton step with quad double arithmetic ...\n"); fail = syscon_read_quaddobl_Laurent_system(); fail = syscon_number_of_quaddobl_Laurentials(&dim); printf("The system container has %d Laurent polynomials.\n",dim); fail = solcon_read_quaddobl_solutions(); fail = solcon_number_of_quaddobl_solutions(&len); printf("The solution container has size %d.\n",len); fail = solcon_dimension_of_quaddobl_solutions(&dim); printf("The solutions in the container have dimension %d.\n",dim); fail = quaddobl_Newton_Laurent_step(); printf("The solutions after the Newton step :\n"); fail = solcon_write_quaddobl_solutions(); }
int Laurent_input_output_on_files ( int precision ) { int fail,rc,nbtasks; if(precision == 0) { fail = syscon_read_standard_Laurent_system(); printf("\nThe system in the container : \n"); fail = syscon_write_standard_Laurent_system(); printf("\nGive the number of tasks : "); scanf("%d",&nbtasks); fail = solve_Laurent_system(&rc,0,nbtasks); /* not silent by default */ printf("\nThe root count : %d\n",rc); printf("\nThe solutions :\n"); fail = solcon_write_standard_solutions(); } else if(precision == 1) { fail = syscon_read_dobldobl_Laurent_system(); printf("\nThe system in the container : \n"); fail = syscon_write_dobldobl_Laurent_system(); printf("\nGive the number of tasks : "); scanf("%d",&nbtasks); fail = solve_dobldobl_Laurent_system(&rc,0,nbtasks); printf("\nThe root count : %d\n",rc); printf("\nThe solutions :\n"); fail = solcon_write_dobldobl_solutions(); } else if(precision == 2) { fail = syscon_read_quaddobl_Laurent_system(); printf("\nThe system in the container : \n"); fail = syscon_write_quaddobl_Laurent_system(); printf("\nGive the number of tasks : "); scanf("%d",&nbtasks); fail = solve_quaddobl_Laurent_system(&rc,0,nbtasks); printf("\nThe root count : %d\n",rc); printf("\nThe solutions :\n"); fail = solcon_write_quaddobl_solutions(); } return fail; }