Mercurial > hg > octave-nkf

comparison liboctave/CMatrix.cc @ 7265:7da4a5262e2e

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2007-12-06 19:16:47 by jwe]
author jwe
date Thu, 06 Dec 2007 19:16:48 +0000
parents d07cb867891b
children f9df7f7520e7 6c36c5d3c38b
comparison
equal deleted inserted replaced
7264:fac10884ddd4 7265:7da4a5262e2e
2863 for (octave_idx_type j = 7; j >= 0; j--) 2863 for (octave_idx_type j = 7; j >= 0; j--)
2864 { 2864 {
2865 for (octave_idx_type i = 0; i < nc; i++) 2865 for (octave_idx_type i = 0; i < nc; i++)
2866 { 2866 {
2867 octave_idx_type k = i * nc + i; 2867 octave_idx_type k = i * nc + i;
2868 pnpp [k] = pnpp [k] + padec [j]; 2868 pnpp[k] += padec[j];
2869 pdpp [k] = pdpp [k] + minus_one_j * padec [j]; 2869 pdpp[k] += minus_one_j * padec[j];
2870 } 2870 }
2871
2871 npp = m * npp; 2872 npp = m * npp;
2873 pnpp = npp.fortran_vec ();
2874
2872 dpp = m * dpp; 2875 dpp = m * dpp;
2876 pdpp = dpp.fortran_vec ();
2877
2873 minus_one_j *= -1; 2878 minus_one_j *= -1;
2874 } 2879 }
2875 2880
2876 // Zero power. 2881 // Zero power.
2877 2882