/usr/include/vtk-5.10/Cosmo/FOFHaloProperties.h is in libvtk5-dev 5.10.1+dfsg-2.1build1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 | /*=========================================================================
Copyright (c) 2007, Los Alamos National Security, LLC
All rights reserved.
Copyright 2007. Los Alamos National Security, LLC.
This software was produced under U.S. Government contract DE-AC52-06NA25396
for Los Alamos National Laboratory (LANL), which is operated by
Los Alamos National Security, LLC for the U.S. Department of Energy.
The U.S. Government has rights to use, reproduce, and distribute this software.
NEITHER THE GOVERNMENT NOR LOS ALAMOS NATIONAL SECURITY, LLC MAKES ANY WARRANTY,
EXPRESS OR IMPLIED, OR ASSUMES ANY LIABILITY FOR THE USE OF THIS SOFTWARE.
If software is modified to produce derivative works, such modified software
should be clearly marked, so as not to confuse it with the version available
from LANL.
Additionally, 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 of Los Alamos National Security, LLC, Los Alamos National
Laboratory, LANL, the U.S. Government, 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 LOS ALAMOS NATIONAL SECURITY, LLC 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 LOS ALAMOS NATIONAL SECURITY, LLC OR
CONTRIBUTORS 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.
=========================================================================*/
// .NAME FOFHaloProperties - calculate properties of all FOF halos
//
// FOFHaloProperties takes data from CosmoHaloFinderP about individual halos
// and data from all particles and calculates properties.
//
#ifndef FOFHaloProperties_h
#define FOFHaloProperties_h
#ifdef USE_VTK_COSMO
#include "CosmoDefinition.h"
#else
#include "Definition.h"
#endif
#include "ChainingMesh.h"
#include <string>
#include <vector>
using namespace std;
#ifdef USE_VTK_COSMO
class COSMO_EXPORT FOFHaloProperties {
#else
class FOFHaloProperties {
#endif
public:
FOFHaloProperties();
~FOFHaloProperties();
// Set parameters for sizes of the dead/alive space
void setParameters(
const string& outName, // Base name of output halo files
POSVEL_T rL, // Box size of the physical problem
POSVEL_T deadSize, // Dead size used to normalize for non periodic
POSVEL_T bb); // Inter particle distance for halos
// Set alive particle vectors which were created elsewhere
void setParticles(
vector<POSVEL_T>* xLoc,
vector<POSVEL_T>* yLoc,
vector<POSVEL_T>* zLoc,
vector<POSVEL_T>* xVel,
vector<POSVEL_T>* yVel,
vector<POSVEL_T>* zVel,
vector<POSVEL_T>* pmass,
vector<POTENTIAL_T>* potential,
vector<ID_T>* id,
vector<MASK_T>* mask,
vector<STATUS_T>* state);
void setParticles(
long count,
POSVEL_T* xLoc,
POSVEL_T* yLoc,
POSVEL_T* zLoc,
POSVEL_T* xVel,
POSVEL_T* yVel,
POSVEL_T* zVel,
POSVEL_T* pmass,
ID_T* id);
// Set the halo information from the FOF halo finder
void setHalos(
int numberOfHalos, // Number of halos found
int* halos, // Index into haloList of first particle
int* haloCount, // Number of particles in the matching halo
int* haloList); // Chain of indices of all particles in halo
// Find the halo centers (minimum potential) finding minimum of array
void FOFHaloCenterMinimumPotential(vector<int>* haloCenter);
// Find the mass of each halo by accumulating individual particle masses
void FOFHaloMass(
vector<POSVEL_T>* haloMass);
// Find the average position of FOF halo particles
void FOFPosition(
vector<POSVEL_T>* xPos,
vector<POSVEL_T>* yPos,
vector<POSVEL_T>* zPos);
// Find the center of mass of FOF halo particles
void FOFCenterOfMass(
vector<POSVEL_T>* xCofMass,
vector<POSVEL_T>* yCofMass,
vector<POSVEL_T>* zCofMass);
// Find the average velocity of FOF halo particles
void FOFVelocity(
vector<POSVEL_T>* xVel,
vector<POSVEL_T>* yVel,
vector<POSVEL_T>* zVel);
// Find the velocity dispersion of FOF halos
void FOFVelocityDispersion(
vector<POSVEL_T>* xVel,
vector<POSVEL_T>* yVel,
vector<POSVEL_T>* zVel,
vector<POSVEL_T>* velDisp);
// Kahan summation of floating point numbers to reduce roundoff error
POSVEL_T KahanSummation(int halo, POSVEL_T* data);
POSVEL_T KahanSummation2(int halo, POSVEL_T* data1, POSVEL_T* data2);
// Dot product
POSVEL_T dotProduct(POSVEL_T x, POSVEL_T y, POSVEL_T z);
// Incremental mean, possibly needed for very large halos
POSVEL_T incrementalMean(int halo, POSVEL_T* data);
// Extract locations and tags for all particles in a halo
void extractLocation(
int halo,
int* actualIndx,
POSVEL_T* xLocHalo,
POSVEL_T* yLocHalo,
POSVEL_T* zLocHalo,
ID_T* tag);
void extractInformation(
int halo,
int* actualIndx,
POSVEL_T* xLocHalo,
POSVEL_T* yLocHalo,
POSVEL_T* zVelHalo,
POSVEL_T* xVelHalo,
POSVEL_T* yVelHalo,
POSVEL_T* zLocHalo,
POSVEL_T* pmass,
ID_T* tag);
#ifndef USE_VTK_COSMO
// Print information about halos for debugging and selection
void FOFHaloCatalog(
vector<int>* haloCenter,
vector<POSVEL_T>* haloMass,
vector<POSVEL_T>* xVel,
vector<POSVEL_T>* yVel,
vector<POSVEL_T>* zVel);
void printHaloSizes(int minSize);
void printLocations(int haloIndex);
void printBoundingBox(int haloIndex);
#endif
private:
int myProc; // My processor number
int numProc; // Total number of processors
string outFile; // File of particles written by this processor
POSVEL_T boxSize; // Physical box size of the data set
POSVEL_T deadSize; // Border size for dead particles
POSVEL_T bb; // Interparticle distance for halos
long particleCount; // Total particles on this processor
POSVEL_T* xx; // X location for particles on this processor
POSVEL_T* yy; // Y location for particles on this processor
POSVEL_T* zz; // Z location for particles on this processor
POSVEL_T* vx; // X velocity for particles on this processor
POSVEL_T* vy; // Y velocity for particles on this processor
POSVEL_T* vz; // Z velocity for particles on this processor
POSVEL_T* mass; // mass of particles on this processor
POTENTIAL_T* pot; // Particle potential
ID_T* tag; // Id tag for particles on this processor
MASK_T* mask; // Particle information
STATUS_T* status; // Particle is ALIVE or labeled with neighbor
// processor index where it is ALIVE
// Information about halos from FOF halo finder
int numberOfHalos; // Number of halos found
int* halos; // First particle index into haloList
int* haloCount; // Size of each halo
int* haloList; // Indices of next particle in halo
};
#endif
|