diff --git a/src/geometry/geometry_map.rs b/src/geometry/geometry_map.rs
index 7aa1f77..06ef6db 100644
--- a/src/geometry/geometry_map.rs
+++ b/src/geometry/geometry_map.rs
@@ -1,5 +1,5 @@
 //! Geometry map
-use crate::{traits::GeometryMap as GeometryMapTrait, types::RealScalar};
+use crate::{traits::GeometryMap as GeometryMapTrait, types::Scalar};
 use ndelement::{reference_cell, traits::MappedFiniteElement, types::ReferenceCellType};
 use rlst::{Array, DynArray, RlstScalar, ValueArrayImpl};
 
@@ -43,7 +43,7 @@ fn cross<T: RlstScalar>(mat: &[T], result: &mut [T]) {
     }
 }
 
-impl<'a, T: RealScalar, B2D: ValueArrayImpl<T, 2>, C2D: ValueArrayImpl<usize, 2>>
+impl<'a, T: Scalar, B2D: ValueArrayImpl<T, 2>, C2D: ValueArrayImpl<usize, 2>>
     GeometryMap<'a, T, B2D, C2D>
 {
     /// Create new
@@ -71,7 +71,7 @@ impl<'a, T: RealScalar, B2D: ValueArrayImpl<T, 2>, C2D: ValueArrayImpl<usize, 2>
     }
 }
 
-impl<T: RealScalar, B2D: ValueArrayImpl<T, 2>, C2D: ValueArrayImpl<usize, 2>> GeometryMapTrait
+impl<T: Scalar, B2D: ValueArrayImpl<T, 2>, C2D: ValueArrayImpl<usize, 2>> GeometryMapTrait
     for GeometryMap<'_, T, B2D, C2D>
 {
     type T = T;
diff --git a/src/geometry/mixed/entity_geometry.rs b/src/geometry/mixed/entity_geometry.rs
index cd0e6ed..c856f80 100644
--- a/src/geometry/mixed/entity_geometry.rs
+++ b/src/geometry/mixed/entity_geometry.rs
@@ -3,13 +3,13 @@ use super::MixedGeometry;
 use crate::{
     geometry::{Point, PointIter},
     traits::Geometry,
-    types::RealScalar,
+    types::Scalar,
 };
 use ndelement::traits::MappedFiniteElement;
 
 /// Geometry of an entity
 #[derive(Debug)]
-pub struct MixedEntityGeometry<'a, T: RealScalar, E: MappedFiniteElement> {
+pub struct MixedEntityGeometry<'a, T: Scalar, E: MappedFiniteElement> {
     geometry: &'a MixedGeometry<T, E>,
     element_index: usize,
     cell_index: usize,
@@ -17,7 +17,7 @@ pub struct MixedEntityGeometry<'a, T: RealScalar, E: MappedFiniteElement> {
     sub_entity_index: usize,
 }
 
