fluid_dynamics.cc

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file fluid_dynamics.cc

/*******************************************************************************
 * Copyright (c) The JETSCAPE Collaboration, 2018
 *
 * Modular, task-based framework for simulating all aspects of heavy-ion collisions
 * 
 * For the list of contributors see AUTHORS.
 *
 * Report issues at https://github.com/JETSCAPE/JETSCAPE/issues
 *
 * or via email to bugs.jetscape@gmail.com
 *
 * Distributed under the GNU General Public License 3.0 (GPLv3 or later).
 * See COPYING for details.
 ******************************************************************************/

#include "FluidDynamics.h"
#include "FluidEvolutionHistory.h"
#include "gtest/gtest.h"

using namespace Jetscape;

void test_not_in_range(EvolutionHistory hist, real tau, real x, real y, real eta) {
    try {
        hist.CheckInRange(tau, x, y, eta);
    } catch (InvalidSpaceTimeRange & e) {
        auto estr = std::string(e.what());
        ASSERT_TRUE(estr.find("not in range") != estr.npos);
    }
}

// test EvolutionHistory Class
TEST(EvolutionHistoryTest, TEST_WRITE){
    auto hist = EvolutionHistory();
    hist.tau_min = 0.6;
    hist.dtau = 0.1;
    hist.x_min = -10;
    hist.y_min = -10;
    hist.eta_min = -10;
    hist.dx = 0.5;
    hist.dy = 0.5;
    hist.deta = 0.5;
    hist.ntau = 11;
    hist.nx = 41;
    hist.ny = 41;
    hist.neta = 41;
    hist.tau_eta_is_tz = false;

    EXPECT_EQ(hist.TauMax(), static_cast<real>(1.6));
    EXPECT_EQ(hist.XMax(), static_cast<real>(10.0));
    EXPECT_EQ(hist.YMax(), static_cast<real>(10.0));
    EXPECT_EQ(hist.EtaMax(), static_cast<real>(10.0));

    // check range test
    test_not_in_range(hist, 0.5, 0.3, 0.3, 0.3);
    test_not_in_range(hist, 20.5, 0.3, 0.3, 0.3);
    test_not_in_range(hist, 10.5, 10.3, 0.3, 0.3);
    test_not_in_range(hist, 10.5, 10., 10.3, 0.3);
    test_not_in_range(hist, 10.5, 10., 10., 10.3);

    real const_ed = 0.8;
    for (int n=0; n != hist.ntau; n++)
        for (int i=0; i != hist.nx; i++)
            for (int j=0; j != hist.ny; j++)
                for (int k=0; k != hist.neta; k++) {
                    auto cell = FluidCellInfo();
                    cell.energy_density = const_ed;
                    hist.data.emplace_back(std::move(cell));
                }
    // check almost equal for two float numbers
    ASSERT_NEAR(hist.get(0.8, 0.0, 0.0, 0.0).energy_density, static_cast<real>(const_ed), 1.0E-6);
}