From 1ee35d53ccab61523ca6ae026d24edd48a9fc1f4 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Fr=C3=A9d=C3=A9ric=20Simonis?= Date: Thu, 11 Jan 2024 14:07:18 +0100 Subject: [PATCH] Fix coordinates in et1d solid rust --- elastic-tube-1d/solid-rust/src/main.rs | 7 +++---- 1 file changed, 3 insertions(+), 4 deletions(-) diff --git a/elastic-tube-1d/solid-rust/src/main.rs b/elastic-tube-1d/solid-rust/src/main.rs index 831ca4071..f85f22f1d 100644 --- a/elastic-tube-1d/solid-rust/src/main.rs +++ b/elastic-tube-1d/solid-rust/src/main.rs @@ -28,6 +28,7 @@ fn main() -> ExitCode { const DOMAIN_SIZE: usize = 100; const CHUNK_SIZE: usize = DOMAIN_SIZE + 1; + const TUBE_LENGTH: f64 = 10.0; let mut participant = precice::Participant::new("Solid", &config, 0, 1); @@ -44,11 +45,9 @@ fn main() -> ExitCode { let grid_size = CHUNK_SIZE * dimensions as usize; let grid: Vec = { let mut v: Vec = vec![0_f64; grid_size]; + const DX : f64 = TUBE_LENGTH / DOMAIN_SIZE as f64; for i in 0..CHUNK_SIZE - 1 { - for j in 0..(dimensions as usize) - 1 { - let idx = i * dimensions as usize + j; - v[idx] = (i * (1 - j)) as f64; - } + v[i * dimensions as usize] = DX * i as f64; } v };