/*
 * Copyright (C) 2003 Joern Thyssen <jth@gnubg.org>
 * Copyright (C) 2003-2021 the AUTHORS
 *
 * This program is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program.  If not, see <https://www.gnu.org/licenses/>.
 */

#include "config.h"
#include "backgammon.h"

#include <glib.h>
#include <string.h>

#include "eval.h"
#include "format.h"

#include "export.h"
#include "positionid.h"


int fOutputMWC = FALSE;
int fOutputWinPC = FALSE;
int fOutputMatchPC = TRUE;

#define OUTPUT_SZ_LENGTH (5 + MAX_OUTPUT_DIGITS)
int fOutputDigits = 3;

float rErrorRateFactor = 1000.0f;

typedef int (*classdumpfunc) (const TanBoard anBoard, char *szOutput, const bgvariation bgv);

extern char *
OutputRolloutResult(const char *szIndent,
                    char asz[][FORMATEDMOVESIZE],
                    float aarOutput[][NUM_ROLLOUT_OUTPUTS],
                    float aarStdDev[][NUM_ROLLOUT_OUTPUTS],
                    const cubeinfo aci[], const int alt, const int cci, const int fCubeful)
{

    static char sz[1024];
    int ici;

    strcpy(sz, "");

    for (ici = 0; ici < cci; ici++) {

        /* header */

        if (asz && *asz[ici]) {
            if (szIndent && *szIndent)
                strcat(sz, szIndent);
            sprintf(strchr(sz, 0), "%s:\n", asz[ici]);
        }

        /* output */

        if (szIndent && *szIndent)
            strcat(sz, szIndent);

        strcat(sz, "  ");
        strcat(sz, OutputPercents(aarOutput[ici], TRUE));
        strcat(sz, " ");
        strcat(sz, _("CL"));
        strcat(sz, " ");
        strcat(sz, OutputEquityScale(aarOutput[ici][OUTPUT_EQUITY], &aci[alt + ici], &aci[0], TRUE));

        if (fCubeful) {
            strcat(sz, " ");
            strcat(sz, _("CF"));
            strcat(sz, " ");
            strcat(sz, OutputMWC(aarOutput[ici][OUTPUT_CUBEFUL_EQUITY], &aci[0], TRUE));
        }

        strcat(sz, "\n");

        /* std dev */

        if (szIndent && *szIndent)
            strcat(sz, szIndent);

        strcat(sz, " [");
        strcat(sz, OutputPercents(aarStdDev[ici], FALSE));
        strcat(sz, " ");
        strcat(sz, _("CL"));
        strcat(sz, " ");
        strcat(sz, OutputEquityScale(aarStdDev[ici][OUTPUT_EQUITY], &aci[alt + ici], &aci[0], FALSE));

        if (fCubeful) {
            strcat(sz, " ");
            strcat(sz, _("CF"));
            strcat(sz, " ");
            strcat(sz, OutputMWC(aarStdDev[ici][OUTPUT_CUBEFUL_EQUITY], &aci[0], FALSE));
        }

        strcat(sz, "]\n");

    }

    return sz;

}


extern char *
OutputEvalContext(const evalcontext * pec, const int fChequer)
{

    static char sz[1024];
    int i;

    sprintf(sz, "%u-%s %s", pec->nPlies, _("ply"), (!fChequer || pec->fCubeful) ? _("cubeful") : _("cubeless"));

    if (pec->fUsePrune) {
        sprintf(strchr(sz, 0), " %s", _("prune"));
    }

    if (fChequer && pec->nPlies) {
        /* FIXME: movefilters!!! */
    }

    if (pec->rNoise > 0.0f)
        sprintf(strchr(sz, 0), ", %s %0.3g (%s)", _("noise"), pec->rNoise, pec->fDeterministic ? _("d") : _("nd"));

    for (i = 0; i < NUM_SETTINGS; i++)

        if (!cmp_evalcontext(&aecSettings[i], pec)) {
            sprintf(strchr(sz, 0), " [%s]", Q_(aszSettings[i]));
            break;
        }

    return sz;

}



/* return -1 if this isn't a predefined setting
 * n = index of predefined setting if it's a match or
 * if the eval context matches and it's 0 ply
 */
