raySharingRadiation.hpp

#ifndef ABLATELIBRARY_RAYSHARINGRADIATION_H
#define ABLATELIBRARY_RAYSHARINGRADIATION_H
 
#include "domain/reverseRange.hpp"
#include "radiation.hpp"
 
namespace ablate::radiation {
 
class RaySharingRadiation : public ablate::radiation::Radiation {
   public:
    RaySharingRadiation(const std::string& solverId, const std::shared_ptr<domain::Region>& region, const PetscInt raynumber,
                        std::shared_ptr<eos::radiationProperties::RadiationModel> radiationModelIn, std::shared_ptr<ablate::monitors::logs::Log> = {});
    ~RaySharingRadiation();
 
    void Setup(const ablate::domain::Range& cellRange, ablate::domain::SubDomain& subDomain) override;
 
    void IdentifyNewRaysOnRank(ablate::domain::SubDomain& subDomain, DM radReturn, PetscInt npoints) override;
 
    void ParticleStep(ablate::domain::SubDomain& subDomain, DM faceDM, const PetscScalar* faceGeomArray, DM radReturn, PetscInt nlocalpoints,
                      PetscInt nglobalpoints) override;  
 
    static inline std::string GetClassType() { return "RaySharingRadiation"; }
 
    void SetBoundary(CellSegment& raySegment, PetscInt index, Identifier identifier) override {
        PetscMPIInt rank;
        MPI_Comm_rank(PETSC_COMM_WORLD, &rank);
        if (identifier.originRank == rank) {
            raySegment.cell = index;
            raySegment.pathLength = -1;
        }
    }
 
   protected:
    ablate::domain::ReverseRange indexLookup;
 
    std::vector<PetscReal> remoteMap;
};
 
}  // namespace ablate::radiation
 
#endif  // ABLATELIBRARY_RAYSHARINGRADIATION_H