hovsep

const (
	restingLungVolume = 700.0 * Milliliter

	ResidualLungVolume = 600.0 * Milliliter
	TotalLungCapacity  = 3000.0 * Milliliter

	defaultLungCompliance   = 100.0 * MlPerCmH2O
	defaultAirwayResistance = 0.002 * CmH2OPerMlPerSecond

	lungVolumeAsymmetry     = 5 * Percent

hovsep

const (
	TidalRespiratoryRate = 12 * PerMinute

	MinRespiratoryRate = 8 * PerMinute
	MaxRespiratoryRate = 30 * PerMinute

	BasePleuralPressure          = -5 * CmH2O // resting pleural pressure at FRC
	InspiratoryPressureAmplitude = 3 * CmH2O  // peak swing during quiet breathing (−5 → −8 cmH₂O)

	inhaleFraction = 1.0 / 3.0

hovsep

// What it checks: run simulation for 100ms and check the ECG, we must see a few R-peaks
func Test_Human(t *testing.T) {
	tests := []struct {
		name       string
		assertions func(t *testing.T, sim *step_sim.Simulation)
	}{
		{
			name: "heart is beating",
			assertions: func(t *testing.T, sim *step_sim.Simulation) {

hovsep

// PackAutonomicTone builds a signal that represents autonomic tone
func PackAutonomicTone(sym, paraSym, noise, gain, cardiacBias, vascularBias, respiratoryBias, giBias float64) *signal.Signal {
	return signal.New(signal.NewGroup().Add(
		NewLevel(sym, common.Sympathetic),
		NewLevel(paraSym, common.Parasympathetic),
		NewLevel(noise, common.Noise),
		NewLevel(gain, common.Gain),
		NewBias(cardiacBias, common.Cardiac),
		NewBias(vascularBias, common.Vascular),
		NewBias(respiratoryBias, common.Respiratory),

hovsep

// WithRunningSimulation is a helper function that runs the simulation and executes a callback after a given duration
func WithRunningSimulation(sim *step_sim.Simulation, duration time.Duration, f func()) {
	go sim.Run()
	defer func() {
		sim.SendCommand(step_sim.Exit)
	}()

	// Let the simulation run for a while
	time.Sleep(duration)

hovsep

func Test_Time(t *testing.T) {
	tests := []struct {
		name       string
		assertions func(t *testing.T, sim *step_sim.Simulation)
	}{
		{
			name: "time advances in timer",
			assertions: func(t *testing.T, sim *step_sim.Simulation) {
				var observedSimWallTime []time.Time

hovsep

// setMeshCommands sets the commands that can be executed on the mesh
func setMeshCommands(mesh *fmesh.FMesh, commands step_sim.MeshCommandMap) {
	timeComponent := mesh.ComponentByName("time")

	// Print current time
	commands["time:now"] = step_sim.NewMeshCommandDescriptor("Print current time", func(_ *fmesh.FMesh) {
		tickCount := timeComponent.State().Get("tick_count")
		simTime := timeComponent.State().Get("sim_time")
		simWallTime := timeComponent.State().Get("sim_wall_time")
		fmt.Println("Current tick count:", tickCount)

hovsep

// I like my main func to be clean and short
func main() {
	simMesh := getSimulationMesh()

	// Run the mesh in a step simulation
	step_sim.NewApp(simMesh, initSim).Run()
}

// Here is our step_sim wrapper, it provides REPL and some built-in commands: pause, resume, exit, help
func NewApp(fm *fmesh.FMesh, simInitFunc SimInitFunc) *Application {

hovsep

// getSimulationMesh returns the main mesh of the simulation
func getSimulationMesh() *fmesh.FMesh {
	// Set up the world
	habitat := env.NewHabitat(component.NewCollection().Add(
		factor.GetTimeComponent(),
		factor.GetAirComponent(),
		factor.GetSunComponent(),
	))

	// Add human beings

hovsep

func main() {
	fm := getMesh()

	// Initialize the mesh: set diagnostic frames to USB port, so the laptop will send them
	laptopInstance.SendDataToUSB(
		diagnostics.FrameGetEngineDTCs,
		diagnostics.FrameGetSpeed,
		diagnostics.FrameGetRPM,
		diagnostics.FrameGetCoolantTemperature,
		diagnostics.FrameGetCalibrationID,