static int
GetPredefinedChequerplaySetting(const evalcontext * pec, const movefilter aamf[MAX_FILTER_PLIES][MAX_FILTER_PLIES])
{

    int nEval;
    int nFilter;
    int nPlies = pec->nPlies;
    int i;
    int fSame;
    int nPreset;
    int Accept;

    if (nPlies > MAX_FILTER_PLIES) {
        return -1;
    }

    for (nEval = 0; nEval < NUM_SETTINGS; ++nEval) {
        if (cmp_evalcontext(aecSettings + nEval, pec) == 0) {
            /* eval matches and it's 0 ply, we have a predefined one */
            if (nPlies == 0)
                return nEval;

            /* see if there's a filter set which matches and uses the
             * same eval context (world class and supremo use the same
             * eval context and different filters) */
            for (nFilter = 0; nFilter < NUM_SETTINGS; ++nFilter) {
                /* see what filter set goes with the predefined settings */
                if ((nPreset = aiSettingsMoveFilter[nFilter]) < 0)
                    continue;

                /* nPreset = 0/tiny, 1/normal, 2/large, 3/huge */
                fSame = 1;
                for (i = 0; i < nPlies - 1; ++i) {
                    if ((Accept = aamf[nPlies - 1][i].Accept) != aaamfMoveFilterSettings[nPreset][nPlies - 1][i].Accept) {
                        fSame = 0;
                        break;
                    }

                    /* if they are both ignore this level, don't check 
                     * the extra and threshold */
                    if (Accept < 0)
                        continue;

                    if (aamf[nPlies - 1][i].Extra != aaamfMoveFilterSettings[nPreset][nPlies - 1][i].Extra) {
                        fSame = 0;
                        break;
                    }

                    if (fabs(aamf[nPlies - 1][i].Threshold - aaamfMoveFilterSettings[nPreset][nPlies - 1][i].Threshold)
                        > 0.1e-6) {
                        fSame = 0;
                        break;
                    }
                }
                /* the filters match (for this nPlies, so we may have a 
                 * preset */
                if (fSame && (nFilter == nEval))
                    return nEval;
            }
        }
    }

    return -1;
}

static char *
OutputMoveFilterPly(const char *szIndent, const int nPlies, const movefilter aamf[MAX_FILTER_PLIES][MAX_FILTER_PLIES])
{

    static char sz[1024];
    int i;

    strcpy(sz, "");

    if (!nPlies)
        return sz;

    for (i = 0; i < nPlies; ++i) {

        const movefilter *pmf = &aamf[nPlies - 1][i];

        if (szIndent && *szIndent)
            strcat(sz, szIndent);

        if (pmf->Accept < 0) {
            sprintf(strchr(sz, 0), _("skip pruning for %d-ply moves"), i);
            strcat(sz, "\n");
            continue;
        }

        if (pmf->Accept == 1)
            sprintf(strchr(sz, 0), _("keep the best %d-ply move"), i);
        else
            sprintf(strchr(sz, 0), _("keep the first %d %d-ply moves"), pmf->Accept, i);

        if (pmf->Extra)
            sprintf(strchr(sz, 0), _(" and up to %d more moves within equity %0.3g"), pmf->Extra, pmf->Threshold);

        strcat(sz, "\n");
    }

    return sz;

}

static void
OutputEvalContextsForRollout(char *sz, const char *szIndent,
                             const evalcontext aecCube[2],
                             const evalcontext aecChequer[2],
                             const movefilter aaamf[2][MAX_FILTER_PLIES][MAX_FILTER_PLIES])
{

    int fCube = !cmp_evalcontext(&aecCube[0], &aecCube[1]);
    int fChequer = !cmp_evalcontext(&aecChequer[0], &aecChequer[1]);
    int fMovefilter = fChequer && (!aecChequer[0].nPlies || equal_movefilter(aecChequer[0].nPlies - 1,
                                                                             aaamf[0][aecChequer[0].nPlies - 1],
                                                                             aaamf[1][aecChequer[0].nPlies - 1]));
    int fIdentical = fCube && fChequer && fMovefilter;
    int i;

    for (i = 0; i < 1 + !fIdentical; i++) {
        int j;

        if (!fIdentical) {

            if (szIndent && *szIndent)
                strcat(sz, szIndent);

            sprintf(strchr(sz, 0), "%s %d:\n", _("Player"), i);
        }

        /* chequer play */

        j = GetPredefinedChequerplaySetting(&aecChequer[i], aaamf[i]);

        if (szIndent && *szIndent)
            strcat(sz, szIndent);

        if (aecChequer[i].nPlies) {

            sprintf(strchr(sz, 0), "%s: %s ", _("Play"), (j < 0) ? "" : Q_(aszSettings[j]));

            strcat(sz, OutputEvalContext(&aecChequer[i], TRUE));
            strcat(sz, "\n");
            strcat(sz, OutputMoveFilterPly(szIndent, aecChequer[i].nPlies, aaamf[i]));

        } else {
            sprintf(strchr(sz, 0), "%s: ", _("Play"));
            strcat(sz, OutputEvalContext(&aecChequer[i], FALSE));
            strcat(sz, "\n");
        }

        if (szIndent && *szIndent)
            strcat(sz, szIndent);

        sprintf(strchr(sz, 0), "%s: ", _("Cube"));
        strcat(sz, OutputEvalContext(&aecCube[i], FALSE));
        strcat(sz, "\n");

    }

}



