Addresses tolerance issues in boundary condition tests

Updates tests in BoundaryConditions.jl to use `isapprox` with a
specified tolerance (rtol=1e-14) instead of direct equality checks.
This resolves potential issues caused by floating-point precision
variations.

Author: jmartfrut

README.md

@@ -144,7 +144,8 @@ src="https://github.com/MultiSimOLab/HyperFEM/blob/main/docs/imgs/sims_.png?raw=
 
 In order to give credit to the HyperFEM contributors, we ask that you please reference the paper:
 
-> In process
+C. Perez‐Garcia, R. Ortigosa, J. Martínez‐Frutos, and D. Garcia‐Gonzalez, **Topology and material optimization in ultra-soft magnetoactive structures: making advantage of residual anisotropies.** Adv. Mater. (2025): e18489. https://https://doi.org/10.1002/adma.202518489
+ 
 
 along with the required citations for [Gridap](https://github.com/gridap/Gridap.jl).
 

test/SimulationsTests/BoundaryConditions.jl

@@ -30,13 +30,13 @@ using Gridap.FESpaces
   Vu = TestFESpace(Ω, reffeu, Du, conformity=:H1)
   Uu = TrialFESpace(Vu, Du, 1.0)
 
-  @test norm(Uu.dirichlet_values) == 0.30990321069650995
-
+@test isapprox(norm(Uu.dirichlet_values) ,  0.30990321069650995,rtol=1e-14)
   TrialFESpace!(Uu, Du, 0.4)
-  @test norm(Uu.dirichlet_values) == 0.12396128427860398
+@test isapprox(norm(Uu.dirichlet_values) ,  0.12396128427860398,rtol=1e-14)
   TrialFESpace!(Uu, Du, 1.0)
-  @test norm(Uu.dirichlet_values) == 0.30990321069650995
+@test isapprox(norm(Uu.dirichlet_values) ,  0.30990321069650995,rtol=1e-14)
 
+ 
 end
 
 
@@ -138,6 +138,7 @@ end
     interpolate_everywhere!(x -> Mapping(x, Δβ * (t)), get_free_dof_values(uh⁻), uh⁻.dirichlet_values, Vu_⁻)
     solve!(comp_model_vacmech; args_vacmech...)
   end
-  @test norm(xh.free_values) == 2.8132015601158087
+ @test isapprox( norm(xh.free_values) ,   2.8132015601158087,rtol=1e-14)
+
 
 end