mercury2pvd.cpp

Go to the documentation of this file.

//Copyright (c) 2015, The MercuryDPM Developers Team. All rights reserved.
//For the list of developers, see <http://www.MercuryDPM.org/Team>.
//
//Redistribution and use in source and binary forms, with or without
//modification, are permitted provided that the following conditions are met:
//  * Redistributions of source code must retain the above copyright
//    notice, this list of conditions and the following disclaimer.
//  * Redistributions in binary form must reproduce the above copyright
//    notice, this list of conditions and the following disclaimer in the
//    documentation and/or other materials provided with the distribution.
//  * Neither the name MercuryDPM nor the
//    names of its contributors may be used to endorse or promote products
//    derived from this software without specific prior written permission.
//
//THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND
//ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
//WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
//DISCLAIMED. IN NO EVENT SHALL THE MERCURYDPM DEVELOPERS TEAM BE LIABLE FOR ANY
//DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
//(INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
//LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND
//ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
//(INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
//SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.

#include "DPMBase.h"

class Mercury2PVD : public DPMBase {

public:

    void setupInitialConditions()
    {
    }
};

int main(int argc UNUSED, char *argv[] UNUSED)
{
    Mercury2PVD m;

    //Check to see if we received one argument
    if (argc == 2) {
        std::string prefix = argv[1];
        std::string::size_type pos = prefix.find(".restart");

        if (pos == std::string::npos) {
            m.setName(argv[1]);
        } else {
            m.setName(prefix.substr(0,pos-1));
            m.restartFile.setName(prefix);
            logger(INFO, "restart file % given; prefix is %",m.restartFile.getName(),m.getName());
        }

    } else {
        m.setName("../Drivers/MercuryCG/NewtonsCradleSelfTest");
        //We didn't. Print a usage and exit the program.
        logger(WARN, "Usage: mercury2pvd [prefix] \n"
         "  This program converts MercuryDPM files into ParaView data files,\n"
         "  which can then be used to visualize your particles.\n"
         "\n"
         "   prefix: Prefix to prepend to the i/o files.\n"
         "\n"
         "   The following files will be generated:\n"
         "     - prefix.pvd\n"
         "     - prefix_0.vtu\n"
         "     - prefix_1.vtu\n"
         "         ( ... )\n"
         "     - prefix_987654321.vtu\n"
         "     depending on the amount of timesteps.", argv[0]);
    }

    //Store path as it is overwritten by readRestartFile
    std::string::size_type pos = m.getName().find_last_of("/");
    std::string path = m.getName().substr(0,pos+1);
    logger(INFO,"Path %",path);

    //Read restart file, and re-add path to data and fstat file name
    logger(INFO,"Reading restart file %",m.restartFile.getName());
    bool success = m.readRestartFile();
    if (!success) logger(FATAL,"Reading restart file % failed",m.restartFile.getName());
    m.dataFile.setName(path+m.dataFile.getName());
    m.fStatFile.setName(path+m.fStatFile.getName());

    //prepare dpm class for reading in
    m.dataFile.open(std::fstream::in);
    m.fStatFile.open(std::fstream::in);
    m.interactionHandler.clear();
    m.dataFile.setCounter(0);
    m.fStatFile.setCounter(0);
    m.interactionHandler.setWriteVTK(FileType::MULTIPLE_FILES);

    //read in
    logger(INFO,"Reading data/fstat files %, %",m.dataFile.getName(), m.fStatFile.getName());
    while (m.readNextDataFile())
    {
        //m.readNextFStatFile();
        logger(INFO,"t % N % C %",m.getTime(),m.particleHandler.getNumberOfObjects(),m.interactionHandler.getNumberOfObjects());
        //m.writeOutputFiles()
        m.particleHandler.writeVTK();
        //m.interactionHandler.writeVTK();
        m.wallHandler.writeVTK();
    }
}