extern char *
OutputRolloutContext(const char *szIndent, const rolloutcontext * prc)
{

    static char sz[1024];

    strcpy(sz, "");

    if (szIndent && *szIndent)
        strcat(sz, szIndent);

    if (prc->nTruncate && prc->fDoTruncate)
        sprintf(strchr(sz, 0),
                prc->fCubeful ?
                _("Truncated cubeful rollout (depth %d)") : _("Truncated cubeless rollout (depth %d)"), prc->nTruncate);
    else
        sprintf(strchr(sz, 0), prc->fCubeful ? _("Full cubeful rollout") : _("Full cubeless rollout"));

    if (prc->fTruncBearoffOS && !prc->fCubeful)
        sprintf(strchr(sz, 0), " (%s)", _("truncated at one-sided bearoff"));
    else if (prc->fTruncBearoff2 && !prc->fCubeful)
        sprintf(strchr(sz, 0), " (%s)", _("truncated at exact bearoff"));

    sprintf(strchr(sz, 0), " %s", prc->fVarRedn ? _("with variance reduction") : _("without variance reduction"));

    strcat(sz, "\n");

    if (szIndent && *szIndent)
        strcat(sz, szIndent);

    sprintf(strchr(sz, 0), _("%u games"), prc->nGamesDone);

    if (prc->fInitial) {
        strcat(sz, ", ");
        strcat(sz, _("rollout as initial position"));
    }

    strcat(sz, ", ");
    if (prc->fRotate)
        sprintf(strchr(sz, 0), _("%s dice gen. with seed %lu and quasi-random dice"), gettext(aszRNG[prc->rngRollout]), prc->nSeed);    /* seed may be unsigned long int */
    else
        sprintf(strchr(sz, 0), _("%s dice generator with seed %lu"), gettext(aszRNG[prc->rngRollout]), prc->nSeed);     /* seed may be unsigned long int */

    strcat(sz, "\n");

    if ((prc->fStopOnJsd || prc->fStopOnSTD)
        && szIndent && *szIndent)
        strcat(sz, szIndent);

    /* stop on std.err */

    if (prc->fStopOnSTD && !prc->fStopOnJsd) {
        sprintf(strchr(sz, 0),
                _("Stop when std.errs. are small enough: limit "
                  "%.4f (min. %u games)"), prc->rStdLimit, prc->nMinimumGames);
        strcat(sz, "\n");
    }

    /* stop on JSD */
    if (prc->fStopOnJsd) {
        sprintf(strchr(sz, 0),
                _("Stop when best play is enough JSDs ahead: limit "
                  "%.4f (min. %u games)"), prc->rJsdLimit, prc->nMinimumJsdGames);
        strcat(sz, "\n");
    }

    /* first play */

    OutputEvalContextsForRollout(sz, szIndent, prc->aecCube, prc->aecChequer, prc->aaamfChequer);

    /* later play */

    if (prc->fLateEvals) {

        if (szIndent && *szIndent)
            strcat(sz, szIndent);

        sprintf(strchr(sz, 0), _("Different evaluations after %d plies:"), prc->nLate);
        strcat(sz, "\n");

        OutputEvalContextsForRollout(sz, szIndent, prc->aecCubeLate, prc->aecChequerLate, prc->aaamfLate);


    }

    return sz;

}

/*
 * Return formatted string with equity or MWC.
 *
 * Input:
 *    r: equity (either money equity for normalised money eq. for match play
 *    pci: cubeinfo
 *    f: indicates equity (TRUE) or std. error (FALSE)
 *    
 *
 * Important: function is not re-entrant. Caller must save output
 * if needed.
 */

extern char *
OutputEquity(const float r, const cubeinfo * pci, const int f)
{

    static char sz[OUTPUT_SZ_LENGTH];

    if (!pci->nMatchTo || !fOutputMWC) {
        if (f)
            snprintf(sz, OUTPUT_SZ_LENGTH, "%+*.*f", fOutputDigits + 3, fOutputDigits, r);
        else
            snprintf(sz, OUTPUT_SZ_LENGTH, "% *.*f", fOutputDigits + 2, fOutputDigits, r);
    } else {
        if (fOutputMatchPC) {
            snprintf(sz, OUTPUT_SZ_LENGTH, "%*.*f%%", fOutputDigits + 3, fOutputDigits > 1 ? fOutputDigits - 1 : 0,
                     100.0f * (f ? eq2mwc(r, pci) : se_eq2mwc(r, pci)));
        } else {
            snprintf(sz, OUTPUT_SZ_LENGTH, "%*.*f", fOutputDigits + 3, fOutputDigits + 1,
                     f ? eq2mwc(r, pci) : se_eq2mwc(r, pci));
        }
    }

    return sz;

}


extern char *
OutputMoneyEquity(const float ar[], const int f)
{

    static char sz[OUTPUT_SZ_LENGTH];
    float eq = 2.0f * ar[OUTPUT_WIN] - 1.0f + ar[OUTPUT_WINGAMMON] + ar[OUTPUT_WINBACKGAMMON] -
        ar[OUTPUT_LOSEGAMMON] - ar[OUTPUT_LOSEBACKGAMMON];

    if (f)
        snprintf(sz, OUTPUT_SZ_LENGTH, "%+*.*f", fOutputDigits + 3, fOutputDigits, eq);
    else
        snprintf(sz, OUTPUT_SZ_LENGTH, "% *.*f", fOutputDigits + 2, fOutputDigits, eq);


    return sz;

}



/*
 * Return formatted string with equity or MWC.
 *
 * Input:
 *    r: equity (either money equity for normalised money eq. for match play
 *    pci: cubeinfo
 *    f: indicates equity (TRUE) or std. error (FALSE)
 *    
 *
 * Important: function is not re-entrant. Caller must save output
 * if needed.
 */

