Вывести все наддиагональные элементы матрицы в порядке, заданном диаграммой

#include <stdio.h>
#include <stdlib.h>
#include <time.h>
 
#define MIN -10
#define MAX 10
#define FIELD_WIDTH 7
#define PREC 2
#define SIZE 8
 
void dump(const double* matrix); /* pass pointer to the first item */
void print(const double* matrix, int offset); /* pass pointer to the first item */
 
int main(void)
{
    double matrix[SIZE][SIZE];
    int i, j;
 
    srand((size_t) time(NULL));
    
    for(i = 0; i < SIZE; ++i)
    for(j = 0; j < SIZE; ++j)
        matrix[i][j] = ((double) rand()) / RAND_MAX * (MAX - MIN + 1) + MIN;
 
    puts("\nGenerated matrix:");
    dump(&matrix[0][0]);
    
    puts("\nElements listed in the given order:");
    print(&matrix[0][0], SIZE);
    
    exit(0);
}
 
void dump(const double* matrix)
{
    int i, j;
 
    for(i = 0; i < SIZE; ++i)
    {
    for(j = 0; j < SIZE; ++j)
        printf("%*.*f", FIELD_WIDTH, PREC, matrix[i * SIZE + j]);
    putchar('\n');
    }
}
 
void print(const double* matrix, int offset)
{
    int i, j;
 
    /* stop recursion */
    if(offset <= 0)
    return;
    
    for(j = 0; j < offset; ++j)
    printf("%*.*f\n", FIELD_WIDTH, PREC, matrix[j]);
 
    for(i = 1, j = offset - 1; i < offset; ++i)
    printf("%*.*f\n", FIELD_WIDTH, PREC, matrix[i * SIZE + j]);
 
    for(i = offset - 2; i > 0; --i)
    printf("%*.*f\n", FIELD_WIDTH, PREC, matrix[i * SIZE + i]);
 
    putchar('\n');
    
    print(&matrix[SIZE + 2], offset - 3);
}