src/preequilibrium/NullPreDynamics.cc

0001 /*******************************************************************************
0002  * Copyright (c) The JETSCAPE Collaboration, 2018
0003  *
0004  * Modular, task-based framework for simulating all aspects of heavy-ion collisions
0005  *
0006  * For the list of contributors see AUTHORS.
0007  *
0008  * Report issues at https://github.com/JETSCAPE/JETSCAPE/issues
0009  *
0010  * or via email to bugs.jetscape@gmail.com
0011  *
0012  * Distributed under the GNU General Public License 3.0 (GPLv3 or later).
0013  * See COPYING for details.
0014  ******************************************************************************/
0015
0016 #include <stdio.h>
0017 #include <sys/stat.h>
0018
0019 #include <cstring>
0020
0021 #include "JetScapeLogger.h"
0022 #include "NullPreDynamics.h"
0023
0024 // Register the module with the base class
0025 RegisterJetScapeModule<NullPreDynamics> NullPreDynamics::reg("NullPreDynamics");
0026
0027 NullPreDynamics::NullPreDynamics() {
0028   preequilibrium_status_ = NOT_STARTED;
0029   SetId("NullPreDynamics");
0030 }
0031
0032 void NullPreDynamics::EvolvePreequilibrium() {
0033   VERBOSE(2) << "Initialize energy density profile in NullPreDynamics ...";
0034   // grab initial energy density from vector from initial state module
0035   std::vector<double> energy_density = ini->GetEntropyDensityDistribution();
0036   preequilibrium_status_ = INIT;
0037   if (preequilibrium_status_ == INIT) {
0038     VERBOSE(2) << "running NullPreDynamics ...";
0039     for (auto const &e_local : energy_density) {
0040       e_.push_back(e_local);
0041       P_.push_back(e_local / 3.);
0042       utau_.push_back(1.);
0043       ux_.push_back(0.);
0044       uy_.push_back(0.);
0045       ueta_.push_back(0.);
0046       pi00_.push_back(0.);
0047       pi01_.push_back(0.);
0048       pi02_.push_back(0.);
0049       pi03_.push_back(0.);
0050       pi11_.push_back(0.);
0051       pi12_.push_back(0.);
0052       pi13_.push_back(0.);
0053       pi22_.push_back(0.);
0054       pi23_.push_back(0.);
0055       pi33_.push_back(0.);
0056       bulk_Pi_.push_back(0.);
0057     }
0058     preequilibrium_status_ = DONE;
0059   }
0060 }