extern char *
OutputEquityScale(const float r, const cubeinfo * pci, const cubeinfo * pciBase, const int f)
{

    static char sz[OUTPUT_SZ_LENGTH];

    if (!pci->nMatchTo) {
        if (f)
            snprintf(sz, OUTPUT_SZ_LENGTH, "%+*.*f", fOutputDigits + 3, fOutputDigits,
                     (float) (pci->nCube / pciBase->nCube) * r);
        else
            snprintf(sz, OUTPUT_SZ_LENGTH, "% *.*f", fOutputDigits + 2, fOutputDigits,
                     (float) (pci->nCube / pciBase->nCube) * r);
    } else {

        if (fOutputMWC) {

            if (fOutputMatchPC) {
                snprintf(sz, OUTPUT_SZ_LENGTH, "%*.*f%%", fOutputDigits + 3, fOutputDigits > 1 ? fOutputDigits - 1 : 0,
                         100.0f * (f ? eq2mwc(r, pci) : se_eq2mwc(r, pci)));
            } else {
                snprintf(sz, OUTPUT_SZ_LENGTH, "%*.*f", fOutputDigits + 3, fOutputDigits + 1,
                         f ? eq2mwc(r, pci) : se_eq2mwc(r, pci));
            }

        } else {
            if (f)
                snprintf(sz, OUTPUT_SZ_LENGTH, "%+*.*f", fOutputDigits + 3, fOutputDigits,
                         mwc2eq(eq2mwc(r, pci), pciBase));
            else
                snprintf(sz, OUTPUT_SZ_LENGTH, "% *.*f", fOutputDigits + 2, fOutputDigits,
                         se_mwc2eq(se_eq2mwc(r, pci), pciBase));
        }


    }

    return sz;

}


/*
 * Return formatted string with equity or MWC for an equity difference.
 *
 * Input:
 *    r1, r2: equities (either money equity for normalised money eq. for match play
 *    pci: cubeinfo
 *
 * Important: function is not re-entrant. Caller must save output
 * if needed.
 */

extern char *
OutputEquityDiff(const float r1, const float r2, const cubeinfo * pci)
{

    static char sz[OUTPUT_SZ_LENGTH];

    if (fOutputDigits > MAX_OUTPUT_DIGITS) {
        g_assert_not_reached();
        return NULL;
    }

    /* If the difference is 0, print it as -0 in the right format.
     * Using IEEE754 negative zero is tricky and would not work
     * if compiled with -ffast-math (which we tend to do for performance
     * reasons). Make a special case and print the minus sign explicitely.
     */

    if (!pci->nMatchTo || !fOutputMWC) {
        if (r1 == r2)
            snprintf(sz, OUTPUT_SZ_LENGTH, "-%*.*f", fOutputDigits + 2, fOutputDigits, 0.0);
        else
            snprintf(sz, OUTPUT_SZ_LENGTH, "%+*.*f", fOutputDigits + 3, fOutputDigits, r1 - r2);
    } else {
        if (fOutputMatchPC) {
            if (r1 == r2)
                snprintf(sz, OUTPUT_SZ_LENGTH, " -%*.*f%%", fOutputDigits + 1,
                         fOutputDigits > 1 ? fOutputDigits - 1 : 0, 0.0);
            else
                snprintf(sz, OUTPUT_SZ_LENGTH, "%+*.*f%%", fOutputDigits + 3, fOutputDigits > 1 ? fOutputDigits - 1 : 0,
                         100.0f * eq2mwc(r1, pci) - 100.0f * eq2mwc(r2, pci));
        } else {
            if (r1 == r2)
                snprintf(sz, OUTPUT_SZ_LENGTH, "-%*.*f", fOutputDigits + 2, fOutputDigits + 1, 0.0);
            else
                snprintf(sz, OUTPUT_SZ_LENGTH, "+%*.*f", fOutputDigits + 3, fOutputDigits + 1,
                         eq2mwc(r1, pci) - eq2mwc(r2, pci));
        }
    }

    return sz;

}

/*
 * Return formatted string with equity or MWC.
 *
 * Input is: equity for money game/MWC for match play.
 *
 * Important: function is not re-entrant. Caller must save output
 * if needed.
 */


extern char *
OutputMWC(const float r, const cubeinfo * pci, const int f)
{

    static char sz[OUTPUT_SZ_LENGTH];

    if (!pci->nMatchTo) {
        if (f)
            snprintf(sz, OUTPUT_SZ_LENGTH, "%+*.*f", fOutputDigits + 3, fOutputDigits, r);
        else
            snprintf(sz, OUTPUT_SZ_LENGTH, "% *.*f", fOutputDigits + 2, fOutputDigits, r);
    } else {

        if (!fOutputMWC) {
            if (f)
                snprintf(sz, OUTPUT_SZ_LENGTH, "%+*.*f", fOutputDigits + 3, fOutputDigits, mwc2eq(r, pci));
            else
                snprintf(sz, OUTPUT_SZ_LENGTH, "% *.*f", fOutputDigits + 3, fOutputDigits, se_mwc2eq(r, pci));
        } else if (fOutputMatchPC) {
            snprintf(sz, OUTPUT_SZ_LENGTH, "%*.*f%%", fOutputDigits + 3, fOutputDigits > 1 ? fOutputDigits - 1 : 0,
                     100.0f * r);
        } else {
            if (f)
                snprintf(sz, OUTPUT_SZ_LENGTH, "%+*.*f", fOutputDigits + 3, fOutputDigits + 1, r);
            else
                snprintf(sz, OUTPUT_SZ_LENGTH, "% *.*f", fOutputDigits + 3, fOutputDigits + 1, r);
        }
    }

    return sz;

}


