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(); }