--- trunk/src/applications/dynamicProps/DynamicProps.cpp 2005/04/15 22:04:00 507 +++ trunk/src/applications/dynamicProps/DynamicProps.cpp 2008/05/27 16:39:06 1245 @@ -52,7 +52,13 @@ #include "applications/dynamicProps/DipoleCorrFunc.hpp" #include "applications/dynamicProps/RCorrFunc.hpp" #include "applications/dynamicProps/VCorrFunc.hpp" +#include "applications/dynamicProps/LegendreCorrFunc.hpp" +#include "applications/dynamicProps/RadialRCorrFunc.hpp" +#include "applications/dynamicProps/ThetaCorrFunc.hpp" +#include "applications/dynamicProps/DirectionalRCorrFunc.hpp" +#include "applications/dynamicProps/EnergyCorrFunc.hpp" + using namespace oopse; int main(int argc, char* argv[]){ @@ -70,9 +76,6 @@ int main(int argc, char* argv[]){ //get the dumpfile name and meta-data file name std::string dumpFileName = args_info.input_arg; - - std::string mdFileName = dumpFileName.substr(0, dumpFileName.rfind(".")) + ".md"; - std::string sele1; std::string sele2; @@ -105,7 +108,7 @@ int main(int argc, char* argv[]){ //parse md file and set up the system SimCreator creator; - SimInfo* info = creator.createSim(mdFileName, false); + SimInfo* info = creator.createSim(dumpFileName, false); TimeCorrFunc* corrFunc; @@ -113,9 +116,29 @@ int main(int argc, char* argv[]){ corrFunc = new DipoleCorrFunc(info, dumpFileName, sele1, sele2); } else if (args_info.rcorr_given) { corrFunc = new RCorrFunc(info, dumpFileName, sele1, sele2); - } - else if (args_info.vcorr_given) { + } else if (args_info.r_rcorr_given) { + corrFunc = new RadialRCorrFunc(info, dumpFileName, sele1, sele2); + } else if (args_info.thetacorr_given) { + corrFunc = new ThetaCorrFunc(info, dumpFileName, sele1, sele2); + } else if (args_info.drcorr_given) { + corrFunc = new DirectionalRCorrFunc(info, dumpFileName, sele1, sele2); + } else if (args_info.vcorr_given) { corrFunc = new VCorrFunc(info, dumpFileName, sele1, sele2); + } else if (args_info.helfandEcorr_given){ + corrFunc = new EnergyCorrFunc(info, dumpFileName, sele1, sele2); + } else if (args_info.lcorr_given) { + int order; + if (args_info.order_given) + order = args_info.order_arg; + else { + sprintf( painCave.errMsg, + "--order must be set if --lcoor is set\n"); + painCave.severity = OOPSE_ERROR; + painCave.isFatal = 1; + simError(); + } + + corrFunc = new LegendreCorrFunc(info, dumpFileName, sele1, sele2, order); } if (args_info.output_given) {