diff --git a/.github/dependabot.yml b/.github/dependabot.yml
index 700707c..1e8a051 100644
--- a/.github/dependabot.yml
+++ b/.github/dependabot.yml
@@ -5,3 +5,6 @@ updates:
     directory: "/" # Location of package manifests
     schedule:
       interval: "weekly"
+    ignore:
+      - dependency-name: "crate-ci/typos"
+        update-types: ["version-update:semver-patch"]
diff --git a/.github/workflows/SpellCheck.yml b/.github/workflows/SpellCheck.yml
new file mode 100644
index 0000000..599253c
--- /dev/null
+++ b/.github/workflows/SpellCheck.yml
@@ -0,0 +1,13 @@
+name: Spell Check
+
+on: [pull_request]
+
+jobs:
+  typos-check:
+    name: Spell Check with Typos
+    runs-on: ubuntu-latest
+    steps:
+      - name: Checkout Actions Repository
+        uses: actions/checkout@v3
+      - name: Check spelling
+        uses: crate-ci/typos@v1.16.23 
\ No newline at end of file
diff --git a/.typos.toml b/.typos.toml
new file mode 100644
index 0000000..1960516
--- /dev/null
+++ b/.typos.toml
@@ -0,0 +1,2 @@
+[default.extend-words]
+AAS = "AAS"
\ No newline at end of file
diff --git a/paper/paper.md b/paper/paper.md
index 673eece..ce631cc 100644
--- a/paper/paper.md
+++ b/paper/paper.md
@@ -31,7 +31,7 @@ aas-journal: Astrophysical Journal <- The name of the AAS journal.
 `HighDimPDE.jl` is a Julia [@Bezanson2017] package that implements solver algorithms to solve highly dimensional non-local non-linear Partial Differential Equations (PDEs). The solver algorithms provided break down the curse of dimensionality, with a computational complexity that only grows polynomially in the number of dimension of the PDE.  It is an open-source project hosted on GitHub and distributed under the MIT license. The package is designed with a user-friendly interface, provides both CPUs and GPUs support, and is integrated within the Sci-ML[@SciML] ecosystem.
 
 # Statement of need
-Non-local nonlinear Partial Differential Equations arise in a variety of scientific domains including physics, engineering, finance and biology. In biology, they are for instance used for modelling the evolution of biological populations that are phenotypically and physically structured. The dimension of the PDEs can be large, corresponding to the number of phenotypic traits and physical dimensions considered. Highly dimensional PDE's cannot be solved with standard numerical methods as their computational cost increases exponentially in the number of dimensions, a problem commonly refered as the curse of dimensionality. 
+Non-local nonlinear Partial Differential Equations arise in a variety of scientific domains including physics, engineering, finance and biology. In biology, they are for instance used for modelling the evolution of biological populations that are phenotypically and physically structured. The dimension of the PDEs can be large, corresponding to the number of phenotypic traits and physical dimensions considered. Highly dimensional PDE's cannot be solved with standard numerical methods as their computational cost increases exponentially in the number of dimensions, a problem commonly referred as the curse of dimensionality. 
 
 # Solving PDEs with HighDimPDE.jl
 HighDimPDE.jl can solve for PDEs of the form
@@ -42,7 +42,7 @@ $$
     & \quad + \big\langle \mu(t,x), ( \nabla_x u )( t,x ) \big\rangle + \tfrac{1}{2} \text{Trace} \big(\sigma(t,x) [ \sigma(t,x) ]^* ( \text{Hess}_x u)(t, x ) \big). \tag{1}
 \end{aligned}
 $$
-where $u \colon [0,T] \times \Omega \to \R$, $\Omega \subset \R^d$ is a function subject to initial conditions $u(0,x) = g(x)$ and Neumann Boudary conditions.
+where $u \colon [0,T] \times \Omega \to \R$, $\Omega \subset \R^d$ is a function subject to initial conditions $u(0,x) = g(x)$ and Neumann Boundary conditions.
 
 `HighDimPDE.jl` currently proposes ? solver algorithms.
 
diff --git a/src/DeepSplitting.jl b/src/DeepSplitting.jl
index 3af3936..bb554ad 100644
--- a/src/DeepSplitting.jl
+++ b/src/DeepSplitting.jl
@@ -189,7 +189,7 @@ function DiffEqBase.solve(prob::PIDEProblem,
 
                 # report on training
                 if epoch % verbose_rate == 1
-                    l = loss(y0, y1, z, t) # explictly computing loss every verbose_rate
+                    l = loss(y0, y1, z, t) # explicitly computing loss every verbose_rate
                     verbose && println("Current loss is: $l")
                     push!(losses[net], l)
                     if l < abstol
diff --git a/src/MCSample.jl b/src/MCSample.jl
index 94d078f..e385ff0 100644
--- a/src/MCSample.jl
+++ b/src/MCSample.jl
@@ -29,7 +29,7 @@ end
 Normal sampling method for the Monte Carlo integration.
 
 # Arguments
-* `σ`: the standard devation of the sampling
+* `σ`: the standard deviation of the sampling
 * `shifted` : if true, the integration is shifted by `x`. Defaults to false.
 """
 struct NormalSampling{T} <: MCSampling{T}