diff --git a/KokkosStream.cpp b/KokkosStream.cpp
index 09c1eaf..00efe92 100644
--- a/KokkosStream.cpp
+++ b/KokkosStream.cpp
@@ -119,6 +119,21 @@ void KokkosStream<T>::triad()
   Kokkos::fence();
 }
 
+template <class T>
+void KokkosStream<T>::nstream()
+{
+  Kokkos::View<T*> a(*d_a);
+  Kokkos::View<T*> b(*d_b);
+  Kokkos::View<T*> c(*d_c);
+
+  const T scalar = startScalar;
+  Kokkos::parallel_for(array_size, KOKKOS_LAMBDA (const long index)
+  {
+    a[index] += b[index] + scalar*c[index];
+  });
+  Kokkos::fence();
+}
+
 template <class T>
 T KokkosStream<T>::dot()
 {
diff --git a/KokkosStream.hpp b/KokkosStream.hpp
index e88622c..3aa7cf5 100644
--- a/KokkosStream.hpp
+++ b/KokkosStream.hpp
@@ -41,6 +41,7 @@ class KokkosStream : public Stream<T>
     virtual void add() override;
     virtual void mul() override;
     virtual void triad() override;
+    virtual void nstream() override;
     virtual T dot() override;
 
     virtual void init_arrays(T initA, T initB, T initC) override;