symmessages.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("   no diagonal compensation\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 PP_PERM
       printf("scaling + ddPP ordering\n");
#elif defined MMD_PP
       printf("scaling + MMD ordering\n");
#elif defined AMF_PP
       printf("scaling + AMF ordering\n");
#elif defined AMD_PP
       printf("scaling + AMD ordering\n");
#elif defined RCM_PP
       printf("scaling + RCM ordering\n");
#elif defined FC_PP
       printf("scaling + fine/coarse grid partitioning\n");
#elif defined METIS_E_PP
       printf("scaling + MeTiS Edge ND ordering\n");
#elif defined METIS_N_PP
       printf("scaling + MeTiS Node ND ordering\n");
#elif defined SMWM_RCM_PP
       printf("scaling + PARDISO sym. max. w. matching + RCM ordering\n");
#elif defined SMWM_MMD_PP
       printf("scaling + PARDISO sym. max. w. matching + MMD ordering\n");
#elif defined SMWM_AMF_PP
       printf("scaling + PARDISO sym. max. w. matching + AMF ordering\n");
#elif defined SMWM_AMD_PP
       printf("scaling + PARDISO sym. max. w. matching + AMD ordering\n");
#elif defined SMWM_METIS_E_PP
       printf("scaling + PARDISO sym. max. w. matching + MeTiS Edge ND ordering\n");
#elif defined SMWM_METIS_N_PP
       printf("scaling + PARDISO sym. max. w. matching + MeTiS Node ND ordering\n");
#elif defined SMC64_RCM_PP
       printf("scaling + MC64 sym. max. w. matching + RCM ordering\n");
#elif defined SMC64_MMD_PP
       printf("scaling + MC64 sym. max. w. matching + MMD ordering\n");
#elif defined SMC64_AMF_PP
       printf("scaling + MC64 sym. max. w. matching + AMF ordering\n");
#elif defined SMC64_AMD_PP
       printf("scaling + MC64 sym. max. w. matching + AMD ordering\n");
#elif defined SMC64_METIS_E_PP
       printf("scaling + MC64 sym. max. w. matching + MeTiS Edge ND ordering\n");
#elif defined SMC64_METIS_N_PP
       printf("scaling + MC64 sym. max. w. matching + MeTiS Node ND ordering\n");
#elif defined FC_SMWM_METIS_E
       printf("scaling + fine/coarse grid partitioning\n");
#elif defined FC_SMWM_METIS_N
       printf("scaling + fine/coarse grid partitioning\n");
#elif defined FC_SMWM_AMD
       printf("scaling + fine/coarse grid partitioning\n");
#elif defined FC_SMWM_AMF
       printf("scaling + fine/coarse grid partitioning\n");
#elif defined FC_SMWM_MMD
       printf("scaling + fine/coarse grid partitioning\n");
#elif defined FC_SMWM_RCM
       printf("scaling + fine/coarse grid partitioning\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 + approximate minimum fill ordering\n");
#elif defined AMD
       printf("scaling + AMD ordering\n");
#elif defined PP_PERM
       printf("scaling + ddPP ordering\n");
#elif defined MMD_PP
       printf("scaling + MMD ordering\n");
#elif defined AMF_PP
       printf("scaling + AMF ordering\n");
#elif defined AMD_PP
       printf("scaling + AMD ordering\n");
#elif defined RCM_PP
       printf("scaling + RCM ordering\n");
#elif defined FC_PP
       printf("scaling + fine/coarse grid partitioning\n");
#elif defined METIS_E_PP
       printf("scaling + MeTiS Edge ND ordering\n");
#elif defined METIS_N_PP
       printf("scaling + MeTiS Node ND ordering\n");
#elif defined SMWM_RCM_PP
       printf("scaling + PARDISO sym. max. w. matching + RCM ordering\n");
#elif defined SMWM_MMD_PP
       printf("scaling + PARDISO sym. max. w. matching + MMD ordering\n");
#elif defined SMWM_AMF_PP
       printf("scaling + PARDISO sym. max. w. matching + AMF ordering\n");
#elif defined SMWM_AMD_PP
       printf("scaling + PARDISO sym. max. w. matching + AMD ordering\n");
#elif defined SMWM_METIS_E_PP
       printf("scaling + PARDISO sym. max. w. matching + MeTiS Edge ND ordering\n");
#elif defined SMWM_METIS_N_PP
       printf("scaling + PARDISO sym. max. w. matching + MeTiS Node ND ordering\n");
#elif defined SMC64_RCM_PP
       printf("scaling + MC64 sym. max. w. matching + RCM ordering\n");
#elif defined SMC64_MMD_PP
       printf("scaling + MC64 sym. max. w. matching + MMD ordering\n");
#elif defined SMC64_AMF_PP
       printf("scaling + MC64 sym. max. w. matching + AMF ordering\n");
#elif defined SMC64_AMD_PP
       printf("scaling + MC64 sym. max. w. matching + AMD ordering\n");
#elif defined SMC64_METIS_E_PP
       printf("scaling + MC64 sym. max. w. matching + MeTiS Edge ND ordering\n");
#elif defined SMC64_METIS_N_PP
       printf("scaling + MC64 sym. max. w. matching + MeTiS Node ND ordering\n");
#elif defined FC_SMWM_METIS_E
       printf("scaling + fine/coarse grid partitioning\n");
#elif defined FC_SMWM_METIS_N
       printf("scaling + fine/coarse grid partitioning\n");
#elif defined FC_SMWM_AMD
       printf("scaling + fine/coarse grid partitioning\n");
#elif defined FC_SMWM_AMF
       printf("scaling + fine/coarse grid partitioning\n");
#elif defined FC_SMWM_MMD
       printf("scaling + fine/coarse grid partitioning\n");
#elif defined FC_SMWM_RCM
       printf("scaling + fine/coarse grid partitioning\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 + approximate minimum fill ordering\n");
#elif defined AMD
       printf("scaling + AMD ordering\n");
#elif defined PP_PERM
       printf("scaling + ddPP ordering\n");
#elif defined MMD_PP
       printf("scaling + ddPP ordering\n");
#elif defined AMF_PP
       printf("scaling + ddPP ordering\n");
#elif defined AMD_PP
       printf("scaling + ddPP ordering\n");
#elif defined RCM_PP
       printf("scaling + ddPP ordering\n");
#elif defined FC_PP
       printf("scaling + ddPP ordering\n");
#elif defined METIS_E_PP
       printf("scaling + ddPP ordering\n");
#elif defined METIS_N_PP
       printf("scaling + ddPP ordering\n");
#elif defined SMWM_RCM_PP
       printf("scaling + ddPP ordering\n");
#elif defined SMWM_MMD_PP
       printf("scaling + ddPP ordering\n");
#elif defined SMWM_AMF_PP
       printf("scaling + ddPP ordering\n");
#elif defined SMWM_AMD_PP
       printf("scaling + ddPP ordering\n");
#elif defined SMWM_METIS_E_PP
       printf("scaling + ddPP ordering\n");
#elif defined SMWM_METIS_N_PP
       printf("scaling + ddPP ordering\n");
#elif defined SMC64_RCM_PP
       printf("scaling + ddPP ordering\n");
#elif defined SMC64_MMD_PP
       printf("scaling + ddPP ordering\n");
#elif defined SMC64_AMF_PP
       printf("scaling + ddPP ordering\n");
#elif defined SMC64_AMD_PP
       printf("scaling + ddPP ordering\n");
#elif defined SMC64_METIS_E_PP
       printf("scaling + ddPP ordering\n");
#elif defined SMC64_METIS_N_PP
       printf("scaling + ddPP ordering\n");
#elif defined FC_SMWM_METIS_E
       printf("scaling + PARDISO sym. max. w. matching + MeTiS (node) ordering\n");
#elif defined FC_SMWM_METIS_N
       printf("scaling + PARDISO sym. max. w. matching + MeTiS (edge) ordering\n");
#elif defined FC_SMWM_AMD
       printf("scaling + PARDISO sym. max. w. matching + AMD ordering\n");
#elif defined FC_SMWM_AMF
       printf("scaling + PARDISO sym. max. w. matching + AMF ordering\n");
#elif defined FC_SMWM_MMD
       printf("scaling + PARDISO sym. max. w. matching + MMD ordering\n");
#elif defined FC_SMWM_RCM
       printf("scaling + PARDISO sym. max. w. matching + RCM 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