Mirrored_Repos/veloren
1
0
Fork 0
mirror of https://gitlab.com/veloren/veloren.git synced 2024-08-30 18:12:32 +00:00
Code Issues Packages Projects Releases Wiki Activity

have a first constant formular for pos but cant get dist from it so far

Browse Source
This commit is contained in:
Marcel Märtens 2022-06-10 13:30:13 +02:00
parent aa68bd6c9f
commit ccf9fdf9a7
1 changed files with 12 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
common/systems/tests/phys/basic.rs
View File

@ -5,7 +5,6 @@ use specs::WorldExt;
use std::error::Error; use std::error::Error;
use utils::{DT, DT_F64, EPSILON}; use utils::{DT, DT_F64, EPSILON};
use vek::{approx, Vec2, Vec3}; use vek::{approx, Vec2, Vec3};
use common::states::utils::unwrap_tool_data;
use veloren_common_systems::add_local_systems; use veloren_common_systems::add_local_systems;
#[test] #[test]
@ -324,29 +323,33 @@ fn physics_theory() -> Result<(), Box<dyn Error>> {
let mass = 1.0; let mass = 1.0;
let air_friction_co = 1.1_f64; let air_friction_co = 0.9_f64;
let air_friction_area = 1.0_f64; let air_friction_area = 0.75_f64;
let air_density = 1.225_f64; let air_density = 1.225_f64;
let c = air_friction_co * air_friction_area * 0.5 * air_density * mass; let c = air_friction_co * air_friction_area * 0.5 * air_density * mass;
//let acc = accel_force / mass * move_dir; //let acc = accel_force / mass * move_dir;
let acc = 9.0_f64; // btw: cant accelerate faster than gravity on foot let acc = 9.2_f64; // btw: cant accelerate faster than gravity on foot
let old_vel = vel; let old_vel = vel;
// controller // controller
// I know what you think, wtf, yep: https://math.stackexchange.com/questions/1929436/line-integral-of-force-of-air-resistanc // I know what you think, wtf, yep: https://math.stackexchange.com/questions/1929436/line-integral-of-force-of-air-resistanc
// basically an integral of the air resistance formula which scales with v^2 transformed with an ODE. // basically an integral of the air resistance formula which scales with v^2 transformed with an ODE.
let vel = (mass * acc / c).sqrt() * ( ( (c / acc / mass ).sqrt() * old_vel ).atanh() + (c * acc / mass).sqrt() * 0.1* dt ).tanh(); let vel = (mass * acc / c).sqrt() * ( ( (c / acc / mass ).sqrt() * old_vel ).atanh() + (c * acc / mass).sqrt() * dt ).tanh();
//physics //physics
let acc2 = (vel - old_vel) / dt; //let acc2 = (vel - old_vel) / dt;
let distance = vel * dt;// 0.5 * vel * dt - 0.5 * acc2 * dt * dt; //let distance = 0.5 * vel * dt - 0.5 * acc2 * dt * dt;
let pos = pos + distance; let distance = mass / c * ( ( ( old_vel * (c/acc / mass).sqrt()).atanh() + dt* (c * acc / mass).sqrt()).cosh() ).ln();
let distance2 = mass / c * ( ( ( vel * (c/acc / mass).sqrt()).atanh() + dt* (c * acc / mass).sqrt()).cosh() ).ln();
let diff = distance2 - distance;
let pos = pos + diff;
//if ((i+1) as f64 *dt * 10.0).round() as i64 % 2 == 0 { //if ((i+1) as f64 *dt * 10.0).round() as i64 % 2 == 0 {
println!("[{:0>2.1}]: move_dir: {:4.1}, acc: {:4.4}, vel: {:4.4}, pos: {:4.4}, c: {:4.4}", (i+1) as f64 *dt, move_dir, acc, vel, pos, c); println!("[{:0>2.1}]: move_dir: {:0>1.1}, acc: {:0>4.4}, vel: {:0>4.4}, dist: {:0>7.4}, dist: {:0>7.4}, pos: {:0>7.4}, c: {:0>4.4}", (i+1) as f64 *dt, move_dir, acc, vel, distance, distance2, pos, c);
//} //}
(acc, vel, pos) (acc, vel, pos)
}; };