Инициализация массива без индексации

double *p_elements = elements; 
//зануление элементов 
for(int i = 0; i < size*size; i++, p_elements++)
 *p_elements = 0; 
//присваивание диагональным
 for(double i = 0; i < size; i++,p_elements = i * (1+col_size)) 
*p_elements= diagonal_value;

#include <iostream>
 
int main()
{
    double d_val = 9;
    
    const int n = 3;
    double x[n*n] = { 0 };
    
    int k = n;
    while(k--) x[k + k*n] = d_val;
    
    for(int i=0, e = n*n; i<e; ++i)
    {
        if(!(i%n)) std::cout << "\n";
        std::cout << x[i] << " ";
    }
    
    return 0;
}

Matrix::Matrix(int size, double diagonal_value)
{
    simple_copy = false;
    
    error_code = 0;
    if (size == 0)
    {
        error_code = -1;
        
        row_size = 0;
        col_size = 0;
        elements = nullptr;
        
        return;
    }
    
    row_size = size;
    col_size = size;
    
    elements = new double[size * size];
    
    double *p_elements = elements;
    
    for(int i = 0; i < size*size; i++, p_elements++)
        *p_elements = 0;
    
 
    for(int i = 0; i < size; i++,p_elements++)
        *p_elements + i * (1 + col_size) = diagonal_value; //вот это нужно исправить 
}

        for (double *p = x, *end = p + n * n; p < end; p += n + 1)
            *p = diagonal_value;

#include <iostream>
 
int main()
{
    double d_val = 6;
    
    const int n = 4;
    double x[n*n] = { 0 };
    
    double* k = x;
    double* e = x + n*n;
    while(k < e)
    {
        *k = d_val;
        k += n + 1; 
    }
    
    for(int i=0, e = n*n; i<e; ++i)
    {
        if(!(i%n)) std::cout << "\n";
        std::cout << x[i] << " ";
    }
    
    return 0;
}