diff --git a/tests/cases/cep/cylinder_purkinje_1d3d/README.md b/tests/cases/cep/cylinder_purkinje_1d3d/README.md
new file mode 100755
index 000000000..a8d3d7722
--- /dev/null
+++ b/tests/cases/cep/cylinder_purkinje_1d3d/README.md
@@ -0,0 +1,19 @@
+
+# **Problem Description**
+
+Simulate the propagation of an electric signal in a cylinder coupled with a network of Purkinje fibers using the ten-Tusscher-Panfilov cell activation model.
+
+Note: Conductivities in the myocardium were selected so that the propagation occurs in the myocardium (if too slow, due to the coarseness of the mesh, the action potential does not propagate beyond the coupled nodes).
+
+<p align="center">
+   <img src="./activation_10fps.gif" width="600">
+</p>
+
+## References
+K. H. W. J. ten Tusscher, D. Noble, P. J. Noble, and A. V. Panfilov. A model for human ventricular tissue. American Journal of Physiology-Heart and Circulatory Physiology,
+286(4):H1573–H1589, apr 2004.
+
+K. H. W. J. ten Tusscher and A. V. Panfilov. Alternans and spiral breakup in a human
+ventricular tissue model. American Journal of Physiology-Heart and Circulatory Physiology,
+291(3):H1088–H1100, sep 2006.
+
diff --git a/tests/cases/cep/cylinder_purkinje_1d3d/activation_10fps.gif b/tests/cases/cep/cylinder_purkinje_1d3d/activation_10fps.gif
new file mode 100644
index 000000000..bdf3734f5
--- /dev/null
+++ b/tests/cases/cep/cylinder_purkinje_1d3d/activation_10fps.gif
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:5fce1e5d81a12dd0e8935b79b61c126299f310259af6422854620b471e8014a8
+size 5013173
diff --git a/tests/cases/cep/cylinder_purkinje_1d3d/mesh/cyl-tree.vtu b/tests/cases/cep/cylinder_purkinje_1d3d/mesh/cyl-tree.vtu
new file mode 100644
index 000000000..a61afb426
--- /dev/null
+++ b/tests/cases/cep/cylinder_purkinje_1d3d/mesh/cyl-tree.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:c246d295275410c214d19ab9f96e813530aaa783ce138cf55dc72bf4afd98b8a
+size 36058
diff --git a/tests/cases/cep/cylinder_purkinje_1d3d/mesh/cyl-tree_domain.dat b/tests/cases/cep/cylinder_purkinje_1d3d/mesh/cyl-tree_domain.dat
new file mode 100644
index 000000000..75a62027e
--- /dev/null
+++ b/tests/cases/cep/cylinder_purkinje_1d3d/mesh/cyl-tree_domain.dat
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:6f0ee50981ad21964c089db40ba5ab3099c018e653dbc36215ca31c0cab3ca4d
+size 890
diff --git a/tests/cases/cep/cylinder_purkinje_1d3d/mesh/cyl-tree_endnodes.txt b/tests/cases/cep/cylinder_purkinje_1d3d/mesh/cyl-tree_endnodes.txt
new file mode 100644
index 000000000..4c3fbce73
--- /dev/null
+++ b/tests/cases/cep/cylinder_purkinje_1d3d/mesh/cyl-tree_endnodes.txt
@@ -0,0 +1,81 @@
+1
+16
+54
+62
+66
+74
+84
+85
+91
+99
+101
+122
+123
+124
+128
+129
+137
+141
+149
+151
+164
+167
+175
+178
+203
+227
+231
+239
+252
+268
+283
+285
+287
+291
+294
+298
+302
+307
+323
+326
+327
+328
+329
+333
+334
+338
+342
+346
+350
+351
+355
+363
+364
+368
+369
+373
+376
+378
+379
+380
+381
+384
+392
+396
+400
+402
+409
+411
+414
+415
+419
+426
+430
+432
+436
+438
+441
+442
+443
+445
+446
diff --git a/tests/cases/cep/cylinder_purkinje_1d3d/mesh/cylinder.vtu b/tests/cases/cep/cylinder_purkinje_1d3d/mesh/cylinder.vtu
new file mode 100644
index 000000000..0297e2499
--- /dev/null
+++ b/tests/cases/cep/cylinder_purkinje_1d3d/mesh/cylinder.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:9195c6ef95a6490b68c40561ee6ec66cd7ac96cd893103e09df8152a7ff7dec1
+size 67480
diff --git a/tests/cases/cep/cylinder_purkinje_1d3d/mesh/cylinder_domain.dat b/tests/cases/cep/cylinder_purkinje_1d3d/mesh/cylinder_domain.dat
new file mode 100644
index 000000000..192f5638b
--- /dev/null
+++ b/tests/cases/cep/cylinder_purkinje_1d3d/mesh/cylinder_domain.dat
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:306178540daa4fb17b2f3abd8059edb4f776a8e5f9a0ffbc3caa3548653dfa32
+size 4560
diff --git a/tests/cases/cep/cylinder_purkinje_1d3d/mesh/cylinder_fiber.vtu b/tests/cases/cep/cylinder_purkinje_1d3d/mesh/cylinder_fiber.vtu
new file mode 100644
index 000000000..e9d9a5968
--- /dev/null
+++ b/tests/cases/cep/cylinder_purkinje_1d3d/mesh/cylinder_fiber.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:248289d115368e3ae26c3661e86b455aabf9a1f886826f8ead4dff3a45d8d3aa
+size 71037
diff --git a/tests/cases/cep/cylinder_purkinje_1d3d/mesh/mesh-surfaces/cylinder_inner.vtp b/tests/cases/cep/cylinder_purkinje_1d3d/mesh/mesh-surfaces/cylinder_inner.vtp
new file mode 100644
index 000000000..7ebc3f425
--- /dev/null
+++ b/tests/cases/cep/cylinder_purkinje_1d3d/mesh/mesh-surfaces/cylinder_inner.vtp
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:b9fc93450da649483d18dd4ce041ee41d646514e53696f3cee36ffd237c84911
+size 22703
diff --git a/tests/cases/cep/cylinder_purkinje_1d3d/result_001.vtu b/tests/cases/cep/cylinder_purkinje_1d3d/result_001.vtu
new file mode 100644
index 000000000..5fdb86311
--- /dev/null
+++ b/tests/cases/cep/cylinder_purkinje_1d3d/result_001.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:23aec8b3fc9a0e69a44f83990e98d8b50e8952b387fc2991098b70440d351caa
+size 91659
diff --git a/tests/cases/cep/cylinder_purkinje_1d3d/solver.xml b/tests/cases/cep/cylinder_purkinje_1d3d/solver.xml
new file mode 100644
index 000000000..35844bcb0
--- /dev/null
+++ b/tests/cases/cep/cylinder_purkinje_1d3d/solver.xml
@@ -0,0 +1,107 @@
+<?xml version="1.0" encoding="UTF-8" ?>
+<svMultiPhysicsFile version="0.1">
+
+<GeneralSimulationParameters>
+  <Continue_previous_simulation> false </Continue_previous_simulation>
+  <Number_of_spatial_dimensions> 3 </Number_of_spatial_dimensions> 
+  <Number_of_time_steps> 1 </Number_of_time_steps> 
+  <Time_step_size> 0.1 </Time_step_size> 
+  <Spectral_radius_of_infinite_time_step> 0.50 </Spectral_radius_of_infinite_time_step> 
+  <Searched_file_name_to_trigger_stop> STOP_SIM </Searched_file_name_to_trigger_stop> 
+
+  <Save_results_to_VTK_format> true </Save_results_to_VTK_format> 
+  <Name_prefix_of_saved_VTK_files> result </Name_prefix_of_saved_VTK_files> 
+  <Increment_in_saving_VTK_files> 1 </Increment_in_saving_VTK_files> 
+  <Start_saving_after_time_step> 1 </Start_saving_after_time_step> 
+
+  <Increment_in_saving_restart_files> 1 </Increment_in_saving_restart_files> 
+  <Convert_BIN_to_VTK_format> 0 </Convert_BIN_to_VTK_format> 
+
+  <Verbose> 1 </Verbose> 
+  <Warning> 0 </Warning> 
+  <Debug> 0 </Debug> 
+</GeneralSimulationParameters>
+
+<Add_mesh name="efib" > 
+  <Set_mesh_as_fibers> true </Set_mesh_as_fibers> 
+  <Mesh_file_path> ./mesh/cyl-tree.vtu  </Mesh_file_path>
+  <Add_face name = "fc_efib">
+      <End_nodes_face_file_path> ./mesh/cyl-tree_endnodes.txt </End_nodes_face_file_path>
+  </Add_face>
+  <Domain_file_path> ./mesh/cyl-tree_domain.dat </Domain_file_path>
+  <Mesh_scale_factor> 10 </Mesh_scale_factor> 
+</Add_mesh>
+
+<Add_mesh name="msh" > 
+
+  <Mesh_file_path> ./mesh/cylinder.vtu   </Mesh_file_path>
+
+  <Add_face name="endo">
+      <Face_file_path> ./mesh/mesh-surfaces/cylinder_inner.vtp </Face_file_path>
+  </Add_face>
+
+  <Domain_file_path> ./mesh/cylinder_domain.dat </Domain_file_path>
+
+   <Fiber_direction_file_path> ./mesh/cylinder_fiber.vtu </Fiber_direction_file_path>
+   <Mesh_scale_factor> 10 </Mesh_scale_factor>
+
+</Add_mesh>
+
+<Add_projection name = "fc_efib">
+   <Project_from_face> endo </Project_from_face>
+   <Projection_tolerance> -1.0 </Projection_tolerance>  <!-- a negative value means using the nearest neighbor, no matter how far it is -->
+</Add_projection>
+
+
+<Add_equation type="CEP" > 
+   <Coupled> true </Coupled>
+   <Min_iterations> 1 </Min_iterations>  
+   <Max_iterations> 3 </Max_iterations> 
+   <Tolerance> 1e-12 </Tolerance> 
+  
+   <Domain id="1" >
+      <Electrophysiology_model> TTP </Electrophysiology_model>
+      <Myocardial_zone> myocardium </Myocardial_zone> 
+      <Isotropic_conductivity> 0.05 </Isotropic_conductivity>
+      <Anisotropic_conductivity> 0.05 </Anisotropic_conductivity>
+      <ODE_solver> RK </ODE_solver>
+   </Domain>
+
+   <Domain id="2" >
+      <Electrophysiology_model> TTP </Electrophysiology_model>
+      <Myocardial_zone> pfib </Myocardial_zone> 
+      <Isotropic_conductivity> 1.1 </Isotropic_conductivity>
+      <ODE_solver> RK </ODE_solver>
+   </Domain> 
+
+   <Domain id="3" >
+      <Electrophysiology_model> TTP </Electrophysiology_model>
+      <Myocardial_zone> pfib </Myocardial_zone> 
+      <Isotropic_conductivity> 1.1 </Isotropic_conductivity>
+      <ODE_solver> RK </ODE_solver>
+      <Stimulus type="Istim" >
+         <Amplitude> -52.0 </Amplitude>
+         <Start_time> 0.0 </Start_time>
+         <Duration> 2.0 </Duration>
+         <Cycle_length> 1000.0 </Cycle_length>
+      </Stimulus>
+   </Domain>
+
+   <Output type="Spatial" >
+      <Action_potential> true </Action_potential>
+   </Output>
+
+   <LS type="GMRES" >
+      <Linear_algebra type="fsils" >
+         <Preconditioner> fsils </Preconditioner>
+      </Linear_algebra>
+      <Max_iterations> 100 </Max_iterations>
+      <Tolerance> 1e-12 </Tolerance>
+      <Krylov_space_dimension> 50 </Krylov_space_dimension>
+   </LS>
+
+</Add_equation>
+
+</svMultiPhysicsFile>
+
+
diff --git a/tests/test_cep.py b/tests/test_cep.py
index 57c318460..6195019b8 100644
--- a/tests/test_cep.py
+++ b/tests/test_cep.py
@@ -32,6 +32,10 @@ def test_purkinje(n_proc):
     run_with_reference(base_folder, test_folder, fields, n_proc)
 
 
+def test_cylinder_purkinje_1d3d(n_proc):
+    test_folder = "cylinder_purkinje_1d3d"
+    run_with_reference(base_folder, test_folder, fields, n_proc)
+
 @pytest.mark.parametrize(
     "confs_ecgs",
     [