caustic 0.0.12

A General-Purpose 6D Collisionless Gravitational Dynamics Solver
Documentation 
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

//! Caustic-specific exit conditions.

use super::super::super::{
    conditions::{ExitCondition, ExitReason},
    diagnostics::GlobalDiagnostics,
};

/// Exit when the first caustic forms (max stream count > 1).
pub struct CausticExitCondition;

impl ExitCondition for CausticExitCondition {
    fn check(&self, _diag: &GlobalDiagnostics, _initial: &GlobalDiagnostics) -> Option<ExitReason> {
        // Stream count requires repr access; not available through GlobalDiagnostics alone
        None
    }
}

/// Exit when virial ratio 2T/|W| stabilises within tolerance for N consecutive steps.
pub struct VirialRelaxedExit {
    pub tolerance: f64,
}

impl ExitCondition for VirialRelaxedExit {
    fn check(&self, diag: &GlobalDiagnostics, _initial: &GlobalDiagnostics) -> Option<ExitReason> {
        if (diag.virial_ratio - 1.0).abs() < self.tolerance {
            Some(ExitReason::VirialRelaxed)
        } else {
            None
        }
    }
}