Merge branch 'gradient_and_user_guides' of https://github.com/tmathis…

…720/HYDRA

Cargo.toml

@@ -14,6 +14,7 @@ rayon = "1.5"
 crossbeam = "0.8.4"
 dashmap = "5.1"
 approx = "0.5"
+lazy_static = "1.5"
 
 [dev-dependencies]
 criterion = "0.4"

src/boundary/bc_handler.rs

@@ -1,5 +1,6 @@
 use dashmap::DashMap;
-use std::sync::Arc;
+use std::sync::{Arc, RwLock};
+use lazy_static::lazy_static;
 use crate::domain::mesh_entity::MeshEntity;
 use crate::boundary::dirichlet::DirichletBC;
 use crate::boundary::neumann::NeumannBC;
@@ -29,6 +30,11 @@ pub struct BoundaryConditionHandler {
     conditions: DashMap<MeshEntity, BoundaryCondition>,
 }
 
+lazy_static! {
+    static ref GLOBAL_BC_HANDLER: Arc<RwLock<BoundaryConditionHandler>> =
+        Arc::new(RwLock::new(BoundaryConditionHandler::new()));
+}
+
 impl BoundaryConditionHandler {
     /// Creates a new BoundaryConditionHandler with an empty map to store boundary conditions.
     pub fn new() -> Self {
@@ -37,6 +43,10 @@ impl BoundaryConditionHandler {
         }
     }
 
+    pub fn global() -> Arc<RwLock<BoundaryConditionHandler>> {
+        GLOBAL_BC_HANDLER.clone()
+    }
+
     /// Sets a boundary condition for a specific mesh entity.
     pub fn set_bc(&self, entity: MeshEntity, condition: BoundaryCondition) {
         self.conditions.insert(entity, condition);
@@ -105,6 +115,7 @@ impl BoundaryConditionHandler {
     }
 }
 
+
 /// The BoundaryConditionApply trait defines the `apply` method, which is used to apply 
 /// a boundary condition to a given mesh entity.
 pub trait BoundaryConditionApply {

src/domain/mesh/geometry.rs

@@ -1,7 +1,9 @@
 use super::Mesh;
+use crate::domain;
 use crate::domain::mesh_entity::MeshEntity;
 use crate::geometry::{Geometry, CellShape, FaceShape};
 use dashmap::DashMap;
+use crate::domain::section::Vector3;
 
 impl Mesh {
     /// Retrieves all the faces of a given cell, filtering only face entities.
@@ -152,7 +154,7 @@ impl Mesh {
         &self,
         face: &MeshEntity,
         reference_cell: Option<&MeshEntity>,
-    ) -> Option<[f64; 3]> {
+    ) -> Option<Vector3> {
         // Retrieve face vertices
         let face_vertices = self.get_face_vertices(face);
         let face_shape = match face_vertices.len() {
@@ -186,10 +188,10 @@ impl Mesh {
 
             if dot_product < 0.0 {
                 // Reverse the normal direction to make it outward-pointing
-                return Some([-normal[0], -normal[1], -normal[2]]);
+                return Some(domain::section::Vector3([-normal[0], -normal[1], -normal[2]]));
             }
         }
 
-        Some(normal)
+        Some(domain::section::Vector3(normal))
     }
 }

src/domain/mesh/mod.rs

@@ -12,6 +12,7 @@ use crate::domain::sieve::Sieve;
 use rustc_hash::{FxHashMap, FxHashSet};
 use std::sync::{Arc, RwLock};
 use crossbeam::channel::{Sender, Receiver};
+use lazy_static::lazy_static;
 
 // Delegate methods to corresponding modules
 
@@ -43,6 +44,10 @@ pub struct Mesh {
     pub boundary_data_receiver: Option<Receiver<FxHashMap<MeshEntity, [f64; 3]>>>,  
 }
 
+lazy_static! {
+    static ref GLOBAL_MESH: Arc<RwLock<Mesh>> = Arc::new(RwLock::new(Mesh::new()));
+}
+
 impl Mesh {
     /// Creates a new instance of the `Mesh` struct with initialized components.  
     /// 
@@ -67,6 +72,10 @@ impl Mesh {
             boundary_data_receiver: Some(receiver),
         }
     }
+
+    pub fn global() -> Arc<RwLock<Mesh>> {
+        GLOBAL_MESH.clone()
+    }
 }
 
 #[cfg(test)]