Valarray C++

/* The following code example is taken from the book
 * "The C++ Standard Library - A Tutorial and Reference"
 * by Nicolai M. Josuttis, Addison-Wesley, 1999
 *
 * (C) Copyright Nicolai M. Josuttis 1999.
 * Permission to copy, use, modify, sell and distribute this software
 * is granted provided this copyright notice appears in all copies.
 * This software is provided "as is" without express or implied
 * warranty, and with no claim as to its suitability for any purpose.
 */
#include 
#include 
using namespace std;
// print three-dimensional valarray line-by-line
template
void printValarray3D (const valarray& va, int dim1, int dim2)
{
    for (int i=0; i        for (int j=0; j            for (int k=0; k                cout << va[i*dim1*dim2+j*dim1+k] << ' ';
            }
            cout << '\n';
        }
        cout << '\n';
    }
    cout << endl;
}
int main()
{
    /* valarray with 24 elements
     * - two groups
     * - four rows
     * - three columns
     */
    valarray va(24);
    // fill valarray with values
    for (int i=0; i<24; i++) {
        va[i] = i;
    }
    // print valarray
    printValarray3D (va, 3, 4);
    // we need two two-dimensional subsets of three times 3 values
    // in two 12-element arrays
    size_t lengthvalues[] = {  2, 3 };
    size_t stridevalues[] = { 12, 3 };
    valarray length(lengthvalues,2);
    valarray stride(stridevalues,2);
    // assign the second column of the first three rows
    // to the first column of the first three rows
    va[gslice(0,length,stride)]
        = valarray(va[gslice(1,length,stride)]);
    // add and assign the third of the first three rows
    // to the first of the first three rows
    va[gslice(0,length,stride)]
        += valarray(va[gslice(2,length,stride)]);
    
    // print valarray
    printValarray3D (va, 3, 4);
}
/* 
0 1 2
3 4 5
6 7 8
9 10 11
12 13 14
15 16 17
18 19 20
21 22 23
3 1 2
9 4 5
15 7 8
9 10 11
27 13 14
33 16 17
39 19 20
21 22 23
 */