diff --git a/level-set/src/geomGamma.cpp b/level-set/src/geomGamma.cpp
index 032f717..962352e 100644
--- a/level-set/src/geomGamma.cpp
+++ b/level-set/src/geomGamma.cpp
@@ -192,7 +192,7 @@ DEFUN_DLD (__levelset_geomGamma, args, nargout,
       const Matrix inout = args(4).matrix_value ();
 
       /* Extract and check the dimensions.  */
-      const unsigned nNodes = phi.nelem ();
+      const unsigned nNodes = phi.numel ();
       const unsigned nElem = getDimension (nodelist, -1, 4);
       const unsigned nBdryEl = getDimension (bdryInd, -1, 1);
       getDimension (edges, nBdryEl, 4);
diff --git a/level-set/src/internal_fastmarching.cpp b/level-set/src/internal_fastmarching.cpp
index 543e6dd..cfccf08 100644
--- a/level-set/src/internal_fastmarching.cpp
+++ b/level-set/src/internal_fastmarching.cpp
@@ -74,7 +74,7 @@ DEFUN_DLD (__levelset_internal_fastmarching, args, nargout,
         {
           const Array<octave_idx_type> idx = getOctaveIdx (c);
           assert (c.size () == D
-                  && static_cast<dimensionT> (idx.length ()) == D);
+                  && static_cast<dimensionT> (idx.numel ()) == D);
 
           if (domain(idx))
             {
@@ -99,7 +99,7 @@ DEFUN_DLD (__levelset_internal_fastmarching, args, nargout,
         {
           const Array<octave_idx_type> idx = getOctaveIdx (c);
           assert (c.size () == D
-                  && static_cast<dimensionT> (idx.length ()) == D);
+                  && static_cast<dimensionT> (idx.numel ()) == D);
 
           const Grid& constGrid(grid);
           const Entry* e = constGrid.get (c);
diff --git a/level-set/src/internal_mesh.cpp b/level-set/src/internal_mesh.cpp
index b772615..acf4cf4 100644
--- a/level-set/src/internal_mesh.cpp
+++ b/level-set/src/internal_mesh.cpp
@@ -300,7 +300,7 @@ getInnerSegment (const octave_scalar_map& segs,
 
   assert (innerPts.empty ());
   const ColumnVector inners = segs.contents ("inners").column_vector_value ();
-  const unsigned nInners = inners.nelem ();
+  const unsigned nInners = inners.numel ();
   for (unsigned i = 0; i < nInners; ++i)
     innerPts.push_back (inners(nInners - i - 1) - 1);
 }
@@ -387,7 +387,7 @@ DEFUN_DLD (__levelset_internal_mesh, args, nargout,
         {
           const unsigned cur = bdryElems(i) - 1;
           const Cell cellSegs = bdryelSegs(i).cell_value ();
-          const unsigned nSegs = cellSegs.nelem ();
+          const unsigned nSegs = cellSegs.numel ();
 
           std::vector<octave_scalar_map> segs;
           indexArr endEdges;