extern char *
OutputPercent(const float r)
{

    static char sz[OUTPUT_SZ_LENGTH];

    if (fOutputWinPC) {
        snprintf(sz, OUTPUT_SZ_LENGTH, "%*.*f", fOutputDigits + 2, fOutputDigits > 2 ? fOutputDigits - 2 : 0,
                 100.0f * r);
    } else {
        snprintf(sz, OUTPUT_SZ_LENGTH, "%*.*f", fOutputDigits + 2, fOutputDigits, r);
    }

    return sz;

}

extern char *
OutputPercents(const float ar[], const int f)
{

    static char sz[80];

    strcpy(sz, "");

    strcat(sz, OutputPercent(ar[OUTPUT_WIN]));
    strcat(sz, " ");
    strcat(sz, OutputPercent(ar[OUTPUT_WINGAMMON]));
    strcat(sz, " ");
    strcat(sz, OutputPercent(ar[OUTPUT_WINBACKGAMMON]));
    strcat(sz, " - ");
    if (f)
        strcat(sz, OutputPercent(1.0f - ar[OUTPUT_WIN]));
    else
        strcat(sz, OutputPercent(ar[OUTPUT_WIN]));
    strcat(sz, " ");
    strcat(sz, OutputPercent(ar[OUTPUT_LOSEGAMMON]));
    strcat(sz, " ");
    strcat(sz, OutputPercent(ar[OUTPUT_LOSEBACKGAMMON]));

    return sz;

}



/*
 * Print cube analysis
 *
 * Input:
 *  pf: output file
 *  arDouble: equitites for cube decisions
 *  fPlayer: player who doubled
 *  esDouble: eval setup
 *  pci: cubeinfo
 *  fDouble: double/no double
 *  fTake: take/drop
 *
 */

extern char *
OutputCubeAnalysisFull(float aarOutput[2][NUM_ROLLOUT_OUTPUTS],
                       float aarStdDev[2][NUM_ROLLOUT_OUTPUTS],
                       const evalsetup * pes, const cubeinfo * pci,
                       int fDouble, int fTake, skilltype stDouble, skilltype stTake)
{

    float r;

    int fMissed;
    int fAnno = FALSE;

    float arDouble[4];

    static char sz[4096];


    strcpy(sz, "");

    /* check if cube analysis should be printed */

    if (pes->et == EVAL_NONE)
        return NULL;            /* no evaluation */

    FindCubeDecision(arDouble, aarOutput, pci);

    /* print alerts */

    fMissed = fDouble > -1 && isMissedDouble(arDouble, aarOutput, fDouble, pci);

    /* print alerts */

    if (fMissed) {

        fAnno = TRUE;

        /* missed double */

        sprintf(strchr(sz, 0), "%s (%s)!\n",
                _("Alert: missed double"),
                OutputEquityDiff(arDouble[OUTPUT_NODOUBLE],
                                 (arDouble[OUTPUT_TAKE] >
                                  arDouble[OUTPUT_DROP]) ? arDouble[OUTPUT_DROP] : arDouble[OUTPUT_TAKE], pci));

        if (badSkill(stDouble))
            sprintf(strchr(sz, 0), " [%s]", gettext(aszSkillType[stDouble]));

    }

    r = arDouble[OUTPUT_TAKE] - arDouble[OUTPUT_DROP];

    if (fTake > 0 && r > 0.0f) {

        fAnno = TRUE;

        /* wrong take */

        sprintf(strchr(sz, 0), "%s (%s)!\n",
                _("Alert: wrong take"), OutputEquityDiff(arDouble[OUTPUT_DROP], arDouble[OUTPUT_TAKE], pci));

        if (badSkill(stTake))
            sprintf(strchr(sz, 0), " [%s]", gettext(aszSkillType[stTake]));

    }

    r = arDouble[OUTPUT_DROP] - arDouble[OUTPUT_TAKE];

    if (fDouble > 0 && !fTake && r > 0.0f) {

        fAnno = TRUE;

        /* wrong pass */

        sprintf(strchr(sz, 0), "%s (%s)!\n",
                _("Alert: wrong pass"), OutputEquityDiff(arDouble[OUTPUT_TAKE], arDouble[OUTPUT_DROP], pci));

        if (badSkill(stTake))
            sprintf(strchr(sz, 0), " [%s]", gettext(aszSkillType[stTake]));

    }


    if (arDouble[OUTPUT_TAKE] > arDouble[OUTPUT_DROP])
        r = arDouble[OUTPUT_NODOUBLE] - arDouble[OUTPUT_DROP];
    else
        r = arDouble[OUTPUT_NODOUBLE] - arDouble[OUTPUT_TAKE];

    if (fDouble > 0 && fTake < 0 && r > 0.0f) {

        fAnno = TRUE;

        /* wrong double */

        sprintf(strchr(sz, 0), "%s (%s)!\n",
                _("Alert: wrong double"),
                OutputEquityDiff((arDouble[OUTPUT_TAKE] >
                                  arDouble[OUTPUT_DROP]) ?
                                 arDouble[OUTPUT_DROP] : arDouble[OUTPUT_TAKE], arDouble[OUTPUT_NODOUBLE], pci));

        if (badSkill(stDouble))
            sprintf(strchr(sz, 0), " [%s]", gettext(aszSkillType[stDouble]));

    }

    if ((badSkill(stDouble) || badSkill(stTake)) && !fAnno) {

        if (badSkill(stDouble)) {
            sprintf(strchr(sz, 0), _("Alert: double decision marked %s"), gettext(aszSkillType[stDouble]));
            strcat(sz, "\n");
        }

        if (badSkill(stTake)) {
            sprintf(strchr(sz, 0), _("Alert: take decision marked %s"), gettext(aszSkillType[stTake]));
            strcat(sz, "\n");
        }

    }

    strcat(sz, OutputCubeAnalysis(aarOutput, aarStdDev, pes, pci, fTake));

    return sz;

}

