field.hpp

#ifndef ABLATELIBRARY_PARTICLEFIELD_HPP
#define ABLATELIBRARY_PARTICLEFIELD_HPP
 
#include <petsc.h>
#include <string>
#include <vector>
#include "domain/field.hpp"
 
namespace ablate::particles {
 
struct Field {
    const std::string name;
 
    const PetscInt numberComponents;
 
    const std::vector<std::string> components;
 
    enum domain::FieldLocation location = domain::FieldLocation::AUX;
 
    const PetscDataType dataType;
 
    const PetscInt offset = 0;
 
    PetscInt dataSize = 0;
 
    template <class IndexType>
    inline IndexType operator[](IndexType particle) const {
        return particle * dataSize + offset;
    }
};
 
}  // namespace ablate::particles
#endif  // ABLATELIBRARY_PARTICLEFIELD_HPP