| 57 |
|
#include "applications/staticProps/GofXyz.hpp" |
| 58 |
|
#include "applications/staticProps/P2OrderParameter.hpp" |
| 59 |
|
#include "applications/staticProps/SCDOrderParameter.hpp" |
| 60 |
+ |
#include "applications/staticProps/DensityPlot.hpp" |
| 61 |
+ |
#include "applications/staticProps/RhoZ.hpp" |
| 62 |
+ |
#if defined(HAVE_FFTW_H) || defined(HAVE_DFFTW_H) || defined(HAVE_FFTW3_H) |
| 63 |
+ |
#include "applications/staticProps/Hxy.hpp" |
| 64 |
+ |
#endif |
| 65 |
|
|
| 66 |
|
using namespace oopse; |
| 67 |
|
|
| 93 |
|
char* sele1Env= getenv("OOPSE_SELE1"); |
| 94 |
|
if (sele1Env) { |
| 95 |
|
sele1 = sele1Env; |
| 96 |
< |
}else { |
| 96 |
> |
}else if (!args_info.scd_given) { |
| 97 |
|
sprintf( painCave.errMsg, |
| 98 |
|
"neither --sele1 option nor $OOPSE_SELE1 is set"); |
| 99 |
|
painCave.severity = OOPSE_ERROR; |
| 108 |
|
char* sele2Env = getenv("OOPSE_SELE2"); |
| 109 |
|
if (sele2Env) { |
| 110 |
|
sele2 = sele2Env; |
| 111 |
< |
} else { |
| 111 |
> |
} else if (args_info.density_given) { |
| 112 |
> |
sele2 = "select all"; |
| 113 |
> |
} else if(!args_info.scd_given && !args_info.density_given && !args_info.slab_density_given) { |
| 114 |
|
sprintf( painCave.errMsg, |
| 115 |
|
"neither --sele2 option nor $OOPSE_SELE2 is set"); |
| 116 |
|
painCave.severity = OOPSE_ERROR; |
| 185 |
|
std::string sele3 = args_info.sele3_arg; |
| 186 |
|
analyser = new SCDOrderParameter(info, dumpFileName, sele1, sele2, sele3); |
| 187 |
|
} |
| 188 |
+ |
}else if (args_info.density_given) { |
| 189 |
+ |
analyser= new DensityPlot(info, dumpFileName, sele1, sele2, maxLen, args_info.nrbins_arg); |
| 190 |
+ |
} else if (args_info.slab_density_given) { |
| 191 |
+ |
Mat3x3d hmat = info->getSnapshotManager()->getCurrentSnapshot()->getHmat(); |
| 192 |
+ |
analyser = new RhoZ(info, dumpFileName, sele1, hmat(2, 2), args_info.nrbins_arg); |
| 193 |
+ |
#if defined(HAVE_FFTW_H) || defined(HAVE_DFFTW_H) || defined(HAVE_FFTW3_H) |
| 194 |
+ |
}else if (args_info.hxy_given) { |
| 195 |
+ |
analyser = new Hxy(info, dumpFileName, sele1, args_info.nbins_x_arg, args_info.nbins_y_arg, args_info.nrbins_arg); |
| 196 |
+ |
#endif |
| 197 |
|
} |
| 198 |
|
|
| 199 |
|
if (args_info.output_given) { |