Mercurial > fem-fenics-eugenio
changeset 74:d123adbe296b
New parameters for the creation of the sparse matrix.
| author | gedeone-octave <marcovass89@hotmail.it> |
|---|---|
| date | Fri, 02 Aug 2013 22:51:34 +0200 |
| parents | 7c74e99fb548 |
| children | d6df5cc8ef53 |
| files | inst/generate_rhs.m |
| diffstat | 1 files changed, 5 insertions(+), 5 deletions(-) [+] |
line wrap: on
line diff
--- a/inst/generate_rhs.m Mon Jul 29 22:55:03 2013 +0200 +++ b/inst/generate_rhs.m Fri Aug 02 22:51:34 2013 +0200 @@ -105,9 +105,9 @@ \n\ std::size_t nr = A.size (0), nc = A.size (1);\n\ // nz shoud be estimated in a better way\n\ - double alpha = 0.005;\n\ - std::size_t nz = alpha * nr * nc;\n\ - SparseMatrix sm (nr, nc, ceil (nz));\n\ + double alpha = 0.0005;\n\ + octave_idx_type nz = alpha * nr * nc;\n\ + SparseMatrix sm (nr, nc, nz);\n\ \n\ std::size_t ii = 0;\n\ sm.cidx (0) = 0;\n\ @@ -120,8 +120,8 @@ {\n\ if (ii == nz)\n\ {\n\ - nz = 1.1 * ((nc * ii) / (j + 1));\n\ - sm.change_capacity (ceil (nz));\n\ + nz = 1.2 * ((nc * ii) / (j + 1));\n\ + sm.change_capacity (nz);\n\ }\n\ sm.data(ii) = tmp;\n\ sm.ridx(ii) = i;\n\