231 lines
11 KiB
C++
231 lines
11 KiB
C++
#include "FluidSim/Scenario.hpp"
|
|
#include "FluidSim/ScenarioCatalog.hpp"
|
|
#include "FluidSim/TwoStrokeScenario.hpp"
|
|
#include <algorithm>
|
|
#include <bit>
|
|
#include <cmath>
|
|
#include <cstdint>
|
|
#include <filesystem>
|
|
#include <fstream>
|
|
#include <iostream>
|
|
#include <memory>
|
|
#include <string>
|
|
#include <type_traits>
|
|
using namespace fluidsim::port;
|
|
namespace {
|
|
class Wav {
|
|
std::ofstream out_;
|
|
void integer(std::uint32_t n, int bytes) {
|
|
for (int i = 0; i < bytes; i++)
|
|
out_.put(char(n >> (8 * i)));
|
|
}
|
|
|
|
public:
|
|
Wav(const std::string &path, int frames) : out_(path, std::ios::binary) {
|
|
if (!out_)
|
|
throw std::runtime_error("cannot open WAV: " + path);
|
|
out_.write("RIFF", 4);
|
|
integer(36 + frames * 8, 4);
|
|
out_.write("WAVEfmt ", 8);
|
|
integer(16, 4);
|
|
integer(3, 2); // IEEE float preserves quiet diagnostic sources.
|
|
integer(2, 2);
|
|
integer(44100, 4);
|
|
integer(44100 * 8, 4);
|
|
integer(8, 2);
|
|
integer(32, 2);
|
|
out_.write("data", 4);
|
|
integer(frames * 8, 4);
|
|
}
|
|
void sample(AudioFrame f) {
|
|
for (float x : {f.left, f.right}) {
|
|
integer(std::bit_cast<std::uint32_t>(x), 4);
|
|
}
|
|
}
|
|
};
|
|
struct Metrics {
|
|
double sum{}, peak{}, delta{};
|
|
std::uint64_t samples{}, clips{};
|
|
AudioFrame previous{};
|
|
void add(AudioFrame f) {
|
|
for (float x : {f.left, f.right}) {
|
|
if (!std::isfinite(x))
|
|
throw std::runtime_error("nonfinite capture");
|
|
sum += x * x;
|
|
peak = std::max(peak, double(std::abs(x)));
|
|
clips += std::abs(x) > 1;
|
|
++samples;
|
|
}
|
|
delta = std::max(
|
|
{delta, double(std::abs(f.left - previous.left)), double(std::abs(f.right - previous.right))});
|
|
previous = f;
|
|
}
|
|
void print(const char *name) const {
|
|
std::cout << name << " peak=" << peak << " rms=" << std::sqrt(sum / samples)
|
|
<< " max_sample_delta=" << delta << " clipping_samples=" << clips << '\n';
|
|
}
|
|
};
|
|
template <class Scenario>
|
|
void capture(const std::string &prefix, int frames, float pedal, OutdoorScene scene, TwoStrokeConfig config, PipeSolver solver, float fixed_rpm, float warmup, const ScenarioExample* example=nullptr, const std::string& solo="all") {
|
|
std::unique_ptr<Scenario> s;
|
|
if constexpr (std::is_same_v<Scenario, TwoStrokeScenario>) s = std::make_unique<Scenario>(44100, config);
|
|
else if constexpr (std::is_same_v<Scenario, FourStrokeScenario>) {
|
|
if(example) s=std::make_unique<Scenario>(44100,*example);
|
|
else s=std::make_unique<Scenario>();
|
|
}
|
|
else s = std::make_unique<Scenario>();
|
|
s->set_pipe_solver(solver);
|
|
s->throttle = pedal;
|
|
s->set_audio_scene(scene);
|
|
bool valid_solo=solo=="all";
|
|
for(std::size_t i=0;i<OutdoorAudio::source_names.size();++i) {
|
|
const bool selected=solo==OutdoorAudio::source_names[i];
|
|
valid_solo|=selected;
|
|
s->set_audio_source_enabled(i,solo=="all" || selected);
|
|
}
|
|
if(!valid_solo) throw std::runtime_error("choose --solo all, intake, exhaust, exhaust-2, or valve-contact");
|
|
for(int i=0; i<int(std::lround(warmup*44100)); ++i) {
|
|
if(fixed_rpm>0) s->crankshaft.omega=fixed_rpm*2*pi/60;
|
|
s->process();
|
|
}
|
|
Wav rendered(prefix + "-rendered.wav", frames), dry(prefix + "-dry.wav", frames),
|
|
raw(prefix + "-unfiltered.wav", frames);
|
|
Metrics m, d, r, tail_metrics;
|
|
std::array<std::unique_ptr<Wav>,4> stems;
|
|
std::array<Metrics,4> source_metrics;
|
|
for(std::size_t i=0;i<stems.size();++i)
|
|
stems[i]=std::make_unique<Wav>(prefix+"-"+OutdoorAudio::source_names[i]+".wav",frames);
|
|
Wav tail(prefix+"-tail.wav",frames);
|
|
int cuts = 0;
|
|
for (int i = 0; i < frames; i++) {
|
|
if(fixed_rpm>0) s->crankshaft.omega=fixed_rpm*2*pi/60;
|
|
auto f = s->process(), df = s->dry_audio(), rf = s->unfiltered_audio();
|
|
for(std::size_t j=0;j<stems.size();++j) {
|
|
const auto stem=s->source_audio()[j].rendered;
|
|
stems[j]->sample(stem); source_metrics[j].add(stem);
|
|
}
|
|
tail.sample(s->tail_audio()); tail_metrics.add(s->tail_audio());
|
|
m.add(f);
|
|
d.add(df);
|
|
r.add(rf);
|
|
rendered.sample(f);
|
|
dry.sample(df);
|
|
raw.sample(rf);
|
|
if (i % 4410 == 0)
|
|
cuts += s->snapshot().fuel_cut;
|
|
}
|
|
for(std::size_t j=0;j<stems.size();++j) source_metrics[j].print(OutdoorAudio::source_names[j]);
|
|
tail_metrics.print("shared tail");
|
|
m.print("rendered");
|
|
d.print("dry filtered");
|
|
r.print("unfiltered spatial");
|
|
std::uint64_t repairs=0, fallbacks=0;
|
|
for(const auto& pipe:s->pipes) { repairs+=pipe.positivity_repairs; fallbacks+=pipe.hllc_fallbacks; }
|
|
if constexpr (std::is_same_v<Scenario, FourStrokeScenario> || std::is_base_of_v<FourStrokeScenario,Scenario>) {
|
|
const auto diagnostics=s->snapshot();
|
|
std::cout<<"junction_model="<<s->junction_config.model<<" junction_fallbacks="<<diagnostics.junction_fallbacks
|
|
<<" junction_limited_fluxes="<<diagnostics.junction_limited_fluxes<<" junction_iterations="<<diagnostics.junction_iterations<<" limited_fluxes="<<diagnostics.boundary_limited_fluxes<<'\n';
|
|
}
|
|
std::cout<<"pipe_positivity_repairs="<<repairs<<" hllc_fallbacks="<<fallbacks<<'\n';
|
|
std::cout << "fuel-cut checkpoints=" << cuts << " final_rpm=" << s->snapshot().rpm
|
|
<< " scene=" << scene_name(scene) << '\n';
|
|
}
|
|
} // namespace
|
|
int main(int argc, char **argv) {
|
|
try {
|
|
std::string example_id, solo="all", junction_model;
|
|
float junction_loss=-1;
|
|
std::string scenario = "twin", prefix = "build/audio/twin", scene = "wall";
|
|
float seconds = 6, pedal = 1, fixed_rpm = 0, warmup = 0;
|
|
TwoStrokeConfig config;
|
|
PipeSolver solver=PipeSolver::Rusanov;
|
|
for (int i = 1; i < argc; i++) {
|
|
std::string key = argv[i];
|
|
if (i + 1 >= argc)
|
|
throw std::runtime_error("missing option value");
|
|
std::string value = argv[++i];
|
|
if(key=="--junctions") {
|
|
if(value!="volume" && value!="characteristic" && value!="legacy") throw std::runtime_error("choose volume, characteristic, or legacy junctions");
|
|
junction_model=value;
|
|
}
|
|
else if(key=="--junction-loss") {
|
|
junction_loss=std::stof(value);
|
|
if(!std::isfinite(junction_loss) || junction_loss<0 || junction_loss>10) throw std::runtime_error("junction loss must be 0..10");
|
|
}
|
|
else if(key=="--solo") solo=value;
|
|
else if(key=="--example") example_id=value;
|
|
else if (key == "--scenario")
|
|
scenario = value;
|
|
else if (key == "--out")
|
|
prefix = value;
|
|
else if (key == "--seconds")
|
|
seconds = std::stof(value);
|
|
else if (key == "--throttle")
|
|
pedal = std::stof(value);
|
|
else if (key == "--warmup") warmup=std::stof(value);
|
|
else if (key == "--rpm") fixed_rpm=std::stof(value);
|
|
else if (key == "--solver") solver=parse_pipe_solver(value);
|
|
else if (key == "--exhaust") {
|
|
if (value == "racing") config.exhaust = TwoStrokeExhaust::Racing;
|
|
else if (value == "straight") config.exhaust = TwoStrokeExhaust::Straight;
|
|
else throw std::runtime_error("choose racing or straight exhaust");
|
|
}
|
|
else if (key == "--scene")
|
|
scene = value;
|
|
else
|
|
throw std::runtime_error("unknown capture option");
|
|
}
|
|
if (!std::isfinite(seconds) || seconds <= 0 || seconds > 600 || !std::isfinite(pedal) || pedal < 0 ||
|
|
pedal > 1)
|
|
throw std::runtime_error("invalid capture duration or throttle");
|
|
if(!std::isfinite(fixed_rpm) || fixed_rpm<0 || fixed_rpm>30000)
|
|
throw std::runtime_error("invalid capture RPM");
|
|
if(!std::isfinite(warmup) || warmup<0 || warmup>600) throw std::runtime_error("invalid warmup");
|
|
OutdoorScene selected;
|
|
if (scene == "dry")
|
|
selected = OutdoorScene::Dry;
|
|
else if (scene == "field")
|
|
selected = OutdoorScene::Field;
|
|
else if (scene == "wall")
|
|
selected = OutdoorScene::Wall;
|
|
else if (scene == "woodland")
|
|
selected = OutdoorScene::Woodland;
|
|
else if (scene == "tunnel")
|
|
selected = OutdoorScene::Tunnel;
|
|
else
|
|
throw std::runtime_error("invalid scene");
|
|
if (!std::filesystem::path(prefix).parent_path().empty())
|
|
std::filesystem::create_directories(std::filesystem::path(prefix).parent_path());
|
|
const int frames = int(std::lround(seconds * 44100));
|
|
auto examples=scenario_examples();
|
|
if(!junction_model.empty() || junction_loss>=0) {
|
|
if(example_id.empty()) example_id=scenario=="v8"?"v8-h":scenario=="b205"?"saab-straight":scenario=="twin"?"twin":"yz-racing";
|
|
for(auto& example:examples) {
|
|
if(!junction_model.empty()) example.junctions.model=junction_model;
|
|
if(junction_loss>=0) example.junctions.turning_loss=junction_loss;
|
|
}
|
|
}
|
|
if(!example_id.empty()) {
|
|
const auto example=std::find_if(examples.begin(),examples.end(),[&](const auto& e){return e.id==example_id;});
|
|
if(example==examples.end()) throw std::runtime_error("unknown example: "+example_id);
|
|
if(example->kind==ScenarioKind::TwoStroke && (!junction_model.empty() || junction_loss>=0))
|
|
throw std::runtime_error("two-stroke exhaust has no configurable pipe junctions");
|
|
if(example->kind==ScenarioKind::TwoStroke)
|
|
capture<TwoStrokeScenario>(prefix,frames,pedal,selected,example->two,solver,fixed_rpm,warmup,nullptr,solo);
|
|
else capture<FourStrokeScenario>(prefix,frames,pedal,selected,config,solver,fixed_rpm,warmup,&*example,solo);
|
|
} else if (scenario == "twin")
|
|
capture<FourStrokeTwinScenario>(prefix, frames, pedal, selected, config, solver, fixed_rpm, warmup,nullptr,solo);
|
|
else if (scenario == "b205")
|
|
capture<SaabB205Scenario>(prefix, frames, pedal, selected, config, solver, fixed_rpm, warmup,nullptr,solo);
|
|
else if (scenario == "two-stroke")
|
|
capture<TwoStrokeScenario>(prefix, frames, pedal, selected, config, solver, fixed_rpm, warmup,nullptr,solo);
|
|
else if (scenario == "v8")
|
|
capture<AmericanV8Scenario>(prefix, frames, pedal, selected, config, solver, fixed_rpm, warmup,nullptr,solo);
|
|
else
|
|
throw std::runtime_error("choose twin, two-stroke, b205 or v8");
|
|
} catch (const std::exception &e) {
|
|
std::cerr << e.what() << '\n';
|
|
return 1;
|
|
}
|
|
}
|