// produces a list of vein materials available on the map. can be run with '-a' modifier to show even unrevealed minerals deep underground
// with -b modifier, it will show base layer materials too

// TODO: use material colors to make the output prettier
// TODO: needs the tiletype filter!
// TODO: tile override materials
// TODO: material types, trees, ice, constructions
// TODO: GUI

#include <iostream>
#include <string.h> // for memset
#include <string>
#include <vector>
#include <map>
#include <stdio.h>
#include <algorithm>
using namespace std;

#include <DFHack.h>
#include <dfhack/DFTileTypes.h>

template<template <typename> class P = std::less >
struct compare_pair_first
{
    template<class T1, class T2>
    bool operator()(const std::pair<T1, T2>& left, const std::pair<T1, T2>& right)
    {
        return P<T1>()(left.first, right.first);
    }
};

template<template <typename> class P = std::less >
struct compare_pair_second
{
    template<class T1, class T2>
    bool operator()(const std::pair<T1, T2>& left, const std::pair<T1, T2>& right)
    {
        return P<T2>()(left.second, right.second);
    }
};

int main (int argc, const char* argv[])
{
    bool showhidden = false;
    bool showbaselayers = false;
    for(int i = 1; i < argc; i++)
    {
        string test = argv[i];
        if(test == "-a")
        {
            showhidden = true;
        }
        else if(test == "-b")
        {
            showbaselayers = true;
        }
        else if(test == "-ab" || test == "-ba")
        {
            showhidden = true;
            showbaselayers = true;
        }
    }
    // let's be more useful when double-clicked on windows
    #ifndef LINUX_BUILD
    showhidden = true;
    #endif
    uint32_t x_max,y_max,z_max;
    /*
    DFHack::tiletypes40d tiletypes;
    DFHack::designations40d designations;
    DFHack::biome_indices40d regionoffsets;
    */
    DFHack::mapblock40d Block;
    map <int16_t, uint32_t> materials;
    materials.clear();
    vector<DFHack::t_feature> global_features;
    std::map <DFHack::planecoord, std::vector<DFHack::t_feature *> > local_features;

    vector< vector <uint16_t> > layerassign;

    DFHack::ContextManager DFMgr("Memory.xml");
    DFHack::Context *DF;
    try
    {
        DF = DFMgr.getSingleContext();
        DF->Attach();
    }
    catch (exception& e)
    {
        cerr << e.what() << endl;
        #ifndef LINUX_BUILD
            cin.ignore();
        #endif
        return 1;
    }

    DFHack::Maps * Maps = DF->getMaps();
    DFHack::Materials * Mats = DF->getMaterials();

    // init the map
    if(!Maps->Start())
    {
        cerr << "Can't init map." << endl;
        #ifndef LINUX_BUILD
            cin.ignore();
        #endif
        return 1;
    }
    Maps->getSize(x_max,y_max,z_max);

    if(!Maps->ReadGlobalFeatures(global_features))
    {
        cerr << "Can't get global features." << endl;
        #ifndef LINUX_BUILD
            cin.ignore();
        #endif
        return 1; 
    }

    if(!Maps->ReadLocalFeatures(local_features))
    {
        cerr << "Can't get local features." << endl;
        #ifndef LINUX_BUILD
            cin.ignore();
        #endif
        return 1; 
    }
    // get stone matgloss mapping
    if(!Mats->ReadInorganicMaterials())
    {
        //DF.DestroyMap();
        cerr << "Can't get the materials." << endl;
        #ifndef LINUX_BUILD
            cin.ignore();
        #endif
        return 1; 
    }

    // get region geology
    if(showbaselayers && !Maps->ReadGeology( layerassign ))
    {
        cerr << "Can't get region geology." << endl;
        #ifndef LINUX_BUILD
            cin.ignore();
        #endif
        return 1; 
    }

    int16_t tempvein [16][16];
    vector <DFHack::t_vein> veins;
    uint32_t maximum_regionoffset = 0;
    uint32_t num_overflows = 0;
    // walk the map!
    for(uint32_t x = 0; x< x_max;x++)
    {
        for(uint32_t y = 0; y< y_max;y++)
        {
            for(uint32_t z = 0; z< z_max;z++)
            {
                if(!Maps->isValidBlock(x,y,z))
                    continue;

                // read data
                Maps->ReadBlock40d(x,y,z, &Block);
                //Maps->ReadTileTypes(x,y,z, &tiletypes);
                //Maps->ReadDesignations(x,y,z, &designations);

                memset(tempvein, -1, sizeof(tempvein));
                veins.clear();
                Maps->ReadVeins(x,y,z,&veins);

                if(showbaselayers)
                {
                    //Maps->ReadRegionOffsets(x,y,z, &regionoffsets);
                    // get the layer materials
                    for(uint32_t xx = 0;xx<16;xx++)
                    {
                        for (uint32_t yy = 0; yy< 16;yy++)
                        {
                            uint8_t test = Block.designation[xx][yy].bits.biome;
                            if(test > maximum_regionoffset)
                                maximum_regionoffset = test;
                            if( test >= sizeof(Block.biome_indices))
                            {
                                num_overflows++;
                                continue;
                            }
                            tempvein[xx][yy] =
                            layerassign
                            [Block.biome_indices[test]]
                            [Block.designation[xx][yy].bits.geolayer_index];
                        }
                    }
                }

                // for each vein
                for(int i = 0; i < (int)veins.size();i++)
                {
                    //iterate through vein rows
                    for(uint32_t j = 0;j<16;j++)
                    {
                        //iterate through the bits
                        for (uint32_t k = 0; k< 16;k++)
                        {
                            // and the bit array with a one-bit mask, check if the bit is set
                            bool set = !!(((1 << k) & veins[i].assignment[j]) >> k);
                            if(set)
                            {
                                // store matgloss
                                tempvein[k][j] = veins[i].type;
                            }
                        }
                    }
                }

                // global feature overrides
                int16_t idx = Block.global_feature;
                if( idx != -1 && uint16_t(idx) < global_features.size() && global_features[idx].type == DFHack::feature_Underworld)
                {
                    for(uint32_t xi = 0 ; xi< 16 ; xi++) for(uint32_t yi = 0 ; yi< 16 ; yi++)
                    {
                        if(Block.designation[xi][yi].bits.feature_global)
                        {
                            if(global_features[idx].main_material == 0) // stone
                            {
                                tempvein[xi][yi] = global_features[idx].sub_material;
                            }
                            else
                            {
                                tempvein[xi][yi] = -1;
                            }
                        }
                    }
                }

                idx = Block.local_feature;
                if( idx != -1 )
                {
                    DFHack::planecoord pc;
                    pc.dim.x = x;
                    pc.dim.y = y;
                    std::map <DFHack::planecoord, std::vector<DFHack::t_feature *> >::iterator it;
                    it = local_features.find(pc);
                    if(it != local_features.end())
                    {
                        std::vector<DFHack::t_feature *>& vectr = (*it).second;
                        if(uint16_t(idx) < vectr.size() && vectr[idx]->type == DFHack::feature_Adamantine_Tube)
                            for(uint32_t xi = 0 ; xi< 16 ; xi++) for(uint32_t yi = 0 ; yi< 16 ; yi++)
                            {
                                if(Block.designation[xi][yi].bits.feature_local && DFHack::isWallTerrain(Block.tiletypes[xi][yi]))
                                {
                                    if(vectr[idx]->main_material == 0) // stone
                                    {
                                        tempvein[xi][yi] = vectr[idx]->sub_material;
                                    }
                                    else
                                    {
                                        tempvein[xi][yi] = -1;
                                    }
                                }
                            }
                    }
                }

                // count the material types
                for(uint32_t xi = 0 ; xi< 16 ; xi++)
                {
                    for(uint32_t yi = 0 ; yi< 16 ; yi++)
                    {
                        // hidden tiles are ignored unless '-a' is provided on the command line
                        // non-wall tiles are ignored
                        if( (Block.designation[xi][yi].bits.hidden && !showhidden) || !DFHack::isWallTerrain(Block.tiletypes[xi][yi]))
                            continue;
                        if(tempvein[xi][yi] < 0)
                            continue;

                        if(materials.count(tempvein[xi][yi]))
                        {
                            materials[tempvein[xi][yi]] += 1;
                        }
                        else
                        {
                            materials[tempvein[xi][yi]] = 1;
                        }
                    }
                }
            }
        }
    }
    // print report
    cout << "Maximal regionoffset seen: " << maximum_regionoffset << ".";
    if(maximum_regionoffset >= sizeof(Block.biome_indices) )
    {
        cout << " This is above the regionoffsets array size!" << endl;
        cout << "Number of overflows: " << num_overflows;
    }
    cout << endl;
    vector <pair <int16_t, uint32_t> > matss;
    map<int16_t, uint32_t>::iterator p;
    for(p = materials.begin(); p != materials.end(); p++)
    {
        if(p->first == -1)
        {
            cout << "Non-stone" << " : " << p->second << endl;
        }
        else
        {
            matss.push_back( pair<int16_t,uint32_t>(p->first, p->second) );
        }
    }
    std::sort(matss.begin(), matss.end(), compare_pair_second<>());
    for(size_t i = 0; i < matss.size();i++)
    {
        if(matss[i].first >= Mats->inorganic.size())
        {
            cerr << "Error, material out of bounds: " << matss[i].first << endl;
            continue;
        }
        cout << Mats->inorganic[matss[i].first].id << " : " << matss[i].second << endl;
    }
    DF->Detach();
    #ifndef LINUX_BUILD
        cout << "Done. Press any key to continue" << endl;
        cin.ignore();
    #endif
    return 0;
}