1 /*
2 Copyright (C) 2018 Fredrik Johansson
3
4 This file is part of Arb.
5
6 Arb is free software: you can redistribute it and/or modify it under
7 the terms of the GNU Lesser General Public License (LGPL) as published
8 by the Free Software Foundation; either version 2.1 of the License, or
9 (at your option) any later version. See <http://www.gnu.org/licenses/>.
10 */
11
12 #include "arb_mat.h"
13
14 void
arb_mat_hilbert(arb_mat_t mat,slong prec)15 arb_mat_hilbert(arb_mat_t mat, slong prec)
16 {
17 slong R, C, i, j;
18
19 R = arb_mat_nrows(mat);
20 C = arb_mat_ncols(mat);
21
22 for (i = 0; i < R; i++)
23 {
24 for (j = 0; j < C; j++)
25 {
26 arb_one(arb_mat_entry(mat, i, j));
27 arb_div_ui(arb_mat_entry(mat, i, j),
28 arb_mat_entry(mat, i, j), i + j + 1, prec);
29 }
30 }
31 }
32
33