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