Mercurial > hg > octave-nkf
diff examples/mysparse.c @ 6580:d2bb3b8a8d20
[project @ 2007-04-25 22:19:03 by dbateman]
author | dbateman |
---|---|
date | Wed, 25 Apr 2007 22:20:28 +0000 |
parents | 11bb9bf343a0 |
children | 2aad75fcc93a |
line wrap: on
line diff
--- a/examples/mysparse.c +++ b/examples/mysparse.c @@ -20,8 +20,8 @@ if (mxIsComplex (prhs[0])) { - mexPrintf ("Matrix is %d-by-%d complex sparse matrix with %d elements\n", - m, n, nz); + mexPrintf ("Matrix is %d-by-%d complex sparse matrix", m, n); + mexPrintf (" with %d elements\n", nz); pr = mxGetPr (prhs[0]); pi = mxGetPi (prhs[0]); @@ -30,8 +30,8 @@ i = n; while (jc[i] == jc[i-1] && i != 0) i--; - mexPrintf ("last non-zero element (%d, %d) = (%g, %g)\n", ir[nz-1]+ 1, - i, pr[nz-1], pi[nz-1]); + mexPrintf ("last non-zero element (%d, %d) = (%g, %g)\n", + ir[nz-1]+ 1, i, pr[nz-1], pi[nz-1]); v = mxCreateSparse (m, n, nz, mxCOMPLEX); pr2 = mxGetPr (v); @@ -40,22 +40,22 @@ jc2 = mxGetJc (v); for (i = 0; i < nz; i++) - { - pr2[i] = 2 * pr[i]; - pi2[i] = 2 * pi[i]; - ir2[i] = ir[i]; - } + { + pr2[i] = 2 * pr[i]; + pi2[i] = 2 * pi[i]; + ir2[i] = ir[i]; + } for (i = 0; i < n + 1; i++) - jc2[i] = jc[i]; + jc2[i] = jc[i]; if (nlhs > 0) - plhs[0] = v; + plhs[0] = v; } else if (mxIsLogical (prhs[0])) { bool *pbr, *pbr2; - mexPrintf ("Matrix is %d-by-%d logical sparse matrix with %d elements\n", - m, n, nz); + mexPrintf ("Matrix is %d-by-%d logical sparse matrix", m, n); + mexPrintf (" with %d elements\n", nz); pbr = mxGetLogicals (prhs[0]); ir = mxGetIr (prhs[0]); @@ -64,7 +64,7 @@ i = n; while (jc[i] == jc[i-1] && i != 0) i--; mexPrintf ("last non-zero element (%d, %d) = %d\n", ir[nz-1]+ 1, - i, pbr[nz-1]); + i, pbr[nz-1]); v = mxCreateSparseLogicalMatrix (m, n, nz); pbr2 = mxGetLogicals (v); @@ -72,21 +72,20 @@ jc2 = mxGetJc (v); for (i = 0; i < nz; i++) - { - pbr2[i] = pbr[i]; - ir2[i] = ir[i]; - } + { + pbr2[i] = pbr[i]; + ir2[i] = ir[i]; + } for (i = 0; i < n + 1; i++) - jc2[i] = jc[i]; + jc2[i] = jc[i]; if (nlhs > 0) - plhs[0] = v; + plhs[0] = v; } else { - - mexPrintf ("Matrix is %d-by-%d real sparse matrix with %d elements\n", - m, n, nz); + mexPrintf ("Matrix is %d-by-%d real sparse matrix", m, n); + mexPrintf (" with %d elements\n", nz); pr = mxGetPr (prhs[0]); ir = mxGetIr (prhs[0]); @@ -95,7 +94,7 @@ i = n; while (jc[i] == jc[i-1] && i != 0) i--; mexPrintf ("last non-zero element (%d, %d) = %g\n", ir[nz-1]+ 1, - i, pr[nz-1]); + i, pr[nz-1]); v = mxCreateSparse (m, n, nz, mxREAL); pr2 = mxGetPr (v); @@ -103,14 +102,14 @@ jc2 = mxGetJc (v); for (i = 0; i < nz; i++) - { - pr2[i] = 2 * pr[i]; - ir2[i] = ir[i]; - } + { + pr2[i] = 2 * pr[i]; + ir2[i] = ir[i]; + } for (i = 0; i < n + 1; i++) - jc2[i] = jc[i]; + jc2[i] = jc[i]; if (nlhs > 0) - plhs[0] = v; + plhs[0] = v; } }