deep_causality_cfd 0.2.0

Counterfactual Computational fluid dynamics solvers and the Flow DSL.
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
32

/*
 * SPDX-License-Identifier: MIT
 * Copyright (c) 2023 - 2026. The DeepCausality Authors and Contributors. All Rights Reserved.
 */

//! [`Marcher`] realization for the DEC solver: one projected RK4 step under the
//! per-step [`Ambient`].

use crate::CfdScalar;
use crate::solvers::dec::DecNsSolver;
use crate::solvers::dec::step_output::StepOutput;
use crate::traits::Marcher;
use crate::types::Ambient;
use deep_causality_physics::{PhysicsError, SolenoidalField};

impl<'m, const D: usize, R: CfdScalar> Marcher<R> for DecNsSolver<'m, D, R> {
    type State = SolenoidalField<R>;
    type Ambient = Ambient<R>;
    type Output = StepOutput<R>;

    /// Advance one projected step. The ambient `ν` is applied to the rate before the
    /// step (between-step write; the four RK4 stages then read it). With a constant
    /// ambient this reproduces the construction-fixed march bit-for-bit.
    fn advance(
        &self,
        state: &SolenoidalField<R>,
        ambient: &Ambient<R>,
    ) -> Result<StepOutput<R>, PhysicsError> {
        self.rate().set_nu(*ambient.nu());
        self.step(state)
    }
}