diff --git a/.gitignore b/.gitignore index 4325418..de35486 100644 --- a/.gitignore +++ b/.gitignore @@ -3,6 +3,7 @@ build/ffi/ *.o *.a target/ +examples/vids/ assets/old/ assets/icons/ Cargo.lock diff --git a/crates/femm-app/src/doc_canvas.rs b/crates/femm-app/src/doc_canvas.rs index 501dc04..abc445f 100644 --- a/crates/femm-app/src/doc_canvas.rs +++ b/crates/femm-app/src/doc_canvas.rs @@ -51,6 +51,7 @@ pub enum Tool { AddNode, AddBlockLabel, AddSegment, + AddProbe, } /// composition mode for a pick action, derived from modifier keys at click/release time. @@ -124,9 +125,10 @@ pub fn view<'a>( view_state: ViewState, show_grid: bool, zoom_window_active: bool, + probes: &'a [crate::probe::Probe], ) -> Element<'a, CanvasMessage> { Canvas::new(DocCanvas { - doc, tool, mesh, solution, render_mode, view_state, show_grid, zoom_window_active, + doc, tool, mesh, solution, render_mode, view_state, show_grid, zoom_window_active, probes, }) .width(Length::Fill) .height(Length::Fill) @@ -142,6 +144,7 @@ struct DocCanvas<'a> { view_state: ViewState, show_grid: bool, zoom_window_active: bool, + probes: &'a [crate::probe::Probe], } impl<'a> canvas::Program for DocCanvas<'a> { @@ -519,6 +522,21 @@ impl<'a> canvas::Program for DocCanvas<'a> { }); } + for probe in self.probes { + let p = view.map(probe.x, probe.y); + let color = Color::from_rgba8(probe.color[0], probe.color[1], probe.color[2], probe.color[3] as f32 / 255.0); + frame.fill(&Path::circle(p, 6.0), color); + frame.stroke(&Path::circle(p, 6.0), + Stroke::default().with_width(1.5).with_color(Color::WHITE)); + frame.fill_text(Text { + content: probe.label.clone(), + position: Point::new(p.x + 9.0, p.y - 9.0), + color, + size: 11.0.into(), + ..Text::default() + }); + } + if let (Some(start), true, Some(now)) = (state.marquee_from, state.dragged, state.cursor_canvas) { @@ -549,7 +567,7 @@ impl<'a> canvas::Program for DocCanvas<'a> { if cursor.position_in(bounds).is_some() { return match self.tool { Tool::Select => mouse::Interaction::Grab, - Tool::AddNode | Tool::AddBlockLabel | Tool::AddSegment => { + Tool::AddNode | Tool::AddBlockLabel | Tool::AddSegment | Tool::AddProbe => { mouse::Interaction::Crosshair } }; diff --git a/crates/femm-app/src/export.rs b/crates/femm-app/src/export.rs index f7356f3..8afd36a 100644 --- a/crates/femm-app/src/export.rs +++ b/crates/femm-app/src/export.rs @@ -29,6 +29,8 @@ pub struct ExportInput { pub height: u32, pub contour: bool, pub crf: u32, + pub probes: Vec, + pub probe_samples: Vec>, pub progress: Arc>, } @@ -86,6 +88,9 @@ pub fn export_webm(input: ExportInput) -> Result { "-pix_fmt", "yuv420p", "-crf", crf_arg.as_str(), "-b:v", "0", + "-r", fps_arg.as_str(), + "-fps_mode", "cfr", + "-vsync", "cfr", "-deadline", "good", "-cpu-used", "2", "-row-mt", "1", @@ -135,7 +140,12 @@ pub fn export_webm(input: ExportInput) -> Result { let t = (f as f64) * input.dt; kinematic::apply_tracks(&mut doc, &input.base_doc, &input.tracks, t); - let pixmap = render_frame(&doc, sol, input.width, input.height, input.contour, xmin, xmax, ymin, ymax); + let pixmap = render_frame( + &doc, sol, input.width, input.height, input.contour, + xmin, xmax, ymin, ymax, + &input.probes, &input.probe_samples, + f, input.buffer_size, + ); let frame_bytes = pixmap.data(); if frame_bytes.len() != expected_bytes_per_frame { drop(stdin); @@ -276,10 +286,16 @@ fn render_frame( height: u32, contour: bool, xmin: f64, xmax: f64, ymin: f64, ymax: f64, + probes: &[crate::probe::Probe], + probe_samples: &[HashMap], + current_frame: i64, + buffer_size: usize, ) -> Pixmap { let mut pixmap = Pixmap::new(width, height).expect("allocate pixmap"); pixmap.fill(Color::from_rgba8(20, 20, 24, 255)); - let view = View::fit(width, height, xmin, xmax, ymin, ymax); + let plot_strip_h = if probes.is_empty() { 0 } else { (height as f32 * 0.33) as u32 }; + let canvas_h = height - plot_strip_h; + let view = View::fit(width, canvas_h, xmin, xmax, ymin, ymax); let (lo, hi) = sol.b_magnitude_range(); let span = if hi > lo { hi - lo } else { 1.0 }; @@ -366,9 +382,163 @@ fn render_frame( } } + for probe in probes { + let (px, py) = view.map(probe.x, probe.y); + paint.set_color(Color::from_rgba8(probe.color[0], probe.color[1], probe.color[2], 255)); + let mut pb = PathBuilder::new(); + pb.push_circle(px, py, 6.0); + if let Some(path) = pb.finish() { + pixmap.fill_path(&path, &paint, FillRule::Winding, Transform::identity(), None); + } + paint.set_color(Color::from_rgba8(255, 255, 255, 255)); + let ring = Stroke { width: 1.6, ..Default::default() }; + let mut pb = PathBuilder::new(); + pb.push_circle(px, py, 6.0); + if let Some(path) = pb.finish() { + pixmap.stroke_path(&path, &paint, &ring, Transform::identity(), None); + } + } + + if plot_strip_h > 0 && !probes.is_empty() { + draw_plot_strip( + &mut pixmap, + 0, canvas_h, width, plot_strip_h, + probes, probe_samples, + current_frame, buffer_size, + ); + } + pixmap } +/// draws one plot pane per probe across a horizontal strip at (x0, y0, w, h), each pane color-matched to its probe and marked with a vertical bar at the current frame. +fn draw_plot_strip( + pixmap: &mut Pixmap, + x0: u32, y0: u32, w: u32, h: u32, + probes: &[crate::probe::Probe], + probe_samples: &[HashMap], + current_frame: i64, + buffer_size: usize, +) { + let n = probes.len().max(1); + let pane_w = (w as f32) / (n as f32); + let pane_x_top_pad = 4.0; + let pane_y_top_pad = 18.0; + let pane_bot_pad = 6.0; + let pane_lr_pad = 4.0; + + let mut paint = Paint::default(); + paint.anti_alias = true; + paint.set_color(Color::from_rgba8(12, 14, 20, 255)); + let strip = tiny_skia::Rect::from_xywh(x0 as f32, y0 as f32, w as f32, h as f32); + if let Some(rect) = strip { + pixmap.fill_rect(rect, &paint, Transform::identity(), None); + } + + for (i, probe) in probes.iter().enumerate() { + let px0 = x0 as f32 + (i as f32) * pane_w + pane_lr_pad; + let py0 = y0 as f32 + pane_x_top_pad; + let pw = pane_w - 2.0 * pane_lr_pad; + let ph = (h as f32) - pane_x_top_pad - pane_bot_pad; + if pw <= 0.0 || ph <= 0.0 { continue; } + + paint.set_color(Color::from_rgba8(18, 20, 26, 255)); + if let Some(r) = tiny_skia::Rect::from_xywh(px0, py0, pw, ph) { + pixmap.fill_rect(r, &paint, Transform::identity(), None); + } + + let plot_x = px0 + 4.0; + let plot_y = py0 + pane_y_top_pad; + let plot_w = pw - 8.0; + let plot_h = ph - pane_y_top_pad - 6.0; + if plot_w <= 0.0 || plot_h <= 0.0 { continue; } + + let probe_color = Color::from_rgba8(probe.color[0], probe.color[1], probe.color[2], 255); + paint.set_color(probe_color); + + let samples = probe_samples.get(i); + let mut points: Vec<(i64, f64)> = Vec::new(); + if let Some(s) = samples { + for (&fi, &(bx, by)) in s { + points.push((fi, probe.mode.extract(bx, by))); + } + } + if points.is_empty() { + continue; + } + let mut ymin = f64::INFINITY; + let mut ymax = f64::NEG_INFINITY; + for &(_, v) in &points { + if v < ymin { ymin = v; } + if v > ymax { ymax = v; } + } + if !ymin.is_finite() || !ymax.is_finite() { ymin = -1.0; ymax = 1.0; } + if (ymax - ymin).abs() < 1e-12 { ymin -= 1.0; ymax += 1.0; } + points.sort_by_key(|p| p.0); + + let buf = buffer_size.max(1) as f64; + let x_for = |fi: i64| -> f32 { + let t = (fi.rem_euclid(buf as i64) as f64) / buf; + plot_x + (t as f32) * plot_w + }; + let y_for = |v: f64| -> f32 { + let t = (v - ymin) / (ymax - ymin); + plot_y + plot_h - (t as f32) * plot_h + }; + + let axis_stroke = Stroke { width: 1.0, ..Default::default() }; + paint.set_color(Color::from_rgba8(60, 65, 75, 255)); + let mut pb = PathBuilder::new(); + pb.move_to(plot_x, plot_y + plot_h); + pb.line_to(plot_x + plot_w, plot_y + plot_h); + if let Some(path) = pb.finish() { + pixmap.stroke_path(&path, &paint, &axis_stroke, Transform::identity(), None); + } + let mid_y = y_for((ymin + ymax) * 0.5); + paint.set_color(Color::from_rgba8(50, 55, 65, 255)); + let mut pb = PathBuilder::new(); + pb.move_to(plot_x, mid_y); + pb.line_to(plot_x + plot_w, mid_y); + if let Some(path) = pb.finish() { + pixmap.stroke_path(&path, &paint, &axis_stroke, Transform::identity(), None); + } + + paint.set_color(probe_color); + let line_stroke = Stroke { width: 1.4, ..Default::default() }; + let mut pb = PathBuilder::new(); + let mut started = false; + for (fi, v) in &points { + let x = x_for(*fi); + let y = y_for(*v); + if !started { pb.move_to(x, y); started = true; } + else { pb.line_to(x, y); } + } + if started { + if let Some(path) = pb.finish() { + pixmap.stroke_path(&path, &paint, &line_stroke, Transform::identity(), None); + } + } + + let cx = x_for(current_frame); + paint.set_color(Color::from_rgba8(255, 255, 255, 200)); + let cur_stroke = Stroke { width: 1.2, ..Default::default() }; + let mut pb = PathBuilder::new(); + pb.move_to(cx, plot_y); + pb.line_to(cx, plot_y + plot_h); + if let Some(path) = pb.finish() { + pixmap.stroke_path(&path, &paint, &cur_stroke, Transform::identity(), None); + } + + let label_dot_r = 4.0; + paint.set_color(probe_color); + let mut pb = PathBuilder::new(); + pb.push_circle(px0 + 8.0, py0 + 9.0, label_dot_r); + if let Some(path) = pb.finish() { + pixmap.fill_path(&path, &paint, FillRule::Winding, Transform::identity(), None); + } + } +} + /// returns a polyline approximating an arc between two endpoints across a CCW degree sweep, capped at max_steps subdivisions. fn arc_polyline(x0: f64, y0: f64, x1: f64, y1: f64, sweep_deg: f64, max_steps: usize) -> Vec<(f64, f64)> { let theta = sweep_deg.to_radians(); diff --git a/crates/femm-app/src/kinematic.rs b/crates/femm-app/src/kinematic.rs index 4eb2d40..3ea09d8 100644 --- a/crates/femm-app/src/kinematic.rs +++ b/crates/femm-app/src/kinematic.rs @@ -86,7 +86,7 @@ impl Expression { } } -/// one continuous-time track: clamped chord, axis, member nodes, the closed-form displacement expression, and the pre-subdivision edge spine used to re-subdivide on demand. +/// one continuous-time track: clamped chord, axis, member nodes, the closed-form displacement expression, the pre-subdivision edge spine, and the per-track subdivision count regenerated against a stashed pristine doc. #[derive(Debug, Clone)] pub struct Track { pub label: String, @@ -97,6 +97,7 @@ pub struct Track { pub expression: Expression, pub edge_endpoints: Vec<(usize, usize)>, pub edge_subdivisions: Vec>, + pub subdivisions: usize, } /// snaps doc nodes to base and layers each track's displacement at simulated time t. @@ -184,148 +185,53 @@ pub fn track_from_selection( expression, edge_endpoints: Vec::new(), edge_subdivisions: Vec::new(), + subdivisions: 0, }) } -/// removes the chain of segments inside each edge and replaces it with a fresh chain at `count` subdivisions, appending new node indices to the doc. orphans the prior subdivision nodes to preserve every other track's index references. -pub fn rebuild_track_subdivisions(track: &mut Track, doc: &mut FemmDoc, count: usize) -> Result<(), String> { - if count < 2 { - return Err(String::from("subdivision count must be at least 2")); - } - if track.edge_endpoints.is_empty() { - let (eps, subs) = infer_edges_from_positions(track, doc); - if eps.is_empty() { - return Err(String::from("track has no recorded edges and inference from positions found none - delete and recreate the track")); +/// resets the live doc to a clone of the pristine doc, then re-subdivides every track's edges at the track's per-track subdivision count, refreshing each track's member_nodes and edge_subdivisions. orphan-free because no prior subdivision nodes survive the reset. +pub fn rebuild_all_tracks_from_pristine( + live_doc: &mut FemmDoc, + pristine: &FemmDoc, + tracks: &mut [Track], +) -> Result<(), String> { + *live_doc = pristine.clone(); + for (ti, track) in tracks.iter_mut().enumerate() { + if track.edge_endpoints.is_empty() { + return Err(format!("track {} has no recorded edges; delete and recreate it", ti + 1)); } - track.edge_endpoints = eps; - track.edge_subdivisions = subs; - } - while track.edge_subdivisions.len() < track.edge_endpoints.len() { - track.edge_subdivisions.push(Vec::new()); - } - - let mut new_subdivisions: Vec> = Vec::with_capacity(track.edge_endpoints.len()); - for (edge_idx, (a, b)) in track.edge_endpoints.iter().copied().enumerate() { - let old_subs = track.edge_subdivisions[edge_idx].clone(); - let edge_nodes: std::collections::HashSet = std::iter::once(a) - .chain(std::iter::once(b)) - .chain(old_subs.iter().copied()) - .collect(); - doc.segments.retain(|s| { - !(edge_nodes.contains(&(s.n0 as usize)) && edge_nodes.contains(&(s.n1 as usize))) - }); - - let pa = match doc.nodes.get(a) { - Some(n) => (n.x, n.y), - None => return Err(format!("edge endpoint {a} out of bounds")), - }; - let pb = match doc.nodes.get(b) { - Some(n) => (n.x, n.y), - None => return Err(format!("edge endpoint {b} out of bounds")), - }; - - let mut new_subs: Vec = Vec::with_capacity(count.saturating_sub(1)); - for i in 1..count { - let t = i as f64 / count as f64; - let x = pa.0 + (pb.0 - pa.0) * t; - let y = pa.1 + (pb.1 - pa.1) * t; - let idx = doc.nodes.len(); - doc.nodes.push(femm_doc_mag::Node { - x, y, - boundary_marker: String::new(), - in_group: 0, - selected: false, + let count = if track.subdivisions >= 2 { track.subdivisions } else { 20 }; + let mut new_edge_subs: Vec> = Vec::with_capacity(track.edge_endpoints.len()); + for &(a, b) in &track.edge_endpoints { + let seg_idx = live_doc.segments.iter().position(|s| { + let n0 = s.n0 as usize; + let n1 = s.n1 as usize; + (n0 == a && n1 == b) || (n0 == b && n1 == a) }); - new_subs.push(idx); - } - - let mut prev = a; - for &mid in &new_subs { - doc.segments.push(femm_doc_mag::Segment { - n0: prev as i32, - n1: mid as i32, - max_side_length: -1.0, - boundary_marker: String::new(), - hidden: false, - in_group: 0, - selected: false, - }); - prev = mid; - } - doc.segments.push(femm_doc_mag::Segment { - n0: prev as i32, - n1: b as i32, - max_side_length: -1.0, - boundary_marker: String::new(), - hidden: false, - in_group: 0, - selected: false, - }); - new_subdivisions.push(new_subs); - } - - let mut corner_set: std::collections::HashSet = std::collections::HashSet::new(); - for &(a, b) in &track.edge_endpoints { - corner_set.insert(a); - corner_set.insert(b); - } - let mut new_members: Vec = Vec::new(); - for &c in &corner_set { - if c != track.anchor_a && c != track.anchor_b { new_members.push(c); } - } - for subs in &new_subdivisions { - for &s in subs { new_members.push(s); } - } - track.member_nodes = new_members; - track.edge_subdivisions = new_subdivisions; - Ok(()) -} - -/// reconstructs a track's edges from its node positions: clusters anchors+members by the lateral coordinate (perpendicular to the displacement axis), and treats each cluster as one edge with extremes as endpoints and middles as subdivisions. -fn infer_edges_from_positions(track: &Track, doc: &FemmDoc) -> (Vec<(usize, usize)>, Vec>) { - let along_x = matches!(track.axis, Axis::PlusY | Axis::MinusY); - let pos = |i: usize| doc.nodes.get(i).map(|n| (n.x, n.y)); - let perp = move |i: usize| pos(i).map(|(x, y)| if along_x { x } else { y }); - let lat = move |i: usize| pos(i).map(|(x, y)| if along_x { y } else { x }); - - let mut all_indices: Vec = Vec::with_capacity(track.member_nodes.len() + 2); - all_indices.push(track.anchor_a); - all_indices.push(track.anchor_b); - all_indices.extend(&track.member_nodes); - - let tol = 0.05_f64; - let mut clusters: Vec<(f64, Vec)> = Vec::new(); - for &idx in &all_indices { - let v = match lat(idx) { Some(v) => v, None => continue }; - let mut placed = false; - for cluster in &mut clusters { - if (v - cluster.0).abs() < tol { - cluster.1.push(idx); - placed = true; - break; + match seg_idx { + Some(i) => { + let subs = live_doc.subdivide_segment(i, count) + .into_iter().map(|n| n as usize).collect(); + new_edge_subs.push(subs); + } + None => { + new_edge_subs.push(Vec::new()); + } } } - if !placed { clusters.push((v, vec![idx])); } - } - - let mut endpoints = Vec::new(); - let mut subdivisions = Vec::new(); - for (_, cluster) in clusters { - if cluster.len() < 2 { continue; } - let mut sorted: Vec<(usize, f64)> = cluster.iter() - .filter_map(|&i| perp(i).map(|p| (i, p))) + let mut corner_set: std::collections::HashSet = std::collections::HashSet::new(); + for &(a, b) in &track.edge_endpoints { + corner_set.insert(a); + corner_set.insert(b); + } + let mut members: Vec = corner_set.into_iter() + .filter(|&c| c != track.anchor_a && c != track.anchor_b) .collect(); - sorted.sort_by(|a, b| a.1.partial_cmp(&b.1).unwrap_or(std::cmp::Ordering::Equal)); - let a = sorted.first().unwrap().0; - let b = sorted.last().unwrap().0; - let subs: Vec = if sorted.len() >= 3 { - sorted[1..sorted.len()-1].iter().map(|(i, _)| *i).collect() - } else { - Vec::new() - }; - endpoints.push((a, b)); - subdivisions.push(subs); + for subs in &new_edge_subs { + for &s in subs { members.push(s); } + } + track.edge_subdivisions = new_edge_subs; + track.member_nodes = members; } - - (endpoints, subdivisions) + Ok(()) } diff --git a/crates/femm-app/src/main.rs b/crates/femm-app/src/main.rs index e325d1a..a1cfc3f 100644 --- a/crates/femm-app/src/main.rs +++ b/crates/femm-app/src/main.rs @@ -3,6 +3,8 @@ mod doc_canvas; mod export; mod kinematic; +mod plot; +mod probe; mod session; mod sim_meta; mod spice; @@ -96,6 +98,9 @@ enum Message { SimExportWebM, SimExportComplete(Result), SimExportPoll, + ProbeSetMode(usize, probe::ProbeMode), + ProbeRemove(usize), + ProbeExportWav(usize), } /// copyable diagnostic shown when a pipeline stage fails. @@ -130,6 +135,8 @@ struct Simulation { fundamental_hz: f64, fundamental_hz_text: String, loop_playback: bool, + probes: Vec, + probe_samples: Vec>, } struct App { @@ -149,6 +156,7 @@ struct App { zoom_window_active: bool, canvas_size: Option<(f32, f32)>, export_progress: Option>>, + pristine_doc: Option, } impl App { @@ -171,6 +179,7 @@ impl App { zoom_window_active: false, canvas_size: None, export_progress: None, + pristine_doc: None, }; (app, Task::none()) } @@ -199,11 +208,80 @@ impl App { } } + /// reads a .femsess from path and replaces the doc, simulation, pristine snapshot, and source path on success. + fn load_session_from_path(&mut self, path: &Path) -> Task { + let text = match std::fs::read_to_string(path) { + Ok(t) => t, + Err(e) => { + self.error = Some(ErrorReport { + title: String::from("read session failed"), + body: format!("{path:?}\n\n{e}"), + }); + return Task::none(); + } + }; + let (base_doc, pristine_doc, meta) = match session::parse(&text) { + Ok(triple) => triple, + Err(e) => { + self.error = Some(ErrorReport { + title: String::from("parse session failed"), + body: format!("{path:?}\n\n{e}"), + }); + return Task::none(); + } + }; + let tracks = match sim_meta::resolve_tracks(&meta, &base_doc) { + Ok(t) => t, + Err(e) => { + self.error = Some(ErrorReport { + title: String::from("session tracks do not match embedded doc"), + body: format!("{path:?}\n\n{e}"), + }); + return Task::none(); + } + }; + self.doc = base_doc.clone(); + self.mesh = None; + self.solution = None; + self.error = None; + self.pristine_doc = pristine_doc; + let mut sim = new_simulation(base_doc, String::new()); + if let Some(secs) = spice::parse_spice(&meta.dt_text) { sim.dt = secs; } + sim.dt_text = meta.dt_text.clone(); + if let Some(secs) = spice::parse_spice(&meta.interval_text) { + sim.wall_interval = Duration::from_secs_f64(secs); + } + sim.interval_text = meta.interval_text.clone(); + sim.subdivisions = meta.subdivisions; + sim.subdivisions_text = meta.subdivisions.to_string(); + sim.buffer_size = meta.buffer_size; + sim.buffer_size_text = meta.buffer_size.to_string(); + sim.match_cycle = meta.match_cycle; + sim.loop_playback = meta.loop_playback; + sim.fundamental_hz = meta.fundamental_hz; + sim.fundamental_hz_text = format!("{}", meta.fundamental_hz); + sim.tracks = tracks; + if !sim.tracks.is_empty() { + sim.selected_track = Some(0); + sim.expression_text = sim.tracks[0].expression.source.clone(); + } + maybe_recompute_buffer(&mut sim); + self.simulation = Some(sim); + self.source_path = meta.source_fem.as_ref().map(PathBuf::from); + self.source_label = path + .file_name() + .and_then(|n| n.to_str()) + .unwrap_or("?") + .to_string(); + self.status = format!("loaded session: {}", path.display()); + Task::none() + } + fn update(&mut self, msg: Message) -> Task { match msg { Message::OpenFem => { let picked = rfd::FileDialog::new() - .add_filter("FEMM magnetostatic", &["fem", "FEM"]) + .add_filter("FEMM magnetostatic or session", &["fem", "FEM", "femsess"]) .add_filter("All files", &["*"]) .pick_file(); @@ -213,11 +291,19 @@ impl App { .and_then(|n| n.to_str()) .unwrap_or("?") .to_string(); + let is_session = path.extension() + .and_then(|e| e.to_str()) + .map(|e| e.eq_ignore_ascii_case("femsess")) + .unwrap_or(false); + if is_session { + return self.load_session_from_path(&path); + } match FemmDoc::open(&path) { Ok(d) => { self.doc = d; self.mesh = None; self.solution = None; + self.pristine_doc = None; self.source_label = label; self.source_path = Some(path); self.status = String::new(); @@ -232,6 +318,7 @@ impl App { self.doc = FemmDoc::default(); self.mesh = None; self.solution = None; + self.pristine_doc = None; self.source_label = String::from("untitled"); self.source_path = None; self.error = None; @@ -435,6 +522,28 @@ impl App { self.solution = None; self.status = format!("block label {idx} at ({:.3}, {:.3})", world.0, world.1); } + Tool::AddProbe => { + if self.simulation.is_none() { + self.simulation = Some(new_simulation(self.doc.clone(), String::new())); + } + let sim = self.simulation.as_mut().unwrap(); + let idx = sim.probes.len(); + let p = probe::Probe { + x: world.0, y: world.1, + label: format!("p{}", idx + 1), + color: probe::palette_color(idx), + mode: probe::ProbeMode::default(), + }; + let mut samples = std::collections::HashMap::new(); + for (&fi, sol) in &sim.lut { + if let Some(bxby) = sol.sample_b_at(p.x, p.y) { + samples.insert(fi, bxby); + } + } + sim.probes.push(p); + sim.probe_samples.push(samples); + self.status = format!("probe p{} placed at ({:.3}, {:.3})", idx + 1, world.0, world.1); + } Tool::Select | Tool::AddSegment => {} } } @@ -495,6 +604,9 @@ impl App { body: String::from("select at least one segment or arc (left-click or marquee). the track moves every node along the selection between its leftmost and rightmost endpoints."), }); } else { + if self.pristine_doc.is_none() { + self.pristine_doc = Some(self.doc.clone()); + } let subdivisions = self.simulation.as_ref() .map(|s| s.subdivisions) .unwrap_or(SIM_DEFAULT_SUBDIVISIONS); @@ -537,6 +649,7 @@ impl App { Ok(mut track) => { track.edge_endpoints = edge_endpoints; track.edge_subdivisions = edge_subdivisions; + track.subdivisions = subdivisions; if self.simulation.is_none() { self.simulation = Some(new_simulation(self.doc.clone(), expr_text)); } @@ -671,20 +784,27 @@ impl App { self.status = format!("subdivisions for new tracks set to {n}"); return Task::none(); }; - let mut track = sim.tracks[selected].clone(); - let rebuild = kinematic::rebuild_track_subdivisions(&mut track, &mut self.doc, n); + sim.tracks[selected].subdivisions = n; + let Some(pristine) = self.pristine_doc.clone() else { + self.error = Some(ErrorReport { + title: String::from("no pristine doc stashed"), + body: String::from("this session was loaded without a pristine snapshot. delete and recreate all tracks (or load from a newer session) to enable subdivision rebuild."), + }); + return Task::none(); + }; + let sim_mut = self.simulation.as_mut().unwrap(); + let rebuild = kinematic::rebuild_all_tracks_from_pristine(&mut self.doc, &pristine, &mut sim_mut.tracks); match rebuild { Ok(()) => { - let sim = self.simulation.as_mut().unwrap(); - sim.tracks[selected] = track; - sim.base_doc = self.doc.clone(); - invalidate_lut(sim); - self.status = format!("rebuilt track {} with {n} subdivisions per edge", selected + 1); + let sim_mut = self.simulation.as_mut().unwrap(); + sim_mut.base_doc = self.doc.clone(); + invalidate_lut(sim_mut); + self.status = format!("rebuilt all tracks from pristine; track {} now at {n} subdivisions", selected + 1); } Err(e) => { self.error = Some(ErrorReport { title: String::from("subdivision rebuild failed"), - body: format!("{e}\n\ndelete and recreate the track to enable rebuild for tracks created before this feature."), + body: e, }); } } @@ -789,6 +909,14 @@ impl App { if frame_idx == sim.frame_idx && sim.current.is_none() { sim.current = Some(sol.clone()); } + while sim.probe_samples.len() < sim.probes.len() { + sim.probe_samples.push(std::collections::HashMap::new()); + } + for (i, p) in sim.probes.iter().enumerate() { + if let Some(bxby) = sol.sample_b_at(p.x, p.y) { + sim.probe_samples[i].insert(frame_idx, bxby); + } + } sim.lut.insert(frame_idx, sol); } Err(report) => { @@ -975,7 +1103,7 @@ impl App { .save_file(); if let Some(path) = picked { let meta = build_sim_meta(sim, self.source_path.as_deref()); - let text = session::serialize(&sim.base_doc, &meta); + let text = session::serialize(&sim.base_doc, self.pristine_doc.as_ref(), &meta); match std::fs::write(&path, text) { Ok(()) => self.status = format!("saved session: {}", path.display()), Err(e) => self.error = Some(ErrorReport { @@ -990,64 +1118,7 @@ impl App { .add_filter("femm42 session", &["femsess"]) .pick_file(); let Some(path) = picked else { return Task::none(); }; - let text = match std::fs::read_to_string(&path) { - Ok(t) => t, - Err(e) => { - self.error = Some(ErrorReport { - title: String::from("read session failed"), - body: format!("{path:?}\n\n{e}"), - }); - return Task::none(); - } - }; - let (base_doc, meta) = match session::parse(&text) { - Ok(pair) => pair, - Err(e) => { - self.error = Some(ErrorReport { - title: String::from("parse session failed"), - body: format!("{path:?}\n\n{e}"), - }); - return Task::none(); - } - }; - let tracks = match sim_meta::resolve_tracks(&meta, &base_doc) { - Ok(t) => t, - Err(e) => { - self.error = Some(ErrorReport { - title: String::from("session tracks do not match embedded doc"), - body: format!("{path:?}\n\n{e}"), - }); - return Task::none(); - } - }; - self.doc = base_doc.clone(); - self.mesh = None; - self.solution = None; - self.error = None; - let mut sim = new_simulation(base_doc, String::new()); - if let Some(secs) = spice::parse_spice(&meta.dt_text) { sim.dt = secs; } - sim.dt_text = meta.dt_text.clone(); - if let Some(secs) = spice::parse_spice(&meta.interval_text) { - sim.wall_interval = Duration::from_secs_f64(secs); - } - sim.interval_text = meta.interval_text.clone(); - sim.subdivisions = meta.subdivisions; - sim.subdivisions_text = meta.subdivisions.to_string(); - sim.buffer_size = meta.buffer_size; - sim.buffer_size_text = meta.buffer_size.to_string(); - sim.match_cycle = meta.match_cycle; - sim.loop_playback = meta.loop_playback; - sim.fundamental_hz = meta.fundamental_hz; - sim.fundamental_hz_text = format!("{}", meta.fundamental_hz); - sim.tracks = tracks; - if !sim.tracks.is_empty() { - sim.selected_track = Some(0); - sim.expression_text = sim.tracks[0].expression.source.clone(); - } - maybe_recompute_buffer(&mut sim); - self.simulation = Some(sim); - self.source_path = meta.source_fem.as_ref().map(PathBuf::from); - self.status = format!("loaded session: {}", path.display()); + return self.load_session_from_path(&path); } Message::SimExportWebM => { let Some(sim) = self.simulation.as_ref() else { @@ -1094,6 +1165,8 @@ impl App { height: 1080, contour, crf: 18, + probes: sim.probes.clone(), + probe_samples: sim.probe_samples.clone(), progress: std::sync::Arc::clone(&progress), }; self.export_progress = Some(progress); @@ -1135,6 +1208,61 @@ impl App { }; } } + Message::ProbeSetMode(i, mode) => { + if let Some(sim) = self.simulation.as_mut() { + if let Some(p) = sim.probes.get_mut(i) { + p.mode = mode; + } + } + } + Message::ProbeRemove(i) => { + if let Some(sim) = self.simulation.as_mut() { + if i < sim.probes.len() { + sim.probes.remove(i); + if i < sim.probe_samples.len() { + sim.probe_samples.remove(i); + } + } + } + } + Message::ProbeExportWav(i) => { + let Some(sim) = self.simulation.as_ref() else { return Task::none(); }; + let Some(p) = sim.probes.get(i) else { return Task::none(); }; + let Some(samples) = sim.probe_samples.get(i) else { return Task::none(); }; + if samples.is_empty() { + self.error = Some(ErrorReport { + title: String::from("no samples"), + body: String::from("fill the buffer (run the sim with loop enabled) before exporting a wav."), + }); + return Task::none(); + } + let default_name = format!("{}_{}.wav", p.label, match p.mode { + probe::ProbeMode::Magnitude => "magB", + probe::ProbeMode::Bx => "Bx", + probe::ProbeMode::By => "By", + probe::ProbeMode::Angle => "ang", + }); + let picked = rfd::FileDialog::new() + .add_filter("WAV audio", &["wav"]) + .set_file_name(&default_name) + .save_file(); + let Some(path) = picked else { return Task::none(); }; + let mut ordered: Vec<(i64, f64)> = samples.iter() + .map(|(&fi, &(bx, by))| (fi, p.mode.extract(bx, by))) + .collect(); + ordered.sort_by_key(|s| s.0); + let values: Vec = ordered.into_iter().map(|(_, v)| v).collect(); + let sample_rate = (1.0 / sim.dt).round().max(1.0) as u32; + match write_wav_mono_16bit_normalized(&path, &values, sample_rate) { + Ok(()) => self.status = format!( + "wrote {}: {} samples at {} Hz", path.display(), values.len(), sample_rate, + ), + Err(e) => self.error = Some(ErrorReport { + title: String::from("wav export failed"), + body: format!("{path:?}\n\n{e}"), + }), + } + } Message::SimToggleLoop(checked) => { if let Some(sim) = self.simulation.as_mut() { sim.loop_playback = checked; @@ -1174,10 +1302,10 @@ impl App { .size(12); let file_group = row![ - icon_button(ICON_FILE_NEW, "New problem", Some(Message::NewDoc)), - icon_button(ICON_FILE_OPEN, "Open .fem", Some(Message::OpenFem)), - icon_button(ICON_FILE_SAVE, "Save", Some(Message::SaveDoc)), - icon_button(ICON_FILE_SAVE_AS, "Save As", Some(Message::SaveDocAs)), + icon_button(ICON_FILE_NEW, "New problem", Some(Message::NewDoc)), + icon_button(ICON_FILE_OPEN, "Open .fem or .femsess", Some(Message::OpenFem)), + icon_button(ICON_FILE_SAVE, "Save", Some(Message::SaveDoc)), + icon_button(ICON_FILE_SAVE_AS, "Save As", Some(Message::SaveDocAs)), ].spacing(2); let tool_group = row![ @@ -1185,6 +1313,7 @@ impl App { tool_icon_button(ICON_TOOL_NODE, "Add Node", Tool::AddNode, self.tool), tool_icon_button(ICON_TOOL_SEGMENT, "Add Segment", Tool::AddSegment, self.tool), tool_icon_button(ICON_TOOL_LABEL, "Add Label", Tool::AddBlockLabel, self.tool), + text_button("Probe", Message::SelectTool(Tool::AddProbe)), ].spacing(2); let edit_group = row![ @@ -1247,9 +1376,12 @@ impl App { .and_then(|s| s.current.as_ref()) .or(self.solution.as_ref()); + let probes_slice: &[probe::Probe] = self.simulation.as_ref() + .map(|s| s.probes.as_slice()) + .unwrap_or(&[]); let canvas = doc_canvas::view( &self.doc, self.tool, self.mesh.as_ref(), active_solution, self.render_mode, - self.view_state, self.show_grid, self.zoom_window_active, + self.view_state, self.show_grid, self.zoom_window_active, probes_slice, ).map(Message::Canvas); let canvas_row = row![canvas, view_strip] @@ -1268,6 +1400,42 @@ impl App { body = body.push(simulation_panel(sim)); } body = body.push(canvas_row); + if let Some(sim) = &self.simulation { + if !sim.probes.is_empty() { + static EMPTY_MAP: std::sync::LazyLock> + = std::sync::LazyLock::new(std::collections::HashMap::new); + let mut pane_row = row![].spacing(8); + for (i, p) in sim.probes.iter().enumerate() { + let samples = sim.probe_samples.get(i).unwrap_or(&*EMPTY_MAP); + let mode_row = row![ + button(text(p.label.clone()).size(11)).style(button::secondary), + button(text("|B|").size(10)).on_press(Message::ProbeSetMode(i, probe::ProbeMode::Magnitude)) + .style(if p.mode == probe::ProbeMode::Magnitude { button::primary } else { button::secondary }), + button(text("Bx").size(10)).on_press(Message::ProbeSetMode(i, probe::ProbeMode::Bx)) + .style(if p.mode == probe::ProbeMode::Bx { button::primary } else { button::secondary }), + button(text("By").size(10)).on_press(Message::ProbeSetMode(i, probe::ProbeMode::By)) + .style(if p.mode == probe::ProbeMode::By { button::primary } else { button::secondary }), + button(text("ang").size(10)).on_press(Message::ProbeSetMode(i, probe::ProbeMode::Angle)) + .style(if p.mode == probe::ProbeMode::Angle { button::primary } else { button::secondary }), + button(text("wav").size(10)).on_press(Message::ProbeExportWav(i)).style(button::secondary), + button(text("x").size(10)).on_press(Message::ProbeRemove(i)).style(button::danger), + ].spacing(2); + let plot_canvas: Element<'_, Message> = + plot::view::(p, samples, sim.frame_idx, sim.buffer_size.max(1)); + let pane = column![mode_row, plot_canvas] + .spacing(2) + .width(Length::FillPortion(1)) + .height(Length::Fill); + pane_row = pane_row.push(pane); + } + body = body.push( + container(pane_row) + .width(Length::Fill) + .height(Length::FillPortion(1)) + .padding(4) + ); + } + } if !self.status.is_empty() { body = body.push(text(&self.status).size(12)); } @@ -1319,6 +1487,8 @@ fn new_simulation(base_doc: FemmDoc, expression_text: String) -> Simulation { fundamental_hz: SIM_DEFAULT_FUNDAMENTAL_HZ, fundamental_hz_text: format!("{SIM_DEFAULT_FUNDAMENTAL_HZ}"), loop_playback: false, + probes: Vec::new(), + probe_samples: Vec::new(), } } @@ -1339,6 +1509,7 @@ fn build_sim_meta(sim: &Simulation, source: Option<&Path>) -> sim_meta::SimMeta axis: t.axis, expression: t.expression.source.clone(), edges, + subdivisions: t.subdivisions, } }).collect(); sim_meta::SimMeta { @@ -1370,6 +1541,46 @@ fn invalidate_lut(sim: &mut Simulation) { sim.pending.clear(); sim.generation = sim.generation.wrapping_add(1); sim.current = None; + for samples in sim.probe_samples.iter_mut() { + samples.clear(); + } +} + +/// writes mono 16-bit PCM WAV with the given sample rate, normalising the input so the peak magnitude maps to i16 full-scale. +fn write_wav_mono_16bit_normalized(path: &Path, values: &[f64], sample_rate: u32) -> std::io::Result<()> { + use std::io::Write; + let mut peak = 0.0_f64; + for &v in values { + let a = v.abs(); + if a.is_finite() && a > peak { peak = a; } + } + let scale = if peak > 0.0 { (i16::MAX as f64 - 1.0) / peak } else { 0.0 }; + let mut pcm: Vec = Vec::with_capacity(values.len()); + for &v in values { + let s = if v.is_finite() { (v * scale).round() } else { 0.0 }; + pcm.push(s.clamp(i16::MIN as f64, i16::MAX as f64) as i16); + } + let data_len = (pcm.len() * 2) as u32; + let chunk_size = 36 + data_len; + let byte_rate = sample_rate * 2; + let mut f = std::fs::File::create(path)?; + f.write_all(b"RIFF")?; + f.write_all(&chunk_size.to_le_bytes())?; + f.write_all(b"WAVE")?; + f.write_all(b"fmt ")?; + f.write_all(&16u32.to_le_bytes())?; + f.write_all(&1u16.to_le_bytes())?; + f.write_all(&1u16.to_le_bytes())?; + f.write_all(&sample_rate.to_le_bytes())?; + f.write_all(&byte_rate.to_le_bytes())?; + f.write_all(&2u16.to_le_bytes())?; + f.write_all(&16u16.to_le_bytes())?; + f.write_all(b"data")?; + f.write_all(&data_len.to_le_bytes())?; + for s in &pcm { + f.write_all(&s.to_le_bytes())?; + } + Ok(()) } /// caps concurrent solve jobs at cores x 4, clamped to [8, 64]. diff --git a/crates/femm-app/src/plot.rs b/crates/femm-app/src/plot.rs new file mode 100644 index 0000000..9e288eb --- /dev/null +++ b/crates/femm-app/src/plot.rs @@ -0,0 +1,161 @@ +//! per-probe waveform plot widget rendered as an iced canvas, color-matched to the probe marker on the main canvas. + +use iced::widget::canvas::{self, Canvas, Frame, Geometry, Path, Stroke, Text}; +use iced::{Color, Element, Length, Point, Rectangle, Renderer, Theme, mouse}; + +use crate::probe::Probe; + +/// renders one probe's plot pane: a labeled waveform of the chosen signal across every frame currently in samples, with a vertical bar marking the current playback frame. +pub fn view<'a, Msg: 'a + Clone>( + probe: &'a Probe, + samples: &'a std::collections::HashMap, + current_frame: i64, + buffer_size: usize, +) -> Element<'a, Msg> { + Canvas::new(PlotProgram { + probe, + samples, + current_frame, + buffer_size, + }) + .width(Length::Fill) + .height(Length::Fill) + .into() +} + +struct PlotProgram<'a> { + probe: &'a Probe, + samples: &'a std::collections::HashMap, + current_frame: i64, + buffer_size: usize, +} + +impl<'a, Msg> canvas::Program for PlotProgram<'a> { + type State = (); + + fn draw( + &self, + _state: &Self::State, + renderer: &Renderer, + _theme: &Theme, + bounds: Rectangle, + _cursor: mouse::Cursor, + ) -> Vec> { + let mut frame = Frame::new(renderer, bounds.size()); + let bg = Color::from_rgba8(18, 20, 26, 1.0); + frame.fill_rectangle(Point::ORIGIN, bounds.size(), bg); + + let pad_l = 8.0_f32; + let pad_r = 8.0_f32; + let pad_t = 18.0_f32; + let pad_b = 16.0_f32; + let plot_w = (bounds.width - pad_l - pad_r).max(1.0); + let plot_h = (bounds.height - pad_t - pad_b).max(1.0); + + let probe_color = Color::from_rgba8( + self.probe.color[0], self.probe.color[1], self.probe.color[2], 1.0, + ); + + frame.fill_text(Text { + content: format!("{} - {}", self.probe.label, self.probe.mode.label()), + position: Point::new(pad_l, 2.0), + color: probe_color, + size: 11.0.into(), + ..Text::default() + }); + + let buf = self.buffer_size.max(1); + let mut points: Vec<(i64, f64)> = Vec::with_capacity(self.samples.len()); + for (&fi, &(bx, by)) in self.samples { + points.push((fi, self.probe.mode.extract(bx, by))); + } + if points.is_empty() { + frame.fill_text(Text { + content: String::from("no samples yet"), + position: Point::new(pad_l + 4.0, pad_t + 4.0), + color: Color::from_rgba8(140, 140, 150, 1.0), + size: 10.0.into(), + ..Text::default() + }); + return vec![frame.into_geometry()]; + } + let mut ymin = f64::INFINITY; + let mut ymax = f64::NEG_INFINITY; + for &(_, v) in &points { + if v < ymin { ymin = v; } + if v > ymax { ymax = v; } + } + if !ymin.is_finite() || !ymax.is_finite() { + ymin = -1.0; ymax = 1.0; + } + if (ymax - ymin).abs() < 1e-12 { + ymin -= 1.0; ymax += 1.0; + } + + points.sort_by_key(|p| p.0); + + let x_for = |fi: i64| -> f32 { + let t = (fi.rem_euclid(buf as i64) as f64) / (buf as f64).max(1.0); + pad_l + (t as f32) * plot_w + }; + let y_for = |v: f64| -> f32 { + let t = (v - ymin) / (ymax - ymin); + pad_t + plot_h - (t as f32) * plot_h + }; + + frame.stroke( + &Path::line( + Point::new(pad_l, pad_t + plot_h), + Point::new(pad_l + plot_w, pad_t + plot_h), + ), + Stroke::default().with_width(0.5).with_color(Color::from_rgba8(80, 80, 90, 1.0)), + ); + let mid_y = y_for((ymin + ymax) * 0.5); + frame.stroke( + &Path::line( + Point::new(pad_l, mid_y), + Point::new(pad_l + plot_w, mid_y), + ), + Stroke::default().with_width(0.5).with_color(Color::from_rgba8(60, 60, 70, 1.0)), + ); + + let mut path = canvas::path::Builder::new(); + let mut started = false; + for (fi, v) in &points { + let x = x_for(*fi); + let y = y_for(*v); + if !started { path.move_to(Point::new(x, y)); started = true; } + else { path.line_to(Point::new(x, y)); } + } + if started { + frame.stroke(&path.build(), + Stroke::default().with_width(1.4).with_color(probe_color)); + } + + let cx = x_for(self.current_frame); + frame.stroke( + &Path::line( + Point::new(cx, pad_t), + Point::new(cx, pad_t + plot_h), + ), + Stroke::default().with_width(1.0).with_color(Color::from_rgba8(255, 255, 255, 0.55)), + ); + + frame.fill_text(Text { + content: format!("{:.3e}", ymax), + position: Point::new(pad_l, pad_t - 12.0), + color: Color::from_rgba8(150, 150, 160, 1.0), + size: 9.0.into(), + ..Text::default() + }); + frame.fill_text(Text { + content: format!("{:.3e}", ymin), + position: Point::new(pad_l, pad_t + plot_h + 2.0), + color: Color::from_rgba8(150, 150, 160, 1.0), + size: 9.0.into(), + ..Text::default() + }); + + vec![frame.into_geometry()] + } +} diff --git a/crates/femm-app/src/probe.rs b/crates/femm-app/src/probe.rs new file mode 100644 index 0000000..7a609f9 --- /dev/null +++ b/crates/femm-app/src/probe.rs @@ -0,0 +1,55 @@ +//! point-sample probes for the magnetic field, time-series plotted per frame in a panel below the canvas. + +/// signal extracted from the (bx, by) sample at a probe location for waveform plotting. +#[derive(Debug, Clone, Copy, Default, PartialEq, Eq)] +pub enum ProbeMode { + #[default] + Magnitude, + Bx, + By, + Angle, +} + +impl ProbeMode { + pub fn label(self) -> &'static str { + match self { + ProbeMode::Magnitude => "|B|", + ProbeMode::Bx => "Bx", + ProbeMode::By => "By", + ProbeMode::Angle => "atan2(By,Bx)", + } + } + pub fn extract(self, bx: f64, by: f64) -> f64 { + match self { + ProbeMode::Magnitude => (bx * bx + by * by).sqrt(), + ProbeMode::Bx => bx, + ProbeMode::By => by, + ProbeMode::Angle => by.atan2(bx), + } + } +} + +/// one user-placed probe with a world position, display color, label, and plot-mode selector. +#[derive(Debug, Clone)] +pub struct Probe { + pub x: f64, + pub y: f64, + pub label: String, + pub color: [u8; 4], + pub mode: ProbeMode, +} + +/// returns a saturated palette color from a small cycling palette suitable for probe markers and matching plot lines. +pub fn palette_color(idx: usize) -> [u8; 4] { + const PALETTE: &[[u8; 4]] = &[ + [231, 76, 60, 255], + [241, 196, 15, 255], + [ 46, 204, 113, 255], + [ 52, 152, 219, 255], + [155, 89, 182, 255], + [230, 126, 34, 255], + [ 26, 188, 156, 255], + [255, 105, 180, 255], + ]; + PALETTE[idx % PALETTE.len()] +} diff --git a/crates/femm-app/src/session.rs b/crates/femm-app/src/session.rs index 41e7193..d1f9b6a 100644 --- a/crates/femm-app/src/session.rs +++ b/crates/femm-app/src/session.rs @@ -4,13 +4,22 @@ use femm_doc_mag::FemmDoc; use crate::sim_meta::{self, SimMeta}; const SESSION_HEADER: &str = "# femm42 session v1\n"; +const PRISTINE_OPEN: &str = "[pristine-doc]\n"; +const PRISTINE_CLOSE: &str = "[/pristine-doc]\n"; const DOC_OPEN: &str = "[base-doc]\n"; const DOC_CLOSE: &str = "[/base-doc]\n"; -/// joins the base .fem text and the sim metadata into one delimited session blob. -pub fn serialize(base_doc: &FemmDoc, meta: &SimMeta) -> String { +/// joins the optional pristine .fem text, the live base .fem text, and the sim metadata into one delimited session blob. +pub fn serialize(base_doc: &FemmDoc, pristine_doc: Option<&FemmDoc>, meta: &SimMeta) -> String { let mut out = String::new(); out.push_str(SESSION_HEADER); + if let Some(p) = pristine_doc { + out.push_str(PRISTINE_OPEN); + let pt = p.write(); + out.push_str(&pt); + if !pt.ends_with('\n') { out.push('\n'); } + out.push_str(PRISTINE_CLOSE); + } out.push_str(DOC_OPEN); let doc_text = base_doc.write(); out.push_str(&doc_text); @@ -20,8 +29,19 @@ pub fn serialize(base_doc: &FemmDoc, meta: &SimMeta) -> String { out } -/// extracts the base FemmDoc and SimMeta from a session blob, failing if the markers or sub-payloads are malformed. -pub fn parse(text: &str) -> Result<(FemmDoc, SimMeta), String> { +/// extracts the optional pristine FemmDoc, the base FemmDoc, and SimMeta from a session blob. +pub fn parse(text: &str) -> Result<(FemmDoc, Option, SimMeta), String> { + let pristine = if let Some(p_start) = text.find(PRISTINE_OPEN) { + let body_start = p_start + PRISTINE_OPEN.len(); + let after_body = text[body_start..] + .find(PRISTINE_CLOSE) + .ok_or_else(|| String::from("missing [/pristine-doc] marker"))?; + let p_text = &text[body_start..body_start + after_body]; + Some(FemmDoc::parse(p_text).map_err(|e| format!("parse [pristine-doc] failed: {e}"))?) + } else { + None + }; + let doc_start = text.find(DOC_OPEN).ok_or_else(|| String::from("missing [base-doc] section"))?; let body_start = doc_start + DOC_OPEN.len(); let after_body = text[body_start..] @@ -33,5 +53,5 @@ pub fn parse(text: &str) -> Result<(FemmDoc, SimMeta), String> { let doc = FemmDoc::parse(doc_text) .map_err(|e| format!("parse [base-doc] failed: {e}"))?; let meta = sim_meta::parse_meta(meta_text)?; - Ok((doc, meta)) + Ok((doc, pristine, meta)) } diff --git a/crates/femm-app/src/sim_meta.rs b/crates/femm-app/src/sim_meta.rs index c645bc3..04730ae 100644 --- a/crates/femm-app/src/sim_meta.rs +++ b/crates/femm-app/src/sim_meta.rs @@ -25,6 +25,7 @@ pub struct TrackMeta { pub axis: Axis, pub expression: String, pub edges: Vec<(usize, usize, Vec)>, + pub subdivisions: usize, } /// renders a SimMeta into the on-disk text format. @@ -51,6 +52,7 @@ pub fn serialize_meta(meta: &SimMeta) -> String { s.push_str(&format!("anchor_b = {}\n", track.anchor_b)); s.push_str(&format!("axis = {}\n", axis_to_str(track.axis))); s.push_str(&format!("expression = {}\n", track.expression)); + s.push_str(&format!("subdivisions = {}\n", track.subdivisions)); s.push_str("members ="); for n in &track.members { s.push_str(&format!(" {n}")); } s.push_str("\n"); @@ -95,6 +97,7 @@ pub fn parse_meta(text: &str) -> Result { axis: Axis::PlusY, expression: String::from("0"), edges: Vec::new(), + subdivisions: 0, }); Section::Track } other => return Err(format!("unknown section {other:?} at line {}", line_no + 1)), }; @@ -134,6 +137,9 @@ pub fn parse_meta(text: &str) -> Result { track.members.push(parse_usize(tok, line_no)?); } } + "subdivisions" => { + track.subdivisions = parse_usize(value, line_no)?; + } "edge" => { let mut iter = value.split_whitespace(); let a = iter.next().ok_or_else(|| format!("edge missing first endpoint at line {}", line_no + 1))?; @@ -203,6 +209,7 @@ pub fn resolve_tracks(meta: &SimMeta, doc: &FemmDoc) -> Result, Strin expression: expr, edge_endpoints, edge_subdivisions, + subdivisions: if tm.subdivisions >= 2 { tm.subdivisions } else { 20 }, }); } Ok(out) diff --git a/crates/femm-doc-mag/src/ans.rs b/crates/femm-doc-mag/src/ans.rs index 8fb26d5..9125d9c 100644 --- a/crates/femm-doc-mag/src/ans.rs +++ b/crates/femm-doc-mag/src/ans.rs @@ -175,6 +175,21 @@ impl MagSolution { + el.b2.re * el.b2.re + el.b2.im * el.b2.im).sqrt() } + /// looks up the mesh triangle containing world point (x, y) and returns its (bx, by) real components, or None when the point lies outside every element. + pub fn sample_b_at(&self, x: f64, y: f64) -> Option<(f64, f64)> { + for el in &self.mesh_elements { + let (Some(a), Some(b), Some(c)) = ( + self.mesh_nodes.get(el.n[0] as usize), + self.mesh_nodes.get(el.n[1] as usize), + self.mesh_nodes.get(el.n[2] as usize), + ) else { continue }; + if point_in_triangle((x, y), (a.x, a.y), (b.x, b.y), (c.x, c.y)) { + return Some((el.b1.re, el.b2.re)); + } + } + None + } + /// (min, max) of |B| across every element, skipping NaN. pub fn b_magnitude_range(&self) -> (f64, f64) { let mut lo = f64::INFINITY; @@ -335,6 +350,19 @@ fn take_i32(s: &str) -> (i32, &str) { (head.parse::().unwrap_or(0), tail) } +/// returns true when point p lies inside (or on the boundary of) the triangle a-b-c, via the half-plane sign test. +fn point_in_triangle(p: (f64, f64), a: (f64, f64), b: (f64, f64), c: (f64, f64)) -> bool { + let sign = |p1: (f64, f64), p2: (f64, f64), p3: (f64, f64)| -> f64 { + (p1.0 - p3.0) * (p2.1 - p3.1) - (p2.0 - p3.0) * (p1.1 - p3.1) + }; + let d1 = sign(p, a, b); + let d2 = sign(p, b, c); + let d3 = sign(p, c, a); + let neg = d1 < 0.0 || d2 < 0.0 || d3 < 0.0; + let pos = d1 > 0.0 || d2 > 0.0 || d3 > 0.0; + !(neg && pos) +} + #[cfg(test)] mod tests { use super::*; diff --git a/examples/guitar_bar_and_discs.femsess b/examples/guitar_bar_and_discs.femsess new file mode 100644 index 0000000..97860a1 --- /dev/null +++ b/examples/guitar_bar_and_discs.femsess @@ -0,0 +1,1134 @@ +# femm42 session v1 +[pristine-doc] +[Format] = 4.00000000000000000e0 +[Frequency] = 0.00000000000000000e0 +[Precision] = 1.00000000000000002e-8 +[MinAngle] = 2.50000000000000000e1 +[DoSmartMesh] = 1 +[Depth] = 5.00000000000000000e0 +[LengthUnits] = millimeters +[ProblemType] = planar +[Coordinates] = cartesian +[ACSolver] = 0 +[PrevType] = 0 +[PrevSoln] = "" +[Comment] = "six plain steel strings (.010-.046 inch gauges) over a 25.5-inch scale, anchored to a NdFeB bar magnet at the left (mag_dir 0, N facing the strings) and a passive steel block at the right. small NdFeB disc magnets sit above and below each string at 5/7 of the way toward the steel block, all magnetised in -y so the disc above each string has its N face downward and the disc below has its S face upward, driving a strong vertical bias field through each string at the probe point. flip every disc mag_dir from 270 to 90 to invert the bias direction." +[PointProps] = 0 +[BdryProps] = 1 + + = "A=0" + = 0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + +[BlockProps] = 3 + + = "Air" + = 1.00000000000000000e0 + = 1.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0 + = 1.00000000000000000e0 + = 0 + = 0.00000000000000000e0 + = 0 + + + = "NdFeB" + = 1.05000000000000004e0 + = 1.05000000000000004e0 + = 9.15000000000000000e5 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 6.67000000000000037e-1 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0 + = 1.00000000000000000e0 + = 0 + = 0.00000000000000000e0 + = 0 + + + = "Steel" + = 2.50000000000000000e3 + = 2.50000000000000000e3 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 5.79999999999999982e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0 + = 1.00000000000000000e0 + = 0 + = 0.00000000000000000e0 + = 0 + +[CircuitProps] = 0 +[NumPoints] = 84 +-1.00000000000000000e2 -1.00000000000000000e2 0 0 +7.50000000000000000e2 -1.00000000000000000e2 0 0 +7.50000000000000000e2 1.00000000000000000e2 0 0 +-1.00000000000000000e2 1.00000000000000000e2 0 0 +-3.00000000000000000e1 -3.00000000000000000e1 0 0 +-3.00000000000000000e1 3.00000000000000000e1 0 0 +0.00000000000000000e0 -3.00000000000000000e1 0 0 +0.00000000000000000e0 -2.68339999999999996e1 0 0 +0.00000000000000000e0 -2.56660000000000004e1 0 0 +0.00000000000000000e0 -1.62070000000000007e1 0 0 +0.00000000000000000e0 -1.52929999999999993e1 0 0 +0.00000000000000000e0 -5.58000000000000007e0 0 0 +0.00000000000000000e0 -4.91999999999999993e0 0 0 +0.00000000000000000e0 5.03399999999999981e0 0 0 +0.00000000000000000e0 5.46600000000000019e0 0 0 +0.00000000000000000e0 1.55850000000000009e1 0 0 +0.00000000000000000e0 1.59149999999999991e1 0 0 +0.00000000000000000e0 2.61230000000000011e1 0 0 +0.00000000000000000e0 2.63769999999999989e1 0 0 +0.00000000000000000e0 3.00000000000000000e1 0 0 +6.47700000000000045e2 -3.00000000000000000e1 0 0 +6.47700000000000045e2 -2.68339999999999996e1 0 0 +6.47700000000000045e2 -2.56660000000000004e1 0 0 +6.47700000000000045e2 -1.62070000000000007e1 0 0 +6.47700000000000045e2 -1.52929999999999993e1 0 0 +6.47700000000000045e2 -5.58000000000000007e0 0 0 +6.47700000000000045e2 -4.91999999999999993e0 0 0 +6.47700000000000045e2 5.03399999999999981e0 0 0 +6.47700000000000045e2 5.46600000000000019e0 0 0 +6.47700000000000045e2 1.55850000000000009e1 0 0 +6.47700000000000045e2 1.59149999999999991e1 0 0 +6.47700000000000045e2 2.61230000000000011e1 0 0 +6.47700000000000045e2 2.63769999999999989e1 0 0 +6.47700000000000045e2 3.00000000000000000e1 0 0 +6.77700000000000045e2 -3.00000000000000000e1 0 0 +6.77700000000000045e2 3.00000000000000000e1 0 0 +4.61139999999999986e2 2.73769999999999989e1 0 0 +4.64139999999999986e2 2.73769999999999989e1 0 0 +4.64139999999999986e2 2.93769999999999989e1 0 0 +4.61139999999999986e2 2.93769999999999989e1 0 0 +4.61139999999999986e2 2.31230000000000011e1 0 0 +4.64139999999999986e2 2.31230000000000011e1 0 0 +4.64139999999999986e2 2.51230000000000011e1 0 0 +4.61139999999999986e2 2.51230000000000011e1 0 0 +4.61139999999999986e2 1.69149999999999991e1 0 0 +4.64139999999999986e2 1.69149999999999991e1 0 0 +4.64139999999999986e2 1.89149999999999991e1 0 0 +4.61139999999999986e2 1.89149999999999991e1 0 0 +4.61139999999999986e2 1.25850000000000009e1 0 0 +4.64139999999999986e2 1.25850000000000009e1 0 0 +4.64139999999999986e2 1.45850000000000009e1 0 0 +4.61139999999999986e2 1.45850000000000009e1 0 0 +4.61139999999999986e2 6.46600000000000019e0 0 0 +4.64139999999999986e2 6.46600000000000019e0 0 0 +4.64139999999999986e2 8.46599999999999930e0 0 0 +4.61139999999999986e2 8.46599999999999930e0 0 0 +4.61139999999999986e2 2.03399999999999981e0 0 0 +4.64139999999999986e2 2.03399999999999981e0 0 0 +4.64139999999999986e2 4.03399999999999981e0 0 0 +4.61139999999999986e2 4.03399999999999981e0 0 0 +4.61139999999999986e2 -3.91999999999999993e0 0 0 +4.64139999999999986e2 -3.91999999999999993e0 0 0 +4.64139999999999986e2 -1.91999999999999993e0 0 0 +4.61139999999999986e2 -1.91999999999999993e0 0 0 +4.61139999999999986e2 -8.58000000000000007e0 0 0 +4.64139999999999986e2 -8.58000000000000007e0 0 0 +4.64139999999999986e2 -6.58000000000000007e0 0 0 +4.61139999999999986e2 -6.58000000000000007e0 0 0 +4.61139999999999986e2 -1.42929999999999993e1 0 0 +4.64139999999999986e2 -1.42929999999999993e1 0 0 +4.64139999999999986e2 -1.22929999999999993e1 0 0 +4.61139999999999986e2 -1.22929999999999993e1 0 0 +4.61139999999999986e2 -1.92070000000000007e1 0 0 +4.64139999999999986e2 -1.92070000000000007e1 0 0 +4.64139999999999986e2 -1.72070000000000007e1 0 0 +4.61139999999999986e2 -1.72070000000000007e1 0 0 +4.61139999999999986e2 -2.46660000000000004e1 0 0 +4.64139999999999986e2 -2.46660000000000004e1 0 0 +4.64139999999999986e2 -2.26660000000000004e1 0 0 +4.61139999999999986e2 -2.26660000000000004e1 0 0 +4.61139999999999986e2 -2.98339999999999996e1 0 0 +4.64139999999999986e2 -2.98339999999999996e1 0 0 +4.64139999999999986e2 -2.78339999999999996e1 0 0 +4.61139999999999986e2 -2.78339999999999996e1 0 0 +[NumSegments] = 96 +0 1 -1 1 0 0 +1 2 -1 1 0 0 +2 3 -1 1 0 0 +3 0 -1 1 0 0 +4 5 -1 0 0 0 +5 19 -1 0 0 0 +6 4 -1 0 0 0 +19 18 -1 0 0 0 +18 17 -1 0 0 0 +17 16 -1 0 0 0 +16 15 -1 0 0 0 +15 14 -1 0 0 0 +14 13 -1 0 0 0 +13 12 -1 0 0 0 +12 11 -1 0 0 0 +11 10 -1 0 0 0 +10 9 -1 0 0 0 +9 8 -1 0 0 0 +8 7 -1 0 0 0 +7 6 -1 0 0 0 +20 34 -1 0 0 0 +34 35 -1 0 0 0 +35 33 -1 0 0 0 +33 32 -1 0 0 0 +32 31 -1 0 0 0 +31 30 -1 0 0 0 +30 29 -1 0 0 0 +29 28 -1 0 0 0 +28 27 -1 0 0 0 +27 26 -1 0 0 0 +26 25 -1 0 0 0 +25 24 -1 0 0 0 +24 23 -1 0 0 0 +23 22 -1 0 0 0 +22 21 -1 0 0 0 +21 20 -1 0 0 0 +8 22 -1 0 0 0 +7 21 -1 0 0 0 +10 24 -1 0 0 0 +9 23 -1 0 0 0 +12 26 -1 0 0 0 +11 25 -1 0 0 0 +14 28 -1 0 0 0 +13 27 -1 0 0 0 +16 30 -1 0 0 0 +15 29 -1 0 0 0 +18 32 -1 0 0 0 +17 31 -1 0 0 0 +36 37 -1 0 0 0 +37 38 -1 0 0 0 +38 39 -1 0 0 0 +39 36 -1 0 0 0 +40 41 -1 0 0 0 +41 42 -1 0 0 0 +42 43 -1 0 0 0 +43 40 -1 0 0 0 +44 45 -1 0 0 0 +45 46 -1 0 0 0 +46 47 -1 0 0 0 +47 44 -1 0 0 0 +48 49 -1 0 0 0 +49 50 -1 0 0 0 +50 51 -1 0 0 0 +51 48 -1 0 0 0 +52 53 -1 0 0 0 +53 54 -1 0 0 0 +54 55 -1 0 0 0 +55 52 -1 0 0 0 +56 57 -1 0 0 0 +57 58 -1 0 0 0 +58 59 -1 0 0 0 +59 56 -1 0 0 0 +60 61 -1 0 0 0 +61 62 -1 0 0 0 +62 63 -1 0 0 0 +63 60 -1 0 0 0 +64 65 -1 0 0 0 +65 66 -1 0 0 0 +66 67 -1 0 0 0 +67 64 -1 0 0 0 +68 69 -1 0 0 0 +69 70 -1 0 0 0 +70 71 -1 0 0 0 +71 68 -1 0 0 0 +72 73 -1 0 0 0 +73 74 -1 0 0 0 +74 75 -1 0 0 0 +75 72 -1 0 0 0 +76 77 -1 0 0 0 +77 78 -1 0 0 0 +78 79 -1 0 0 0 +79 76 -1 0 0 0 +80 81 -1 0 0 0 +81 82 -1 0 0 0 +82 83 -1 0 0 0 +83 80 -1 0 0 0 +[NumArcSegments] = 0 +[NumHoles] = 0 +[NumBlockLabels] = 21 +-5.00000000000000000e1 -5.00000000000000000e1 1 -1 0 0.00000000000000000e0 0 1 2 +-1.50000000000000000e1 0.00000000000000000e0 2 5.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +6.62700000000000045e2 0.00000000000000000e0 3 5.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +3.20000000000000000e2 -2.62500000000000000e1 3 1.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +3.20000000000000000e2 -1.57500000000000000e1 3 1.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +3.20000000000000000e2 -5.25000000000000000e0 3 1.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +3.20000000000000000e2 5.25000000000000000e0 3 1.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +3.20000000000000000e2 1.57500000000000000e1 3 1.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +3.20000000000000000e2 2.62500000000000000e1 3 1.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +4.62639999999999986e2 2.83769999999999989e1 2 5.00000000000000000e-1 0 2.70000000000000000e2 0 1 0 +4.62639999999999986e2 2.41230000000000011e1 2 5.00000000000000000e-1 0 2.70000000000000000e2 0 1 0 +4.62639999999999986e2 1.79149999999999991e1 2 5.00000000000000000e-1 0 2.70000000000000000e2 0 1 0 +4.62639999999999986e2 1.35850000000000009e1 2 5.00000000000000000e-1 0 2.70000000000000000e2 0 1 0 +4.62639999999999986e2 7.46600000000000019e0 2 5.00000000000000000e-1 0 2.70000000000000000e2 0 1 0 +4.62639999999999986e2 3.03399999999999981e0 2 5.00000000000000000e-1 0 2.70000000000000000e2 0 1 0 +4.62639999999999986e2 -2.91999999999999993e0 2 5.00000000000000000e-1 0 2.70000000000000000e2 0 1 0 +4.62639999999999986e2 -7.58000000000000007e0 2 5.00000000000000000e-1 0 2.70000000000000000e2 0 1 0 +4.62639999999999986e2 -1.32929999999999993e1 2 5.00000000000000000e-1 0 2.70000000000000000e2 0 1 0 +4.62639999999999986e2 -1.82070000000000007e1 2 5.00000000000000000e-1 0 2.70000000000000000e2 0 1 0 +4.62639999999999986e2 -2.36660000000000004e1 2 5.00000000000000000e-1 0 2.70000000000000000e2 0 1 0 +4.62639999999999986e2 -2.88339999999999996e1 2 5.00000000000000000e-1 0 2.70000000000000000e2 0 1 0 +[/pristine-doc] +[base-doc] +[Format] = 4.00000000000000000e0 +[Frequency] = 0.00000000000000000e0 +[Precision] = 1.00000000000000002e-8 +[MinAngle] = 2.50000000000000000e1 +[DoSmartMesh] = 1 +[Depth] = 5.00000000000000000e0 +[LengthUnits] = millimeters +[ProblemType] = planar +[Coordinates] = cartesian +[ACSolver] = 0 +[PrevType] = 0 +[PrevSoln] = "" +[Comment] = "six plain steel strings (.010-.046 inch gauges) over a 25.5-inch scale, anchored to a NdFeB bar magnet at the left (mag_dir 0, N facing the strings) and a passive steel block at the right. small NdFeB disc magnets sit above and below each string at 5/7 of the way toward the steel block, all magnetised in -y so the disc above each string has its N face downward and the disc below has its S face upward, driving a strong vertical bias field through each string at the probe point. flip every disc mag_dir from 270 to 90 to invert the bias direction." +[PointProps] = 0 +[BdryProps] = 1 + + = "A=0" + = 0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + +[BlockProps] = 3 + + = "Air" + = 1.00000000000000000e0 + = 1.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0 + = 1.00000000000000000e0 + = 0 + = 0.00000000000000000e0 + = 0 + + + = "NdFeB" + = 1.05000000000000004e0 + = 1.05000000000000004e0 + = 9.15000000000000000e5 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 6.67000000000000037e-1 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0 + = 1.00000000000000000e0 + = 0 + = 0.00000000000000000e0 + = 0 + + + = "Steel" + = 2.50000000000000000e3 + = 2.50000000000000000e3 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 5.79999999999999982e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0 + = 1.00000000000000000e0 + = 0 + = 0.00000000000000000e0 + = 0 + +[CircuitProps] = 0 +[NumPoints] = 312 +-1.00000000000000000e2 -1.00000000000000000e2 0 0 +7.50000000000000000e2 -1.00000000000000000e2 0 0 +7.50000000000000000e2 1.00000000000000000e2 0 0 +-1.00000000000000000e2 1.00000000000000000e2 0 0 +-3.00000000000000000e1 -3.00000000000000000e1 0 0 +-3.00000000000000000e1 3.00000000000000000e1 0 0 +0.00000000000000000e0 -3.00000000000000000e1 0 0 +0.00000000000000000e0 -2.68339999999999996e1 0 0 +0.00000000000000000e0 -2.56660000000000004e1 0 0 +0.00000000000000000e0 -1.62070000000000007e1 0 0 +0.00000000000000000e0 -1.52929999999999993e1 0 0 +0.00000000000000000e0 -5.58000000000000007e0 0 0 +0.00000000000000000e0 -4.91999999999999993e0 0 0 +0.00000000000000000e0 5.03399999999999981e0 0 0 +0.00000000000000000e0 5.46600000000000019e0 0 0 +0.00000000000000000e0 1.55850000000000009e1 0 0 +0.00000000000000000e0 1.59149999999999991e1 0 0 +0.00000000000000000e0 2.61230000000000011e1 0 0 +0.00000000000000000e0 2.63769999999999989e1 0 0 +0.00000000000000000e0 3.00000000000000000e1 0 0 +6.47700000000000045e2 -3.00000000000000000e1 0 0 +6.47700000000000045e2 -2.68339999999999996e1 0 0 +6.47700000000000045e2 -2.56660000000000004e1 0 0 +6.47700000000000045e2 -1.62070000000000007e1 0 0 +6.47700000000000045e2 -1.52929999999999993e1 0 0 +6.47700000000000045e2 -5.58000000000000007e0 0 0 +6.47700000000000045e2 -4.91999999999999993e0 0 0 +6.47700000000000045e2 5.03399999999999981e0 0 0 +6.47700000000000045e2 5.46600000000000019e0 0 0 +6.47700000000000045e2 1.55850000000000009e1 0 0 +6.47700000000000045e2 1.59149999999999991e1 0 0 +6.47700000000000045e2 2.61230000000000011e1 0 0 +6.47700000000000045e2 2.63769999999999989e1 0 0 +6.47700000000000045e2 3.00000000000000000e1 0 0 +6.77700000000000045e2 -3.00000000000000000e1 0 0 +6.77700000000000045e2 3.00000000000000000e1 0 0 +4.61139999999999986e2 2.73769999999999989e1 0 0 +4.64139999999999986e2 2.73769999999999989e1 0 0 +4.64139999999999986e2 2.93769999999999989e1 0 0 +4.61139999999999986e2 2.93769999999999989e1 0 0 +4.61139999999999986e2 2.31230000000000011e1 0 0 +4.64139999999999986e2 2.31230000000000011e1 0 0 +4.64139999999999986e2 2.51230000000000011e1 0 0 +4.61139999999999986e2 2.51230000000000011e1 0 0 +4.61139999999999986e2 1.69149999999999991e1 0 0 +4.64139999999999986e2 1.69149999999999991e1 0 0 +4.64139999999999986e2 1.89149999999999991e1 0 0 +4.61139999999999986e2 1.89149999999999991e1 0 0 +4.61139999999999986e2 1.25850000000000009e1 0 0 +4.64139999999999986e2 1.25850000000000009e1 0 0 +4.64139999999999986e2 1.45850000000000009e1 0 0 +4.61139999999999986e2 1.45850000000000009e1 0 0 +4.61139999999999986e2 6.46600000000000019e0 0 0 +4.64139999999999986e2 6.46600000000000019e0 0 0 +4.64139999999999986e2 8.46599999999999930e0 0 0 +4.61139999999999986e2 8.46599999999999930e0 0 0 +4.61139999999999986e2 2.03399999999999981e0 0 0 +4.64139999999999986e2 2.03399999999999981e0 0 0 +4.64139999999999986e2 4.03399999999999981e0 0 0 +4.61139999999999986e2 4.03399999999999981e0 0 0 +4.61139999999999986e2 -3.91999999999999993e0 0 0 +4.64139999999999986e2 -3.91999999999999993e0 0 0 +4.64139999999999986e2 -1.91999999999999993e0 0 0 +4.61139999999999986e2 -1.91999999999999993e0 0 0 +4.61139999999999986e2 -8.58000000000000007e0 0 0 +4.64139999999999986e2 -8.58000000000000007e0 0 0 +4.64139999999999986e2 -6.58000000000000007e0 0 0 +4.61139999999999986e2 -6.58000000000000007e0 0 0 +4.61139999999999986e2 -1.42929999999999993e1 0 0 +4.64139999999999986e2 -1.42929999999999993e1 0 0 +4.64139999999999986e2 -1.22929999999999993e1 0 0 +4.61139999999999986e2 -1.22929999999999993e1 0 0 +4.61139999999999986e2 -1.92070000000000007e1 0 0 +4.64139999999999986e2 -1.92070000000000007e1 0 0 +4.64139999999999986e2 -1.72070000000000007e1 0 0 +4.61139999999999986e2 -1.72070000000000007e1 0 0 +4.61139999999999986e2 -2.46660000000000004e1 0 0 +4.64139999999999986e2 -2.46660000000000004e1 0 0 +4.64139999999999986e2 -2.26660000000000004e1 0 0 +4.61139999999999986e2 -2.26660000000000004e1 0 0 +4.61139999999999986e2 -2.98339999999999996e1 0 0 +4.64139999999999986e2 -2.98339999999999996e1 0 0 +4.64139999999999986e2 -2.78339999999999996e1 0 0 +4.61139999999999986e2 -2.78339999999999996e1 0 0 +3.23850000000000051e1 -2.68339999999999996e1 0 0 +6.47700000000000102e1 -2.68339999999999996e1 0 0 +9.71550000000000153e1 -2.68339999999999996e1 0 0 +1.29540000000000020e2 -2.68339999999999996e1 0 0 +1.61925000000000011e2 -2.68339999999999996e1 0 0 +1.94310000000000031e2 -2.68339999999999996e1 0 0 +2.26695000000000050e2 -2.68339999999999996e1 0 0 +2.59080000000000041e2 -2.68339999999999996e1 0 0 +2.91465000000000032e2 -2.68339999999999996e1 0 0 +3.23850000000000023e2 -2.68339999999999996e1 0 0 +3.56235000000000070e2 -2.68339999999999996e1 0 0 +3.88620000000000061e2 -2.68339999999999996e1 0 0 +4.21005000000000052e2 -2.68339999999999996e1 0 0 +4.53390000000000100e2 -2.68339999999999996e1 0 0 +4.85775000000000091e2 -2.68339999999999996e1 0 0 +5.18160000000000082e2 -2.68339999999999996e1 0 0 +5.50545000000000073e2 -2.68339999999999996e1 0 0 +5.82930000000000064e2 -2.68339999999999996e1 0 0 +6.15315000000000055e2 -2.68339999999999996e1 0 0 +3.23850000000000051e1 -2.56660000000000004e1 0 0 +6.47700000000000102e1 -2.56660000000000004e1 0 0 +9.71550000000000153e1 -2.56660000000000004e1 0 0 +1.29540000000000020e2 -2.56660000000000004e1 0 0 +1.61925000000000011e2 -2.56660000000000004e1 0 0 +1.94310000000000031e2 -2.56660000000000004e1 0 0 +2.26695000000000050e2 -2.56660000000000004e1 0 0 +2.59080000000000041e2 -2.56660000000000004e1 0 0 +2.91465000000000032e2 -2.56660000000000004e1 0 0 +3.23850000000000023e2 -2.56660000000000004e1 0 0 +3.56235000000000070e2 -2.56660000000000004e1 0 0 +3.88620000000000061e2 -2.56660000000000004e1 0 0 +4.21005000000000052e2 -2.56660000000000004e1 0 0 +4.53390000000000100e2 -2.56660000000000004e1 0 0 +4.85775000000000091e2 -2.56660000000000004e1 0 0 +5.18160000000000082e2 -2.56660000000000004e1 0 0 +5.50545000000000073e2 -2.56660000000000004e1 0 0 +5.82930000000000064e2 -2.56660000000000004e1 0 0 +6.15315000000000055e2 -2.56660000000000004e1 0 0 +3.23850000000000051e1 -1.62070000000000007e1 0 0 +6.47700000000000102e1 -1.62070000000000007e1 0 0 +9.71550000000000153e1 -1.62070000000000007e1 0 0 +1.29540000000000020e2 -1.62070000000000007e1 0 0 +1.61925000000000011e2 -1.62070000000000007e1 0 0 +1.94310000000000031e2 -1.62070000000000007e1 0 0 +2.26695000000000050e2 -1.62070000000000007e1 0 0 +2.59080000000000041e2 -1.62070000000000007e1 0 0 +2.91465000000000032e2 -1.62070000000000007e1 0 0 +3.23850000000000023e2 -1.62070000000000007e1 0 0 +3.56235000000000070e2 -1.62070000000000007e1 0 0 +3.88620000000000061e2 -1.62070000000000007e1 0 0 +4.21005000000000052e2 -1.62070000000000007e1 0 0 +4.53390000000000100e2 -1.62070000000000007e1 0 0 +4.85775000000000091e2 -1.62070000000000007e1 0 0 +5.18160000000000082e2 -1.62070000000000007e1 0 0 +5.50545000000000073e2 -1.62070000000000007e1 0 0 +5.82930000000000064e2 -1.62070000000000007e1 0 0 +6.15315000000000055e2 -1.62070000000000007e1 0 0 +3.23850000000000051e1 -1.52929999999999993e1 0 0 +6.47700000000000102e1 -1.52929999999999993e1 0 0 +9.71550000000000153e1 -1.52929999999999993e1 0 0 +1.29540000000000020e2 -1.52929999999999993e1 0 0 +1.61925000000000011e2 -1.52929999999999993e1 0 0 +1.94310000000000031e2 -1.52929999999999993e1 0 0 +2.26695000000000050e2 -1.52929999999999993e1 0 0 +2.59080000000000041e2 -1.52929999999999993e1 0 0 +2.91465000000000032e2 -1.52929999999999993e1 0 0 +3.23850000000000023e2 -1.52929999999999993e1 0 0 +3.56235000000000070e2 -1.52929999999999993e1 0 0 +3.88620000000000061e2 -1.52929999999999993e1 0 0 +4.21005000000000052e2 -1.52929999999999993e1 0 0 +4.53390000000000100e2 -1.52929999999999993e1 0 0 +4.85775000000000091e2 -1.52929999999999993e1 0 0 +5.18160000000000082e2 -1.52929999999999993e1 0 0 +5.50545000000000073e2 -1.52929999999999993e1 0 0 +5.82930000000000064e2 -1.52929999999999993e1 0 0 +6.15315000000000055e2 -1.52929999999999993e1 0 0 +3.23850000000000051e1 -5.58000000000000007e0 0 0 +6.47700000000000102e1 -5.58000000000000007e0 0 0 +9.71550000000000153e1 -5.58000000000000007e0 0 0 +1.29540000000000020e2 -5.58000000000000007e0 0 0 +1.61925000000000011e2 -5.58000000000000007e0 0 0 +1.94310000000000031e2 -5.58000000000000007e0 0 0 +2.26695000000000050e2 -5.58000000000000007e0 0 0 +2.59080000000000041e2 -5.58000000000000007e0 0 0 +2.91465000000000032e2 -5.58000000000000007e0 0 0 +3.23850000000000023e2 -5.58000000000000007e0 0 0 +3.56235000000000070e2 -5.58000000000000007e0 0 0 +3.88620000000000061e2 -5.58000000000000007e0 0 0 +4.21005000000000052e2 -5.58000000000000007e0 0 0 +4.53390000000000100e2 -5.58000000000000007e0 0 0 +4.85775000000000091e2 -5.58000000000000007e0 0 0 +5.18160000000000082e2 -5.58000000000000007e0 0 0 +5.50545000000000073e2 -5.58000000000000007e0 0 0 +5.82930000000000064e2 -5.58000000000000007e0 0 0 +6.15315000000000055e2 -5.58000000000000007e0 0 0 +3.23850000000000051e1 -4.91999999999999993e0 0 0 +6.47700000000000102e1 -4.91999999999999993e0 0 0 +9.71550000000000153e1 -4.91999999999999993e0 0 0 +1.29540000000000020e2 -4.91999999999999993e0 0 0 +1.61925000000000011e2 -4.91999999999999993e0 0 0 +1.94310000000000031e2 -4.91999999999999993e0 0 0 +2.26695000000000050e2 -4.91999999999999993e0 0 0 +2.59080000000000041e2 -4.91999999999999993e0 0 0 +2.91465000000000032e2 -4.91999999999999993e0 0 0 +3.23850000000000023e2 -4.91999999999999993e0 0 0 +3.56235000000000070e2 -4.91999999999999993e0 0 0 +3.88620000000000061e2 -4.91999999999999993e0 0 0 +4.21005000000000052e2 -4.91999999999999993e0 0 0 +4.53390000000000100e2 -4.91999999999999993e0 0 0 +4.85775000000000091e2 -4.91999999999999993e0 0 0 +5.18160000000000082e2 -4.91999999999999993e0 0 0 +5.50545000000000073e2 -4.91999999999999993e0 0 0 +5.82930000000000064e2 -4.91999999999999993e0 0 0 +6.15315000000000055e2 -4.91999999999999993e0 0 0 +3.23850000000000051e1 5.03399999999999981e0 0 0 +6.47700000000000102e1 5.03399999999999981e0 0 0 +9.71550000000000153e1 5.03399999999999981e0 0 0 +1.29540000000000020e2 5.03399999999999981e0 0 0 +1.61925000000000011e2 5.03399999999999981e0 0 0 +1.94310000000000031e2 5.03399999999999981e0 0 0 +2.26695000000000050e2 5.03399999999999981e0 0 0 +2.59080000000000041e2 5.03399999999999981e0 0 0 +2.91465000000000032e2 5.03399999999999981e0 0 0 +3.23850000000000023e2 5.03399999999999981e0 0 0 +3.56235000000000070e2 5.03399999999999981e0 0 0 +3.88620000000000061e2 5.03399999999999981e0 0 0 +4.21005000000000052e2 5.03399999999999981e0 0 0 +4.53390000000000100e2 5.03399999999999981e0 0 0 +4.85775000000000091e2 5.03399999999999981e0 0 0 +5.18160000000000082e2 5.03399999999999981e0 0 0 +5.50545000000000073e2 5.03399999999999981e0 0 0 +5.82930000000000064e2 5.03399999999999981e0 0 0 +6.15315000000000055e2 5.03399999999999981e0 0 0 +3.23850000000000051e1 5.46600000000000019e0 0 0 +6.47700000000000102e1 5.46600000000000019e0 0 0 +9.71550000000000153e1 5.46600000000000019e0 0 0 +1.29540000000000020e2 5.46600000000000019e0 0 0 +1.61925000000000011e2 5.46600000000000019e0 0 0 +1.94310000000000031e2 5.46600000000000019e0 0 0 +2.26695000000000050e2 5.46600000000000019e0 0 0 +2.59080000000000041e2 5.46600000000000019e0 0 0 +2.91465000000000032e2 5.46600000000000019e0 0 0 +3.23850000000000023e2 5.46600000000000019e0 0 0 +3.56235000000000070e2 5.46600000000000019e0 0 0 +3.88620000000000061e2 5.46600000000000019e0 0 0 +4.21005000000000052e2 5.46600000000000019e0 0 0 +4.53390000000000100e2 5.46600000000000019e0 0 0 +4.85775000000000091e2 5.46600000000000019e0 0 0 +5.18160000000000082e2 5.46600000000000019e0 0 0 +5.50545000000000073e2 5.46600000000000019e0 0 0 +5.82930000000000064e2 5.46600000000000019e0 0 0 +6.15315000000000055e2 5.46600000000000019e0 0 0 +3.23850000000000051e1 1.55850000000000009e1 0 0 +6.47700000000000102e1 1.55850000000000009e1 0 0 +9.71550000000000153e1 1.55850000000000009e1 0 0 +1.29540000000000020e2 1.55850000000000009e1 0 0 +1.61925000000000011e2 1.55850000000000009e1 0 0 +1.94310000000000031e2 1.55850000000000009e1 0 0 +2.26695000000000050e2 1.55850000000000009e1 0 0 +2.59080000000000041e2 1.55850000000000009e1 0 0 +2.91465000000000032e2 1.55850000000000009e1 0 0 +3.23850000000000023e2 1.55850000000000009e1 0 0 +3.56235000000000070e2 1.55850000000000009e1 0 0 +3.88620000000000061e2 1.55850000000000009e1 0 0 +4.21005000000000052e2 1.55850000000000009e1 0 0 +4.53390000000000100e2 1.55850000000000009e1 0 0 +4.85775000000000091e2 1.55850000000000009e1 0 0 +5.18160000000000082e2 1.55850000000000009e1 0 0 +5.50545000000000073e2 1.55850000000000009e1 0 0 +5.82930000000000064e2 1.55850000000000009e1 0 0 +6.15315000000000055e2 1.55850000000000009e1 0 0 +3.23850000000000051e1 1.59149999999999991e1 0 0 +6.47700000000000102e1 1.59149999999999991e1 0 0 +9.71550000000000153e1 1.59149999999999991e1 0 0 +1.29540000000000020e2 1.59149999999999991e1 0 0 +1.61925000000000011e2 1.59149999999999991e1 0 0 +1.94310000000000031e2 1.59149999999999991e1 0 0 +2.26695000000000050e2 1.59149999999999991e1 0 0 +2.59080000000000041e2 1.59149999999999991e1 0 0 +2.91465000000000032e2 1.59149999999999991e1 0 0 +3.23850000000000023e2 1.59149999999999991e1 0 0 +3.56235000000000070e2 1.59149999999999991e1 0 0 +3.88620000000000061e2 1.59149999999999991e1 0 0 +4.21005000000000052e2 1.59149999999999991e1 0 0 +4.53390000000000100e2 1.59149999999999991e1 0 0 +4.85775000000000091e2 1.59149999999999991e1 0 0 +5.18160000000000082e2 1.59149999999999991e1 0 0 +5.50545000000000073e2 1.59149999999999991e1 0 0 +5.82930000000000064e2 1.59149999999999991e1 0 0 +6.15315000000000055e2 1.59149999999999991e1 0 0 +3.23850000000000051e1 2.61230000000000011e1 0 0 +6.47700000000000102e1 2.61230000000000011e1 0 0 +9.71550000000000153e1 2.61230000000000011e1 0 0 +1.29540000000000020e2 2.61230000000000011e1 0 0 +1.61925000000000011e2 2.61230000000000011e1 0 0 +1.94310000000000031e2 2.61230000000000011e1 0 0 +2.26695000000000050e2 2.61230000000000011e1 0 0 +2.59080000000000041e2 2.61230000000000011e1 0 0 +2.91465000000000032e2 2.61230000000000011e1 0 0 +3.23850000000000023e2 2.61230000000000011e1 0 0 +3.56235000000000070e2 2.61230000000000011e1 0 0 +3.88620000000000061e2 2.61230000000000011e1 0 0 +4.21005000000000052e2 2.61230000000000011e1 0 0 +4.53390000000000100e2 2.61230000000000011e1 0 0 +4.85775000000000091e2 2.61230000000000011e1 0 0 +5.18160000000000082e2 2.61230000000000011e1 0 0 +5.50545000000000073e2 2.61230000000000011e1 0 0 +5.82930000000000064e2 2.61230000000000011e1 0 0 +6.15315000000000055e2 2.61230000000000011e1 0 0 +3.23850000000000051e1 2.63769999999999989e1 0 0 +6.47700000000000102e1 2.63769999999999989e1 0 0 +9.71550000000000153e1 2.63769999999999989e1 0 0 +1.29540000000000020e2 2.63769999999999989e1 0 0 +1.61925000000000011e2 2.63769999999999989e1 0 0 +1.94310000000000031e2 2.63769999999999989e1 0 0 +2.26695000000000050e2 2.63769999999999989e1 0 0 +2.59080000000000041e2 2.63769999999999989e1 0 0 +2.91465000000000032e2 2.63769999999999989e1 0 0 +3.23850000000000023e2 2.63769999999999989e1 0 0 +3.56235000000000070e2 2.63769999999999989e1 0 0 +3.88620000000000061e2 2.63769999999999989e1 0 0 +4.21005000000000052e2 2.63769999999999989e1 0 0 +4.53390000000000100e2 2.63769999999999989e1 0 0 +4.85775000000000091e2 2.63769999999999989e1 0 0 +5.18160000000000082e2 2.63769999999999989e1 0 0 +5.50545000000000073e2 2.63769999999999989e1 0 0 +5.82930000000000064e2 2.63769999999999989e1 0 0 +6.15315000000000055e2 2.63769999999999989e1 0 0 +[NumSegments] = 324 +0 1 -1 1 0 0 +1 2 -1 1 0 0 +2 3 -1 1 0 0 +3 0 -1 1 0 0 +4 5 -1 0 0 0 +5 19 -1 0 0 0 +6 4 -1 0 0 0 +19 18 -1 0 0 0 +18 17 -1 0 0 0 +17 16 -1 0 0 0 +16 15 -1 0 0 0 +15 14 -1 0 0 0 +14 13 -1 0 0 0 +13 12 -1 0 0 0 +12 11 -1 0 0 0 +11 10 -1 0 0 0 +10 9 -1 0 0 0 +9 8 -1 0 0 0 +8 7 -1 0 0 0 +7 6 -1 0 0 0 +20 34 -1 0 0 0 +34 35 -1 0 0 0 +35 33 -1 0 0 0 +33 32 -1 0 0 0 +32 31 -1 0 0 0 +31 30 -1 0 0 0 +30 29 -1 0 0 0 +29 28 -1 0 0 0 +28 27 -1 0 0 0 +27 26 -1 0 0 0 +26 25 -1 0 0 0 +25 24 -1 0 0 0 +24 23 -1 0 0 0 +23 22 -1 0 0 0 +22 21 -1 0 0 0 +21 20 -1 0 0 0 +36 37 -1 0 0 0 +37 38 -1 0 0 0 +38 39 -1 0 0 0 +39 36 -1 0 0 0 +40 41 -1 0 0 0 +41 42 -1 0 0 0 +42 43 -1 0 0 0 +43 40 -1 0 0 0 +44 45 -1 0 0 0 +45 46 -1 0 0 0 +46 47 -1 0 0 0 +47 44 -1 0 0 0 +48 49 -1 0 0 0 +49 50 -1 0 0 0 +50 51 -1 0 0 0 +51 48 -1 0 0 0 +52 53 -1 0 0 0 +53 54 -1 0 0 0 +54 55 -1 0 0 0 +55 52 -1 0 0 0 +56 57 -1 0 0 0 +57 58 -1 0 0 0 +58 59 -1 0 0 0 +59 56 -1 0 0 0 +60 61 -1 0 0 0 +61 62 -1 0 0 0 +62 63 -1 0 0 0 +63 60 -1 0 0 0 +64 65 -1 0 0 0 +65 66 -1 0 0 0 +66 67 -1 0 0 0 +67 64 -1 0 0 0 +68 69 -1 0 0 0 +69 70 -1 0 0 0 +70 71 -1 0 0 0 +71 68 -1 0 0 0 +72 73 -1 0 0 0 +73 74 -1 0 0 0 +74 75 -1 0 0 0 +75 72 -1 0 0 0 +76 77 -1 0 0 0 +77 78 -1 0 0 0 +78 79 -1 0 0 0 +79 76 -1 0 0 0 +80 81 -1 0 0 0 +81 82 -1 0 0 0 +82 83 -1 0 0 0 +83 80 -1 0 0 0 +7 84 -1 0 0 0 +84 85 -1 0 0 0 +85 86 -1 0 0 0 +86 87 -1 0 0 0 +87 88 -1 0 0 0 +88 89 -1 0 0 0 +89 90 -1 0 0 0 +90 91 -1 0 0 0 +91 92 -1 0 0 0 +92 93 -1 0 0 0 +93 94 -1 0 0 0 +94 95 -1 0 0 0 +95 96 -1 0 0 0 +96 97 -1 0 0 0 +97 98 -1 0 0 0 +98 99 -1 0 0 0 +99 100 -1 0 0 0 +100 101 -1 0 0 0 +101 102 -1 0 0 0 +102 21 -1 0 0 0 +8 103 -1 0 0 0 +103 104 -1 0 0 0 +104 105 -1 0 0 0 +105 106 -1 0 0 0 +106 107 -1 0 0 0 +107 108 -1 0 0 0 +108 109 -1 0 0 0 +109 110 -1 0 0 0 +110 111 -1 0 0 0 +111 112 -1 0 0 0 +112 113 -1 0 0 0 +113 114 -1 0 0 0 +114 115 -1 0 0 0 +115 116 -1 0 0 0 +116 117 -1 0 0 0 +117 118 -1 0 0 0 +118 119 -1 0 0 0 +119 120 -1 0 0 0 +120 121 -1 0 0 0 +121 22 -1 0 0 0 +9 122 -1 0 0 0 +122 123 -1 0 0 0 +123 124 -1 0 0 0 +124 125 -1 0 0 0 +125 126 -1 0 0 0 +126 127 -1 0 0 0 +127 128 -1 0 0 0 +128 129 -1 0 0 0 +129 130 -1 0 0 0 +130 131 -1 0 0 0 +131 132 -1 0 0 0 +132 133 -1 0 0 0 +133 134 -1 0 0 0 +134 135 -1 0 0 0 +135 136 -1 0 0 0 +136 137 -1 0 0 0 +137 138 -1 0 0 0 +138 139 -1 0 0 0 +139 140 -1 0 0 0 +140 23 -1 0 0 0 +10 141 -1 0 0 0 +141 142 -1 0 0 0 +142 143 -1 0 0 0 +143 144 -1 0 0 0 +144 145 -1 0 0 0 +145 146 -1 0 0 0 +146 147 -1 0 0 0 +147 148 -1 0 0 0 +148 149 -1 0 0 0 +149 150 -1 0 0 0 +150 151 -1 0 0 0 +151 152 -1 0 0 0 +152 153 -1 0 0 0 +153 154 -1 0 0 0 +154 155 -1 0 0 0 +155 156 -1 0 0 0 +156 157 -1 0 0 0 +157 158 -1 0 0 0 +158 159 -1 0 0 0 +159 24 -1 0 0 0 +11 160 -1 0 0 0 +160 161 -1 0 0 0 +161 162 -1 0 0 0 +162 163 -1 0 0 0 +163 164 -1 0 0 0 +164 165 -1 0 0 0 +165 166 -1 0 0 0 +166 167 -1 0 0 0 +167 168 -1 0 0 0 +168 169 -1 0 0 0 +169 170 -1 0 0 0 +170 171 -1 0 0 0 +171 172 -1 0 0 0 +172 173 -1 0 0 0 +173 174 -1 0 0 0 +174 175 -1 0 0 0 +175 176 -1 0 0 0 +176 177 -1 0 0 0 +177 178 -1 0 0 0 +178 25 -1 0 0 0 +12 179 -1 0 0 0 +179 180 -1 0 0 0 +180 181 -1 0 0 0 +181 182 -1 0 0 0 +182 183 -1 0 0 0 +183 184 -1 0 0 0 +184 185 -1 0 0 0 +185 186 -1 0 0 0 +186 187 -1 0 0 0 +187 188 -1 0 0 0 +188 189 -1 0 0 0 +189 190 -1 0 0 0 +190 191 -1 0 0 0 +191 192 -1 0 0 0 +192 193 -1 0 0 0 +193 194 -1 0 0 0 +194 195 -1 0 0 0 +195 196 -1 0 0 0 +196 197 -1 0 0 0 +197 26 -1 0 0 0 +13 198 -1 0 0 0 +198 199 -1 0 0 0 +199 200 -1 0 0 0 +200 201 -1 0 0 0 +201 202 -1 0 0 0 +202 203 -1 0 0 0 +203 204 -1 0 0 0 +204 205 -1 0 0 0 +205 206 -1 0 0 0 +206 207 -1 0 0 0 +207 208 -1 0 0 0 +208 209 -1 0 0 0 +209 210 -1 0 0 0 +210 211 -1 0 0 0 +211 212 -1 0 0 0 +212 213 -1 0 0 0 +213 214 -1 0 0 0 +214 215 -1 0 0 0 +215 216 -1 0 0 0 +216 27 -1 0 0 0 +14 217 -1 0 0 0 +217 218 -1 0 0 0 +218 219 -1 0 0 0 +219 220 -1 0 0 0 +220 221 -1 0 0 0 +221 222 -1 0 0 0 +222 223 -1 0 0 0 +223 224 -1 0 0 0 +224 225 -1 0 0 0 +225 226 -1 0 0 0 +226 227 -1 0 0 0 +227 228 -1 0 0 0 +228 229 -1 0 0 0 +229 230 -1 0 0 0 +230 231 -1 0 0 0 +231 232 -1 0 0 0 +232 233 -1 0 0 0 +233 234 -1 0 0 0 +234 235 -1 0 0 0 +235 28 -1 0 0 0 +15 236 -1 0 0 0 +236 237 -1 0 0 0 +237 238 -1 0 0 0 +238 239 -1 0 0 0 +239 240 -1 0 0 0 +240 241 -1 0 0 0 +241 242 -1 0 0 0 +242 243 -1 0 0 0 +243 244 -1 0 0 0 +244 245 -1 0 0 0 +245 246 -1 0 0 0 +246 247 -1 0 0 0 +247 248 -1 0 0 0 +248 249 -1 0 0 0 +249 250 -1 0 0 0 +250 251 -1 0 0 0 +251 252 -1 0 0 0 +252 253 -1 0 0 0 +253 254 -1 0 0 0 +254 29 -1 0 0 0 +16 255 -1 0 0 0 +255 256 -1 0 0 0 +256 257 -1 0 0 0 +257 258 -1 0 0 0 +258 259 -1 0 0 0 +259 260 -1 0 0 0 +260 261 -1 0 0 0 +261 262 -1 0 0 0 +262 263 -1 0 0 0 +263 264 -1 0 0 0 +264 265 -1 0 0 0 +265 266 -1 0 0 0 +266 267 -1 0 0 0 +267 268 -1 0 0 0 +268 269 -1 0 0 0 +269 270 -1 0 0 0 +270 271 -1 0 0 0 +271 272 -1 0 0 0 +272 273 -1 0 0 0 +273 30 -1 0 0 0 +17 274 -1 0 0 0 +274 275 -1 0 0 0 +275 276 -1 0 0 0 +276 277 -1 0 0 0 +277 278 -1 0 0 0 +278 279 -1 0 0 0 +279 280 -1 0 0 0 +280 281 -1 0 0 0 +281 282 -1 0 0 0 +282 283 -1 0 0 0 +283 284 -1 0 0 0 +284 285 -1 0 0 0 +285 286 -1 0 0 0 +286 287 -1 0 0 0 +287 288 -1 0 0 0 +288 289 -1 0 0 0 +289 290 -1 0 0 0 +290 291 -1 0 0 0 +291 292 -1 0 0 0 +292 31 -1 0 0 0 +18 293 -1 0 0 0 +293 294 -1 0 0 0 +294 295 -1 0 0 0 +295 296 -1 0 0 0 +296 297 -1 0 0 0 +297 298 -1 0 0 0 +298 299 -1 0 0 0 +299 300 -1 0 0 0 +300 301 -1 0 0 0 +301 302 -1 0 0 0 +302 303 -1 0 0 0 +303 304 -1 0 0 0 +304 305 -1 0 0 0 +305 306 -1 0 0 0 +306 307 -1 0 0 0 +307 308 -1 0 0 0 +308 309 -1 0 0 0 +309 310 -1 0 0 0 +310 311 -1 0 0 0 +311 32 -1 0 0 0 +[NumArcSegments] = 0 +[NumHoles] = 0 +[NumBlockLabels] = 21 +-5.00000000000000000e1 -5.00000000000000000e1 1 -1 0 0.00000000000000000e0 0 1 2 +-1.50000000000000000e1 0.00000000000000000e0 2 5.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +6.62700000000000045e2 0.00000000000000000e0 3 5.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +3.20000000000000000e2 -2.62500000000000000e1 3 1.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +3.20000000000000000e2 -1.57500000000000000e1 3 1.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +3.20000000000000000e2 -5.25000000000000000e0 3 1.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +3.20000000000000000e2 5.25000000000000000e0 3 1.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +3.20000000000000000e2 1.57500000000000000e1 3 1.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +3.20000000000000000e2 2.62500000000000000e1 3 1.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +4.62639999999999986e2 2.83769999999999989e1 2 5.00000000000000000e-1 0 2.70000000000000000e2 0 1 0 +4.62639999999999986e2 2.41230000000000011e1 2 5.00000000000000000e-1 0 2.70000000000000000e2 0 1 0 +4.62639999999999986e2 1.79149999999999991e1 2 5.00000000000000000e-1 0 2.70000000000000000e2 0 1 0 +4.62639999999999986e2 1.35850000000000009e1 2 5.00000000000000000e-1 0 2.70000000000000000e2 0 1 0 +4.62639999999999986e2 7.46600000000000019e0 2 5.00000000000000000e-1 0 2.70000000000000000e2 0 1 0 +4.62639999999999986e2 3.03399999999999981e0 2 5.00000000000000000e-1 0 2.70000000000000000e2 0 1 0 +4.62639999999999986e2 -2.91999999999999993e0 2 5.00000000000000000e-1 0 2.70000000000000000e2 0 1 0 +4.62639999999999986e2 -7.58000000000000007e0 2 5.00000000000000000e-1 0 2.70000000000000000e2 0 1 0 +4.62639999999999986e2 -1.32929999999999993e1 2 5.00000000000000000e-1 0 2.70000000000000000e2 0 1 0 +4.62639999999999986e2 -1.82070000000000007e1 2 5.00000000000000000e-1 0 2.70000000000000000e2 0 1 0 +4.62639999999999986e2 -2.36660000000000004e1 2 5.00000000000000000e-1 0 2.70000000000000000e2 0 1 0 +4.62639999999999986e2 -2.88339999999999996e1 2 5.00000000000000000e-1 0 2.70000000000000000e2 0 1 0 +[/base-doc] +# femm42 simulation metadata +# track node indices reference the doc loaded above; editing the doc may invalidate them. + +[source] +fem = /Volumes/External/Repositories/femm42src/examples/guitar_bar_and_discs.fem + +[sim] +dt = 278.300us +interval = 50ms +subdivisions = 20 +buffer_size = 60 +match_cycle = true +loop_playback = true +fundamental_hz = 60 + +[track] +label = track 1 +anchor_a = 8 +anchor_b = 22 +axis = +y +expression = series(s, t, 60, 5) +subdivisions = 20 +members = 7 21 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 +edge = 8 22 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 +edge = 7 21 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 + +[track] +label = track 2 +anchor_a = 10 +anchor_b = 24 +axis = +y +expression = series(s, t, 90, 5) +subdivisions = 20 +members = 9 23 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 +edge = 10 24 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 +edge = 9 23 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 + +[track] +label = track 3 +anchor_a = 12 +anchor_b = 26 +axis = +y +expression = series(s, t, 120, 5) +subdivisions = 20 +members = 11 25 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 +edge = 12 26 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 +edge = 11 25 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 + +[track] +label = track 4 +anchor_a = 14 +anchor_b = 28 +axis = +y +expression = series(s, t, 180, 5) +subdivisions = 20 +members = 13 27 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 +edge = 14 28 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 +edge = 13 27 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 + +[track] +label = track 5 +anchor_a = 16 +anchor_b = 30 +axis = +y +expression = series(s, t, 200, 5) +subdivisions = 20 +members = 15 29 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 +edge = 16 30 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 +edge = 15 29 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 + +[track] +label = track 6 +anchor_a = 18 +anchor_b = 32 +axis = +y +expression = series(s, t, 240, 5) +subdivisions = 20 +members = 17 31 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 +edge = 18 32 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 +edge = 17 31 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 + diff --git a/examples/magnettes_singles_guitar_qtrlength.fem b/examples/magnettes_singles_guitar_qtrlength.fem new file mode 100644 index 0000000..6b9db2d --- /dev/null +++ b/examples/magnettes_singles_guitar_qtrlength.fem @@ -0,0 +1,233 @@ +[Format] = 4.0 +[Frequency] = 0 +[Precision] = 1e-008 +[MinAngle] = 25 +[Depth] = 5 +[LengthUnits] = millimeters +[ProblemType] = planar +[Coordinates] = cartesian +[Comment] = "quarter-scale fretboard: six plain steel strings (.009/.011/.016/.026/.036/.046) spanning 162 mm. each string is anchored on the left to its own isolated NdFeB pole-piece magnet (mag_dir 0, N face into the string) and on the right to a shared steel block that ties all six strings together. a bar magnet on the right face of the steel block (mag_dir 0) drives flux through the steel and out into the air return." +[PointProps] = 0 +[BdryProps] = 1 + + = "A=0" + = 0 + = 0 + = 0 + = 0 + = 0 + = 0 + = 0 + = 0 + = 0 + = 0 + = 0 + +[BlockProps] = 3 + + = "Air" + = 1 + = 1 + = 0 + = 0 + = 0 + = 0 + = 0 + = 0 + = 0 + = 0 + = 0 + = 0 + = 1 + = 0 + = 0 + = 0 + + + = "NdFeB" + = 1.05 + = 1.05 + = 915000 + = 0 + = 0 + = 0 + = 0.667 + = 0 + = 0 + = 0 + = 0 + = 0 + = 1 + = 0 + = 0 + = 0 + + + = "Steel" + = 2500 + = 2500 + = 0 + = 0 + = 0 + = 0 + = 5.8 + = 0 + = 0 + = 0 + = 0 + = 0 + = 1 + = 0 + = 0 + = 0 + +[CircuitProps] = 0 +[NumPoints] = 58 +-100 -100 0 0 +300 -100 0 0 +300 100 0 0 +-100 100 0 0 +-10 -30.25 0 0 +0 -30.25 0 0 +0 -26.8342 0 0 +0 -25.6658 0 0 +0 -22.25 0 0 +-10 -22.25 0 0 +-10 -19.75 0 0 +0 -19.75 0 0 +0 -16.2072 0 0 +0 -15.2928 0 0 +0 -11.75 0 0 +-10 -11.75 0 0 +-10 -9.25 0 0 +0 -9.25 0 0 +0 -5.5802 0 0 +0 -4.9198 0 0 +0 -1.25 0 0 +-10 -1.25 0 0 +-10 1.25 0 0 +0 1.25 0 0 +0 5.0468 0 0 +0 5.4532 0 0 +0 9.25 0 0 +-10 9.25 0 0 +-10 11.75 0 0 +0 11.75 0 0 +0 15.6103 0 0 +0 15.8897 0 0 +0 19.75 0 0 +-10 19.75 0 0 +-10 22.25 0 0 +0 22.25 0 0 +0 26.1357 0 0 +0 26.3643 0 0 +0 30.25 0 0 +-10 30.25 0 0 +162 -30 0 0 +162 -26.8342 0 0 +162 -25.6658 0 0 +162 -16.2072 0 0 +162 -15.2928 0 0 +162 -5.5802 0 0 +162 -4.9198 0 0 +162 5.0468 0 0 +162 5.4532 0 0 +162 15.6103 0 0 +162 15.8897 0 0 +162 26.1357 0 0 +162 26.3643 0 0 +162 30 0 0 +192 -30 0 0 +192 30 0 0 +222 -30 0 0 +222 30 0 0 +[NumSegments] = 71 +0 1 -1 1 0 0 +1 2 -1 1 0 0 +2 3 -1 1 0 0 +3 0 -1 1 0 0 +4 9 -1 0 0 0 +9 8 -1 0 0 0 +8 7 -1 0 0 0 +7 6 -1 0 0 0 +6 5 -1 0 0 0 +5 4 -1 0 0 0 +10 15 -1 0 0 0 +15 14 -1 0 0 0 +14 13 -1 0 0 0 +13 12 -1 0 0 0 +12 11 -1 0 0 0 +11 10 -1 0 0 0 +16 21 -1 0 0 0 +21 20 -1 0 0 0 +20 19 -1 0 0 0 +19 18 -1 0 0 0 +18 17 -1 0 0 0 +17 16 -1 0 0 0 +22 27 -1 0 0 0 +27 26 -1 0 0 0 +26 25 -1 0 0 0 +25 24 -1 0 0 0 +24 23 -1 0 0 0 +23 22 -1 0 0 0 +28 33 -1 0 0 0 +33 32 -1 0 0 0 +32 31 -1 0 0 0 +31 30 -1 0 0 0 +30 29 -1 0 0 0 +29 28 -1 0 0 0 +34 39 -1 0 0 0 +39 38 -1 0 0 0 +38 37 -1 0 0 0 +37 36 -1 0 0 0 +36 35 -1 0 0 0 +35 34 -1 0 0 0 +7 42 -1 0 0 0 +6 41 -1 0 0 0 +13 44 -1 0 0 0 +12 43 -1 0 0 0 +19 46 -1 0 0 0 +18 45 -1 0 0 0 +25 48 -1 0 0 0 +24 47 -1 0 0 0 +31 50 -1 0 0 0 +30 49 -1 0 0 0 +37 52 -1 0 0 0 +36 51 -1 0 0 0 +40 41 -1 0 0 0 +41 42 -1 0 0 0 +42 43 -1 0 0 0 +43 44 -1 0 0 0 +44 45 -1 0 0 0 +45 46 -1 0 0 0 +46 47 -1 0 0 0 +47 48 -1 0 0 0 +48 49 -1 0 0 0 +49 50 -1 0 0 0 +50 51 -1 0 0 0 +51 52 -1 0 0 0 +52 53 -1 0 0 0 +54 55 -1 0 0 0 +53 55 -1 0 0 0 +54 40 -1 0 0 0 +56 57 -1 0 0 0 +55 57 -1 0 0 0 +56 54 -1 0 0 0 +[NumArcSegments] = 0 +[NumHoles] = 0 +[NumBlockLabels] = 15 +-50 0 1 -1 0 0 0 1 2 +-5 -26.25 2 5 0 0 0 1 0 +-5 -15.75 2 5 0 0 0 1 0 +-5 -5.25 2 5 0 0 0 1 0 +-5 5.25 2 5 0 0 0 1 0 +-5 15.75 2 5 0 0 0 1 0 +-5 26.25 2 5 0 0 0 1 0 +81 -26.25 3 1 0 0 0 1 0 +81 -15.75 3 1 0 0 0 1 0 +81 -5.25 3 1 0 0 0 1 0 +81 5.25 3 1 0 0 0 1 0 +81 15.75 3 1 0 0 0 1 0 +81 26.25 3 1 0 0 0 1 0 +177 0 3 5 0 0 0 1 0 +207 0 2 5 0 0 0 1 0 diff --git a/examples/magnettes_singles_guitar_qtrlength.femsess b/examples/magnettes_singles_guitar_qtrlength.femsess new file mode 100644 index 0000000..a6fc6b7 --- /dev/null +++ b/examples/magnettes_singles_guitar_qtrlength.femsess @@ -0,0 +1,1020 @@ +# femm42 session v1 +[pristine-doc] +[Format] = 4.00000000000000000e0 +[Frequency] = 0.00000000000000000e0 +[Precision] = 1.00000000000000002e-8 +[MinAngle] = 2.50000000000000000e1 +[DoSmartMesh] = 1 +[Depth] = 5.00000000000000000e0 +[LengthUnits] = millimeters +[ProblemType] = planar +[Coordinates] = cartesian +[ACSolver] = 0 +[PrevType] = 0 +[PrevSoln] = "" +[Comment] = "quarter-scale fretboard: six plain steel strings (.009/.011/.016/.026/.036/.046) spanning 162 mm. each string is anchored on the left to its own isolated NdFeB pole-piece magnet (mag_dir 0, N face into the string) and on the right to a shared steel block that ties all six strings together. a bar magnet on the right face of the steel block (mag_dir 0) drives flux through the steel and out into the air return." +[PointProps] = 0 +[BdryProps] = 1 + + = "A=0" + = 0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + +[BlockProps] = 3 + + = "Air" + = 1.00000000000000000e0 + = 1.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0 + = 1.00000000000000000e0 + = 0 + = 0.00000000000000000e0 + = 0 + + + = "NdFeB" + = 1.05000000000000004e0 + = 1.05000000000000004e0 + = 9.15000000000000000e5 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 6.67000000000000037e-1 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0 + = 1.00000000000000000e0 + = 0 + = 0.00000000000000000e0 + = 0 + + + = "Steel" + = 2.50000000000000000e3 + = 2.50000000000000000e3 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 5.79999999999999982e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0 + = 1.00000000000000000e0 + = 0 + = 0.00000000000000000e0 + = 0 + +[CircuitProps] = 0 +[NumPoints] = 58 +-1.00000000000000000e2 -1.00000000000000000e2 0 0 +3.00000000000000000e2 -1.00000000000000000e2 0 0 +3.00000000000000000e2 1.00000000000000000e2 0 0 +-1.00000000000000000e2 1.00000000000000000e2 0 0 +-1.00000000000000000e1 -3.02500000000000000e1 0 0 +0.00000000000000000e0 -3.02500000000000000e1 0 0 +0.00000000000000000e0 -2.68341999999999992e1 0 0 +0.00000000000000000e0 -2.56658000000000008e1 0 0 +0.00000000000000000e0 -2.22500000000000000e1 0 0 +-1.00000000000000000e1 -2.22500000000000000e1 0 0 +-1.00000000000000000e1 -1.97500000000000000e1 0 0 +0.00000000000000000e0 -1.97500000000000000e1 0 0 +0.00000000000000000e0 -1.62072000000000003e1 0 0 +0.00000000000000000e0 -1.52927999999999997e1 0 0 +0.00000000000000000e0 -1.17500000000000000e1 0 0 +-1.00000000000000000e1 -1.17500000000000000e1 0 0 +-1.00000000000000000e1 -9.25000000000000000e0 0 0 +0.00000000000000000e0 -9.25000000000000000e0 0 0 +0.00000000000000000e0 -5.58019999999999960e0 0 0 +0.00000000000000000e0 -4.91980000000000040e0 0 0 +0.00000000000000000e0 -1.25000000000000000e0 0 0 +-1.00000000000000000e1 -1.25000000000000000e0 0 0 +-1.00000000000000000e1 1.25000000000000000e0 0 0 +0.00000000000000000e0 1.25000000000000000e0 0 0 +0.00000000000000000e0 5.04680000000000017e0 0 0 +0.00000000000000000e0 5.45319999999999983e0 0 0 +0.00000000000000000e0 9.25000000000000000e0 0 0 +-1.00000000000000000e1 9.25000000000000000e0 0 0 +-1.00000000000000000e1 1.17500000000000000e1 0 0 +0.00000000000000000e0 1.17500000000000000e1 0 0 +0.00000000000000000e0 1.56103000000000005e1 0 0 +0.00000000000000000e0 1.58896999999999995e1 0 0 +0.00000000000000000e0 1.97500000000000000e1 0 0 +-1.00000000000000000e1 1.97500000000000000e1 0 0 +-1.00000000000000000e1 2.22500000000000000e1 0 0 +0.00000000000000000e0 2.22500000000000000e1 0 0 +0.00000000000000000e0 2.61356999999999999e1 0 0 +0.00000000000000000e0 2.63643000000000001e1 0 0 +0.00000000000000000e0 3.02500000000000000e1 0 0 +-1.00000000000000000e1 3.02500000000000000e1 0 0 +1.62000000000000000e2 -3.00000000000000000e1 0 0 +1.62000000000000000e2 -2.68341999999999992e1 0 0 +1.62000000000000000e2 -2.56658000000000008e1 0 0 +1.62000000000000000e2 -1.62072000000000003e1 0 0 +1.62000000000000000e2 -1.52927999999999997e1 0 0 +1.62000000000000000e2 -5.58019999999999960e0 0 0 +1.62000000000000000e2 -4.91980000000000040e0 0 0 +1.62000000000000000e2 5.04680000000000017e0 0 0 +1.62000000000000000e2 5.45319999999999983e0 0 0 +1.62000000000000000e2 1.56103000000000005e1 0 0 +1.62000000000000000e2 1.58896999999999995e1 0 0 +1.62000000000000000e2 2.61356999999999999e1 0 0 +1.62000000000000000e2 2.63643000000000001e1 0 0 +1.62000000000000000e2 3.00000000000000000e1 0 0 +1.92000000000000000e2 -3.00000000000000000e1 0 0 +1.92000000000000000e2 3.00000000000000000e1 0 0 +2.22000000000000000e2 -3.00000000000000000e1 0 0 +2.22000000000000000e2 3.00000000000000000e1 0 0 +[NumSegments] = 71 +0 1 -1 1 0 0 +1 2 -1 1 0 0 +2 3 -1 1 0 0 +3 0 -1 1 0 0 +4 9 -1 0 0 0 +9 8 -1 0 0 0 +8 7 -1 0 0 0 +7 6 -1 0 0 0 +6 5 -1 0 0 0 +5 4 -1 0 0 0 +10 15 -1 0 0 0 +15 14 -1 0 0 0 +14 13 -1 0 0 0 +13 12 -1 0 0 0 +12 11 -1 0 0 0 +11 10 -1 0 0 0 +16 21 -1 0 0 0 +21 20 -1 0 0 0 +20 19 -1 0 0 0 +19 18 -1 0 0 0 +18 17 -1 0 0 0 +17 16 -1 0 0 0 +22 27 -1 0 0 0 +27 26 -1 0 0 0 +26 25 -1 0 0 0 +25 24 -1 0 0 0 +24 23 -1 0 0 0 +23 22 -1 0 0 0 +28 33 -1 0 0 0 +33 32 -1 0 0 0 +32 31 -1 0 0 0 +31 30 -1 0 0 0 +30 29 -1 0 0 0 +29 28 -1 0 0 0 +34 39 -1 0 0 0 +39 38 -1 0 0 0 +38 37 -1 0 0 0 +37 36 -1 0 0 0 +36 35 -1 0 0 0 +35 34 -1 0 0 0 +7 42 -1 0 0 0 +6 41 -1 0 0 0 +13 44 -1 0 0 0 +12 43 -1 0 0 0 +19 46 -1 0 0 0 +18 45 -1 0 0 0 +25 48 -1 0 0 0 +24 47 -1 0 0 0 +31 50 -1 0 0 0 +30 49 -1 0 0 0 +37 52 -1 0 0 0 +36 51 -1 0 0 0 +40 41 -1 0 0 0 +41 42 -1 0 0 0 +42 43 -1 0 0 0 +43 44 -1 0 0 0 +44 45 -1 0 0 0 +45 46 -1 0 0 0 +46 47 -1 0 0 0 +47 48 -1 0 0 0 +48 49 -1 0 0 0 +49 50 -1 0 0 0 +50 51 -1 0 0 0 +51 52 -1 0 0 0 +52 53 -1 0 0 0 +54 55 -1 0 0 0 +53 55 -1 0 0 0 +54 40 -1 0 0 0 +56 57 -1 0 0 0 +55 57 -1 0 0 0 +56 54 -1 0 0 0 +[NumArcSegments] = 0 +[NumHoles] = 0 +[NumBlockLabels] = 15 +-5.00000000000000000e1 0.00000000000000000e0 1 -1 0 0.00000000000000000e0 0 1 2 +-5.00000000000000000e0 -2.62500000000000000e1 2 5.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +-5.00000000000000000e0 -1.57500000000000000e1 2 5.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +-5.00000000000000000e0 -5.25000000000000000e0 2 5.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +-5.00000000000000000e0 5.25000000000000000e0 2 5.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +-5.00000000000000000e0 1.57500000000000000e1 2 5.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +-5.00000000000000000e0 2.62500000000000000e1 2 5.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +8.10000000000000000e1 -2.62500000000000000e1 3 1.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +8.10000000000000000e1 -1.57500000000000000e1 3 1.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +8.10000000000000000e1 -5.25000000000000000e0 3 1.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +8.10000000000000000e1 5.25000000000000000e0 3 1.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +8.10000000000000000e1 1.57500000000000000e1 3 1.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +8.10000000000000000e1 2.62500000000000000e1 3 1.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +1.77000000000000000e2 0.00000000000000000e0 3 5.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +2.07000000000000000e2 0.00000000000000000e0 2 5.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +[/pristine-doc] +[base-doc] +[Format] = 4.00000000000000000e0 +[Frequency] = 0.00000000000000000e0 +[Precision] = 1.00000000000000002e-8 +[MinAngle] = 2.50000000000000000e1 +[DoSmartMesh] = 1 +[Depth] = 5.00000000000000000e0 +[LengthUnits] = millimeters +[ProblemType] = planar +[Coordinates] = cartesian +[ACSolver] = 0 +[PrevType] = 0 +[PrevSoln] = "" +[Comment] = "quarter-scale fretboard: six plain steel strings (.009/.011/.016/.026/.036/.046) spanning 162 mm. each string is anchored on the left to its own isolated NdFeB pole-piece magnet (mag_dir 0, N face into the string) and on the right to a shared steel block that ties all six strings together. a bar magnet on the right face of the steel block (mag_dir 0) drives flux through the steel and out into the air return." +[PointProps] = 0 +[BdryProps] = 1 + + = "A=0" + = 0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + +[BlockProps] = 3 + + = "Air" + = 1.00000000000000000e0 + = 1.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0 + = 1.00000000000000000e0 + = 0 + = 0.00000000000000000e0 + = 0 + + + = "NdFeB" + = 1.05000000000000004e0 + = 1.05000000000000004e0 + = 9.15000000000000000e5 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 6.67000000000000037e-1 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0 + = 1.00000000000000000e0 + = 0 + = 0.00000000000000000e0 + = 0 + + + = "Steel" + = 2.50000000000000000e3 + = 2.50000000000000000e3 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 5.79999999999999982e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0.00000000000000000e0 + = 0 + = 1.00000000000000000e0 + = 0 + = 0.00000000000000000e0 + = 0 + +[CircuitProps] = 0 +[NumPoints] = 286 +-1.00000000000000000e2 -1.00000000000000000e2 0 0 +3.00000000000000000e2 -1.00000000000000000e2 0 0 +3.00000000000000000e2 1.00000000000000000e2 0 0 +-1.00000000000000000e2 1.00000000000000000e2 0 0 +-1.00000000000000000e1 -3.02500000000000000e1 0 0 +0.00000000000000000e0 -3.02500000000000000e1 0 0 +0.00000000000000000e0 -2.68341999999999992e1 0 0 +0.00000000000000000e0 -2.56658000000000008e1 0 0 +0.00000000000000000e0 -2.22500000000000000e1 0 0 +-1.00000000000000000e1 -2.22500000000000000e1 0 0 +-1.00000000000000000e1 -1.97500000000000000e1 0 0 +0.00000000000000000e0 -1.97500000000000000e1 0 0 +0.00000000000000000e0 -1.62072000000000003e1 0 0 +0.00000000000000000e0 -1.52927999999999997e1 0 0 +0.00000000000000000e0 -1.17500000000000000e1 0 0 +-1.00000000000000000e1 -1.17500000000000000e1 0 0 +-1.00000000000000000e1 -9.25000000000000000e0 0 0 +0.00000000000000000e0 -9.25000000000000000e0 0 0 +0.00000000000000000e0 -5.58019999999999960e0 0 0 +0.00000000000000000e0 -4.91980000000000040e0 0 0 +0.00000000000000000e0 -1.25000000000000000e0 0 0 +-1.00000000000000000e1 -1.25000000000000000e0 0 0 +-1.00000000000000000e1 1.25000000000000000e0 0 0 +0.00000000000000000e0 1.25000000000000000e0 0 0 +0.00000000000000000e0 5.04680000000000017e0 0 0 +0.00000000000000000e0 5.45319999999999983e0 0 0 +0.00000000000000000e0 9.25000000000000000e0 0 0 +-1.00000000000000000e1 9.25000000000000000e0 0 0 +-1.00000000000000000e1 1.17500000000000000e1 0 0 +0.00000000000000000e0 1.17500000000000000e1 0 0 +0.00000000000000000e0 1.56103000000000005e1 0 0 +0.00000000000000000e0 1.58896999999999995e1 0 0 +0.00000000000000000e0 1.97500000000000000e1 0 0 +-1.00000000000000000e1 1.97500000000000000e1 0 0 +-1.00000000000000000e1 2.22500000000000000e1 0 0 +0.00000000000000000e0 2.22500000000000000e1 0 0 +0.00000000000000000e0 2.61356999999999999e1 0 0 +0.00000000000000000e0 2.63643000000000001e1 0 0 +0.00000000000000000e0 3.02500000000000000e1 0 0 +-1.00000000000000000e1 3.02500000000000000e1 0 0 +1.62000000000000000e2 -3.00000000000000000e1 0 0 +1.62000000000000000e2 -2.68341999999999992e1 0 0 +1.62000000000000000e2 -2.56658000000000008e1 0 0 +1.62000000000000000e2 -1.62072000000000003e1 0 0 +1.62000000000000000e2 -1.52927999999999997e1 0 0 +1.62000000000000000e2 -5.58019999999999960e0 0 0 +1.62000000000000000e2 -4.91980000000000040e0 0 0 +1.62000000000000000e2 5.04680000000000017e0 0 0 +1.62000000000000000e2 5.45319999999999983e0 0 0 +1.62000000000000000e2 1.56103000000000005e1 0 0 +1.62000000000000000e2 1.58896999999999995e1 0 0 +1.62000000000000000e2 2.61356999999999999e1 0 0 +1.62000000000000000e2 2.63643000000000001e1 0 0 +1.62000000000000000e2 3.00000000000000000e1 0 0 +1.92000000000000000e2 -3.00000000000000000e1 0 0 +1.92000000000000000e2 3.00000000000000000e1 0 0 +2.22000000000000000e2 -3.00000000000000000e1 0 0 +2.22000000000000000e2 3.00000000000000000e1 0 0 +8.09999999999999964e0 -2.68341999999999992e1 0 0 +1.61999999999999993e1 -2.68341999999999992e1 0 0 +2.42999999999999972e1 -2.68341999999999992e1 0 0 +3.23999999999999986e1 -2.68341999999999992e1 0 0 +4.05000000000000000e1 -2.68341999999999992e1 0 0 +4.85999999999999943e1 -2.68341999999999992e1 0 0 +5.66999999999999957e1 -2.68341999999999992e1 0 0 +6.47999999999999972e1 -2.68341999999999992e1 0 0 +7.28999999999999915e1 -2.68341999999999992e1 0 0 +8.10000000000000000e1 -2.68341999999999992e1 0 0 +8.90999999999999943e1 -2.68341999999999992e1 0 0 +9.71999999999999886e1 -2.68341999999999992e1 0 0 +1.05299999999999997e2 -2.68341999999999992e1 0 0 +1.13399999999999991e2 -2.68341999999999992e1 0 0 +1.21500000000000000e2 -2.68341999999999992e1 0 0 +1.29599999999999994e2 -2.68341999999999992e1 0 0 +1.37699999999999989e2 -2.68341999999999992e1 0 0 +1.45799999999999983e2 -2.68341999999999992e1 0 0 +1.53900000000000006e2 -2.68341999999999992e1 0 0 +8.09999999999999964e0 -2.56658000000000008e1 0 0 +1.61999999999999993e1 -2.56658000000000008e1 0 0 +2.42999999999999972e1 -2.56658000000000008e1 0 0 +3.23999999999999986e1 -2.56658000000000008e1 0 0 +4.05000000000000000e1 -2.56658000000000008e1 0 0 +4.85999999999999943e1 -2.56658000000000008e1 0 0 +5.66999999999999957e1 -2.56658000000000008e1 0 0 +6.47999999999999972e1 -2.56658000000000008e1 0 0 +7.28999999999999915e1 -2.56658000000000008e1 0 0 +8.10000000000000000e1 -2.56658000000000008e1 0 0 +8.90999999999999943e1 -2.56658000000000008e1 0 0 +9.71999999999999886e1 -2.56658000000000008e1 0 0 +1.05299999999999997e2 -2.56658000000000008e1 0 0 +1.13399999999999991e2 -2.56658000000000008e1 0 0 +1.21500000000000000e2 -2.56658000000000008e1 0 0 +1.29599999999999994e2 -2.56658000000000008e1 0 0 +1.37699999999999989e2 -2.56658000000000008e1 0 0 +1.45799999999999983e2 -2.56658000000000008e1 0 0 +1.53900000000000006e2 -2.56658000000000008e1 0 0 +8.09999999999999964e0 -1.62072000000000003e1 0 0 +1.61999999999999993e1 -1.62072000000000003e1 0 0 +2.42999999999999972e1 -1.62072000000000003e1 0 0 +3.23999999999999986e1 -1.62072000000000003e1 0 0 +4.05000000000000000e1 -1.62072000000000003e1 0 0 +4.85999999999999943e1 -1.62072000000000003e1 0 0 +5.66999999999999957e1 -1.62072000000000003e1 0 0 +6.47999999999999972e1 -1.62072000000000003e1 0 0 +7.28999999999999915e1 -1.62072000000000003e1 0 0 +8.10000000000000000e1 -1.62072000000000003e1 0 0 +8.90999999999999943e1 -1.62072000000000003e1 0 0 +9.71999999999999886e1 -1.62072000000000003e1 0 0 +1.05299999999999997e2 -1.62072000000000003e1 0 0 +1.13399999999999991e2 -1.62072000000000003e1 0 0 +1.21500000000000000e2 -1.62072000000000003e1 0 0 +1.29599999999999994e2 -1.62072000000000003e1 0 0 +1.37699999999999989e2 -1.62072000000000003e1 0 0 +1.45799999999999983e2 -1.62072000000000003e1 0 0 +1.53900000000000006e2 -1.62072000000000003e1 0 0 +8.09999999999999964e0 -1.52927999999999997e1 0 0 +1.61999999999999993e1 -1.52927999999999997e1 0 0 +2.42999999999999972e1 -1.52927999999999997e1 0 0 +3.23999999999999986e1 -1.52927999999999997e1 0 0 +4.05000000000000000e1 -1.52927999999999997e1 0 0 +4.85999999999999943e1 -1.52927999999999997e1 0 0 +5.66999999999999957e1 -1.52927999999999997e1 0 0 +6.47999999999999972e1 -1.52927999999999997e1 0 0 +7.28999999999999915e1 -1.52927999999999997e1 0 0 +8.10000000000000000e1 -1.52927999999999997e1 0 0 +8.90999999999999943e1 -1.52927999999999997e1 0 0 +9.71999999999999886e1 -1.52927999999999997e1 0 0 +1.05299999999999997e2 -1.52927999999999997e1 0 0 +1.13399999999999991e2 -1.52927999999999997e1 0 0 +1.21500000000000000e2 -1.52927999999999997e1 0 0 +1.29599999999999994e2 -1.52927999999999997e1 0 0 +1.37699999999999989e2 -1.52927999999999997e1 0 0 +1.45799999999999983e2 -1.52927999999999997e1 0 0 +1.53900000000000006e2 -1.52927999999999997e1 0 0 +8.09999999999999964e0 -5.58019999999999960e0 0 0 +1.61999999999999993e1 -5.58019999999999960e0 0 0 +2.42999999999999972e1 -5.58019999999999960e0 0 0 +3.23999999999999986e1 -5.58019999999999960e0 0 0 +4.05000000000000000e1 -5.58019999999999960e0 0 0 +4.85999999999999943e1 -5.58019999999999960e0 0 0 +5.66999999999999957e1 -5.58019999999999960e0 0 0 +6.47999999999999972e1 -5.58019999999999960e0 0 0 +7.28999999999999915e1 -5.58019999999999960e0 0 0 +8.10000000000000000e1 -5.58019999999999960e0 0 0 +8.90999999999999943e1 -5.58019999999999960e0 0 0 +9.71999999999999886e1 -5.58019999999999960e0 0 0 +1.05299999999999997e2 -5.58019999999999960e0 0 0 +1.13399999999999991e2 -5.58019999999999960e0 0 0 +1.21500000000000000e2 -5.58019999999999960e0 0 0 +1.29599999999999994e2 -5.58019999999999960e0 0 0 +1.37699999999999989e2 -5.58019999999999960e0 0 0 +1.45799999999999983e2 -5.58019999999999960e0 0 0 +1.53900000000000006e2 -5.58019999999999960e0 0 0 +8.09999999999999964e0 -4.91980000000000040e0 0 0 +1.61999999999999993e1 -4.91980000000000040e0 0 0 +2.42999999999999972e1 -4.91980000000000040e0 0 0 +3.23999999999999986e1 -4.91980000000000040e0 0 0 +4.05000000000000000e1 -4.91980000000000040e0 0 0 +4.85999999999999943e1 -4.91980000000000040e0 0 0 +5.66999999999999957e1 -4.91980000000000040e0 0 0 +6.47999999999999972e1 -4.91980000000000040e0 0 0 +7.28999999999999915e1 -4.91980000000000040e0 0 0 +8.10000000000000000e1 -4.91980000000000040e0 0 0 +8.90999999999999943e1 -4.91980000000000040e0 0 0 +9.71999999999999886e1 -4.91980000000000040e0 0 0 +1.05299999999999997e2 -4.91980000000000040e0 0 0 +1.13399999999999991e2 -4.91980000000000040e0 0 0 +1.21500000000000000e2 -4.91980000000000040e0 0 0 +1.29599999999999994e2 -4.91980000000000040e0 0 0 +1.37699999999999989e2 -4.91980000000000040e0 0 0 +1.45799999999999983e2 -4.91980000000000040e0 0 0 +1.53900000000000006e2 -4.91980000000000040e0 0 0 +8.09999999999999964e0 5.04680000000000017e0 0 0 +1.61999999999999993e1 5.04680000000000017e0 0 0 +2.42999999999999972e1 5.04680000000000017e0 0 0 +3.23999999999999986e1 5.04680000000000017e0 0 0 +4.05000000000000000e1 5.04680000000000017e0 0 0 +4.85999999999999943e1 5.04680000000000017e0 0 0 +5.66999999999999957e1 5.04680000000000017e0 0 0 +6.47999999999999972e1 5.04680000000000017e0 0 0 +7.28999999999999915e1 5.04680000000000017e0 0 0 +8.10000000000000000e1 5.04680000000000017e0 0 0 +8.90999999999999943e1 5.04680000000000017e0 0 0 +9.71999999999999886e1 5.04680000000000017e0 0 0 +1.05299999999999997e2 5.04680000000000017e0 0 0 +1.13399999999999991e2 5.04680000000000017e0 0 0 +1.21500000000000000e2 5.04680000000000017e0 0 0 +1.29599999999999994e2 5.04680000000000017e0 0 0 +1.37699999999999989e2 5.04680000000000017e0 0 0 +1.45799999999999983e2 5.04680000000000017e0 0 0 +1.53900000000000006e2 5.04680000000000017e0 0 0 +8.09999999999999964e0 5.45319999999999983e0 0 0 +1.61999999999999993e1 5.45319999999999983e0 0 0 +2.42999999999999972e1 5.45319999999999983e0 0 0 +3.23999999999999986e1 5.45319999999999983e0 0 0 +4.05000000000000000e1 5.45319999999999983e0 0 0 +4.85999999999999943e1 5.45319999999999983e0 0 0 +5.66999999999999957e1 5.45319999999999983e0 0 0 +6.47999999999999972e1 5.45319999999999983e0 0 0 +7.28999999999999915e1 5.45319999999999983e0 0 0 +8.10000000000000000e1 5.45319999999999983e0 0 0 +8.90999999999999943e1 5.45319999999999983e0 0 0 +9.71999999999999886e1 5.45319999999999983e0 0 0 +1.05299999999999997e2 5.45319999999999983e0 0 0 +1.13399999999999991e2 5.45319999999999983e0 0 0 +1.21500000000000000e2 5.45319999999999983e0 0 0 +1.29599999999999994e2 5.45319999999999983e0 0 0 +1.37699999999999989e2 5.45319999999999983e0 0 0 +1.45799999999999983e2 5.45319999999999983e0 0 0 +1.53900000000000006e2 5.45319999999999983e0 0 0 +8.09999999999999964e0 1.56103000000000005e1 0 0 +1.61999999999999993e1 1.56103000000000005e1 0 0 +2.42999999999999972e1 1.56103000000000005e1 0 0 +3.23999999999999986e1 1.56103000000000005e1 0 0 +4.05000000000000000e1 1.56103000000000005e1 0 0 +4.85999999999999943e1 1.56103000000000005e1 0 0 +5.66999999999999957e1 1.56103000000000005e1 0 0 +6.47999999999999972e1 1.56103000000000005e1 0 0 +7.28999999999999915e1 1.56103000000000005e1 0 0 +8.10000000000000000e1 1.56103000000000005e1 0 0 +8.90999999999999943e1 1.56103000000000005e1 0 0 +9.71999999999999886e1 1.56103000000000005e1 0 0 +1.05299999999999997e2 1.56103000000000005e1 0 0 +1.13399999999999991e2 1.56103000000000005e1 0 0 +1.21500000000000000e2 1.56103000000000005e1 0 0 +1.29599999999999994e2 1.56103000000000005e1 0 0 +1.37699999999999989e2 1.56103000000000005e1 0 0 +1.45799999999999983e2 1.56103000000000005e1 0 0 +1.53900000000000006e2 1.56103000000000005e1 0 0 +8.09999999999999964e0 1.58896999999999995e1 0 0 +1.61999999999999993e1 1.58896999999999995e1 0 0 +2.42999999999999972e1 1.58896999999999995e1 0 0 +3.23999999999999986e1 1.58896999999999995e1 0 0 +4.05000000000000000e1 1.58896999999999995e1 0 0 +4.85999999999999943e1 1.58896999999999995e1 0 0 +5.66999999999999957e1 1.58896999999999995e1 0 0 +6.47999999999999972e1 1.58896999999999995e1 0 0 +7.28999999999999915e1 1.58896999999999995e1 0 0 +8.10000000000000000e1 1.58896999999999995e1 0 0 +8.90999999999999943e1 1.58896999999999995e1 0 0 +9.71999999999999886e1 1.58896999999999995e1 0 0 +1.05299999999999997e2 1.58896999999999995e1 0 0 +1.13399999999999991e2 1.58896999999999995e1 0 0 +1.21500000000000000e2 1.58896999999999995e1 0 0 +1.29599999999999994e2 1.58896999999999995e1 0 0 +1.37699999999999989e2 1.58896999999999995e1 0 0 +1.45799999999999983e2 1.58896999999999995e1 0 0 +1.53900000000000006e2 1.58896999999999995e1 0 0 +8.09999999999999964e0 2.61356999999999999e1 0 0 +1.61999999999999993e1 2.61356999999999999e1 0 0 +2.42999999999999972e1 2.61356999999999999e1 0 0 +3.23999999999999986e1 2.61356999999999999e1 0 0 +4.05000000000000000e1 2.61356999999999999e1 0 0 +4.85999999999999943e1 2.61356999999999999e1 0 0 +5.66999999999999957e1 2.61356999999999999e1 0 0 +6.47999999999999972e1 2.61356999999999999e1 0 0 +7.28999999999999915e1 2.61356999999999999e1 0 0 +8.10000000000000000e1 2.61356999999999999e1 0 0 +8.90999999999999943e1 2.61356999999999999e1 0 0 +9.71999999999999886e1 2.61356999999999999e1 0 0 +1.05299999999999997e2 2.61356999999999999e1 0 0 +1.13399999999999991e2 2.61356999999999999e1 0 0 +1.21500000000000000e2 2.61356999999999999e1 0 0 +1.29599999999999994e2 2.61356999999999999e1 0 0 +1.37699999999999989e2 2.61356999999999999e1 0 0 +1.45799999999999983e2 2.61356999999999999e1 0 0 +1.53900000000000006e2 2.61356999999999999e1 0 0 +8.09999999999999964e0 2.63643000000000001e1 0 0 +1.61999999999999993e1 2.63643000000000001e1 0 0 +2.42999999999999972e1 2.63643000000000001e1 0 0 +3.23999999999999986e1 2.63643000000000001e1 0 0 +4.05000000000000000e1 2.63643000000000001e1 0 0 +4.85999999999999943e1 2.63643000000000001e1 0 0 +5.66999999999999957e1 2.63643000000000001e1 0 0 +6.47999999999999972e1 2.63643000000000001e1 0 0 +7.28999999999999915e1 2.63643000000000001e1 0 0 +8.10000000000000000e1 2.63643000000000001e1 0 0 +8.90999999999999943e1 2.63643000000000001e1 0 0 +9.71999999999999886e1 2.63643000000000001e1 0 0 +1.05299999999999997e2 2.63643000000000001e1 0 0 +1.13399999999999991e2 2.63643000000000001e1 0 0 +1.21500000000000000e2 2.63643000000000001e1 0 0 +1.29599999999999994e2 2.63643000000000001e1 0 0 +1.37699999999999989e2 2.63643000000000001e1 0 0 +1.45799999999999983e2 2.63643000000000001e1 0 0 +1.53900000000000006e2 2.63643000000000001e1 0 0 +[NumSegments] = 299 +0 1 -1 1 0 0 +1 2 -1 1 0 0 +2 3 -1 1 0 0 +3 0 -1 1 0 0 +4 9 -1 0 0 0 +9 8 -1 0 0 0 +8 7 -1 0 0 0 +7 6 -1 0 0 0 +6 5 -1 0 0 0 +5 4 -1 0 0 0 +10 15 -1 0 0 0 +15 14 -1 0 0 0 +14 13 -1 0 0 0 +13 12 -1 0 0 0 +12 11 -1 0 0 0 +11 10 -1 0 0 0 +16 21 -1 0 0 0 +21 20 -1 0 0 0 +20 19 -1 0 0 0 +19 18 -1 0 0 0 +18 17 -1 0 0 0 +17 16 -1 0 0 0 +22 27 -1 0 0 0 +27 26 -1 0 0 0 +26 25 -1 0 0 0 +25 24 -1 0 0 0 +24 23 -1 0 0 0 +23 22 -1 0 0 0 +28 33 -1 0 0 0 +33 32 -1 0 0 0 +32 31 -1 0 0 0 +31 30 -1 0 0 0 +30 29 -1 0 0 0 +29 28 -1 0 0 0 +34 39 -1 0 0 0 +39 38 -1 0 0 0 +38 37 -1 0 0 0 +37 36 -1 0 0 0 +36 35 -1 0 0 0 +35 34 -1 0 0 0 +40 41 -1 0 0 0 +41 42 -1 0 0 0 +42 43 -1 0 0 0 +43 44 -1 0 0 0 +44 45 -1 0 0 0 +45 46 -1 0 0 0 +46 47 -1 0 0 0 +47 48 -1 0 0 0 +48 49 -1 0 0 0 +49 50 -1 0 0 0 +50 51 -1 0 0 0 +51 52 -1 0 0 0 +52 53 -1 0 0 0 +54 55 -1 0 0 0 +53 55 -1 0 0 0 +54 40 -1 0 0 0 +56 57 -1 0 0 0 +55 57 -1 0 0 0 +56 54 -1 0 0 0 +6 58 -1 0 0 0 +58 59 -1 0 0 0 +59 60 -1 0 0 0 +60 61 -1 0 0 0 +61 62 -1 0 0 0 +62 63 -1 0 0 0 +63 64 -1 0 0 0 +64 65 -1 0 0 0 +65 66 -1 0 0 0 +66 67 -1 0 0 0 +67 68 -1 0 0 0 +68 69 -1 0 0 0 +69 70 -1 0 0 0 +70 71 -1 0 0 0 +71 72 -1 0 0 0 +72 73 -1 0 0 0 +73 74 -1 0 0 0 +74 75 -1 0 0 0 +75 76 -1 0 0 0 +76 41 -1 0 0 0 +7 77 -1 0 0 0 +77 78 -1 0 0 0 +78 79 -1 0 0 0 +79 80 -1 0 0 0 +80 81 -1 0 0 0 +81 82 -1 0 0 0 +82 83 -1 0 0 0 +83 84 -1 0 0 0 +84 85 -1 0 0 0 +85 86 -1 0 0 0 +86 87 -1 0 0 0 +87 88 -1 0 0 0 +88 89 -1 0 0 0 +89 90 -1 0 0 0 +90 91 -1 0 0 0 +91 92 -1 0 0 0 +92 93 -1 0 0 0 +93 94 -1 0 0 0 +94 95 -1 0 0 0 +95 42 -1 0 0 0 +12 96 -1 0 0 0 +96 97 -1 0 0 0 +97 98 -1 0 0 0 +98 99 -1 0 0 0 +99 100 -1 0 0 0 +100 101 -1 0 0 0 +101 102 -1 0 0 0 +102 103 -1 0 0 0 +103 104 -1 0 0 0 +104 105 -1 0 0 0 +105 106 -1 0 0 0 +106 107 -1 0 0 0 +107 108 -1 0 0 0 +108 109 -1 0 0 0 +109 110 -1 0 0 0 +110 111 -1 0 0 0 +111 112 -1 0 0 0 +112 113 -1 0 0 0 +113 114 -1 0 0 0 +114 43 -1 0 0 0 +13 115 -1 0 0 0 +115 116 -1 0 0 0 +116 117 -1 0 0 0 +117 118 -1 0 0 0 +118 119 -1 0 0 0 +119 120 -1 0 0 0 +120 121 -1 0 0 0 +121 122 -1 0 0 0 +122 123 -1 0 0 0 +123 124 -1 0 0 0 +124 125 -1 0 0 0 +125 126 -1 0 0 0 +126 127 -1 0 0 0 +127 128 -1 0 0 0 +128 129 -1 0 0 0 +129 130 -1 0 0 0 +130 131 -1 0 0 0 +131 132 -1 0 0 0 +132 133 -1 0 0 0 +133 44 -1 0 0 0 +18 134 -1 0 0 0 +134 135 -1 0 0 0 +135 136 -1 0 0 0 +136 137 -1 0 0 0 +137 138 -1 0 0 0 +138 139 -1 0 0 0 +139 140 -1 0 0 0 +140 141 -1 0 0 0 +141 142 -1 0 0 0 +142 143 -1 0 0 0 +143 144 -1 0 0 0 +144 145 -1 0 0 0 +145 146 -1 0 0 0 +146 147 -1 0 0 0 +147 148 -1 0 0 0 +148 149 -1 0 0 0 +149 150 -1 0 0 0 +150 151 -1 0 0 0 +151 152 -1 0 0 0 +152 45 -1 0 0 0 +19 153 -1 0 0 0 +153 154 -1 0 0 0 +154 155 -1 0 0 0 +155 156 -1 0 0 0 +156 157 -1 0 0 0 +157 158 -1 0 0 0 +158 159 -1 0 0 0 +159 160 -1 0 0 0 +160 161 -1 0 0 0 +161 162 -1 0 0 0 +162 163 -1 0 0 0 +163 164 -1 0 0 0 +164 165 -1 0 0 0 +165 166 -1 0 0 0 +166 167 -1 0 0 0 +167 168 -1 0 0 0 +168 169 -1 0 0 0 +169 170 -1 0 0 0 +170 171 -1 0 0 0 +171 46 -1 0 0 0 +24 172 -1 0 0 0 +172 173 -1 0 0 0 +173 174 -1 0 0 0 +174 175 -1 0 0 0 +175 176 -1 0 0 0 +176 177 -1 0 0 0 +177 178 -1 0 0 0 +178 179 -1 0 0 0 +179 180 -1 0 0 0 +180 181 -1 0 0 0 +181 182 -1 0 0 0 +182 183 -1 0 0 0 +183 184 -1 0 0 0 +184 185 -1 0 0 0 +185 186 -1 0 0 0 +186 187 -1 0 0 0 +187 188 -1 0 0 0 +188 189 -1 0 0 0 +189 190 -1 0 0 0 +190 47 -1 0 0 0 +25 191 -1 0 0 0 +191 192 -1 0 0 0 +192 193 -1 0 0 0 +193 194 -1 0 0 0 +194 195 -1 0 0 0 +195 196 -1 0 0 0 +196 197 -1 0 0 0 +197 198 -1 0 0 0 +198 199 -1 0 0 0 +199 200 -1 0 0 0 +200 201 -1 0 0 0 +201 202 -1 0 0 0 +202 203 -1 0 0 0 +203 204 -1 0 0 0 +204 205 -1 0 0 0 +205 206 -1 0 0 0 +206 207 -1 0 0 0 +207 208 -1 0 0 0 +208 209 -1 0 0 0 +209 48 -1 0 0 0 +30 210 -1 0 0 0 +210 211 -1 0 0 0 +211 212 -1 0 0 0 +212 213 -1 0 0 0 +213 214 -1 0 0 0 +214 215 -1 0 0 0 +215 216 -1 0 0 0 +216 217 -1 0 0 0 +217 218 -1 0 0 0 +218 219 -1 0 0 0 +219 220 -1 0 0 0 +220 221 -1 0 0 0 +221 222 -1 0 0 0 +222 223 -1 0 0 0 +223 224 -1 0 0 0 +224 225 -1 0 0 0 +225 226 -1 0 0 0 +226 227 -1 0 0 0 +227 228 -1 0 0 0 +228 49 -1 0 0 0 +31 229 -1 0 0 0 +229 230 -1 0 0 0 +230 231 -1 0 0 0 +231 232 -1 0 0 0 +232 233 -1 0 0 0 +233 234 -1 0 0 0 +234 235 -1 0 0 0 +235 236 -1 0 0 0 +236 237 -1 0 0 0 +237 238 -1 0 0 0 +238 239 -1 0 0 0 +239 240 -1 0 0 0 +240 241 -1 0 0 0 +241 242 -1 0 0 0 +242 243 -1 0 0 0 +243 244 -1 0 0 0 +244 245 -1 0 0 0 +245 246 -1 0 0 0 +246 247 -1 0 0 0 +247 50 -1 0 0 0 +36 248 -1 0 0 0 +248 249 -1 0 0 0 +249 250 -1 0 0 0 +250 251 -1 0 0 0 +251 252 -1 0 0 0 +252 253 -1 0 0 0 +253 254 -1 0 0 0 +254 255 -1 0 0 0 +255 256 -1 0 0 0 +256 257 -1 0 0 0 +257 258 -1 0 0 0 +258 259 -1 0 0 0 +259 260 -1 0 0 0 +260 261 -1 0 0 0 +261 262 -1 0 0 0 +262 263 -1 0 0 0 +263 264 -1 0 0 0 +264 265 -1 0 0 0 +265 266 -1 0 0 0 +266 51 -1 0 0 0 +37 267 -1 0 0 0 +267 268 -1 0 0 0 +268 269 -1 0 0 0 +269 270 -1 0 0 0 +270 271 -1 0 0 0 +271 272 -1 0 0 0 +272 273 -1 0 0 0 +273 274 -1 0 0 0 +274 275 -1 0 0 0 +275 276 -1 0 0 0 +276 277 -1 0 0 0 +277 278 -1 0 0 0 +278 279 -1 0 0 0 +279 280 -1 0 0 0 +280 281 -1 0 0 0 +281 282 -1 0 0 0 +282 283 -1 0 0 0 +283 284 -1 0 0 0 +284 285 -1 0 0 0 +285 52 -1 0 0 0 +[NumArcSegments] = 0 +[NumHoles] = 0 +[NumBlockLabels] = 15 +-5.00000000000000000e1 0.00000000000000000e0 1 -1 0 0.00000000000000000e0 0 1 2 +-5.00000000000000000e0 -2.62500000000000000e1 2 5.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +-5.00000000000000000e0 -1.57500000000000000e1 2 5.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +-5.00000000000000000e0 -5.25000000000000000e0 2 5.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +-5.00000000000000000e0 5.25000000000000000e0 2 5.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +-5.00000000000000000e0 1.57500000000000000e1 2 5.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +-5.00000000000000000e0 2.62500000000000000e1 2 5.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +8.10000000000000000e1 -2.62500000000000000e1 3 1.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +8.10000000000000000e1 -1.57500000000000000e1 3 1.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +8.10000000000000000e1 -5.25000000000000000e0 3 1.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +8.10000000000000000e1 5.25000000000000000e0 3 1.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +8.10000000000000000e1 1.57500000000000000e1 3 1.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +8.10000000000000000e1 2.62500000000000000e1 3 1.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +1.77000000000000000e2 0.00000000000000000e0 3 5.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +2.07000000000000000e2 0.00000000000000000e0 2 5.00000000000000000e0 0 0.00000000000000000e0 0 1 0 +[/base-doc] +# femm42 simulation metadata +# track node indices reference the doc loaded above; editing the doc may invalidate them. + +[source] +fem = /Volumes/External/Repositories/femm42src/examples/magnettes_singles_guitar_qtrlength.fem + +[sim] +dt = 100us +interval = 50ms +subdivisions = 20 +buffer_size = 42 +match_cycle = true +loop_playback = true +fundamental_hz = 240 + +[track] +label = track 1 +anchor_a = 7 +anchor_b = 42 +axis = +y +expression = series(s, t, 60, 5) +subdivisions = 20 +members = 6 41 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 +edge = 7 42 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 +edge = 6 41 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 + +[track] +label = track 2 +anchor_a = 13 +anchor_b = 44 +axis = +y +expression = series(s, t, 90, 5) +subdivisions = 20 +members = 12 43 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 +edge = 13 44 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 +edge = 12 43 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 + +[track] +label = track 3 +anchor_a = 19 +anchor_b = 46 +axis = +y +expression = series(s, t, 120, 5) +subdivisions = 20 +members = 18 45 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 +edge = 19 46 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 +edge = 18 45 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 + +[track] +label = track 4 +anchor_a = 25 +anchor_b = 48 +axis = +y +expression = series(s, t, 240, 5) +subdivisions = 20 +members = 24 47 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 +edge = 25 48 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 +edge = 24 47 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 + +[track] +label = track 5 +anchor_a = 31 +anchor_b = 50 +axis = +y +expression = series(s, t, 300, 5) +subdivisions = 20 +members = 30 49 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 +edge = 31 50 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 +edge = 30 49 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 + +[track] +label = track 6 +anchor_a = 37 +anchor_b = 52 +axis = +y +expression = series(s, t, 360, 5) +subdivisions = 20 +members = 36 51 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 +edge = 37 52 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 +edge = 36 51 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 +