extern char *
OutputCubeAnalysis(float aarOutput[2][NUM_ROLLOUT_OUTPUTS],
                   float aarStdDev[2][NUM_ROLLOUT_OUTPUTS], const evalsetup * pes, const cubeinfo * pci, int fTake)
{

    static char sz[4096];
    int i;
    float arDouble[4];
    const char *aszCube[] = {
        NULL,
        N_("No double"),
        N_("Double, take"),
        N_("Double, pass")
    };

    int ai[3];
    cubedecision cd;
    float r;

    FindCubeDecision(arDouble, aarOutput, pci);

    /* header */

    sprintf(sz, "\n%s\n", _("Cube analysis"));

    /* ply & cubeless equity */

    switch (pes->et) {
    case EVAL_NONE:
        strcat(sz, _("n/a"));
        break;
    case EVAL_EVAL:
        sprintf(strchr(sz, 0), "%u-%s", pes->ec.nPlies, _("ply"));
        break;
    case EVAL_ROLLOUT:
        strcat(sz, _("Rollout"));
        break;
    }

    if (fTake >= 0) { /* take or drop */
        InvertEvaluationR(aarOutput[0], pci);
        InvertEvaluationR(aarOutput[1], pci);
        if (pes->et == EVAL_ROLLOUT) {
            InvertStdDev(aarStdDev[0]);
            InvertStdDev(aarStdDev[1]);
        }
    }

    if (pci->nMatchTo)
        sprintf(strchr(sz, 0), " %s %s (%s: %s)\n",
                fOutputMWC ? _("cubeless MWC") : _("cubeless equity"),
                OutputEquity(aarOutput[0][OUTPUT_EQUITY], pci, TRUE),
                _("Money"), OutputMoneyEquity(aarOutput[0], TRUE));
    else
        sprintf(strchr(sz, 0), " %s %s\n", _("cubeless equity"), OutputMoneyEquity(aarOutput[0], TRUE));




    /* Output percentags for evaluations */

    if (exsExport.fCubeDetailProb && pes->et == EVAL_EVAL) {

        strcat(sz, "  ");
        strcat(sz, OutputPercents(aarOutput[0], TRUE));

    }

    strcat(sz, "\n");

    /* equities */

    strcat(sz, _("Cubeful equities"));
    strcat(sz, ":\n");
    if (pes->et == EVAL_EVAL && exsExport.afCubeParameters[0]) {
        strcat(sz, "  ");
        strcat(sz, OutputEvalContext(&pes->ec, FALSE));
        strcat(sz, "\n");
    }

    getCubeDecisionOrdering(ai, arDouble, aarOutput, pci);

    for (i = 0; i < 3; i++) {

        sprintf(strchr(sz, 0), "%d. %-20s", i + 1, gettext(aszCube[ai[i]]));

        if (fTake == -1) { /* double */
            strcat(sz, OutputEquity(arDouble[ai[i]], pci, TRUE));

            if (i)
                sprintf(strchr(sz, 0), "  (%s)", OutputEquityDiff(arDouble[ai[i]], arDouble[OUTPUT_OPTIMAL], pci));
        } else { /* take or drop */
            strcat(sz, OutputEquity(-arDouble[ai[i]], pci, TRUE));

            if (i)
                sprintf(strchr(sz, 0), "  (%s)", OutputEquityDiff(-arDouble[ai[i]], -arDouble[OUTPUT_OPTIMAL], pci));
        }
        strcat(sz, "\n");

    }

    /* cube decision */

    cd = FindBestCubeDecision(arDouble, aarOutput, pci);

    sprintf(strchr(sz, 0), "%s: %s", _("Proper cube action"), GetCubeRecommendation(cd));

    if ((r = getPercent(cd, arDouble)) >= 0.0f)
        sprintf(strchr(sz, 0), " (%.1f%%)", 100.0f * r);

    strcat(sz, "\n");

    /* dump rollout */

    if (pes->et == EVAL_ROLLOUT && exsExport.fCubeDetailProb) {

        char asz[2][FORMATEDMOVESIZE];
        cubeinfo aci[2];

        for (i = 0; i < 2; i++) {

            memcpy(&aci[i], pci, sizeof(cubeinfo));

            if (i) {
                aci[i].fCubeOwner = !pci->fMove;
                aci[i].nCube *= 2;
            }

            FormatCubePosition(asz[i], &aci[i]);

        }

        sprintf(strchr(sz, 0), "\n%s:\n", _("Rollout details"));

        strcat(strchr(sz, 0), OutputRolloutResult(NULL, asz, aarOutput, aarStdDev, aci, 0, 2, pes->rc.fCubeful));


    }

    if (fTake >= 0) { /* take or drop */
        InvertEvaluationR(aarOutput[0], pci);
        InvertEvaluationR(aarOutput[1], pci);
        if (pes->et == EVAL_ROLLOUT) {
            InvertStdDev(aarStdDev[0]);
            InvertStdDev(aarStdDev[1]);
        }
    }

    if (pes->et == EVAL_ROLLOUT && exsExport.afCubeParameters[1])
        strcat(strchr(sz, 0), OutputRolloutContext(NULL, &pes->rc));

    return sz;
}


