messages.c

数学算法的实现库。可以实现常见的线性计算。

#ifdef USE_MPILUC
    printf("M-ILUPACK PARAMETERS:\n");
#else
    printf("ILUPACK PARAMETERS:\n");
#endif
    printf("   droptol=%8.1le\n",droptols[1]);
    printf("   condest=%8.1le\n",condest);
    printf("   elbow space factor=%4d\n",  ELBOW);
    if (flags & TISMENETSKY_SC)
       printf("   Tismenetsky-like Schur complement\n");
    else
       if (flags & SIMPLE_SC)
	 printf("   simple Schur complement\n");
       else
	 printf("   medium Schur complement\n");
    if (flags & DIAGONAL_COMPENSATION)
       printf("   diagonal compensation\n");
    else
       printf("   diagonal entries unmodified\n");
    if (flags & AGGRESSIVE_DROPPING)
       printf("   aggressive dropping applied\n");

    if (flags & DISCARD_MATRIX)
       printf("   coarse grid matrix is discarded at the earliest stage\n");


    if (flags & PREPROCESS_INITIAL_SYSTEM) {
       printf("   preproc. initial system: ");
#ifdef MINIMUM_DEGREE
       printf("scaling + MMD ordering\n");
#elif defined REVERSE_CM
       printf("scaling + RCM ordering\n");
#elif defined NESTED_DISSECTION
       printf("scaling + ND ordering\n");
#elif defined IND_SET
       printf("scaling + independent set ordering\n");
#elif defined AMF
       printf("scaling + approximate minimum fill ordering\n");
#elif defined AMD
       printf("scaling + AMD ordering\n");
#elif defined MC64_RCM_PQ
       printf("scaling + MC64 max. weight matching + RCM ordering\n");
#elif defined MC64_MMD_PQ
       printf("scaling + MC64 max. weight matching + MMD ordering\n");
#elif defined MC64_AMF_PQ
       printf("scaling + MC64 max. weight matching + AMF ordering\n");
#elif defined MC64_AMD_PQ
       printf("scaling + MC64 max. weight matching + AMD ordering\n");
#elif defined MC64_METIS_E_PQ
       printf("scaling + MC64 max. weight matching + MeTis Edge ND ordering\n");
#elif defined MC64_METIS_N_PQ
       printf("scaling + MC64 max. weight matching + MeTis Node ND ordering\n");
#elif defined PQ_PERM
       printf("scaling + ddPQ ordering\n");
#elif defined MMD_PQ
       printf("scaling + MMD ordering\n");
#elif defined AMF_PQ
       printf("scaling + AMF ordering\n");
#elif defined AMD_PQ
       printf("scaling + AMD ordering\n");
#elif defined RCM_PQ
       printf("scaling + RCM ordering\n");
#elif defined FC_PQ
       printf("scaling + fine/coarse grid partitioning\n");
#elif defined METIS_E_PQ
       printf("scaling + MeTiS Edge ND ordering\n");
#elif defined METIS_N_PQ
       printf("scaling + MeTiS Node ND ordering\n");
#elif defined MWM_RCM_PQ
       printf("scaling + PARDISO max. weight matching + RCM ordering\n");
#elif defined MWM_MMD_PQ
       printf("scaling + PARDISO max. weight matching + MMD ordering\n");
#elif defined MWM_AMF_PQ
       printf("scaling + PARDISO max. weight matching + AMF ordering\n");
#elif defined MWM_AMD_PQ
       printf("scaling + PARDISO max. weight matching + AMD ordering\n");
#elif defined MWM_METIS_E_PQ
       printf("scaling + PARDISO max. weight matching + MeTiS Edge ND ordering\n");
#elif defined MWM_METIS_N_PQ
       printf("scaling + PARDISO max. weight matching + MeTiS Node ND ordering\n");
#endif
    }
    else
       printf("   initial system left unchanged\n");

    if (flags & PREPROCESS_SUBSYSTEMS) {
       printf("   reordering subsystems:   ");
#ifdef MINIMUM_DEGREE
       printf("scaling + MMD ordering\n");
#elif defined REVERSE_CM
       printf("scaling + RCM ordering\n");
#elif defined NESTED_DISSECTION
       printf("scaling + ND ordering\n");
#elif defined IND_SET
       printf("scaling + independent set ordering\n");
#elif defined AMF
       printf("scaling + AMF ordering\n");
#elif defined AMD
       printf("scaling + AMD ordering\n");
#elif defined MC64_RCM_PQ
       printf("scaling + RCM ordering\n");
#elif defined MC64_MMD_PQ
       printf("scaling + MMD ordering\n");
#elif defined MC64_AMF_PQ
       printf("scaling + AMF ordering\n");
#elif defined MC64_AMD_PQ
       printf("scaling + AMD ordering\n");
#elif defined MC64_METIS_E_PQ
       printf("scaling + MeTiS Edge ND ordering\n");
#elif defined MC64_METIS_N_PQ
       printf("scaling + MeTiS Node ND ordering\n");
#elif defined PQ_PERM
       printf("scaling + ddPQ ordering\n");
#elif defined MMD_PQ
       printf("scaling + MMD ordering\n");
#elif defined AMF_PQ
       printf("scaling + AMF ordering\n");
#elif defined AMD_PQ
       printf("scaling + AMD ordering\n");
#elif defined RCM_PQ
       printf("scaling + RCM ordering\n");
#elif defined FC_PQ
       printf("scaling + fine/coarse grid partitioning\n");
#elif defined METIS_E_PQ
       printf("scaling + MeTiS Edge ND ordering\n");
#elif defined METIS_N_PQ
       printf("scaling + MeTiS Node ND ordering\n");
#elif defined MWM_RCM_PQ
       printf("scaling + RCM ordering\n");
#elif defined MWM_MMD_PQ
       printf("scaling + MMD ordering\n");
#elif defined MWM_AMF_PQ
       printf("scaling + AMF ordering\n");
#elif defined MWM_AMD_PQ
       printf("scaling + AMD ordering\n");
#elif defined MWM_METIS_E_PQ
       printf("scaling + MeTiS Edge ND ordering\n");
#elif defined MWM_METIS_N_PQ
       printf("scaling + MeTiS Node ND ordering\n");
#endif
    }
    else
       printf("   remaining subsystems left unchanged\n");
    fflush(stdout);

    if (flags & FINAL_PIVOTING) {
       printf("   pivoting final system:   ");
#ifdef MINIMUM_DEGREE
       printf("scaling + MMD ordering\n");
#elif defined REVERSE_CM
       printf("scaling + RCM ordering\n");
#elif defined NESTED_DISSECTION
       printf("scaling + ND ordering\n");
#elif defined IND_SET
       printf("scaling + independent set ordering\n");
#elif defined AMF
       printf("scaling + AMF ordering\n");
#elif defined AMD
       printf("scaling + AMD ordering\n");
#elif defined MC64_RCM_PQ
       printf("scaling + ddPQ ordering\n");
#elif defined MC64_MMD_PQ
       printf("scaling + ddPQ ordering\n");
#elif defined MC64_AMF_PQ
       printf("scaling + ddPQ ordering\n");
#elif defined MC64_AMD_PQ
       printf("scaling + ddPQ ordering\n");
#elif defined MC64_METIS_E_PQ
       printf("scaling + ddPQ ordering\n");
#elif defined MC64_METIS_N_PQ
       printf("scaling + ddPQ ordering\n");
#elif defined PQ_PERM
       printf("scaling + ddPQ ordering\n");
#elif defined MMD_PQ
       printf("scaling + ddPQ ordering\n");
#elif defined AMF_PQ
       printf("scaling + ddPQ ordering\n");
#elif defined AMD_PQ
       printf("scaling + ddPQ ordering\n");
#elif defined RCM_PQ
       printf("scaling + ddPQ ordering\n");
#elif defined FC_PQ
       printf("scaling + ddPQ ordering\n");
#elif defined METIS_E_PQ
       printf("scaling + ddPQ ordering\n");
#elif defined METIS_N_PQ
       printf("scaling + ddPQ ordering\n");
#elif defined MWM_RCM_PQ
       printf("scaling + ddPQ ordering\n");
#elif defined MWM_MMD_PQ
       printf("scaling + ddPQ ordering\n");
#elif defined MWM_AMF_PQ
       printf("scaling + ddPQ ordering\n");
#elif defined MWM_AMD_PQ
       printf("scaling + ddPQ ordering\n");
#elif defined MWM_METIS_E_PQ
       printf("scaling + ddPQ ordering\n");
#elif defined MWM_METIS_N_PQ
       printf("scaling + ddPQ ordering\n");
#endif
    }
    else
       printf("   final subsystem(s) left unchanged\n");
    fflush(stdout);

#ifdef USE_MPILUC
    if (flags & TISMENETSKY_SC)
      fprintf(fo,"M-ILPCK(T)|");
    else
       if (flags & SIMPLE_SC)
	 fprintf(fo,"M-ILPCK(S)|");
       else
	 fprintf(fo,"M-ILPCK(M)|");
#else
    if (flags & TISMENETSKY_SC)
      fprintf(fo,"ILUPACK(T)|");
    else
       if (flags & SIMPLE_SC)
	 fprintf(fo,"ILUPACK(S)|");
       else
	 fprintf(fo,"ILUPACK(M)|");
#endif