-impl<'a, T: RealScalar, E: MappedFiniteElement> MixedEntityGeometry<'a, T, E> {
+impl<'a, T: Scalar, E: MappedFiniteElement> MixedEntityGeometry<'a, T, E> {
     /// Create new
     pub fn new(
         geometry: &'a MixedGeometry<T, E>,
@@ -36,7 +36,7 @@ impl<'a, T: RealScalar, E: MappedFiniteElement> MixedEntityGeometry<'a, T, E> {
     }
 }
 
-impl<T: RealScalar, E: MappedFiniteElement> Geometry for MixedEntityGeometry<'_, T, E> {
+impl<T: Scalar, E: MappedFiniteElement> Geometry for MixedEntityGeometry<'_, T, E> {
     type T = T;
     type Point<'a>
         = Point<'a, T>
diff --git a/src/geometry/mixed/geometry.rs b/src/geometry/mixed/geometry.rs
index d349ce9..0702646 100644
--- a/src/geometry/mixed/geometry.rs
+++ b/src/geometry/mixed/geometry.rs
@@ -1,7 +1,7 @@
 //! Geometry where each entity of a given dimension is represented by the same element
 #[cfg(feature = "serde")]
 use crate::traits::ConvertToSerializable;
-use crate::types::RealScalar;
+use crate::types::Scalar;
 use itertools::izip;
 #[cfg(feature = "serde")]
 use ndelement::{
@@ -19,7 +19,7 @@ use std::collections::HashMap;
 use std::fmt::{Debug, Formatter};
 
 /// Single element geometry
-pub struct MixedGeometry<T: RealScalar, E: MappedFiniteElement> {
+pub struct MixedGeometry<T: Scalar, E: MappedFiniteElement> {
     points: DynArray<T, 2>,
     cells: Vec<DynArray<usize, 2>>,
     elements: Vec<HashMap<ReferenceCellType, E>>,
@@ -27,7 +27,7 @@ pub struct MixedGeometry<T: RealScalar, E: MappedFiniteElement> {
     pub(crate) insertion_indices_to_cell_indices: Vec<usize>,
 }
 
-impl<T: RealScalar, E: MappedFiniteElement> Debug for MixedGeometry<T, E> {
+impl<T: Scalar, E: MappedFiniteElement> Debug for MixedGeometry<T, E> {
     fn fmt(&self, f: &mut Formatter<'_>) -> Result<(), std::fmt::Error> {
         f.debug_struct("MixedGeometry")
             .field("points", &self.points)
@@ -39,7 +39,7 @@ impl<T: RealScalar, E: MappedFiniteElement> Debug for MixedGeometry<T, E> {
 #[cfg(feature = "serde")]
 #[derive(serde::Serialize, Debug, serde::Deserialize)]
 #[serde(bound = "for<'de2> T: serde::Deserialize<'de2>")]
-pub struct SerializableGeometry<T: RealScalar + serde::Serialize>
+pub struct SerializableGeometry<T: Scalar + serde::Serialize>
 where
     for<'de2> T: serde::Deserialize<'de2>,
 {
@@ -51,8 +51,8 @@ where
 }
 
 #[cfg(feature = "serde")]
-impl<T: RealScalar + serde::Serialize> ConvertToSerializable
-    for MixedGeometry<T, CiarletElement<T, IdentityMap>>
+impl<T: Scalar + serde::Serialize> ConvertToSerializable
+    for MixedGeometry<T, CiarletElement<T, IdentityMap, T>>
 where
     for<'de2> T: serde::Deserialize<'de2>,
 {
@@ -115,7 +115,7 @@ where
     }
 }
 
-impl<T: RealScalar, E: MappedFiniteElement> MixedGeometry<T, E> {
+impl<T: Scalar, E: MappedFiniteElement> MixedGeometry<T, E> {
     /// Create single element geometry
     pub fn new(
         cell_types_in: &[ReferenceCellType],
diff --git a/src/geometry/point.rs b/src/geometry/point.rs
index cf10589..9ef2ad1 100644
--- a/src/geometry/point.rs
+++ b/src/geometry/point.rs
@@ -1,20 +1,20 @@
 //! Point
-use crate::{traits::Point as PointTrait, types::RealScalar};
+use crate::{traits::Point as PointTrait, types::Scalar};
 
 /// A points
 #[derive(Debug, Clone, Copy)]
-pub struct Point<'a, T: RealScalar> {
+pub struct Point<'a, T: Scalar> {
     index: usize,
     coordinates: &'a [T],
 }
 
-impl<'a, T: RealScalar> Point<'a, T> {
+impl<'a, T: Scalar> Point<'a, T> {
     /// Create new
     pub fn new(index: usize, coordinates: &'a [T]) -> Self {
         Self { index, coordinates }
     }
 }
-impl<T: RealScalar> PointTrait for Point<'_, T> {
+impl<T: Scalar> PointTrait for Point<'_, T> {
     type T = T;
 
     fn index(&self) -> usize {
@@ -32,17 +32,17 @@ impl<T: RealScalar> PointTrait for Point<'_, T> {
 
 /// Iterator over points
 #[derive(Debug)]
-pub struct PointIter<'a, T: RealScalar> {
+pub struct PointIter<'a, T: Scalar> {
     points: Vec<(usize, &'a [T])>,
     index: usize,
 }
-impl<'a, T: RealScalar> PointIter<'a, T> {
+impl<'a, T: Scalar> PointIter<'a, T> {
     /// Create new
     pub fn new(points: Vec<(usize, &'a [T])>) -> Self {
         Self { points, index: 0 }
     }
 }
-impl<'a, T: RealScalar> Iterator for PointIter<'a, T> {
+impl<'a, T: Scalar> Iterator for PointIter<'a, T> {
     type Item = Point<'a, T>;
 
     fn next(&mut self) -> Option<Point<'a, T>> {
diff --git a/src/geometry/single_element/entity_geometry.rs b/src/geometry/single_element/entity_geometry.rs
index 4bcfca5..2438d84 100644
--- a/src/geometry/single_element/entity_geometry.rs
+++ b/src/geometry/single_element/entity_geometry.rs
@@ -3,20 +3,20 @@ use super::SingleElementGeometry;
 use crate::{
     geometry::{Point, PointIter},
     traits::Geometry,
-    types::RealScalar,
+    types::Scalar,
 };
 use ndelement::traits::MappedFiniteElement;
 
 /// Geometry of an entity
 #[derive(Debug)]
-pub struct SingleElementEntityGeometry<'a, T: RealScalar, E: MappedFiniteElement> {
+pub struct SingleElementEntityGeometry<'a, T: Scalar, E: MappedFiniteElement> {
     geometry: &'a SingleElementGeometry<T, E>,
     cell_index: usize,
     sub_entity_dimension: usize,
     sub_entity_index: usize,
 }
 
-impl<'a, T: RealScalar, E: MappedFiniteElement> SingleElementEntityGeometry<'a, T, E> {
+impl<'a, T: Scalar, E: MappedFiniteElement> SingleElementEntityGeometry<'a, T, E> {
     /// Create new
     pub fn new(
         geometry: &'a SingleElementGeometry<T, E>,
@@ -33,7 +33,7 @@ impl<'a, T: RealScalar, E: MappedFiniteElement> SingleElementEntityGeometry<'a,
     }
 }
 
-impl<T: RealScalar, E: MappedFiniteElement> Geometry for SingleElementEntityGeometry<'_, T, E> {
+impl<T: Scalar, E: MappedFiniteElement> Geometry for SingleElementEntityGeometry<'_, T, E> {
     type T = T;
     type Point<'a>
         = Point<'a, T>
diff --git a/src/geometry/single_element/geometry.rs b/src/geometry/single_element/geometry.rs
index 72c7ddf..f8fad06 100644
--- a/src/geometry/single_element/geometry.rs
+++ b/src/geometry/single_element/geometry.rs
@@ -1,7 +1,7 @@
 //! Geometry where each entity of a given dimension is represented by the same element
 #[cfg(feature = "serde")]
 use crate::traits::ConvertToSerializable;
-use crate::types::RealScalar;
+use crate::types::Scalar;
 #[cfg(feature = "serde")]
 use ndelement::{
     ciarlet::{CiarletElement, lagrange},
@@ -18,13 +18,13 @@ use rlst::{DynArray, rlst_dynamic_array};
 use std::fmt::{Debug, Formatter};
 
 /// Single element geometry
-pub struct SingleElementGeometry<T: RealScalar, E: MappedFiniteElement> {
+pub struct SingleElementGeometry<T: Scalar, E: MappedFiniteElement> {
     points: DynArray<T, 2>,
     cells: DynArray<usize, 2>,
     elements: Vec<E>,
 }
 
-impl<T: RealScalar, E: MappedFiniteElement> Debug for SingleElementGeometry<T, E> {
+impl<T: Scalar, E: MappedFiniteElement> Debug for SingleElementGeometry<T, E> {
     fn fmt(&self, f: &mut Formatter<'_>) -> Result<(), std::fmt::Error> {
         f.debug_struct("SingleElementGeometry")
             .field("points", &self.points)
@@ -36,7 +36,7 @@ impl<T: RealScalar, E: MappedFiniteElement> Debug for SingleElementGeometry<T, E
 #[cfg(feature = "serde")]
 #[derive(serde::Serialize, Debug, serde::Deserialize)]
 #[serde(bound = "for<'de2> T: serde::Deserialize<'de2>")]
-pub struct SerializableGeometry<T: RealScalar + serde::Serialize>
+pub struct SerializableGeometry<T: Scalar + serde::Serialize>
 where
     for<'de2> T: serde::Deserialize<'de2>,
 {
@@ -46,8 +46,8 @@ where
 }
 
 #[cfg(feature = "serde")]
-impl<T: RealScalar + serde::Serialize> ConvertToSerializable
-    for SingleElementGeometry<T, CiarletElement<T, IdentityMap>>
+impl<T: Scalar + serde::Serialize> ConvertToSerializable
+    for SingleElementGeometry<T, CiarletElement<T, IdentityMap, T>>
 where
     for<'de2> T: serde::Deserialize<'de2>,
 {
@@ -86,7 +86,7 @@ where
     }
 }
 
-impl<T: RealScalar, E: MappedFiniteElement> SingleElementGeometry<T, E> {
+impl<T: Scalar, E: MappedFiniteElement> SingleElementGeometry<T, E> {
     /// Create single element geometry
     pub fn new(
         cell_type: ReferenceCellType,
diff --git a/src/grid/local_grid/mixed/builder.rs b/src/grid/local_grid/mixed/builder.rs
index 76ccf58..14651fa 100644
--- a/src/grid/local_grid/mixed/builder.rs
+++ b/src/grid/local_grid/mixed/builder.rs
@@ -5,7 +5,7 @@ use crate::{
     geometry::MixedGeometry,
     topology::mixed::MixedTopology,
     traits::{Builder, GeometryBuilder, GridBuilder, TopologyBuilder},
-    types::RealScalar,
+    types::Scalar,
 };
 use itertools::izip;
 use ndelement::{
@@ -20,10 +20,10 @@ use std::collections::{HashMap, HashSet};
 
 /// Grid builder for a grid with a mixture of element types
 #[derive(Debug)]
-pub struct MixedGridBuilder<T: RealScalar> {
+pub struct MixedGridBuilder<T: Scalar> {
     gdim: usize,
     element_indices: HashMap<(ReferenceCellType, usize), usize>,
-    elements: Vec<CiarletElement<T, IdentityMap>>,
+    elements: Vec<CiarletElement<T, IdentityMap, T>>,
     points_per_cell: Vec<usize>,
     pub(crate) points: Vec<T>,
     cells: Vec<usize>,
@@ -39,7 +39,7 @@ pub struct MixedGridBuilder<T: RealScalar> {
     cell_indices: HashSet<usize>,
 }
 
-impl<T: RealScalar> MixedGridBuilder<T> {
+impl<T: Scalar> MixedGridBuilder<T> {
     /// Create a new grid builder
     pub fn new(gdim: usize) -> Self {
         Self {
@@ -63,8 +63,8 @@ impl<T: RealScalar> MixedGridBuilder<T> {
     }
 }
 
-impl<T: RealScalar> Builder for MixedGridBuilder<T> {
-    type Grid = MixedGrid<T, CiarletElement<T, IdentityMap>>;
+impl<T: Scalar> Builder for MixedGridBuilder<T> {
+    type Grid = MixedGrid<T, CiarletElement<T, IdentityMap, T>>;
     type T = T;
     type CellData<'a> = (ReferenceCellType, usize, &'a [usize]);
     type EntityDescriptor = ReferenceCellType;
@@ -89,7 +89,7 @@ impl<T: RealScalar> Builder for MixedGridBuilder<T> {
             .entry((cell_data.0, cell_data.1))
             .or_insert_with(|| {
                 let n = self.cell_indices_to_ids.len();
-                self.elements.push(lagrange::create::<T>(
+                self.elements.push(lagrange::create::<T, T>(
                     cell_data.0,
                     cell_data.1,
                     Continuity::Standard,
@@ -127,7 +127,7 @@ impl<T: RealScalar> Builder for MixedGridBuilder<T> {
         self.add_cell(id, (cell_type, cell_degree, nodes));
     }
 
-    fn create_grid(&self) -> MixedGrid<T, CiarletElement<T, IdentityMap>> {
+    fn create_grid(&self) -> MixedGrid<T, CiarletElement<T, IdentityMap, T>> {
         let cell_vertices =
             self.extract_vertices(&self.cells, &self.cell_types, &self.cell_degrees);
 
@@ -209,8 +209,8 @@ impl<T: RealScalar> Builder for MixedGridBuilder<T> {
     }
 }
 
-impl<T: RealScalar> GeometryBuilder for MixedGridBuilder<T> {
-    type GridGeometry = MixedGeometry<T, CiarletElement<T, IdentityMap>>;
+impl<T: Scalar> GeometryBuilder for MixedGridBuilder<T> {
+    type GridGeometry = MixedGeometry<T, CiarletElement<T, IdentityMap, T>>;
     fn create_geometry(
         &self,
         point_ids: &[usize],
@@ -218,7 +218,7 @@ impl<T: RealScalar> GeometryBuilder for MixedGridBuilder<T> {
         cell_points: &[usize],
         cell_types: &[ReferenceCellType],
         cell_degrees: &[usize],
-    ) -> MixedGeometry<T, CiarletElement<T, IdentityMap>> {
+    ) -> MixedGeometry<T, CiarletElement<T, IdentityMap, T>> {
         let npts = point_ids.len();
         let mut points = rlst_dynamic_array!(T, [self.gdim(), npts]);
         points.data_mut().unwrap().copy_from_slice(coordinates);
@@ -229,7 +229,7 @@ impl<T: RealScalar> GeometryBuilder for MixedGridBuilder<T> {
         for degree in cell_degrees {
             cell_families.push(*ef_indices.entry(*degree).or_insert_with(|| {
                 let n = element_families.len();
-                element_families.push(LagrangeElementFamily::<T>::new(
+                element_families.push(LagrangeElementFamily::<T, T>::new(
                     *degree,
                     Continuity::Standard,
                 ));
@@ -237,7 +237,7 @@ impl<T: RealScalar> GeometryBuilder for MixedGridBuilder<T> {
             }))
         }
 
-        MixedGeometry::<T, CiarletElement<T, IdentityMap>>::new(
+        MixedGeometry::<T, CiarletElement<T, IdentityMap, T>>::new(
             cell_types,
             points,
             cell_points,
@@ -247,7 +247,7 @@ impl<T: RealScalar> GeometryBuilder for MixedGridBuilder<T> {
     }
 }
 
-impl<T: RealScalar> TopologyBuilder for MixedGridBuilder<T> {
+impl<T: Scalar> TopologyBuilder for MixedGridBuilder<T> {
     type GridTopology = MixedTopology;
     fn create_topology(
         &self,
@@ -297,7 +297,7 @@ impl<T: RealScalar> TopologyBuilder for MixedGridBuilder<T> {
     }
 }
 
-impl<T: RealScalar> GridBuilder for MixedGridBuilder<T> {
+impl<T: Scalar> GridBuilder for MixedGridBuilder<T> {
     fn create_grid_from_topology_geometry(
         &self,
         topology: <Self as TopologyBuilder>::GridTopology,
diff --git a/src/grid/local_grid/mixed/grid.rs b/src/grid/local_grid/mixed/grid.rs
index 793c294..8cf4d66 100644
--- a/src/grid/local_grid/mixed/grid.rs
+++ b/src/grid/local_grid/mixed/grid.rs
@@ -16,7 +16,7 @@ use crate::{
     geometry::{GeometryMap, MixedEntityGeometry, MixedGeometry},
     topology::mixed::{MixedEntityTopology, MixedTopology},
     traits::{Entity, Grid},
-    types::{Ownership, RealScalar},
+    types::{Ownership, Scalar},
 };
 use itertools::izip;
 #[cfg(feature = "mpi")]
@@ -36,7 +36,7 @@ use std::collections::HashMap;
 #[derive(Debug)]
 pub struct MixedGridEntity<
     'a,
-    T: RealScalar,
+    T: Scalar,
     E: MappedFiniteElement<CellType = ReferenceCellType, T = T>,
 > {
     grid: &'a MixedGrid<T, E>,
@@ -48,7 +48,7 @@ pub struct MixedGridEntity<
     geometry_cell_index: usize,
 }
 
-impl<'e, T: RealScalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>>
+impl<'e, T: Scalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>>
     MixedGridEntity<'e, T, E>
 {
     /// Create new
@@ -72,7 +72,7 @@ impl<'e, T: RealScalar, E: MappedFiniteElement<CellType = ReferenceCellType, T =
         }
     }
 }
-impl<T: RealScalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>> Entity
+impl<T: Scalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>> Entity
     for MixedGridEntity<'_, T, E>
 {
     type T = T;
@@ -125,7 +125,7 @@ impl<T: RealScalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>>
 #[derive(Debug)]
 pub struct MixedGridEntityIter<
     'a,
-    T: RealScalar,
+    T: Scalar,
     E: MappedFiniteElement<CellType = ReferenceCellType, T = T>,
 > {
     grid: &'a MixedGrid<T, E>,
@@ -133,7 +133,7 @@ pub struct MixedGridEntityIter<
     index: usize,
 }
 
-impl<'a, T: RealScalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>>
+impl<'a, T: Scalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>>
     MixedGridEntityIter<'a, T, E>
 {
     /// Create new
@@ -145,7 +145,7 @@ impl<'a, T: RealScalar, E: MappedFiniteElement<CellType = ReferenceCellType, T =
         }
     }
 }
-impl<'a, T: RealScalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>> Iterator
+impl<'a, T: Scalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>> Iterator
     for MixedGridEntityIter<'a, T, E>
 {
     type Item = MixedGridEntity<'a, T, E>;
@@ -158,7 +158,7 @@ impl<'a, T: RealScalar, E: MappedFiniteElement<CellType = ReferenceCellType, T =
 
 /// Serial mixed element grid
 #[derive(Debug)]
-pub struct MixedGrid<T: RealScalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>> {
+pub struct MixedGrid<T: Scalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>> {
     topology: MixedTopology,
     geometry: MixedGeometry<T, E>,
 }
@@ -166,7 +166,7 @@ pub struct MixedGrid<T: RealScalar, E: MappedFiniteElement<CellType = ReferenceC
 #[cfg(feature = "serde")]
 #[derive(serde::Serialize, Debug, serde::Deserialize)]
 #[serde(bound = "for<'de2> T: serde::Deserialize<'de2>")]
-pub struct SerializableGrid<T: RealScalar + serde::Serialize>
+pub struct SerializableGrid<T: Scalar + serde::Serialize>
 where
     for<'de2> T: serde::Deserialize<'de2>,
 {
@@ -175,8 +175,8 @@ where
 }
 
 #[cfg(feature = "serde")]
-impl<T: RealScalar + serde::Serialize> ConvertToSerializable
-    for MixedGrid<T, CiarletElement<T, IdentityMap>>
+impl<T: Scalar + serde::Serialize> ConvertToSerializable
+    for MixedGrid<T, CiarletElement<T, IdentityMap, T>>
 {
     type SerializableType = SerializableGrid<T>;
     fn to_serializable(&self) -> SerializableGrid<T> {
@@ -193,14 +193,14 @@ impl<T: RealScalar + serde::Serialize> ConvertToSerializable
     }
 }
 
-impl<T: RealScalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>> MixedGrid<T, E> {
+impl<T: Scalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>> MixedGrid<T, E> {
     /// Create new
     pub fn new(topology: MixedTopology, geometry: MixedGeometry<T, E>) -> Self {
         Self { topology, geometry }
     }
 }
 
-impl<T: RealScalar> MixedGrid<T, CiarletElement<T, IdentityMap>> {
+impl<T: Scalar> MixedGrid<T, CiarletElement<T, IdentityMap, T>> {
     /// Create new from raw data
     pub fn new_from_raw_data(
         coordinates: &[T],
@@ -222,13 +222,13 @@ impl<T: RealScalar> MixedGrid<T, CiarletElement<T, IdentityMap>> {
                 *element_family_indices.entry(*d).or_insert_with(|| {
                     let index = element_families.len();
                     element_families
-                        .push(LagrangeElementFamily::<T>::new(*d, Continuity::Standard));
+                        .push(LagrangeElementFamily::<T, T>::new(*d, Continuity::Standard));
                     index
                 })
             })
             .collect::<Vec<_>>();
 
-        let geometry = MixedGeometry::<T, CiarletElement<T, IdentityMap>>::new(
+        let geometry = MixedGeometry::<T, CiarletElement<T, IdentityMap, T>>::new(
             cell_types,
             points,
             cells,
@@ -252,7 +252,7 @@ impl<T: RealScalar> MixedGrid<T, CiarletElement<T, IdentityMap>> {
     }
 }
 
-impl<T: RealScalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>> Grid
+impl<T: Scalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>> Grid
     for MixedGrid<T, E>
 {
     type T = T;
@@ -361,11 +361,11 @@ impl<T: RealScalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>>
 }
 
 #[cfg(feature = "mpi")]
-impl<T: RealScalar + Equivalence, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>>
+impl<T: Scalar + Equivalence, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>>
     DistributableGrid for MixedGrid<T, E>
 {
     type ParallelGrid<'a, C: Communicator + 'a> =
-        ParallelGridImpl<'a, C, MixedGrid<T, CiarletElement<T, IdentityMap>>>;
+        ParallelGridImpl<'a, C, MixedGrid<T, CiarletElement<T, IdentityMap, T>>>;
 
     fn distribute<'a, C: Communicator>(
         &self,
diff --git a/src/grid/local_grid/single_element/builder.rs b/src/grid/local_grid/single_element/builder.rs
index 6e74f3f..593efef 100644
--- a/src/grid/local_grid/single_element/builder.rs
+++ b/src/grid/local_grid/single_element/builder.rs
@@ -5,7 +5,7 @@ use crate::{
     geometry::SingleElementGeometry,
     topology::single_type::SingleTypeTopology,
     traits::{Builder, GeometryBuilder, GridBuilder, TopologyBuilder},
-    types::RealScalar,
+    types::Scalar,
 };
 use ndelement::{
     ciarlet::{CiarletElement, LagrangeElementFamily, lagrange},
@@ -41,10 +41,10 @@ use std::collections::{HashMap, HashSet};
 /// let grid = builder.create_grid();
 /// ```
 #[derive(Debug)]
-pub struct SingleElementGridBuilder<T: RealScalar> {
+pub struct SingleElementGridBuilder<T: Scalar> {
     gdim: usize,
     element_data: (ReferenceCellType, usize),
-    element: CiarletElement<T, IdentityMap>,
+    element: CiarletElement<T, IdentityMap, T>,
     points_per_cell: usize,
     pub(crate) points: Vec<T>,
     cells: Vec<usize>,
@@ -56,7 +56,7 @@ pub struct SingleElementGridBuilder<T: RealScalar> {
     cell_indices: HashSet<usize>,
 }
 
-impl<T: RealScalar> SingleElementGridBuilder<T> {
+impl<T: Scalar> SingleElementGridBuilder<T> {
     /// Create a new grid builder
     pub fn new(gdim: usize, data: (ReferenceCellType, usize)) -> Self {
         Self::new_with_capacity(gdim, 0, 0, data)
@@ -69,7 +69,7 @@ impl<T: RealScalar> SingleElementGridBuilder<T> {
         ncells: usize,
         data: (ReferenceCellType, usize),
     ) -> Self {
-        let element = lagrange::create::<T>(data.0, data.1, Continuity::Standard);
+        let element = lagrange::create::<T, T>(data.0, data.1, Continuity::Standard);
         let points_per_cell = element.dim();
         Self {
             gdim,
@@ -88,8 +88,8 @@ impl<T: RealScalar> SingleElementGridBuilder<T> {
     }
 }
 
-impl<T: RealScalar> Builder for SingleElementGridBuilder<T> {
-    type Grid = SingleElementGrid<T, CiarletElement<T, IdentityMap>>;
+impl<T: Scalar> Builder for SingleElementGridBuilder<T> {
+    type Grid = SingleElementGrid<T, CiarletElement<T, IdentityMap, T>>;
     type T = T;
     type CellData<'a> = &'a [usize];
     type EntityDescriptor = ReferenceCellType;
@@ -135,7 +135,7 @@ impl<T: RealScalar> Builder for SingleElementGridBuilder<T> {
         self.add_cell(id, nodes);
     }
 
-    fn create_grid(&self) -> SingleElementGrid<T, CiarletElement<T, IdentityMap>> {
+    fn create_grid(&self) -> SingleElementGrid<T, CiarletElement<T, IdentityMap, T>> {
         let cell_vertices =
             self.extract_vertices(&self.cells, &[self.element_data.0], &[self.element_data.1]);
 
@@ -216,8 +216,8 @@ impl<T: RealScalar> Builder for SingleElementGridBuilder<T> {
     }
 }
 
-impl<T: RealScalar> GeometryBuilder for SingleElementGridBuilder<T> {
-    type GridGeometry = SingleElementGeometry<T, CiarletElement<T, IdentityMap>>;
+impl<T: Scalar> GeometryBuilder for SingleElementGridBuilder<T> {
+    type GridGeometry = SingleElementGeometry<T, CiarletElement<T, IdentityMap, T>>;
     fn create_geometry(
         &self,
         point_ids: &[usize],
@@ -225,14 +225,14 @@ impl<T: RealScalar> GeometryBuilder for SingleElementGridBuilder<T> {
         cell_points: &[usize],
         _cell_types: &[ReferenceCellType],
         _cell_degrees: &[usize],
-    ) -> SingleElementGeometry<T, CiarletElement<T, IdentityMap>> {
+    ) -> SingleElementGeometry<T, CiarletElement<T, IdentityMap, T>> {
         let npts = point_ids.len();
         let mut points = rlst_dynamic_array!(T, [self.gdim(), npts]);
         points.data_mut().unwrap().copy_from_slice(coordinates);
 
-        let family = LagrangeElementFamily::<T>::new(self.element_data.1, Continuity::Standard);
+        let family = LagrangeElementFamily::<T, T>::new(self.element_data.1, Continuity::Standard);
 
-        SingleElementGeometry::<T, CiarletElement<T, IdentityMap>>::new(
+        SingleElementGeometry::<T, CiarletElement<T, IdentityMap, T>>::new(
             self.element_data.0,
             points,
             cell_points,
@@ -241,7 +241,7 @@ impl<T: RealScalar> GeometryBuilder for SingleElementGridBuilder<T> {
     }
 }
 
-impl<T: RealScalar> TopologyBuilder for SingleElementGridBuilder<T> {
+impl<T: Scalar> TopologyBuilder for SingleElementGridBuilder<T> {
     type GridTopology = SingleTypeTopology;
     fn create_topology(
         &self,
@@ -301,7 +301,7 @@ impl<T: RealScalar> TopologyBuilder for SingleElementGridBuilder<T> {
     }
 }
 
-impl<T: RealScalar> GridBuilder for SingleElementGridBuilder<T> {
+impl<T: Scalar> GridBuilder for SingleElementGridBuilder<T> {
     fn create_grid_from_topology_geometry(
         &self,
         topology: <Self as TopologyBuilder>::GridTopology,
diff --git a/src/grid/local_grid/single_element/grid.rs b/src/grid/local_grid/single_element/grid.rs
index bf9294b..b1b2cdf 100644
--- a/src/grid/local_grid/single_element/grid.rs
+++ b/src/grid/local_grid/single_element/grid.rs
@@ -16,7 +16,7 @@ use crate::{
     geometry::{GeometryMap, SingleElementEntityGeometry, SingleElementGeometry},
     topology::single_type::{SingleTypeEntityTopology, SingleTypeTopology},
     traits::{Entity, Grid},
-    types::{Ownership, RealScalar},
+    types::{Ownership, Scalar},
 };
 #[cfg(feature = "mpi")]
 use mpi::traits::{Communicator, Equivalence};
@@ -34,7 +34,7 @@ use rlst::{SliceArray, rlst_dynamic_array};
 #[derive(Debug)]
 pub struct SingleElementGridEntity<
     'a,
-    T: RealScalar,
+    T: Scalar,
     E: MappedFiniteElement<CellType = ReferenceCellType, T = T>,
 > {
     grid: &'a SingleElementGrid<T, E>,
@@ -43,7 +43,7 @@ pub struct SingleElementGridEntity<
     entity_index: usize,
 }
 
-impl<'e, T: RealScalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>>
+impl<'e, T: Scalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>>
     SingleElementGridEntity<'e, T, E>
 {
     /// Create new
@@ -61,7 +61,7 @@ impl<'e, T: RealScalar, E: MappedFiniteElement<CellType = ReferenceCellType, T =
         }
     }
 }
-impl<T: RealScalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>> Entity
+impl<T: Scalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>> Entity
     for SingleElementGridEntity<'_, T, E>
 {
     type T = T;
@@ -110,7 +110,7 @@ impl<T: RealScalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>>
 #[derive(Debug)]
 pub struct SingleElementGridEntityIter<
     'a,
-    T: RealScalar,
+    T: Scalar,
     E: MappedFiniteElement<CellType = ReferenceCellType, T = T>,
 > {
     grid: &'a SingleElementGrid<T, E>,
@@ -118,7 +118,7 @@ pub struct SingleElementGridEntityIter<
     index: usize,
 }
 
-impl<'a, T: RealScalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>>
+impl<'a, T: Scalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>>
     SingleElementGridEntityIter<'a, T, E>
 {
     /// Create new
@@ -130,7 +130,7 @@ impl<'a, T: RealScalar, E: MappedFiniteElement<CellType = ReferenceCellType, T =
         }
     }
 }
-impl<'a, T: RealScalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>> Iterator
+impl<'a, T: Scalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>> Iterator
     for SingleElementGridEntityIter<'a, T, E>
 {
     type Item = SingleElementGridEntity<'a, T, E>;
@@ -143,10 +143,8 @@ impl<'a, T: RealScalar, E: MappedFiniteElement<CellType = ReferenceCellType, T =
 
 /// Serial single element grid
 #[derive(Debug)]
-pub struct SingleElementGrid<
-    T: RealScalar,
-    E: MappedFiniteElement<CellType = ReferenceCellType, T = T>,
-> {
+pub struct SingleElementGrid<T: Scalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>>
+{
     topology: SingleTypeTopology,
     geometry: SingleElementGeometry<T, E>,
 }
@@ -154,7 +152,7 @@ pub struct SingleElementGrid<
 #[cfg(feature = "serde")]
 #[derive(serde::Serialize, Debug, serde::Deserialize)]
 #[serde(bound = "for<'de2> T: serde::Deserialize<'de2>")]
-pub struct SerializableGrid<T: RealScalar + serde::Serialize>
+pub struct SerializableGrid<T: Scalar + serde::Serialize>
 where
     for<'de2> T: serde::Deserialize<'de2>,
 {
@@ -163,8 +161,8 @@ where
 }
 
 #[cfg(feature = "serde")]
-impl<T: RealScalar + serde::Serialize> ConvertToSerializable
-    for SingleElementGrid<T, CiarletElement<T, IdentityMap>>
+impl<T: Scalar + serde::Serialize> ConvertToSerializable
+    for SingleElementGrid<T, CiarletElement<T, IdentityMap, T>>
 {
     type SerializableType = SerializableGrid<T>;
     fn to_serializable(&self) -> SerializableGrid<T> {
@@ -181,7 +179,7 @@ impl<T: RealScalar + serde::Serialize> ConvertToSerializable
     }
 }
 
-impl<T: RealScalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>>
+impl<T: Scalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>>
     SingleElementGrid<T, E>
 {
     /// Create new
@@ -190,7 +188,7 @@ impl<T: RealScalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>>
     }
 }
 
-impl<T: RealScalar> SingleElementGrid<T, CiarletElement<T, IdentityMap>> {
+impl<T: Scalar> SingleElementGrid<T, CiarletElement<T, IdentityMap, T>> {
     /// Create new from raw data
     pub fn new_from_raw_data(
         coordinates: &[T],
@@ -203,9 +201,9 @@ impl<T: RealScalar> SingleElementGrid<T, CiarletElement<T, IdentityMap>> {
         let mut points = rlst_dynamic_array!(T, [gdim, npts]);
         points.data_mut().unwrap().copy_from_slice(coordinates);
 
-        let family = LagrangeElementFamily::<T>::new(geometry_degree, Continuity::Standard);
+        let family = LagrangeElementFamily::<T, T>::new(geometry_degree, Continuity::Standard);
 
-        let geometry = SingleElementGeometry::<T, CiarletElement<T, IdentityMap>>::new(
+        let geometry = SingleElementGeometry::<T, CiarletElement<T, IdentityMap, T>>::new(
             cell_type, points, cells, &family,
         );
 
@@ -226,7 +224,7 @@ impl<T: RealScalar> SingleElementGrid<T, CiarletElement<T, IdentityMap>> {
     }
 }
 
-impl<T: RealScalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>> Grid
+impl<T: Scalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>> Grid
     for SingleElementGrid<T, E>
 {
     type T = T;
@@ -325,11 +323,11 @@ impl<T: RealScalar, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>>
 }
 
 #[cfg(feature = "mpi")]
-impl<T: RealScalar + Equivalence, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>>
+impl<T: Scalar + Equivalence, E: MappedFiniteElement<CellType = ReferenceCellType, T = T>>
     DistributableGrid for SingleElementGrid<T, E>
 {
     type ParallelGrid<'a, C: Communicator + 'a> =
-        ParallelGridImpl<'a, C, SingleElementGrid<T, CiarletElement<T, IdentityMap>>>;
+        ParallelGridImpl<'a, C, SingleElementGrid<T, CiarletElement<T, IdentityMap, T>>>;
 
     fn distribute<'a, C: Communicator>(
         &self,
diff --git a/src/grid/parallel_grid.rs b/src/grid/parallel_grid.rs
index e66a49e..3a380ee 100644
--- a/src/grid/parallel_grid.rs
+++ b/src/grid/parallel_grid.rs
@@ -4,7 +4,7 @@ use crate::traits::{ConvertToSerializable, RONImportParallel};
 use crate::{
     grid::local_grid::LocalGrid,
     traits::{Grid, ParallelGrid},
-    types::{Ownership, RealScalar},
+    types::{Ownership, Scalar},
 };
 use mpi::traits::Communicator;
 use rlst::distributed_tools::IndexLayout;
@@ -60,7 +60,7 @@ where
     }
 }
 
-impl<T: RealScalar, C: Communicator, G: Grid<T = T> + Sync> ParallelGrid
+impl<T: Scalar, C: Communicator, G: Grid<T = T> + Sync> ParallelGrid
     for ParallelGridImpl<'_, C, G>
 {
     type LocalGrid = LocalGrid<G>;
diff --git a/src/shapes/cube.rs b/src/shapes/cube.rs
index d19ac5a..4aad4d9 100644
--- a/src/shapes/cube.rs
+++ b/src/shapes/cube.rs
@@ -3,16 +3,16 @@
 use crate::{
     grid::local_grid::{SingleElementGrid, SingleElementGridBuilder},
     traits::Builder,
-    types::RealScalar,
+    types::Scalar,
 };
 use ndelement::{ciarlet::CiarletElement, map::IdentityMap, types::ReferenceCellType};
 
 /// Create a unit interval grid
 ///
 /// The unit interval is the interval between (0,) and (1,)
-pub fn unit_interval<T: RealScalar>(
+pub fn unit_interval<T: Scalar>(
     nx: usize,
-) -> SingleElementGrid<T, CiarletElement<T, IdentityMap>> {
+) -> SingleElementGrid<T, CiarletElement<T, IdentityMap, T>> {
     let mut b = SingleElementGridBuilder::new_with_capacity(
         1,
         nx + 1,
@@ -33,11 +33,11 @@ pub fn unit_interval<T: RealScalar>(
 /// Create a unit square grid
 ///
 /// The unit square is the square with corners at (0,0), (1,0), (0,1) and (1,1)
-pub fn unit_square<T: RealScalar>(
+pub fn unit_square<T: Scalar>(
     nx: usize,
     ny: usize,
     cell_type: ReferenceCellType,
-) -> SingleElementGrid<T, CiarletElement<T, IdentityMap>> {
+) -> SingleElementGrid<T, CiarletElement<T, IdentityMap, T>> {
     let mut b = SingleElementGridBuilder::new_with_capacity(
         2,
         (nx + 1) * (ny + 1),
@@ -99,10 +99,10 @@ pub fn unit_square<T: RealScalar>(
 /// Create a grid of the boundary of a unit square
 ///
 /// The unit square is the square with corners at (0,0), (1,0), (0,1) and (1,1)
-pub fn unit_square_boundary<T: RealScalar>(
+pub fn unit_square_boundary<T: Scalar>(
     nx: usize,
     ny: usize,
-) -> SingleElementGrid<T, CiarletElement<T, IdentityMap>> {
+) -> SingleElementGrid<T, CiarletElement<T, IdentityMap, T>> {
     let mut b = SingleElementGridBuilder::new_with_capacity(
         2,
         2 * (nx + ny),
@@ -153,12 +153,12 @@ pub fn unit_square_boundary<T: RealScalar>(
 ///
 /// The unit cube is the cube with corners at (0,0,0), (1,0,0), (0,1,0), (1,1,0), (0,0,1),
 /// (1,0,1), (0,1,1) and (1,1,1)
-pub fn unit_cube<T: RealScalar>(
+pub fn unit_cube<T: Scalar>(
     nx: usize,
     ny: usize,
     nz: usize,
     cell_type: ReferenceCellType,
-) -> SingleElementGrid<T, CiarletElement<T, IdentityMap>> {
+) -> SingleElementGrid<T, CiarletElement<T, IdentityMap, T>> {
     let mut b = SingleElementGridBuilder::new_with_capacity(
         3,
         (nx + 1) * (ny + 1) * (nz + 1),
@@ -257,12 +257,12 @@ pub fn unit_cube<T: RealScalar>(
 ///
 /// The unit cube is the cube with corners at (0,0,0), (1,0,0), (0,1,0), (1,1,0), (0,0,1),
 /// (1,0,1), (0,1,1) and (1,1,1)
-pub fn unit_cube_boundary<T: RealScalar>(
+pub fn unit_cube_boundary<T: Scalar>(
     nx: usize,
     ny: usize,
     nz: usize,
     cell_type: ReferenceCellType,
-) -> SingleElementGrid<T, CiarletElement<T, IdentityMap>> {
+) -> SingleElementGrid<T, CiarletElement<T, IdentityMap, T>> {
     let mut b = SingleElementGridBuilder::new_with_capacity(
         3,
         (nx + 1) * (ny + 1) * (nz + 1) - (nx - 1) * (ny - 1) * (nz - 1),
@@ -394,11 +394,11 @@ pub fn unit_cube_boundary<T: RealScalar>(
 ///
 /// The unit cube is the cube with corners at (0,0,0), (1,0,0), (0,1,0), (1,1,0), (0,0,1),
 /// (1,0,1), (0,1,1) and (1,1,1)
-pub fn unit_cube_edges<T: RealScalar>(
+pub fn unit_cube_edges<T: Scalar>(
     nx: usize,
     ny: usize,
     nz: usize,
-) -> SingleElementGrid<T, CiarletElement<T, IdentityMap>> {
+) -> SingleElementGrid<T, CiarletElement<T, IdentityMap, T>> {
     let mut b = SingleElementGridBuilder::new_with_capacity(
         3,
         4 * (nx + ny + nz + 1),
diff --git a/src/shapes/regular_sphere.rs b/src/shapes/regular_sphere.rs
index a86d9dd..6dfc939 100644
--- a/src/shapes/regular_sphere.rs
+++ b/src/shapes/regular_sphere.rs
@@ -3,10 +3,9 @@
 use crate::{
     grid::local_grid::{SingleElementGrid, SingleElementGridBuilder},
     traits::Builder,
-    types::RealScalar,
+    types::Scalar,
 };
 use ndelement::{ciarlet::CiarletElement, map::IdentityMap, types::ReferenceCellType};
-use num::Float;
 use std::collections::{HashMap, hash_map::Entry::Vacant};
 
 /// Create a surface grid of a regular sphere
@@ -14,9 +13,9 @@ use std::collections::{HashMap, hash_map::Entry::Vacant};
 /// A regular sphere is created by starting with a regular octahedron. The shape is then refined `refinement_level` times.
 /// Each time the grid is refined, each triangle is split into four triangles (by adding lines connecting the midpoints of
 /// each edge). The new points are then scaled so that they are a distance of 1 from the origin.
-pub fn regular_sphere<T: RealScalar>(
+pub fn regular_sphere<T: Scalar>(
     refinement_level: u32,
-) -> SingleElementGrid<T, CiarletElement<T, IdentityMap>> {
+) -> SingleElementGrid<T, CiarletElement<T, IdentityMap, T>> {
     let mut b = SingleElementGridBuilder::new_with_capacity(
         3,
         2 + usize::pow(4, refinement_level + 1),
@@ -71,8 +70,7 @@ pub fn regular_sphere<T: RealScalar>(
                             half * (v_i[1] + v_j[1]),
                             half * (v_i[2] + v_j[2]),
                         ];
-                        let size =
-                            Float::sqrt(new_pt.iter().map(|x| Float::powi(*x, 2)).sum::<T>());
+                        let size = (new_pt.iter().map(|x| x.powi(2)).sum::<T>()).sqrt();
                         for i in new_pt.iter_mut() {
                             *i /= size;
                         }
diff --git a/src/shapes/screen.rs b/src/shapes/screen.rs
index 9726510..2cb25d1 100644
--- a/src/shapes/screen.rs
+++ b/src/shapes/screen.rs
@@ -3,16 +3,16 @@
 use crate::{
     grid::local_grid::{SingleElementGrid, SingleElementGridBuilder},
     traits::Builder,
-    types::RealScalar,
+    types::Scalar,
 };
 use ndelement::{ciarlet::CiarletElement, map::IdentityMap, types::ReferenceCellType};
 /// Create a square grid with triangle cells
 ///
 /// Create a grid of the square \[0,1\]^2 with triangle cells. The input ncells is the number of cells
 /// along each side of the square.
-pub fn screen_triangles<T: RealScalar>(
+pub fn screen_triangles<T: Scalar>(
     ncells: usize,
-) -> SingleElementGrid<T, CiarletElement<T, IdentityMap>> {
+) -> SingleElementGrid<T, CiarletElement<T, IdentityMap, T>> {
     if ncells == 0 {
         panic!("Cannot create a grid with 0 cells");
     }
@@ -61,9 +61,9 @@ pub fn screen_triangles<T: RealScalar>(
 ///
 /// Create a grid of the square \[0,1\]^2 with quadrilateral cells. The input ncells is the number of
 /// cells along each side of the square.
-pub fn screen_quadrilaterals<T: RealScalar>(
+pub fn screen_quadrilaterals<T: Scalar>(
     ncells: usize,
-) -> SingleElementGrid<T, CiarletElement<T, IdentityMap>> {
+) -> SingleElementGrid<T, CiarletElement<T, IdentityMap, T>> {
     if ncells == 0 {
         panic!("Cannot create a grid with 0 cells");
     }
@@ -107,7 +107,7 @@ pub fn screen_quadrilaterals<T: RealScalar>(
 /// Create a grid of the square \[0,2\]x\[0,1\] with triangle cells on the left half and quadrilateral
 /// cells on the right half. The input ncells is the number of cells along each side of the unit
 /// square.
-pub fn screen_mixed<T: Float + RealScalar>(ncells: usize) -> MixedGrid<T>
+pub fn screen_mixed<T: Float + Scalar>(ncells: usize) -> MixedGrid<T>
 {
     if ncells == 0 {
         panic!("Cannot create a grid with 0 cells");
diff --git a/src/traits/builder.rs b/src/traits/builder.rs
index 9e10a9b..3a6125b 100644
--- a/src/traits/builder.rs
+++ b/src/traits/builder.rs
@@ -1,5 +1,5 @@
 //! Grid builder
-use crate::{traits::Grid, types::RealScalar};
+use crate::{traits::Grid, types::Scalar};
 #[cfg(feature = "mpi")]
 use crate::{traits::ParallelGrid, types::GraphPartitioner};
 #[cfg(feature = "mpi")]
@@ -18,7 +18,7 @@ pub trait Builder {
     /// The type of the grid that the builder creates
     type Grid: Grid<EntityDescriptor = Self::EntityDescriptor>;
     /// The floating point type used for coordinates
-    type T: RealScalar;
+    type T: Scalar;
     /// The type of the data that is input to add a cell
     type CellData<'a>;
 
diff --git a/src/traits/entity.rs b/src/traits/entity.rs
index 0c27dba..fbd6791 100644
--- a/src/traits/entity.rs
+++ b/src/traits/entity.rs
@@ -1,6 +1,6 @@
 //! Traits for a mesh entity
 use crate::traits::{Geometry, Topology};
-use crate::types::{Ownership, RealScalar};
+use crate::types::{Ownership, Scalar};
 use std::fmt::Debug;
 use std::hash::Hash;
 
@@ -10,7 +10,7 @@ use std::hash::Hash;
 /// provides a unified interface to any type of entity.
 pub trait Entity {
     /// Scalar type
-    type T: RealScalar;
+    type T: Scalar;
 
     /// Type used as identifier of different entity types.
     /// In most cases this is given by [ReferenceCellType](ndelement::types::ReferenceCellType).
diff --git a/src/traits/geometry.rs b/src/traits/geometry.rs
index 2f54043..738a855 100644
--- a/src/traits/geometry.rs
+++ b/src/traits/geometry.rs
@@ -1,10 +1,10 @@
 //! Entity geometry
-use crate::types::RealScalar;
+use crate::types::Scalar;
 
 /// A point
 pub trait Point {
     /// Scalar type
-    type T: RealScalar;
+    type T: Scalar;
 
     /// The point's index
     fn index(&self) -> usize;
@@ -21,7 +21,7 @@ pub trait Point {
 /// The geometry contains information about all points that make up the entity.
 pub trait Geometry {
     /// Scalar type
-    type T: RealScalar;
+    type T: Scalar;
 
     /// Point Type
     type Point<'a>: Point<T = Self::T>
diff --git a/src/traits/geometry_map.rs b/src/traits/geometry_map.rs
index 24ee078..5b90c50 100644
--- a/src/traits/geometry_map.rs
+++ b/src/traits/geometry_map.rs
@@ -1,6 +1,6 @@
 //! Map from reference to physical space.
 
-use crate::types::RealScalar;
+use crate::types::Scalar;
 
 /// A geometry map allows the computation of maps from reference to physical space and their derivatives.
 ///
@@ -11,7 +11,7 @@ use crate::types::RealScalar;
 /// - The jacobians, transformation determinants and the normals of the physical entity.
 pub trait GeometryMap {
     /// Scalar type
-    type T: RealScalar;
+    type T: Scalar;
 
     /// The topoloical dimension of the entity being mapped.
     ///
diff --git a/src/traits/grid.rs b/src/traits/grid.rs
index dea05ff..7a1e96f 100644
--- a/src/traits/grid.rs
+++ b/src/traits/grid.rs
@@ -2,7 +2,7 @@
 use super::{Entity, GeometryMap};
 #[cfg(feature = "mpi")]
 use crate::types::GraphPartitioner;
-use crate::types::{Ownership, RealScalar};
+use crate::types::{Ownership, Scalar};
 #[cfg(feature = "mpi")]
 use mpi::traits::Communicator;
 #[cfg(feature = "mpi")]
@@ -16,7 +16,7 @@ use std::rc::Rc;
 /// A grid provides access to entities, their geometrical and their topological properties.
 pub trait Grid {
     /// Scalar type
-    type T: RealScalar;
+    type T: Scalar;
 
     /// Type used as identifier of different entity types
     type Entity<'a>: Entity<EntityDescriptor = Self::EntityDescriptor, T = Self::T>
@@ -109,7 +109,7 @@ pub trait Grid {
 #[cfg(feature = "mpi")]
 pub trait ParallelGrid {
     /// Type of the Grid
-    type T: RealScalar;
+    type T: Scalar;
 
     /// Local grid type
     type LocalGrid: Grid<T = Self::T>;
diff --git a/src/types.rs b/src/types.rs
index d674bc1..98f73d8 100644
--- a/src/types.rs
+++ b/src/types.rs
@@ -1,17 +1,15 @@
 //! Types
 #[cfg(feature = "mpi")]
 use mpi::traits::Equivalence;
-use num::Float;
 use rlst::{
     RlstScalar,
     dense::linalg::lapack::interface::{getrf::Getrf, getri::Getri},
 };
 
-/// A real scalar
-pub trait RealScalar: Float + RlstScalar<Real = Self> + Getrf + Getri {}
+/// A scalar.
+pub trait Scalar: RlstScalar + Getrf + Getri {}
 
-impl RealScalar for f32 {}
-impl RealScalar for f64 {}
+impl<T: RlstScalar + Getrf + Getri> Scalar for T {}
 
 /// A (cell, local index) pair
 ///