extern void
FormatCubePositions(const cubeinfo * pci, char asz[2][FORMATEDMOVESIZE])
{

    cubeinfo aci[2];

    SetCubeInfo(&aci[0], pci->nCube, pci->fCubeOwner, pci->fMove,
                pci->nMatchTo, pci->anScore, pci->fCrawford, pci->fJacoby, pci->fBeavers, pci->bgv);

    FormatCubePosition(asz[0], &aci[0]);

    SetCubeInfo(&aci[1], 2 * pci->nCube, !pci->fMove, pci->fMove,
                pci->nMatchTo, pci->anScore, pci->fCrawford, pci->fJacoby, pci->fBeavers, pci->bgv);

    FormatCubePosition(asz[1], &aci[1]);

}

extern char *
FormatCubePosition(char *sz, cubeinfo * pci)
{
    /* FIXME: functions showing rollout progress can call either
     * FormatMove() from drawboard.c where we know that sz is 
     * shorter than FORMATEDMOVESIZE, or this one. Moreover they
     * depend on it to be this short.
     * Truncate the result if necessary to avoid crashes with long
     * player names or some translations. Remove the initial
     * "Player " to save a few characters.
     */
    if (pci->fCubeOwner == -1)
        snprintf(sz, FORMATEDMOVESIZE, _("Centered %d-cube"), pci->nCube);
    else
        snprintf(sz, FORMATEDMOVESIZE, _("%s owns %d-cube"), ap[pci->fCubeOwner].szName, pci->nCube);

    return sz;

}

static int
DumpOver(const TanBoard anBoard, char *pchOutput, const bgvariation bgv)
{

    float ar[NUM_OUTPUTS] = { 0, 0, 0, 0, 0 };  /* NUM_OUTPUTS is 5 */

    if (EvalOver(anBoard, ar, bgv, NULL))
        return -1;

    if (ar[OUTPUT_WIN] > 0.0f)
        strcpy(pchOutput, _("Win"));
    else
        strcpy(pchOutput, _("Loss"));
    strcat(pchOutput, " ");

    if (ar[OUTPUT_WINBACKGAMMON] > 0.0f || ar[OUTPUT_LOSEBACKGAMMON] > 0.0f)
        sprintf(pchOutput, "(%s)\n", _("backgammon"));
    else if (ar[OUTPUT_WINGAMMON] > 0.0f || ar[OUTPUT_LOSEGAMMON] > 0.0f)
        sprintf(pchOutput, "(%s)\n", _("gammon"));
    else
        sprintf(pchOutput, "(%s)\n", _("single"));

    return 0;

}


static int
DumpBearoff1(const TanBoard anBoard, char *szOutput, const bgvariation UNUSED(bgv))
{

    g_assert(pbc1);
    return BearoffDump(pbc1, anBoard, szOutput);

}

static int
DumpBearoff2(const TanBoard anBoard, char *szOutput, const bgvariation UNUSED(bgv))
{

    g_assert(pbc2);

    if (BearoffDump(pbc2, anBoard, szOutput))
        return -1;

    if (pbc1)
        if (BearoffDump(pbc1, anBoard, szOutput))
            return -1;

    return 0;

}


static int
DumpBearoffOS(const TanBoard anBoard, char *szOutput, const bgvariation UNUSED(bgv))
{

    g_assert(pbcOS);
    return BearoffDump(pbcOS, anBoard, szOutput);

}


static int
DumpBearoffTS(const TanBoard anBoard, char *szOutput, const bgvariation UNUSED(bgv))
{

    g_assert(pbcTS);
    return BearoffDump(pbcTS, anBoard, szOutput);

}


static int
DumpRace(const TanBoard UNUSED(anBoard), char *UNUSED(szOutput), const bgvariation UNUSED(bgv))
{

    /* no-op -- nothing much we can say, really (pip count?) */
    return 0;

}

static void
DumpAnyContact(const TanBoard UNUSED(anBoard), char *UNUSED(szOutput),
               const bgvariation UNUSED(bgv), int UNUSED(isCrashed))
{
    return;
}

static int
DumpContact(const TanBoard anBoard, char *szOutput, const bgvariation bgv)
{
    DumpAnyContact(anBoard, szOutput, bgv, 0);
    return 0;
}

static int
DumpCrashed(const TanBoard anBoard, char *szOutput, const bgvariation bgv)
{
    DumpAnyContact(anBoard, szOutput, bgv, 1);
    return 0;
}

static int
DumpHypergammon1(const TanBoard anBoard, char *szOutput, const bgvariation UNUSED(bgv))
{

    g_assert(apbcHyper[0]);
    return BearoffDump(apbcHyper[0], anBoard, szOutput);

}

static int
DumpHypergammon2(const TanBoard anBoard, char *szOutput, const bgvariation UNUSED(bgv))
{

    g_assert(apbcHyper[1]);
    return BearoffDump(apbcHyper[1], anBoard, szOutput);

}

static int
DumpHypergammon3(const TanBoard anBoard, char *szOutput, const bgvariation UNUSED(bgv))
{

    g_assert(apbcHyper[2]);
    return BearoffDump(apbcHyper[2], anBoard, szOutput);

}

static classdumpfunc acdf[N_CLASSES] = {
    DumpOver,
    DumpHypergammon1, DumpHypergammon2, DumpHypergammon3,
    DumpBearoff2, DumpBearoffTS,
    DumpBearoff1, DumpBearoffOS,
    DumpRace, DumpCrashed, DumpContact
};

extern int
DumpPosition(const TanBoard anBoard, char *szOutput,
             const evalcontext * pec, cubeinfo * pci, int fOutputMWC,
             int UNUSED(fOutputWinPC), int fOutputInvert, const char *szMatchID)
{

    float aarOutput[2][NUM_ROLLOUT_OUTPUTS];
    positionclass pc = ClassifyPosition(anBoard, pci->bgv);
    int i, nPlies;
    int j;
    evalcontext ec;
    static const char *aszEvaluator[] = {
        N_("Over"),
        N_("Hypergammon-1"),
        N_("Hypergammon-2"),
        N_("Hypergammon-3"),
        N_("Bearoff2"),
        N_("Bearoff-TS"),
        N_("Bearoff1"),
        N_("Bearoff-OS"),
        N_("Race"),
        N_("Crashed"),
        N_("Contact")
    };

    strcpy(szOutput, "");

    sprintf(strchr(szOutput, 0), "%s:\t", _("Position ID"));
    strcat(szOutput, PositionID(anBoard));
    strcat(szOutput, "\n");
    if (szMatchID) {
        sprintf(strchr(szOutput, 0), "%s:\t", _("Match ID"));
        strcat(szOutput, szMatchID);
        strcat(szOutput, "\n");
    }
    strcat(szOutput, "\n");

    sprintf(strchr(szOutput, 0), "%s: \t", _("Evaluator"));
    strcat(szOutput, gettext(aszEvaluator[pc]));
    strcat(szOutput, "\n\n");
    acdf[pc] (anBoard, strchr(szOutput, 0), pci->bgv);
    szOutput = strchr(szOutput, 0);

    sprintf(strchr(szOutput, 0),
            "\n"
            "        %-7s %-7s %-7s %-7s %-7s %-9s %-9s\n",
            _("Win"), _("W(g)"), _("W(bg)"), _("L(g)"), _("L(bg)"),
            (!pci->nMatchTo || !fOutputMWC) ? _("Equity") : _("MWC"), _("Cubeful"));

    nPlies = pec->nPlies > 9 ? 9 : pec->nPlies;

    memcpy(&ec, pec, sizeof(evalcontext));

    for (i = 0; i <= nPlies; i++) {
        szOutput = strchr(szOutput, 0);


        ec.nPlies = i;

        if (GeneralCubeDecisionE(aarOutput, anBoard, pci, &ec, 0) < 0)
            return -1;

        if (!i)
            strcpy(szOutput, _("static"));
        else
            sprintf(szOutput, "%2d %s", i, _("ply"));

        szOutput = strchr(szOutput, 0);

        if (fOutputInvert) {
            InvertEvaluationR(aarOutput[0], pci);
            InvertEvaluationR(aarOutput[1], pci);
            pci->fMove = !pci->fMove;
        }

        /* Print %'s and equities */

        strcat(szOutput, ": ");

        for (j = 0; j < 5; ++j) {
            sprintf(strchr(szOutput, 0), "%-7s ", OutputPercent(aarOutput[0][j]));
        }

        if (pci->nMatchTo)
            sprintf(strchr(szOutput, 0), "%-9s ", OutputEquity(Utility(aarOutput[0], pci), pci, TRUE));
        else
            sprintf(strchr(szOutput, 0), "%-9s ", OutputMoneyEquity(aarOutput[0], TRUE));

        sprintf(strchr(szOutput, 0), "%-9s ", OutputMWC(aarOutput[0][6], pci, TRUE));

        strcat(szOutput, "\n");

        if (fOutputInvert) {
            pci->fMove = !pci->fMove;
        }
    }

    /* if cube is available, output cube action */
    if (GetDPEq(NULL, NULL, pci)) {

        evalsetup es;

        es.et = EVAL_EVAL;
        es.ec = *pec;

        strcat(szOutput, "\n\n");
        strcat(szOutput, OutputCubeAnalysis(aarOutput, NULL, &es, pci, -1));

    }

    return 0;
}