Diff of the two buildlogs: -- --- b1/build.log 2025-10-23 05:38:47.583316868 +0000 +++ b2/build.log 2025-10-23 05:42:05.499542717 +0000 @@ -1,6 +1,6 @@ I: pbuilder: network access will be disabled during build -I: Current time: Wed Oct 22 17:33:20 -12 2025 -I: pbuilder-time-stamp: 1761197600 +I: Current time: Thu Nov 26 02:01:49 +14 2026 +I: pbuilder-time-stamp: 1795608109 I: Building the build Environment I: extracting base tarball [/var/cache/pbuilder/forky-reproducible-base.tgz] I: copying local configuration @@ -22,52 +22,84 @@ dpkg-source: info: unpacking octave-linear-algebra_2.2.3-5.debian.tar.xz I: using fakeroot in build. I: Installing the build-deps -I: user script /srv/workspace/pbuilder/994964/tmp/hooks/D02_print_environment starting +I: user script /srv/workspace/pbuilder/842317/tmp/hooks/D01_modify_environment starting +debug: Running on infom02-amd64. +I: Changing host+domainname to test build reproducibility +I: Adding a custom variable just for the fun of it... +I: Changing /bin/sh to bash +'/bin/sh' -> '/bin/bash' +lrwxrwxrwx 1 root root 9 Nov 25 12:02 /bin/sh -> /bin/bash +I: Setting pbuilder2's login shell to /bin/bash +I: Setting pbuilder2's GECOS to second user,second room,second work-phone,second home-phone,second other +I: user script /srv/workspace/pbuilder/842317/tmp/hooks/D01_modify_environment finished +I: user script /srv/workspace/pbuilder/842317/tmp/hooks/D02_print_environment starting I: set - BUILDDIR='/build/reproducible-path' - BUILDUSERGECOS='first user,first room,first work-phone,first home-phone,first other' - BUILDUSERNAME='pbuilder1' - BUILD_ARCH='amd64' - DEBIAN_FRONTEND='noninteractive' - DEB_BUILD_OPTIONS='buildinfo=+all reproducible=+all parallel=12 ' - DISTRIBUTION='forky' - HOME='/root' - HOST_ARCH='amd64' + BASH=/bin/sh + BASHOPTS=checkwinsize:cmdhist:complete_fullquote:extquote:force_fignore:globasciiranges:globskipdots:hostcomplete:interactive_comments:patsub_replacement:progcomp:promptvars:sourcepath + BASH_ALIASES=() + BASH_ARGC=() + BASH_ARGV=() + BASH_CMDS=() + BASH_LINENO=([0]="12" [1]="0") + BASH_LOADABLES_PATH=/usr/local/lib/bash:/usr/lib/bash:/opt/local/lib/bash:/usr/pkg/lib/bash:/opt/pkg/lib/bash:. + BASH_SOURCE=([0]="/tmp/hooks/D02_print_environment" [1]="/tmp/hooks/D02_print_environment") + BASH_VERSINFO=([0]="5" [1]="3" [2]="3" [3]="1" [4]="release" [5]="x86_64-pc-linux-gnu") + BASH_VERSION='5.3.3(1)-release' + BUILDDIR=/build/reproducible-path + BUILDUSERGECOS='second user,second room,second work-phone,second home-phone,second other' + BUILDUSERNAME=pbuilder2 + BUILD_ARCH=amd64 + DEBIAN_FRONTEND=noninteractive + DEB_BUILD_OPTIONS='buildinfo=+all reproducible=+all parallel=16 ' + DIRSTACK=() + DISTRIBUTION=forky + EUID=0 + FUNCNAME=([0]="Echo" [1]="main") + GROUPS=() + HOME=/root + HOSTNAME=i-capture-the-hostname + HOSTTYPE=x86_64 + HOST_ARCH=amd64 IFS=' ' - INVOCATION_ID='10664bdb00e1404997221f966c17c74b' - LANG='C' - LANGUAGE='en_US:en' - LC_ALL='C' - MAIL='/var/mail/root' - OPTIND='1' - PATH='/usr/sbin:/usr/bin:/sbin:/bin:/usr/games' - PBCURRENTCOMMANDLINEOPERATION='build' - PBUILDER_OPERATION='build' - PBUILDER_PKGDATADIR='/usr/share/pbuilder' - PBUILDER_PKGLIBDIR='/usr/lib/pbuilder' - PBUILDER_SYSCONFDIR='/etc' - PPID='994964' - PS1='# ' - PS2='> ' + INVOCATION_ID=138e5584070d44708ac44077b3696e0a + LANG=C + LANGUAGE=et_EE:et + LC_ALL=C + MACHTYPE=x86_64-pc-linux-gnu + MAIL=/var/mail/root + OPTERR=1 + OPTIND=1 + OSTYPE=linux-gnu + PATH=/usr/sbin:/usr/bin:/sbin:/bin:/usr/games:/i/capture/the/path + PBCURRENTCOMMANDLINEOPERATION=build + PBUILDER_OPERATION=build + PBUILDER_PKGDATADIR=/usr/share/pbuilder + PBUILDER_PKGLIBDIR=/usr/lib/pbuilder + PBUILDER_SYSCONFDIR=/etc + PIPESTATUS=([0]="0") + POSIXLY_CORRECT=y + PPID=842317 PS4='+ ' - PWD='/' - SHELL='/bin/bash' - SHLVL='2' - SUDO_COMMAND='/usr/bin/timeout -k 18.1h 18h /usr/bin/ionice -c 3 /usr/bin/nice /usr/sbin/pbuilder --build --configfile /srv/reproducible-results/rbuild-debian/r-b-build.Sst4VMmV/pbuilderrc_5Ixl --distribution forky --hookdir /etc/pbuilder/first-build-hooks --debbuildopts -b --basetgz /var/cache/pbuilder/forky-reproducible-base.tgz --buildresult /srv/reproducible-results/rbuild-debian/r-b-build.Sst4VMmV/b1 --logfile b1/build.log octave-linear-algebra_2.2.3-5.dsc' - SUDO_GID='109' - SUDO_HOME='/var/lib/jenkins' - SUDO_UID='104' - SUDO_USER='jenkins' - TERM='unknown' - TZ='/usr/share/zoneinfo/Etc/GMT+12' - USER='root' - _='/usr/bin/systemd-run' + PWD=/ + SHELL=/bin/bash + SHELLOPTS=braceexpand:errexit:hashall:interactive-comments:posix + SHLVL=3 + SUDO_COMMAND='/usr/bin/timeout -k 24.1h 24h /usr/bin/ionice -c 3 /usr/bin/nice -n 11 /usr/bin/unshare --uts -- /usr/sbin/pbuilder --build --configfile /srv/reproducible-results/rbuild-debian/r-b-build.Sst4VMmV/pbuilderrc_5fyA --distribution forky --hookdir /etc/pbuilder/rebuild-hooks --debbuildopts -b --basetgz /var/cache/pbuilder/forky-reproducible-base.tgz --buildresult /srv/reproducible-results/rbuild-debian/r-b-build.Sst4VMmV/b2 --logfile b2/build.log octave-linear-algebra_2.2.3-5.dsc' + SUDO_GID=109 + SUDO_HOME=/var/lib/jenkins + SUDO_UID=104 + SUDO_USER=jenkins + TERM=unknown + TZ=/usr/share/zoneinfo/Etc/GMT-14 + UID=0 + USER=root + _='I: set' I: uname -a - Linux infom01-amd64 6.12.48+deb13-cloud-amd64 #1 SMP PREEMPT_DYNAMIC Debian 6.12.48-1 (2025-09-20) x86_64 GNU/Linux + Linux i-capture-the-hostname 6.12.48+deb13-amd64 #1 SMP PREEMPT_DYNAMIC Debian 6.12.48-1 (2025-09-20) x86_64 GNU/Linux I: ls -l /bin - lrwxrwxrwx 1 root root 7 Aug 10 12:30 /bin -> usr/bin -I: user script /srv/workspace/pbuilder/994964/tmp/hooks/D02_print_environment finished + lrwxrwxrwx 1 root root 7 Aug 10 2025 /bin -> usr/bin +I: user script /srv/workspace/pbuilder/842317/tmp/hooks/D02_print_environment finished -> Attempting to satisfy build-dependencies -> Creating pbuilder-satisfydepends-dummy package Package: pbuilder-satisfydepends-dummy @@ -642,7 +674,7 @@ Get: 535 http://deb.debian.org/debian forky/main amd64 gfortran amd64 4:15.2.0-4 [1436 B] Get: 536 http://deb.debian.org/debian forky/main amd64 octave-dev amd64 10.3.0-2 [1072 kB] Get: 537 http://deb.debian.org/debian forky/main amd64 dh-octave all 1.11.4 [24.5 kB] -Fetched 282 MB in 4s (72.6 MB/s) +Fetched 282 MB in 4s (66.1 MB/s) Preconfiguring packages ... Selecting previously unselected package libexpat1:amd64. (Reading database ... (Reading database ... 5% (Reading database ... 10% (Reading database ... 15% (Reading database ... 20% (Reading database ... 25% (Reading database ... 30% (Reading database ... 35% (Reading database ... 40% (Reading database ... 45% (Reading database ... 50% (Reading database ... 55% (Reading database ... 60% (Reading database ... 65% (Reading database ... 70% (Reading database ... 75% (Reading database ... 80% (Reading database ... 85% (Reading database ... 90% (Reading database ... 95% (Reading database ... 100% (Reading database ... 19898 files and directories currently installed.) @@ -2848,7 +2880,11 @@ fakeroot is already the newest version (1.37.1.2-1). 0 upgraded, 0 newly installed, 0 to remove and 0 not upgraded. I: Building the package -I: Running cd /build/reproducible-path/octave-linear-algebra-2.2.3/ && env PATH="/usr/sbin:/usr/bin:/sbin:/bin:/usr/games" HOME="/nonexistent/first-build" dpkg-buildpackage -us -uc -b && env PATH="/usr/sbin:/usr/bin:/sbin:/bin:/usr/games" HOME="/nonexistent/first-build" dpkg-genchanges -S > ../octave-linear-algebra_2.2.3-5_source.changes +I: user script /srv/workspace/pbuilder/842317/tmp/hooks/A99_set_merged_usr starting +Not re-configuring usrmerge for forky +I: user script /srv/workspace/pbuilder/842317/tmp/hooks/A99_set_merged_usr finished +hostname: Name or service not known +I: Running cd /build/reproducible-path/octave-linear-algebra-2.2.3/ && env PATH="/usr/sbin:/usr/bin:/sbin:/bin:/usr/games:/i/capture/the/path" HOME="/nonexistent/second-build" dpkg-buildpackage -us -uc -b && env PATH="/usr/sbin:/usr/bin:/sbin:/bin:/usr/games:/i/capture/the/path" HOME="/nonexistent/second-build" dpkg-genchanges -S > ../octave-linear-algebra_2.2.3-5_source.changes dpkg-buildpackage: info: source package octave-linear-algebra dpkg-buildpackage: info: source version 2.2.3-5 dpkg-buildpackage: info: source distribution unstable @@ -2889,15 +2925,23 @@ Checking package... Run the unit tests... Checking m files ... -[inst/cartprod.m] ->>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/cartprod.m -***** assert(cartprod(1:2,0:1),[1 0; 2 0; 1 1; 2 1]) +[inst/vec_projection.m] +>>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/vec_projection.m +***** test + assert (vec_projection ([1,0,0], [0.5,0.5,0]), 0.70711,1e-5); + assert (vec_projection ([1,2000,0], [0.5,15,0]), 1998.9, 1e-1); + assert (vec_projection ([1,-2000,0], [0.5,15,0]), -1998.9, 1e-1); + assert (vec_projection ([7,7,0], [15,0,0]), 7.000, 1e-10); + assert (vec_projection ([1,1,0], [1.05,0.94,0]), 1.4121, 1e-4); + assert (vec_projection ([1,1.1,0], [1.05,0.94,0]), 1.4788, 1e-4); 1 test, 1 passed, 0 known failure, 0 skipped -[inst/circulant_inv.m] ->>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/circulant_inv.m -***** shared v +[inst/circulant_eig.m] +>>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/circulant_eig.m +***** shared v,C,vs,lambda v = [1 2 3]'; -***** assert (gallery ("circul", circulant_inv (v)), inv (gallery ("circul", v)), 10*eps); + C = gallery("circul", v)'; + [vs lambda] = circulant_eig(v); +***** assert (vs*lambda, C*vs, 100*eps); 1 test, 1 passed, 0 known failure, 0 skipped [inst/circulant_matrix_vector_product.m] >>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/circulant_matrix_vector_product.m @@ -2909,19 +2953,23 @@ circulant_make_matrix at line 43 column 5 __test__ at line 2 column 1 test at line 685 column 11 - /tmp/tmp.5XEBms9iOM at line 30 column 2 + /tmp/tmp.hZ0ZyEmVbp at line 30 column 2 1 test, 1 passed, 0 known failure, 0 skipped -[inst/smwsolve.m] ->>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/smwsolve.m +[inst/funm.m] +>>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/funm.m +***** function b = fsin (a) + b = sin (a); +***** endfunction ***** test - A = 2.1*eye (10); - u = rand (10, 2); u /= diag (norm (u, "cols")); - v = rand (10, 2); v /= diag (norm (v, "cols")); - b = rand (10, 2); - x1 = (A + u*v') \ b; - x2 = smwsolve (A, u, v, b); - assert (x1, x2, 1e-13); + mtx = randn (100); + assert (funm (mtx, "fsin"), thfm (mtx, "sin"), 1e-9) +1 test, 1 passed, 0 known failure, 0 skipped +[inst/circulant_inv.m] +>>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/circulant_inv.m +***** shared v + v = [1 2 3]'; +***** assert (gallery ("circul", circulant_inv (v)), inv (gallery ("circul", v)), 10*eps); 1 test, 1 passed, 0 known failure, 0 skipped [inst/cod.m] >>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/cod.m @@ -2938,12 +2986,15 @@ [q, r, z, p] = cod (a); assert (norm (q*r*z' - a*p) / norm (a) < 1e-10); 3 tests, 3 passed, 0 known failure, 0 skipped -[inst/@blksparse/ctranspose.m] ->>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/@blksparse/ctranspose.m -***** test - r = blksparse ([1,2],[1,2],cat(3,eye(2),[1 2; -2 1])); - rt = r'; - assert (full(rt'),full(r)); +[inst/circulant_make_matrix.m] +>>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/circulant_make_matrix.m +***** shared v,C + v = [1 2 3]'; C = [1 3 2; 2 1 3; 3 2 1]; +***** assert (circulant_make_matrix(v), C); +1 test, 1 passed, 0 known failure, 0 skipped +[inst/cartprod.m] +>>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/cartprod.m +***** assert(cartprod(1:2,0:1),[1 0; 2 0; 1 1; 2 1]) 1 test, 1 passed, 0 known failure, 0 skipped [inst/nmf_bpas.m] >>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/nmf_bpas.m @@ -2957,12 +3008,12 @@ Stop: tolerance reached. ***** test [W,H,iter,HIS]=nmf_bpas(A,k,'verbose'); - initGrNormW = 347.82 - initGrNormH = 420.35 - initNorm = 14.016 - SC1 = 1.0912 - SC2 = 545.59 - SC3 = 0.4840 + initGrNormW = 446.42 + initGrNormH = 524.75 + initNorm = 14.236 + SC1 = 1.3779 + SC2 = 688.95 + SC3 = 0.5238 W_density = 1 H_density = 1 par = @@ -2973,8 +3024,8 @@ n = 20 type = regularized nnls_solver = bp - alpha = 0.4979 - beta = 0.4979 + alpha = 0.5052 + beta = 0.5052 max_iter = 100 min_iter = 20 max_time = 1000 @@ -2986,438 +3037,438 @@ scalar structure containing the fields: - initGrNormW = 347.82 - initGrNormH = 420.35 - initNorm = 14.016 - SC1 = 0.029083 - SC2 = 0.028502 - SC3 = 2.0715 - W_density = 0.9700 - H_density = 0.6900 + initGrNormW = 446.42 + initGrNormH = 524.75 + initNorm = 14.236 + SC1 = 0.026889 + SC2 = 0.026512 + SC3 = 2.2711 + W_density = 0.9833 + H_density = 0.6300 iter = 1 - elapsed = 4.7410e-03 - tTotal = 4.7410e-03 - subIterW = 7 - subIterH = 29 - relError = 0.4228 + elapsed = 5.1410e-03 + tTotal = 5.1410e-03 + subIterW = 3 + subIterH = 31 + relError = 0.4343 toshow = scalar structure containing the fields: - initGrNormW = 347.82 - initGrNormH = 420.35 - initNorm = 14.016 - SC1 = 0.016316 - SC2 = 0.015500 - SC3 = 1.1181 - W_density = 0.9200 - H_density = 0.5700 + initGrNormW = 446.42 + initGrNormH = 524.75 + initNorm = 14.236 + SC1 = 0.014096 + SC2 = 0.013561 + SC3 = 1.1859 + W_density = 0.9433 + H_density = 0.5400 iter = 2 - elapsed = 4.5500e-03 - tTotal = 9.2910e-03 - subIterW = 11 - subIterH = 28 - relError = 0.3291 + elapsed = 4.9060e-03 + tTotal = 0.010047 + subIterW = 8 + subIterH = 23 + relError = 0.3465 toshow = scalar structure containing the fields: - initGrNormW = 347.82 - initGrNormH = 420.35 - initNorm = 14.016 - SC1 = 0.010846 - SC2 = 0.010021 - SC3 = 0.7366 - W_density = 0.8933 - H_density = 0.6150 + initGrNormW = 446.42 + initGrNormH = 524.75 + initNorm = 14.236 + SC1 = 9.2548e-03 + SC2 = 8.6440e-03 + SC3 = 0.7577 + W_density = 0.9033 + H_density = 0.6100 iter = 3 - elapsed = 2.9670e-03 - tTotal = 0.012258 - subIterW = 9 - subIterH = 14 - relError = 0.2970 + elapsed = 4.4670e-03 + tTotal = 0.014514 + subIterW = 11 + subIterH = 17 + relError = 0.3163 toshow = scalar structure containing the fields: - initGrNormW = 347.82 - initGrNormH = 420.35 - initNorm = 14.016 - SC1 = 7.7322e-03 - SC2 = 7.0363e-03 - SC3 = 0.5219 - W_density = 0.8800 - H_density = 0.6300 + initGrNormW = 446.42 + initGrNormH = 524.75 + initNorm = 14.236 + SC1 = 6.6260e-03 + SC2 = 6.0032e-03 + SC3 = 0.5295 + W_density = 0.8733 + H_density = 0.6350 iter = 4 - elapsed = 2.3690e-03 - tTotal = 0.014627 - subIterW = 8 - subIterH = 8 - relError = 0.2801 + elapsed = 4.5570e-03 + tTotal = 0.019071 + subIterW = 7 + subIterH = 13 + relError = 0.3020 toshow = scalar structure containing the fields: - initGrNormW = 347.82 - initGrNormH = 420.35 - initNorm = 14.016 - SC1 = 5.8675e-03 - SC2 = 5.1751e-03 - SC3 = 0.3967 - W_density = 0.8600 - H_density = 0.6400 + initGrNormW = 446.42 + initGrNormH = 524.75 + initNorm = 14.236 + SC1 = 5.0197e-03 + SC2 = 4.3972e-03 + SC3 = 0.3979 + W_density = 0.8500 + H_density = 0.6600 iter = 5 - elapsed = 2.0510e-03 - tTotal = 0.016678 + elapsed = 3.6450e-03 + tTotal = 0.022716 subIterW = 7 subIterH = 8 - relError = 0.2714 + relError = 0.2922 toshow = scalar structure containing the fields: - initGrNormW = 347.82 - initGrNormH = 420.35 - initNorm = 14.016 - SC1 = 4.5356e-03 - SC2 = 3.9097e-03 - SC3 = 0.3032 - W_density = 0.8400 - H_density = 0.6450 + initGrNormW = 446.42 + initGrNormH = 524.75 + initNorm = 14.236 + SC1 = 3.9105e-03 + SC2 = 3.3396e-03 + SC3 = 0.3042 + W_density = 0.8267 + H_density = 0.6550 iter = 6 - elapsed = 2.0340e-03 - tTotal = 0.018712 + elapsed = 3.8660e-03 + tTotal = 0.026582 subIterW = 6 - subIterH = 8 - relError = 0.2662 + subIterH = 9 + relError = 0.2852 toshow = scalar structure containing the fields: - initGrNormW = 347.82 - initGrNormH = 420.35 - initNorm = 14.016 - SC1 = 3.6126e-03 - SC2 = 3.0418e-03 - SC3 = 0.2393 - W_density = 0.8233 - H_density = 0.6700 + initGrNormW = 446.42 + initGrNormH = 524.75 + initNorm = 14.236 + SC1 = 3.0787e-03 + SC2 = 2.5738e-03 + SC3 = 0.2352 + W_density = 0.8100 + H_density = 0.6750 iter = 7 - elapsed = 2.2030e-03 - tTotal = 0.020915 - subIterW = 4 - subIterH = 5 - relError = 0.2626 + elapsed = 3.7070e-03 + tTotal = 0.030289 + subIterW = 5 + subIterH = 8 + relError = 0.2808 toshow = scalar structure containing the fields: - initGrNormW = 347.82 - initGrNormH = 420.35 - initNorm = 14.016 - SC1 = 2.9021e-03 - SC2 = 2.3971e-03 - SC3 = 0.1921 - W_density = 0.8133 - H_density = 0.6900 + initGrNormW = 446.42 + initGrNormH = 524.75 + initNorm = 14.236 + SC1 = 2.5085e-03 + SC2 = 2.0369e-03 + SC3 = 0.1932 + W_density = 0.8000 + H_density = 0.6800 iter = 8 - elapsed = 2.0540e-03 - tTotal = 0.022969 - subIterW = 5 + elapsed = 3.5310e-03 + tTotal = 0.033820 + subIterW = 2 subIterH = 3 - relError = 0.2600 + relError = 0.2774 toshow = scalar structure containing the fields: - initGrNormW = 347.82 - initGrNormH = 420.35 - initNorm = 14.016 - SC1 = 2.3191e-03 - SC2 = 1.9063e-03 - SC3 = 0.1528 - W_density = 0.8100 - H_density = 0.7050 + initGrNormW = 446.42 + initGrNormH = 524.75 + initNorm = 14.236 + SC1 = 2.0517e-03 + SC2 = 1.6373e-03 + SC3 = 0.1582 + W_density = 0.7933 + H_density = 0.6850 iter = 9 - elapsed = 1.8910e-03 - tTotal = 0.024860 - subIterW = 2 - subIterH = 4 - relError = 0.2582 + elapsed = 3.6310e-03 + tTotal = 0.037451 + subIterW = 7 + subIterH = 0 + relError = 0.2752 toshow = scalar structure containing the fields: - initGrNormW = 347.82 - initGrNormH = 420.35 - initNorm = 14.016 - SC1 = 1.9070e-03 - SC2 = 1.5447e-03 - SC3 = 0.1262 - W_density = 0.8067 - H_density = 0.7150 + initGrNormW = 446.42 + initGrNormH = 524.75 + initNorm = 14.236 + SC1 = 1.6901e-03 + SC2 = 1.3318e-03 + SC3 = 0.1294 + W_density = 0.7867 + H_density = 0.7300 iter = 10 - elapsed = 1.8260e-03 - tTotal = 0.026686 - subIterW = 0 - subIterH = 4 - relError = 0.2569 + elapsed = 4.3090e-03 + tTotal = 0.041760 + subIterW = 4 + subIterH = 6 + relError = 0.2734 toshow = scalar structure containing the fields: - initGrNormW = 347.82 - initGrNormH = 420.35 - initNorm = 14.016 - SC1 = 1.5874e-03 - SC2 = 1.2763e-03 - SC3 = 0.1043 - W_density = 0.8033 - H_density = 0.7400 + initGrNormW = 446.42 + initGrNormH = 524.75 + initNorm = 14.236 + SC1 = 1.3928e-03 + SC2 = 1.0947e-03 + SC3 = 0.1045 + W_density = 0.7800 + H_density = 0.7350 iter = 11 - elapsed = 1.9710e-03 - tTotal = 0.028657 - subIterW = 3 - subIterH = 4 - relError = 0.2560 + elapsed = 3.5140e-03 + tTotal = 0.045274 + subIterW = 4 + subIterH = 0 + relError = 0.2720 toshow = scalar structure containing the fields: - initGrNormW = 347.82 - initGrNormH = 420.35 - initNorm = 14.016 - SC1 = 1.3459e-03 - SC2 = 1.0714e-03 - SC3 = 0.087699 - W_density = 0.8000 - H_density = 0.7400 + initGrNormW = 446.42 + initGrNormH = 524.75 + initNorm = 14.236 + SC1 = 1.1622e-03 + SC2 = 9.1113e-04 + SC3 = 0.085666 + W_density = 0.7800 + H_density = 0.7250 iter = 12 - elapsed = 1.8870e-03 - tTotal = 0.030544 - subIterW = 5 + elapsed = 3.5000e-03 + tTotal = 0.048774 + subIterW = 2 subIterH = 2 - relError = 0.2553 + relError = 0.2707 toshow = scalar structure containing the fields: - initGrNormW = 347.82 - initGrNormH = 420.35 - initNorm = 14.016 - SC1 = 1.1378e-03 - SC2 = 8.9887e-04 - SC3 = 0.073048 - W_density = 0.7933 - H_density = 0.7350 + initGrNormW = 446.42 + initGrNormH = 524.75 + initNorm = 14.236 + SC1 = 9.8729e-04 + SC2 = 7.7009e-04 + SC3 = 0.072096 + W_density = 0.7833 + H_density = 0.7300 iter = 13 - elapsed = 1.8530e-03 - tTotal = 0.032397 - subIterW = 2 - subIterH = 0 - relError = 0.2547 + elapsed = 3.7480e-03 + tTotal = 0.052522 + subIterW = 5 + subIterH = 4 + relError = 0.2696 toshow = scalar structure containing the fields: - initGrNormW = 347.82 - initGrNormH = 420.35 - initNorm = 14.016 - SC1 = 9.7249e-04 - SC2 = 7.6438e-04 - SC3 = 0.061417 - W_density = 0.7900 - H_density = 0.7450 + initGrNormW = 446.42 + initGrNormH = 524.75 + initNorm = 14.236 + SC1 = 8.5454e-04 + SC2 = 6.5800e-04 + SC3 = 0.061723 + W_density = 0.7767 + H_density = 0.7300 iter = 14 - elapsed = 1.8410e-03 - tTotal = 0.034238 - subIterW = 0 + elapsed = 3.5670e-03 + tTotal = 0.056089 + subIterW = 4 subIterH = 2 - relError = 0.2541 + relError = 0.2686 toshow = scalar structure containing the fields: - initGrNormW = 347.82 - initGrNormH = 420.35 - initNorm = 14.016 - SC1 = 8.3523e-04 - SC2 = 6.5482e-04 - SC3 = 0.051935 - W_density = 0.7900 - H_density = 0.7500 + initGrNormW = 446.42 + initGrNormH = 524.75 + initNorm = 14.236 + SC1 = 7.3982e-04 + SC2 = 5.6966e-04 + SC3 = 0.052655 + W_density = 0.7767 + H_density = 0.7400 iter = 15 - elapsed = 1.7990e-03 - tTotal = 0.036037 + elapsed = 3.4720e-03 + tTotal = 0.059561 subIterW = 2 - subIterH = 0 - relError = 0.2536 + subIterH = 2 + relError = 0.2678 toshow = scalar structure containing the fields: - initGrNormW = 347.82 - initGrNormH = 420.35 - initNorm = 14.016 - SC1 = 7.1435e-04 - SC2 = 5.5862e-04 - SC3 = 0.043666 - W_density = 0.7867 - H_density = 0.7450 + initGrNormW = 446.42 + initGrNormH = 524.75 + initNorm = 14.236 + SC1 = 6.5056e-04 + SC2 = 4.9572e-04 + SC3 = 0.045182 + W_density = 0.7633 + H_density = 0.7400 iter = 16 - elapsed = 1.7530e-03 - tTotal = 0.037790 - subIterW = 0 + elapsed = 3.1440e-03 + tTotal = 0.062705 + subIterW = 6 subIterH = 0 - relError = 0.2532 + relError = 0.2670 toshow = scalar structure containing the fields: - initGrNormW = 347.82 - initGrNormH = 420.35 - initNorm = 14.016 - SC1 = 6.1119e-04 - SC2 = 4.7795e-04 - SC3 = 0.037060 - W_density = 0.7867 - H_density = 0.7500 + initGrNormW = 446.42 + initGrNormH = 524.75 + initNorm = 14.236 + SC1 = 5.7024e-04 + SC2 = 4.3452e-04 + SC3 = 0.022305 + W_density = 0.7633 + H_density = 0.7350 iter = 17 - elapsed = 1.7890e-03 - tTotal = 0.039579 - subIterW = 2 - subIterH = 0 - relError = 0.2528 + elapsed = 2.8830e-03 + tTotal = 0.065588 + subIterW = 0 + subIterH = 3 + relError = 0.2662 toshow = scalar structure containing the fields: - initGrNormW = 347.82 - initGrNormH = 420.35 - initNorm = 14.016 - SC1 = 5.2594e-04 - SC2 = 4.0708e-04 - SC3 = 0.031557 - W_density = 0.7767 - H_density = 0.7450 + initGrNormW = 446.42 + initGrNormH = 524.75 + initNorm = 14.236 + SC1 = 4.9659e-04 + SC2 = 3.7940e-04 + SC3 = 0.033641 + W_density = 0.7600 + H_density = 0.7400 iter = 18 - elapsed = 1.8680e-03 - tTotal = 0.041447 - subIterW = 3 - subIterH = 0 - relError = 0.2525 + elapsed = 3.5710e-03 + tTotal = 0.069159 + subIterW = 0 + subIterH = 5 + relError = 0.2655 toshow = scalar structure containing the fields: - initGrNormW = 347.82 - initGrNormH = 420.35 - initNorm = 14.016 - SC1 = 4.5668e-04 - SC2 = 3.4891e-04 - SC3 = 0.027423 - W_density = 0.7700 - H_density = 0.7350 + initGrNormW = 446.42 + initGrNormH = 524.75 + initNorm = 14.236 + SC1 = 4.5144e-04 + SC2 = 3.4400e-04 + SC3 = 0.030174 + W_density = 0.7600 + H_density = 0.7450 iter = 19 - elapsed = 1.8410e-03 - tTotal = 0.043288 + elapsed = 3.5890e-03 + tTotal = 0.072748 subIterW = 2 - subIterH = 2 - relError = 0.2523 + subIterH = 3 + relError = 0.2649 toshow = scalar structure containing the fields: - initGrNormW = 347.82 - initGrNormH = 420.35 - initNorm = 14.016 - SC1 = 3.9644e-04 - SC2 = 3.0288e-04 - SC3 = 0.013966 - W_density = 0.7700 - H_density = 0.7400 + initGrNormW = 446.42 + initGrNormH = 524.75 + initNorm = 14.236 + SC1 = 4.1322e-04 + SC2 = 3.1157e-04 + SC3 = 0.026756 + W_density = 0.7567 + H_density = 0.7350 iter = 20 - elapsed = 1.4590e-03 - tTotal = 0.044747 - subIterW = 0 - subIterH = 0 - relError = 0.2521 + elapsed = 3.6720e-03 + tTotal = 0.076420 + subIterW = 5 + subIterH = 5 + relError = 0.2645 toshow = scalar structure containing the fields: - initGrNormW = 347.82 - initGrNormH = 420.35 - initNorm = 14.016 - SC1 = 3.4890e-04 - SC2 = 2.6517e-04 - SC3 = 0.020643 - W_density = 0.7633 - H_density = 0.7400 + initGrNormW = 446.42 + initGrNormH = 524.75 + initNorm = 14.236 + SC1 = 3.6127e-04 + SC2 = 2.7312e-04 + SC3 = 0.023599 + W_density = 0.7667 + H_density = 0.7250 iter = 21 - elapsed = 1.5050e-03 - tTotal = 0.046252 - subIterW = 2 - subIterH = 0 - relError = 0.2520 + elapsed = 3.4980e-03 + tTotal = 0.079918 + subIterW = 3 + subIterH = 2 + relError = 0.2641 toshow = scalar structure containing the fields: - initGrNormW = 347.82 - initGrNormH = 420.35 - initNorm = 14.016 - SC1 = 3.0110e-04 - SC2 = 2.2703e-04 - SC3 = 0.017888 - W_density = 0.7600 - H_density = 0.7350 + initGrNormW = 446.42 + initGrNormH = 524.75 + initNorm = 14.236 + SC1 = 3.2743e-04 + SC2 = 2.4819e-04 + SC3 = 0.021205 + W_density = 0.7700 + H_density = 0.7300 iter = 22 - elapsed = 1.9900e-03 - tTotal = 0.048242 - subIterW = 2 - subIterH = 3 - relError = 0.2519 + elapsed = 3.2540e-03 + tTotal = 0.083172 + subIterW = 0 + subIterH = 0 + relError = 0.2639 toshow = scalar structure containing the fields: - initGrNormW = 347.82 - initGrNormH = 420.35 - initNorm = 14.016 - SC1 = 2.6308e-04 - SC2 = 1.9784e-04 - SC3 = 0.015560 - W_density = 0.7600 - H_density = 0.7350 + initGrNormW = 446.42 + initGrNormH = 524.75 + initNorm = 14.236 + SC1 = 3.0573e-04 + SC2 = 2.3113e-04 + SC3 = 0.019545 + W_density = 0.7667 + H_density = 0.7300 iter = 23 - elapsed = 1.6040e-03 - tTotal = 0.049846 - subIterW = 2 - subIterH = 0 - relError = 0.2519 + elapsed = 3.3610e-03 + tTotal = 0.086533 + subIterW = 0 + subIterH = 2 + relError = 0.2636 Stop: tolerance reached. final = @@ -3425,19 +3476,19 @@ scalar structure containing the fields: iterations = 23 - elapsed_total = 0.049846 - relative_error = 0.2519 - W_density = 0.7600 - H_density = 0.7350 + elapsed_total = 0.086533 + relative_error = 0.2636 + W_density = 0.7667 + H_density = 0.7300 ***** test [W,H,iter,HIS]=nmf_bpas(A,k,'verbose','nnlssolver','as'); - initGrNormW = 418.38 - initGrNormH = 521.76 - initNorm = 14.016 - SC1 = 1.3376 - SC2 = 668.78 - SC3 = 0.5167 + initGrNormW = 367.18 + initGrNormH = 477.66 + initNorm = 14.236 + SC1 = 1.2050 + SC2 = 602.48 + SC3 = 0.5051 W_density = 1 H_density = 1 par = @@ -3448,8 +3499,8 @@ n = 20 type = regularized nnls_solver = as - alpha = 0.4979 - beta = 0.4979 + alpha = 0.5052 + beta = 0.5052 max_iter = 100 min_iter = 20 max_time = 1000 @@ -3461,438 +3512,438 @@ scalar structure containing the fields: - initGrNormW = 418.38 - initGrNormH = 521.76 - initNorm = 14.016 - SC1 = 0.028999 - SC2 = 0.028477 - SC3 = 1.3111 - W_density = 0.9733 - H_density = 0.6750 + initGrNormW = 367.18 + initGrNormH = 477.66 + initNorm = 14.236 + SC1 = 0.031303 + SC2 = 0.030864 + SC3 = 1.3116 + W_density = 0.9767 + H_density = 0.6700 iter = 1 - elapsed = 2.1050e-03 - tTotal = 2.1050e-03 - subIterW = 6 - subIterH = 29 - relError = 0.4351 + elapsed = 4.2930e-03 + tTotal = 4.2930e-03 + subIterW = 5 + subIterH = 40 + relError = 0.4452 toshow = scalar structure containing the fields: - initGrNormW = 418.38 - initGrNormH = 521.76 - initNorm = 14.016 - SC1 = 0.015405 - SC2 = 0.014604 - SC3 = 0.6872 + initGrNormW = 367.18 + initGrNormH = 477.66 + initNorm = 14.236 + SC1 = 0.016626 + SC2 = 0.015795 + SC3 = 0.6780 W_density = 0.9233 H_density = 0.5350 iter = 2 - elapsed = 2.2870e-03 - tTotal = 4.3920e-03 - subIterW = 15 + elapsed = 4.0780e-03 + tTotal = 8.3710e-03 + subIterW = 12 subIterH = 40 - relError = 0.3442 + relError = 0.3461 toshow = scalar structure containing the fields: - initGrNormW = 418.38 - initGrNormH = 521.76 - initNorm = 14.016 - SC1 = 9.8651e-03 - SC2 = 9.1548e-03 - SC3 = 0.4437 - W_density = 0.9033 - H_density = 0.5650 + initGrNormW = 367.18 + initGrNormH = 477.66 + initNorm = 14.236 + SC1 = 0.010732 + SC2 = 0.010088 + SC3 = 0.4316 + W_density = 0.9133 + H_density = 0.5850 iter = 3 - elapsed = 2.8570e-03 - tTotal = 7.2490e-03 - subIterW = 18 - subIterH = 42 - relError = 0.3101 + elapsed = 4.1850e-03 + tTotal = 0.012556 + subIterW = 14 + subIterH = 36 + relError = 0.3139 toshow = scalar structure containing the fields: - initGrNormW = 418.38 - initGrNormH = 521.76 - initNorm = 14.016 - SC1 = 7.0829e-03 - SC2 = 6.3888e-03 - SC3 = 0.3224 - W_density = 0.8833 - H_density = 0.6300 + initGrNormW = 367.18 + initGrNormH = 477.66 + initNorm = 14.236 + SC1 = 7.5496e-03 + SC2 = 6.9607e-03 + SC3 = 0.3057 + W_density = 0.8933 + H_density = 0.6350 iter = 4 - elapsed = 2.4610e-03 - tTotal = 9.7100e-03 - subIterW = 20 - subIterH = 35 - relError = 0.2939 + elapsed = 3.8520e-03 + tTotal = 0.016408 + subIterW = 18 + subIterH = 31 + relError = 0.2977 toshow = scalar structure containing the fields: - initGrNormW = 418.38 - initGrNormH = 521.76 - initNorm = 14.016 - SC1 = 5.2566e-03 - SC2 = 4.6679e-03 - SC3 = 0.2362 - W_density = 0.8667 - H_density = 0.6650 + initGrNormW = 367.18 + initGrNormH = 477.66 + initNorm = 14.236 + SC1 = 5.6333e-03 + SC2 = 5.0249e-03 + SC3 = 0.2384 + W_density = 0.8633 + H_density = 0.6600 iter = 5 - elapsed = 2.4050e-03 - tTotal = 0.012115 - subIterW = 21 + elapsed = 4.5460e-03 + tTotal = 0.020954 + subIterW = 23 subIterH = 32 - relError = 0.2833 + relError = 0.2883 toshow = scalar structure containing the fields: - initGrNormW = 418.38 - initGrNormH = 521.76 - initNorm = 14.016 - SC1 = 4.0455e-03 - SC2 = 3.4953e-03 - SC3 = 0.1782 - W_density = 0.8467 - H_density = 0.6800 + initGrNormW = 367.18 + initGrNormH = 477.66 + initNorm = 14.236 + SC1 = 4.3579e-03 + SC2 = 3.7827e-03 + SC3 = 0.1754 + W_density = 0.8367 + H_density = 0.6700 iter = 6 - elapsed = 2.3740e-03 - tTotal = 0.014489 - subIterW = 25 - subIterH = 33 - relError = 0.2758 + elapsed = 4.9440e-03 + tTotal = 0.025898 + subIterW = 26 + subIterH = 31 + relError = 0.2818 toshow = scalar structure containing the fields: - initGrNormW = 418.38 - initGrNormH = 521.76 - initNorm = 14.016 - SC1 = 3.1456e-03 - SC2 = 2.6611e-03 - SC3 = 0.1369 - W_density = 0.8367 - H_density = 0.7050 + initGrNormW = 367.18 + initGrNormH = 477.66 + initNorm = 14.236 + SC1 = 3.4717e-03 + SC2 = 2.9231e-03 + SC3 = 0.1409 + W_density = 0.8067 + H_density = 0.7150 iter = 7 - elapsed = 2.2480e-03 - tTotal = 0.016737 - subIterW = 25 - subIterH = 31 - relError = 0.2704 + elapsed = 4.5750e-03 + tTotal = 0.030473 + subIterW = 30 + subIterH = 26 + relError = 0.2773 toshow = scalar structure containing the fields: - initGrNormW = 418.38 - initGrNormH = 521.76 - initNorm = 14.016 - SC1 = 2.5100e-03 - SC2 = 2.0883e-03 - SC3 = 0.1145 - W_density = 0.8200 - H_density = 0.7300 + initGrNormW = 367.18 + initGrNormH = 477.66 + initNorm = 14.236 + SC1 = 2.7981e-03 + SC2 = 2.2889e-03 + SC3 = 0.1131 + W_density = 0.7833 + H_density = 0.7200 iter = 8 - elapsed = 2.4360e-03 - tTotal = 0.019173 - subIterW = 26 - subIterH = 25 - relError = 0.2667 + elapsed = 5.0340e-03 + tTotal = 0.035507 + subIterW = 34 + subIterH = 29 + relError = 0.2743 toshow = scalar structure containing the fields: - initGrNormW = 418.38 - initGrNormH = 521.76 - initNorm = 14.016 - SC1 = 2.0484e-03 - SC2 = 1.6838e-03 - SC3 = 0.092192 - W_density = 0.8067 + initGrNormW = 367.18 + initGrNormH = 477.66 + initNorm = 14.236 + SC1 = 2.2892e-03 + SC2 = 1.8268e-03 + SC3 = 0.090920 + W_density = 0.7767 H_density = 0.7450 iter = 9 - elapsed = 2.4390e-03 - tTotal = 0.021612 - subIterW = 28 + elapsed = 4.7900e-03 + tTotal = 0.040297 + subIterW = 32 subIterH = 25 - relError = 0.2642 + relError = 0.2723 toshow = scalar structure containing the fields: - initGrNormW = 418.38 - initGrNormH = 521.76 - initNorm = 14.016 - SC1 = 1.7074e-03 - SC2 = 1.3830e-03 - SC3 = 0.073033 - W_density = 0.8000 - H_density = 0.7300 + initGrNormW = 367.18 + initGrNormH = 477.66 + initNorm = 14.236 + SC1 = 1.8917e-03 + SC2 = 1.4982e-03 + SC3 = 0.075599 + W_density = 0.7767 + H_density = 0.7350 iter = 10 - elapsed = 2.1980e-03 - tTotal = 0.023810 - subIterW = 30 - subIterH = 28 - relError = 0.2623 + elapsed = 4.3730e-03 + tTotal = 0.044670 + subIterW = 33 + subIterH = 27 + relError = 0.2707 toshow = scalar structure containing the fields: - initGrNormW = 418.38 - initGrNormH = 521.76 - initNorm = 14.016 - SC1 = 1.4364e-03 - SC2 = 1.1433e-03 - SC3 = 0.064257 - W_density = 0.7867 - H_density = 0.7250 + initGrNormW = 367.18 + initGrNormH = 477.66 + initNorm = 14.236 + SC1 = 1.5941e-03 + SC2 = 1.2562e-03 + SC3 = 0.065448 + W_density = 0.7767 + H_density = 0.7350 iter = 11 - elapsed = 2.8120e-03 - tTotal = 0.026622 - subIterW = 34 - subIterH = 26 - relError = 0.2610 + elapsed = 4.7430e-03 + tTotal = 0.049413 + subIterW = 33 + subIterH = 28 + relError = 0.2694 toshow = scalar structure containing the fields: - initGrNormW = 418.38 - initGrNormH = 521.76 - initNorm = 14.016 - SC1 = 1.1886e-03 - SC2 = 9.5091e-04 - SC3 = 0.050525 - W_density = 0.7900 - H_density = 0.7250 + initGrNormW = 367.18 + initGrNormH = 477.66 + initNorm = 14.236 + SC1 = 1.3704e-03 + SC2 = 1.0635e-03 + SC3 = 0.053249 + W_density = 0.7600 + H_density = 0.7350 iter = 12 - elapsed = 2.7750e-03 - tTotal = 0.029397 - subIterW = 34 - subIterH = 28 - relError = 0.2600 + elapsed = 4.8770e-03 + tTotal = 0.054290 + subIterW = 37 + subIterH = 25 + relError = 0.2683 toshow = scalar structure containing the fields: - initGrNormW = 418.38 - initGrNormH = 521.76 - initNorm = 14.016 - SC1 = 1.0056e-03 - SC2 = 8.0243e-04 - SC3 = 0.041371 - W_density = 0.7933 - H_density = 0.7350 + initGrNormW = 367.18 + initGrNormH = 477.66 + initNorm = 14.236 + SC1 = 1.1597e-03 + SC2 = 8.9529e-04 + SC3 = 0.044969 + W_density = 0.7600 + H_density = 0.7400 iter = 13 - elapsed = 2.4560e-03 - tTotal = 0.031853 - subIterW = 35 - subIterH = 25 - relError = 0.2591 + elapsed = 5.4610e-03 + tTotal = 0.059751 + subIterW = 34 + subIterH = 28 + relError = 0.2674 toshow = scalar structure containing the fields: - initGrNormW = 418.38 - initGrNormH = 521.76 - initNorm = 14.016 - SC1 = 8.5988e-04 - SC2 = 6.7587e-04 - SC3 = 0.037870 - W_density = 0.7800 - H_density = 0.7400 + initGrNormW = 367.18 + initGrNormH = 477.66 + initNorm = 14.236 + SC1 = 9.7287e-04 + SC2 = 7.5495e-04 + SC3 = 0.037819 + W_density = 0.7633 + H_density = 0.7450 iter = 14 - elapsed = 2.3600e-03 - tTotal = 0.034213 - subIterW = 37 - subIterH = 26 - relError = 0.2584 + elapsed = 4.2660e-03 + tTotal = 0.064017 + subIterW = 32 + subIterH = 25 + relError = 0.2667 toshow = scalar structure containing the fields: - initGrNormW = 418.38 - initGrNormH = 521.76 - initNorm = 14.016 - SC1 = 7.2386e-04 - SC2 = 5.6172e-04 - SC3 = 0.030755 - W_density = 0.7700 + initGrNormW = 367.18 + initGrNormH = 477.66 + initNorm = 14.236 + SC1 = 8.3661e-04 + SC2 = 6.3917e-04 + SC3 = 0.031807 + W_density = 0.7567 H_density = 0.7450 iter = 15 - elapsed = 2.5130e-03 - tTotal = 0.036726 - subIterW = 40 + elapsed = 5.1830e-03 + tTotal = 0.069200 + subIterW = 34 subIterH = 25 - relError = 0.2580 + relError = 0.2661 toshow = scalar structure containing the fields: - initGrNormW = 418.38 - initGrNormH = 521.76 - initNorm = 14.016 - SC1 = 6.0696e-04 - SC2 = 4.6736e-04 - SC3 = 0.024831 - W_density = 0.7633 + initGrNormW = 367.18 + initGrNormH = 477.66 + initNorm = 14.236 + SC1 = 7.0152e-04 + SC2 = 5.3175e-04 + SC3 = 0.027307 + W_density = 0.7500 H_density = 0.7400 iter = 16 - elapsed = 2.1240e-03 - tTotal = 0.038850 + elapsed = 4.9650e-03 + tTotal = 0.074165 subIterW = 35 - subIterH = 28 - relError = 0.2577 + subIterH = 25 + relError = 0.2657 toshow = scalar structure containing the fields: - initGrNormW = 418.38 - initGrNormH = 521.76 - initNorm = 14.016 - SC1 = 5.1250e-04 - SC2 = 3.9462e-04 - SC3 = 0.021138 - W_density = 0.7600 - H_density = 0.7450 + initGrNormW = 367.18 + initGrNormH = 477.66 + initNorm = 14.236 + SC1 = 5.9014e-04 + SC2 = 4.4614e-04 + SC3 = 0.021505 + W_density = 0.7500 + H_density = 0.7400 iter = 17 - elapsed = 2.3720e-03 - tTotal = 0.041222 + elapsed = 4.7650e-03 + tTotal = 0.078930 subIterW = 35 - subIterH = 26 - relError = 0.2575 + subIterH = 25 + relError = 0.2654 toshow = scalar structure containing the fields: - initGrNormW = 418.38 - initGrNormH = 521.76 - initNorm = 14.016 - SC1 = 4.4182e-04 - SC2 = 3.4020e-04 - SC3 = 0.017401 - W_density = 0.7633 - H_density = 0.7500 + initGrNormW = 367.18 + initGrNormH = 477.66 + initNorm = 14.236 + SC1 = 5.1160e-04 + SC2 = 3.8473e-04 + SC3 = 0.018907 + W_density = 0.7500 + H_density = 0.7400 iter = 18 - elapsed = 2.2940e-03 - tTotal = 0.043516 - subIterW = 34 - subIterH = 24 - relError = 0.2573 + elapsed = 4.1320e-03 + tTotal = 0.083062 + subIterW = 36 + subIterH = 25 + relError = 0.2652 toshow = scalar structure containing the fields: - initGrNormW = 418.38 - initGrNormH = 521.76 - initNorm = 14.016 - SC1 = 3.9291e-04 - SC2 = 3.0018e-04 - SC3 = 0.015186 - W_density = 0.7600 - H_density = 0.7550 + initGrNormW = 367.18 + initGrNormH = 477.66 + initNorm = 14.236 + SC1 = 4.5032e-04 + SC2 = 3.3774e-04 + SC3 = 0.016188 + W_density = 0.7467 + H_density = 0.7400 iter = 19 - elapsed = 2.1990e-03 - tTotal = 0.045715 - subIterW = 35 - subIterH = 23 - relError = 0.2571 + elapsed = 4.3220e-03 + tTotal = 0.087384 + subIterW = 37 + subIterH = 26 + relError = 0.2649 toshow = scalar structure containing the fields: - initGrNormW = 418.38 - initGrNormH = 521.76 - initNorm = 14.016 - SC1 = 3.5502e-04 - SC2 = 2.6697e-04 - SC3 = 0.013753 - W_density = 0.7433 - H_density = 0.7550 + initGrNormW = 367.18 + initGrNormH = 477.66 + initNorm = 14.236 + SC1 = 4.0435e-04 + SC2 = 3.0246e-04 + SC3 = 0.013869 + W_density = 0.7467 + H_density = 0.7400 iter = 20 - elapsed = 2.2690e-03 - tTotal = 0.047984 - subIterW = 38 - subIterH = 22 - relError = 0.2570 + elapsed = 4.4880e-03 + tTotal = 0.091872 + subIterW = 36 + subIterH = 27 + relError = 0.2647 toshow = scalar structure containing the fields: - initGrNormW = 418.38 - initGrNormH = 521.76 - initNorm = 14.016 - SC1 = 3.2140e-04 - SC2 = 2.3976e-04 - SC3 = 0.012512 - W_density = 0.7400 - H_density = 0.7550 + initGrNormW = 367.18 + initGrNormH = 477.66 + initNorm = 14.236 + SC1 = 3.6972e-04 + SC2 = 2.7655e-04 + SC3 = 0.013085 + W_density = 0.7467 + H_density = 0.7350 iter = 21 - elapsed = 2.4140e-03 - tTotal = 0.050398 - subIterW = 38 - subIterH = 21 - relError = 0.2569 + elapsed = 4.5950e-03 + tTotal = 0.096467 + subIterW = 37 + subIterH = 28 + relError = 0.2645 toshow = scalar structure containing the fields: - initGrNormW = 418.38 - initGrNormH = 521.76 - initNorm = 14.016 - SC1 = 2.9325e-04 - SC2 = 2.1876e-04 - SC3 = 0.010666 - W_density = 0.7433 - H_density = 0.7450 + initGrNormW = 367.18 + initGrNormH = 477.66 + initNorm = 14.236 + SC1 = 3.4478e-04 + SC2 = 2.5721e-04 + SC3 = 0.012543 + W_density = 0.7500 + H_density = 0.7300 iter = 22 - elapsed = 2.1210e-03 - tTotal = 0.052519 - subIterW = 36 - subIterH = 24 - relError = 0.2568 + elapsed = 4.4670e-03 + tTotal = 0.1009 + subIterW = 38 + subIterH = 29 + relError = 0.2642 toshow = scalar structure containing the fields: - initGrNormW = 418.38 - initGrNormH = 521.76 - initNorm = 14.016 - SC1 = 2.7487e-04 - SC2 = 2.0450e-04 - SC3 = 0.010906 - W_density = 0.7433 - H_density = 0.7400 + initGrNormW = 367.18 + initGrNormH = 477.66 + initNorm = 14.236 + SC1 = 3.2792e-04 + SC2 = 2.4397e-04 + SC3 = 0.010898 + W_density = 0.7500 + H_density = 0.7300 iter = 23 - elapsed = 1.9880e-03 - tTotal = 0.054507 - subIterW = 34 - subIterH = 24 - relError = 0.2568 + elapsed = 4.7650e-03 + tTotal = 0.1057 + subIterW = 39 + subIterH = 30 + relError = 0.2640 Stop: tolerance reached. final = @@ -3900,19 +3951,19 @@ scalar structure containing the fields: iterations = 23 - elapsed_total = 0.054507 - relative_error = 0.2568 - W_density = 0.7433 - H_density = 0.7400 + elapsed_total = 0.1057 + relative_error = 0.2640 + W_density = 0.7500 + H_density = 0.7300 ***** test [W,H,iter,HIS]=nmf_bpas(A,k,'verbose','type','sparse'); - initGrNormW = 406.79 - initGrNormH = 524.35 - initNorm = 14.016 - SC1 = 1.3273 - SC2 = 663.64 - SC3 = 0.5116 + initGrNormW = 336.80 + initGrNormH = 432.64 + initNorm = 14.236 + SC1 = 1.0966 + SC2 = 548.28 + SC3 = 0.4811 W_density = 1 H_density = 1 par = @@ -3923,8 +3974,8 @@ n = 20 type = sparse nnls_solver = bp - alpha = 0.4979 - beta = 0.4979 + alpha = 0.5052 + beta = 0.5052 max_iter = 100 min_iter = 20 max_time = 1000 @@ -3936,438 +3987,438 @@ scalar structure containing the fields: - initGrNormW = 406.79 - initGrNormH = 524.35 - initNorm = 14.016 - SC1 = 0.020834 - SC2 = 0.018792 - SC3 = 1.5166 - W_density = 0.9300 - H_density = 0.5850 + initGrNormW = 336.80 + initGrNormH = 432.64 + initNorm = 14.236 + SC1 = 0.024469 + SC2 = 0.022805 + SC3 = 1.7041 + W_density = 0.9667 + H_density = 0.5550 iter = 1 - elapsed = 3.0620e-03 - tTotal = 3.0620e-03 - subIterW = 10 - subIterH = 31 - relError = 0.4295 + elapsed = 4.7870e-03 + tTotal = 4.7870e-03 + subIterW = 6 + subIterH = 35 + relError = 0.4270 toshow = scalar structure containing the fields: - initGrNormW = 406.79 - initGrNormH = 524.35 - initNorm = 14.016 - SC1 = 0.010485 - SC2 = 8.4303e-03 - SC3 = 0.7733 - W_density = 0.9067 - H_density = 0.3600 + initGrNormW = 336.80 + initGrNormH = 432.64 + initNorm = 14.236 + SC1 = 0.010919 + SC2 = 9.3245e-03 + SC3 = 0.7306 + W_density = 0.9433 + H_density = 0.3500 iter = 2 - elapsed = 3.2300e-03 - tTotal = 6.2920e-03 - subIterW = 8 - subIterH = 34 - relError = 0.3565 + elapsed = 5.4660e-03 + tTotal = 0.010253 + subIterW = 7 + subIterH = 33 + relError = 0.3505 toshow = scalar structure containing the fields: - initGrNormW = 406.79 - initGrNormH = 524.35 - initNorm = 14.016 - SC1 = 6.5342e-03 - SC2 = 4.8092e-03 - SC3 = 0.5095 - W_density = 0.8967 - H_density = 0.3000 + initGrNormW = 336.80 + initGrNormH = 432.64 + initNorm = 14.236 + SC1 = 6.4026e-03 + SC2 = 5.0324e-03 + SC3 = 0.4350 + W_density = 0.9400 + H_density = 0.3300 iter = 3 - elapsed = 2.6680e-03 - tTotal = 8.9600e-03 - subIterW = 9 - subIterH = 13 - relError = 0.3332 + elapsed = 4.2090e-03 + tTotal = 0.014462 + subIterW = 6 + subIterH = 15 + relError = 0.3245 toshow = scalar structure containing the fields: - initGrNormW = 406.79 - initGrNormH = 524.35 - initNorm = 14.016 - SC1 = 4.4536e-03 - SC2 = 3.0374e-03 - SC3 = 0.3830 - W_density = 0.8833 - H_density = 0.2850 + initGrNormW = 336.80 + initGrNormH = 432.64 + initNorm = 14.236 + SC1 = 4.0978e-03 + SC2 = 2.9586e-03 + SC3 = 0.3084 + W_density = 0.9333 + H_density = 0.2900 iter = 4 - elapsed = 1.9420e-03 - tTotal = 0.010902 - subIterW = 7 - subIterH = 9 - relError = 0.3236 + elapsed = 3.5240e-03 + tTotal = 0.017986 + subIterW = 4 + subIterH = 12 + relError = 0.3159 toshow = scalar structure containing the fields: - initGrNormW = 406.79 - initGrNormH = 524.35 - initNorm = 14.016 - SC1 = 3.1471e-03 - SC2 = 2.0897e-03 - SC3 = 0.2723 - W_density = 0.8900 - H_density = 0.2700 + initGrNormW = 336.80 + initGrNormH = 432.64 + initNorm = 14.236 + SC1 = 2.7271e-03 + SC2 = 1.8763e-03 + SC3 = 0.2119 + W_density = 0.9200 + H_density = 0.2750 iter = 5 - elapsed = 1.7930e-03 - tTotal = 0.012695 + elapsed = 3.4570e-03 + tTotal = 0.021443 subIterW = 4 subIterH = 5 - relError = 0.3207 + relError = 0.3123 toshow = scalar structure containing the fields: - initGrNormW = 406.79 - initGrNormH = 524.35 - initNorm = 14.016 - SC1 = 2.3539e-03 - SC2 = 1.4971e-03 - SC3 = 0.2063 - W_density = 0.8700 - H_density = 0.2550 + initGrNormW = 336.80 + initGrNormH = 432.64 + initNorm = 14.236 + SC1 = 1.9769e-03 + SC2 = 1.3364e-03 + SC3 = 0.1468 + W_density = 0.9133 + H_density = 0.2700 iter = 6 - elapsed = 2.2660e-03 - tTotal = 0.014961 - subIterW = 6 - subIterH = 6 - relError = 0.3203 + elapsed = 3.2040e-03 + tTotal = 0.024647 + subIterW = 5 + subIterH = 4 + relError = 0.3114 toshow = scalar structure containing the fields: - initGrNormW = 406.79 - initGrNormH = 524.35 - initNorm = 14.016 - SC1 = 1.8277e-03 - SC2 = 1.1368e-03 - SC3 = 0.1547 - W_density = 0.8633 - H_density = 0.2400 + initGrNormW = 336.80 + initGrNormH = 432.64 + initNorm = 14.236 + SC1 = 1.6154e-03 + SC2 = 1.0597e-03 + SC3 = 0.1244 + W_density = 0.9033 + H_density = 0.2450 iter = 7 - elapsed = 1.6910e-03 - tTotal = 0.016652 - subIterW = 2 - subIterH = 3 - relError = 0.3206 + elapsed = 2.8700e-03 + tTotal = 0.027517 + subIterW = 3 + subIterH = 5 + relError = 0.3122 toshow = scalar structure containing the fields: - initGrNormW = 406.79 - initGrNormH = 524.35 - initNorm = 14.016 - SC1 = 1.5260e-03 - SC2 = 9.3394e-04 - SC3 = 0.1298 - W_density = 0.8567 - H_density = 0.2350 + initGrNormW = 336.80 + initGrNormH = 432.64 + initNorm = 14.236 + SC1 = 1.4854e-03 + SC2 = 9.5957e-04 + SC3 = 0.1148 + W_density = 0.8933 + H_density = 0.2450 iter = 8 - elapsed = 1.6030e-03 - tTotal = 0.018255 - subIterW = 0 - subIterH = 0 - relError = 0.3210 + elapsed = 2.8540e-03 + tTotal = 0.030371 + subIterW = 3 + subIterH = 2 + relError = 0.3136 toshow = scalar structure containing the fields: - initGrNormW = 406.79 - initGrNormH = 524.35 - initNorm = 14.016 - SC1 = 1.3717e-03 - SC2 = 8.3402e-04 - SC3 = 0.1172 - W_density = 0.8533 + initGrNormW = 336.80 + initGrNormH = 432.64 + initNorm = 14.236 + SC1 = 1.6509e-03 + SC2 = 1.0533e-03 + SC3 = 0.1232 + W_density = 0.8900 H_density = 0.2400 iter = 9 - elapsed = 1.5970e-03 - tTotal = 0.019852 - subIterW = 0 + elapsed = 2.7260e-03 + tTotal = 0.033097 + subIterW = 3 subIterH = 0 - relError = 0.3216 + relError = 0.3151 toshow = scalar structure containing the fields: - initGrNormW = 406.79 - initGrNormH = 524.35 - initNorm = 14.016 - SC1 = 1.3106e-03 - SC2 = 7.8375e-04 - SC3 = 0.1103 - W_density = 0.8400 - H_density = 0.2350 + initGrNormW = 336.80 + initGrNormH = 432.64 + initNorm = 14.236 + SC1 = 1.6693e-03 + SC2 = 1.0517e-03 + SC3 = 0.1042 + W_density = 0.8800 + H_density = 0.2400 iter = 10 - elapsed = 2.0400e-03 - tTotal = 0.021892 - subIterW = 3 - subIterH = 0 - relError = 0.3224 + elapsed = 2.8780e-03 + tTotal = 0.035975 + subIterW = 5 + subIterH = 3 + relError = 0.3164 toshow = scalar structure containing the fields: - initGrNormW = 406.79 - initGrNormH = 524.35 - initNorm = 14.016 - SC1 = 1.2589e-03 - SC2 = 7.5032e-04 - SC3 = 0.1067 - W_density = 0.8367 - H_density = 0.2300 + initGrNormW = 336.80 + initGrNormH = 432.64 + initNorm = 14.236 + SC1 = 1.1832e-03 + SC2 = 7.1940e-04 + SC3 = 0.083314 + W_density = 0.8533 + H_density = 0.2350 iter = 11 - elapsed = 1.6590e-03 - tTotal = 0.023551 - subIterW = 0 - subIterH = 0 - relError = 0.3234 + elapsed = 2.9350e-03 + tTotal = 0.038910 + subIterW = 7 + subIterH = 3 + relError = 0.3177 toshow = scalar structure containing the fields: - initGrNormW = 406.79 - initGrNormH = 524.35 - initNorm = 14.016 - SC1 = 1.2964e-03 - SC2 = 7.7526e-04 - SC3 = 0.1094 - W_density = 0.8400 - H_density = 0.2300 + initGrNormW = 336.80 + initGrNormH = 432.64 + initNorm = 14.236 + SC1 = 8.5756e-04 + SC2 = 5.2311e-04 + SC3 = 0.061302 + W_density = 0.8500 + H_density = 0.2400 iter = 12 - elapsed = 1.2220e-03 - tTotal = 0.024773 + elapsed = 2.7730e-03 + tTotal = 0.041683 subIterW = 0 subIterH = 0 - relError = 0.3244 + relError = 0.3188 toshow = scalar structure containing the fields: - initGrNormW = 406.79 - initGrNormH = 524.35 - initNorm = 14.016 - SC1 = 1.4127e-03 - SC2 = 8.4480e-04 - SC3 = 0.027995 - W_density = 0.8400 - H_density = 0.2300 + initGrNormW = 336.80 + initGrNormH = 432.64 + initNorm = 14.236 + SC1 = 7.0218e-04 + SC2 = 4.2833e-04 + SC3 = 0.053497 + W_density = 0.8533 + H_density = 0.2400 iter = 13 - elapsed = 9.3700e-04 - tTotal = 0.025710 + elapsed = 2.0350e-03 + tTotal = 0.043718 subIterW = 0 subIterH = 0 - relError = 0.3255 + relError = 0.3199 toshow = scalar structure containing the fields: - initGrNormW = 406.79 - initGrNormH = 524.35 - initNorm = 14.016 - SC1 = 1.5036e-03 - SC2 = 8.8414e-04 - SC3 = 0.1087 - W_density = 0.8300 - H_density = 0.2250 + initGrNormW = 336.80 + initGrNormH = 432.64 + initNorm = 14.236 + SC1 = 6.2648e-04 + SC2 = 3.7965e-04 + SC3 = 0.047664 + W_density = 0.8433 + H_density = 0.2400 iter = 14 - elapsed = 1.7780e-03 - tTotal = 0.027488 + elapsed = 2.2170e-03 + tTotal = 0.045935 subIterW = 3 subIterH = 0 - relError = 0.3267 + relError = 0.3210 toshow = scalar structure containing the fields: - initGrNormW = 406.79 - initGrNormH = 524.35 - initNorm = 14.016 - SC1 = 1.3794e-03 - SC2 = 7.9180e-04 - SC3 = 0.1038 - W_density = 0.8100 - H_density = 0.2200 + initGrNormW = 336.80 + initGrNormH = 432.64 + initNorm = 14.236 + SC1 = 5.6190e-04 + SC2 = 3.3939e-04 + SC3 = 0.011996 + W_density = 0.8433 + H_density = 0.2400 iter = 15 - elapsed = 1.7580e-03 - tTotal = 0.029246 - subIterW = 6 - subIterH = 0 - relError = 0.3279 + elapsed = 2.2430e-03 + tTotal = 0.048178 + subIterW = 0 + subIterH = 2 + relError = 0.3219 toshow = scalar structure containing the fields: - initGrNormW = 406.79 - initGrNormH = 524.35 - initNorm = 14.016 - SC1 = 1.2575e-03 - SC2 = 7.1932e-04 - SC3 = 0.1038 - W_density = 0.8133 - H_density = 0.2100 + initGrNormW = 336.80 + initGrNormH = 432.64 + initNorm = 14.236 + SC1 = 5.3347e-04 + SC2 = 3.2115e-04 + SC3 = 0.041257 + W_density = 0.8400 + H_density = 0.2450 iter = 16 - elapsed = 1.6190e-03 - tTotal = 0.030865 + elapsed = 2.6110e-03 + tTotal = 0.050789 subIterW = 0 - subIterH = 2 - relError = 0.3289 + subIterH = 0 + relError = 0.3228 toshow = scalar structure containing the fields: - initGrNormW = 406.79 - initGrNormH = 524.35 - initNorm = 14.016 - SC1 = 1.1678e-03 - SC2 = 6.6562e-04 - SC3 = 0.094757 - W_density = 0.8067 - H_density = 0.2100 + initGrNormW = 336.80 + initGrNormH = 432.64 + initNorm = 14.236 + SC1 = 5.3415e-04 + SC2 = 3.1942e-04 + SC3 = 0.041178 + W_density = 0.8333 + H_density = 0.2450 iter = 17 - elapsed = 1.4430e-03 - tTotal = 0.032308 - subIterW = 2 + elapsed = 2.1280e-03 + tTotal = 0.052917 + subIterW = 1 subIterH = 0 - relError = 0.3301 + relError = 0.3236 toshow = scalar structure containing the fields: - initGrNormW = 406.79 - initGrNormH = 524.35 - initNorm = 14.016 - SC1 = 1.0588e-03 - SC2 = 6.0138e-04 - SC3 = 0.086351 - W_density = 0.8033 - H_density = 0.2100 + initGrNormW = 336.80 + initGrNormH = 432.64 + initNorm = 14.236 + SC1 = 5.6735e-04 + SC2 = 3.3814e-04 + SC3 = 0.043302 + W_density = 0.8300 + H_density = 0.2450 iter = 18 - elapsed = 1.2770e-03 - tTotal = 0.033585 + elapsed = 2.0490e-03 + tTotal = 0.054966 subIterW = 0 subIterH = 0 - relError = 0.3314 + relError = 0.3245 toshow = scalar structure containing the fields: - initGrNormW = 406.79 - initGrNormH = 524.35 - initNorm = 14.016 - SC1 = 9.6638e-04 - SC2 = 5.4697e-04 - SC3 = 0.079821 - W_density = 0.8000 - H_density = 0.2100 + initGrNormW = 336.80 + initGrNormH = 432.64 + initNorm = 14.236 + SC1 = 6.3445e-04 + SC2 = 3.7940e-04 + SC3 = 0.047360 + W_density = 0.8300 + H_density = 0.2500 iter = 19 - elapsed = 1.3850e-03 - tTotal = 0.034970 - subIterW = 3 + elapsed = 2.8310e-03 + tTotal = 0.057797 + subIterW = 2 subIterH = 0 - relError = 0.3327 + relError = 0.3254 toshow = scalar structure containing the fields: - initGrNormW = 406.79 - initGrNormH = 524.35 - initNorm = 14.016 - SC1 = 8.7237e-04 - SC2 = 4.9202e-04 - SC3 = 0.016631 - W_density = 0.8000 - H_density = 0.2050 + initGrNormW = 336.80 + initGrNormH = 432.64 + initNorm = 14.236 + SC1 = 7.5629e-04 + SC2 = 4.4923e-04 + SC3 = 0.054162 + W_density = 0.8167 + H_density = 0.2500 iter = 20 - elapsed = 1.2630e-03 - tTotal = 0.036233 - subIterW = 0 + elapsed = 2.3150e-03 + tTotal = 0.060112 + subIterW = 4 subIterH = 0 - relError = 0.3339 + relError = 0.3263 toshow = scalar structure containing the fields: - initGrNormW = 406.79 - initGrNormH = 524.35 - initNorm = 14.016 - SC1 = 8.1075e-04 - SC2 = 4.5564e-04 - SC3 = 0.070481 - W_density = 0.7967 - H_density = 0.2100 + initGrNormW = 336.80 + initGrNormH = 432.64 + initNorm = 14.236 + SC1 = 9.8296e-04 + SC2 = 5.7994e-04 + SC3 = 0.062276 + W_density = 0.8100 + H_density = 0.2550 iter = 21 - elapsed = 1.6480e-03 - tTotal = 0.037881 - subIterW = 0 + elapsed = 2.7580e-03 + tTotal = 0.062870 + subIterW = 2 subIterH = 0 - relError = 0.3350 + relError = 0.3276 toshow = scalar structure containing the fields: - initGrNormW = 406.79 - initGrNormH = 524.35 - initNorm = 14.016 - SC1 = 7.6379e-04 - SC2 = 4.2467e-04 - SC3 = 0.065213 - W_density = 0.7867 - H_density = 0.2100 + initGrNormW = 336.80 + initGrNormH = 432.64 + initNorm = 14.236 + SC1 = 8.4207e-04 + SC2 = 4.9177e-04 + SC3 = 0.057272 + W_density = 0.8067 + H_density = 0.2450 iter = 22 - elapsed = 1.4250e-03 - tTotal = 0.039306 + elapsed = 3.0360e-03 + tTotal = 0.065906 subIterW = 3 - subIterH = 0 - relError = 0.3360 + subIterH = 2 + relError = 0.3284 toshow = scalar structure containing the fields: - initGrNormW = 406.79 - initGrNormH = 524.35 - initNorm = 14.016 - SC1 = 6.9574e-04 - SC2 = 3.8544e-04 - SC3 = 0.062853 - W_density = 0.7867 - H_density = 0.2050 + initGrNormW = 336.80 + initGrNormH = 432.64 + initNorm = 14.236 + SC1 = 7.7695e-04 + SC2 = 4.5374e-04 + SC3 = 0.054236 + W_density = 0.8100 + H_density = 0.2450 iter = 23 - elapsed = 1.8270e-03 - tTotal = 0.041133 - subIterW = 2 - subIterH = 0 - relError = 0.3368 + elapsed = 2.9660e-03 + tTotal = 0.068872 + subIterW = 0 + subIterH = 2 + relError = 0.3290 Stop: tolerance reached. final = @@ -4375,19 +4426,19 @@ scalar structure containing the fields: iterations = 23 - elapsed_total = 0.041133 - relative_error = 0.3368 - W_density = 0.7867 - H_density = 0.2050 + elapsed_total = 0.068872 + relative_error = 0.3290 + W_density = 0.8100 + H_density = 0.2450 ***** test [W,H,iter,HIS]=nmf_bpas(A,k,'verbose','type','sparse','nnlssolver','bp','alpha',1.1,'beta',1.3); - initGrNormW = 356.97 - initGrNormH = 500.18 - initNorm = 14.016 - SC1 = 1.2290 - SC2 = 614.50 - SC3 = 0.4849 + initGrNormW = 458.61 + initGrNormH = 637.45 + initNorm = 14.236 + SC1 = 1.5706 + SC2 = 785.28 + SC3 = 0.5287 W_density = 1 H_density = 1 par = @@ -4411,438 +4462,438 @@ scalar structure containing the fields: - initGrNormW = 356.97 - initGrNormH = 500.18 - initNorm = 14.016 - SC1 = 0.024599 - SC2 = 0.024501 - SC3 = 1.0442 + initGrNormW = 458.61 + initGrNormH = 637.45 + initNorm = 14.236 + SC1 = 0.022463 + SC2 = 0.022373 + SC3 = 1.0658 W_density = 1 - H_density = 0.4800 + H_density = 0.5500 iter = 1 - elapsed = 2.7440e-03 - tTotal = 2.7440e-03 + elapsed = 4.2990e-03 + tTotal = 4.2990e-03 subIterW = 0 - subIterH = 37 - relError = 0.5473 + subIterH = 41 + relError = 0.5691 toshow = scalar structure containing the fields: - initGrNormW = 356.97 - initGrNormH = 500.18 - initNorm = 14.016 - SC1 = 8.2229e-03 - SC2 = 5.3449e-03 - SC3 = 0.5920 - W_density = 0.8400 - H_density = 0.2200 + initGrNormW = 458.61 + initGrNormH = 637.45 + initNorm = 14.236 + SC1 = 0.011332 + SC2 = 5.6206e-03 + SC3 = 0.6842 + W_density = 0.6300 + H_density = 0.2250 iter = 2 - elapsed = 2.8770e-03 - tTotal = 5.6210e-03 - subIterW = 5 - subIterH = 35 - relError = 0.4257 + elapsed = 4.5930e-03 + tTotal = 8.8920e-03 + subIterW = 4 + subIterH = 34 + relError = 0.4403 toshow = scalar structure containing the fields: - initGrNormW = 356.97 - initGrNormH = 500.18 - initNorm = 14.016 - SC1 = 6.4035e-03 - SC2 = 3.3682e-03 - SC3 = 0.3947 - W_density = 0.6700 - H_density = 0.2300 + initGrNormW = 458.61 + initGrNormH = 637.45 + initNorm = 14.236 + SC1 = 7.7237e-03 + SC2 = 3.1204e-03 + SC3 = 0.4327 + W_density = 0.5167 + H_density = 0.2000 iter = 3 - elapsed = 2.5790e-03 - tTotal = 8.2000e-03 + elapsed = 3.6820e-03 + tTotal = 0.012574 subIterW = 5 - subIterH = 13 - relError = 0.4196 + subIterH = 9 + relError = 0.4393 toshow = scalar structure containing the fields: - initGrNormW = 356.97 - initGrNormH = 500.18 - initNorm = 14.016 - SC1 = 6.2455e-03 - SC2 = 2.9604e-03 - SC3 = 0.3158 - W_density = 0.6200 - H_density = 0.1900 + initGrNormW = 458.61 + initGrNormH = 637.45 + initNorm = 14.236 + SC1 = 7.9624e-03 + SC2 = 2.4365e-03 + SC3 = 0.3493 + W_density = 0.3800 + H_density = 0.1750 iter = 4 - elapsed = 1.9120e-03 - tTotal = 0.010112 - subIterW = 4 - subIterH = 6 - relError = 0.4216 + elapsed = 4.1280e-03 + tTotal = 0.016702 + subIterW = 5 + subIterH = 8 + relError = 0.4443 toshow = scalar structure containing the fields: - initGrNormW = 356.97 - initGrNormH = 500.18 - initNorm = 14.016 - SC1 = 5.9976e-03 - SC2 = 2.2671e-03 - SC3 = 0.3336 - W_density = 0.4833 - H_density = 0.2000 + initGrNormW = 458.61 + initGrNormH = 637.45 + initNorm = 14.236 + SC1 = 6.5834e-03 + SC2 = 1.9092e-03 + SC3 = 0.3137 + W_density = 0.3600 + H_density = 0.1700 iter = 5 - elapsed = 2.0030e-03 - tTotal = 0.012115 - subIterW = 3 - subIterH = 7 - relError = 0.4244 + elapsed = 2.2240e-03 + tTotal = 0.018926 + subIterW = 2 + subIterH = 3 + relError = 0.4491 toshow = scalar structure containing the fields: - initGrNormW = 356.97 - initGrNormH = 500.18 - initNorm = 14.016 - SC1 = 9.2842e-03 - SC2 = 3.4723e-03 - SC3 = 0.5179 - W_density = 0.4800 - H_density = 0.1900 + initGrNormW = 458.61 + initGrNormH = 637.45 + initNorm = 14.236 + SC1 = 5.4242e-03 + SC2 = 1.2910e-03 + SC3 = 0.2382 + W_density = 0.2833 + H_density = 0.1650 iter = 6 - elapsed = 1.9260e-03 - tTotal = 0.014041 - subIterW = 0 - subIterH = 2 - relError = 0.4299 + elapsed = 2.7360e-03 + tTotal = 0.021662 + subIterW = 2 + subIterH = 0 + relError = 0.4526 toshow = scalar structure containing the fields: - initGrNormW = 356.97 - initGrNormH = 500.18 - initNorm = 14.016 - SC1 = 0.012646 - SC2 = 4.3250e-03 - SC3 = 0.6754 - W_density = 0.4400 - H_density = 0.1750 + initGrNormW = 458.61 + initGrNormH = 637.45 + initNorm = 14.236 + SC1 = 5.5547e-03 + SC2 = 1.3220e-03 + SC3 = 0.2428 + W_density = 0.2800 + H_density = 0.1650 iter = 7 - elapsed = 1.5340e-03 - tTotal = 0.015575 - subIterW = 5 - subIterH = 3 - relError = 0.4417 + elapsed = 1.4640e-03 + tTotal = 0.023126 + subIterW = 0 + subIterH = 0 + relError = 0.4565 toshow = scalar structure containing the fields: - initGrNormW = 356.97 - initGrNormH = 500.18 - initNorm = 14.016 - SC1 = 0.013090 - SC2 = 3.5866e-03 - SC3 = 0.4368 - W_density = 0.3400 - H_density = 0.1600 + initGrNormW = 458.61 + initGrNormH = 637.45 + initNorm = 14.236 + SC1 = 7.0669e-03 + SC2 = 1.6112e-03 + SC3 = 0.2959 + W_density = 0.2667 + H_density = 0.1650 iter = 8 - elapsed = 1.8380e-03 - tTotal = 0.017413 - subIterW = 4 - subIterH = 3 - relError = 0.4526 + elapsed = 2.1370e-03 + tTotal = 0.025263 + subIterW = 1 + subIterH = 2 + relError = 0.4623 toshow = scalar structure containing the fields: - initGrNormW = 356.97 - initGrNormH = 500.18 - initNorm = 14.016 - SC1 = 8.6612e-03 - SC2 = 1.8882e-03 - SC3 = 0.2148 - W_density = 0.2600 - H_density = 0.1450 + initGrNormW = 458.61 + initGrNormH = 637.45 + initNorm = 14.236 + SC1 = 7.9710e-03 + SC2 = 1.6899e-03 + SC3 = 0.3520 + W_density = 0.2500 + H_density = 0.1550 iter = 9 - elapsed = 1.6690e-03 - tTotal = 0.019082 + elapsed = 2.0470e-03 + tTotal = 0.027310 subIterW = 1 - subIterH = 2 - relError = 0.4565 + subIterH = 1 + relError = 0.4679 toshow = scalar structure containing the fields: - initGrNormW = 356.97 - initGrNormH = 500.18 - initNorm = 14.016 - SC1 = 6.1400e-03 - SC2 = 1.0929e-03 - SC3 = 0.1253 + initGrNormW = 458.61 + initGrNormH = 637.45 + initNorm = 14.236 + SC1 = 0.011095 + SC2 = 1.9527e-03 + SC3 = 0.2155 W_density = 0.2000 H_density = 0.1400 iter = 10 - elapsed = 8.7400e-04 - tTotal = 0.019956 + elapsed = 1.5060e-03 + tTotal = 0.028816 subIterW = 0 - subIterH = 2 - relError = 0.4570 + subIterH = 1 + relError = 0.4706 toshow = scalar structure containing the fields: - initGrNormW = 356.97 - initGrNormH = 500.18 - initNorm = 14.016 - SC1 = 8.1539e-03 - SC2 = 1.4514e-03 - SC3 = 0.1296 + initGrNormW = 458.61 + initGrNormH = 637.45 + initNorm = 14.236 + SC1 = 0.014171 + SC2 = 2.4657e-03 + SC3 = 0.2898 W_density = 0.2000 - H_density = 0.1450 + H_density = 0.1350 iter = 11 - elapsed = 7.8900e-04 - tTotal = 0.020745 + elapsed = 1.3270e-03 + tTotal = 0.030143 subIterW = 0 subIterH = 0 - relError = 0.4576 + relError = 0.4733 toshow = scalar structure containing the fields: - initGrNormW = 356.97 - initGrNormH = 500.18 - initNorm = 14.016 - SC1 = 0.010765 - SC2 = 1.9161e-03 - SC3 = 0.1759 - W_density = 0.2000 - H_density = 0.1450 + initGrNormW = 458.61 + initGrNormH = 637.45 + initNorm = 14.236 + SC1 = 0.015218 + SC2 = 2.4044e-03 + SC3 = 0.5114 + W_density = 0.1800 + H_density = 0.1250 iter = 12 - elapsed = 4.4800e-04 - tTotal = 0.021193 - subIterW = 0 - subIterH = 0 - relError = 0.4587 + elapsed = 1.9450e-03 + tTotal = 0.032088 + subIterW = 1 + subIterH = 1 + relError = 0.4769 toshow = scalar structure containing the fields: - initGrNormW = 356.97 - initGrNormH = 500.18 - initNorm = 14.016 - SC1 = 0.014020 - SC2 = 2.4956e-03 - SC3 = 0.2004 - W_density = 0.2000 - H_density = 0.1450 + initGrNormW = 458.61 + initGrNormH = 637.45 + initNorm = 14.236 + SC1 = 0.011698 + SC2 = 1.6846e-03 + SC3 = 0.2385 + W_density = 0.1633 + H_density = 0.1150 iter = 13 - elapsed = 4.3100e-04 - tTotal = 0.021624 - subIterW = 0 - subIterH = 0 - relError = 0.4607 + elapsed = 2.0250e-03 + tTotal = 0.034113 + subIterW = 1 + subIterH = 1 + relError = 0.4818 toshow = scalar structure containing the fields: - initGrNormW = 356.97 - initGrNormH = 500.18 - initNorm = 14.016 - SC1 = 0.016186 - SC2 = 2.6545e-03 - SC3 = 0.5097 - W_density = 0.1867 - H_density = 0.1300 + initGrNormW = 458.61 + initGrNormH = 637.45 + initNorm = 14.236 + SC1 = 8.5883e-03 + SC2 = 1.1852e-03 + SC3 = 0.085820 + W_density = 0.1600 + H_density = 0.1050 iter = 14 - elapsed = 1.2900e-03 - tTotal = 0.022914 - subIterW = 1 + elapsed = 1.8850e-03 + tTotal = 0.035998 + subIterW = 0 subIterH = 1 - relError = 0.4642 + relError = 0.4860 toshow = scalar structure containing the fields: - initGrNormW = 356.97 - initGrNormH = 500.18 - initNorm = 14.016 - SC1 = 0.014814 - SC2 = 2.3702e-03 - SC3 = 0.4594 - W_density = 0.1833 - H_density = 0.1250 + initGrNormW = 458.61 + initGrNormH = 637.45 + initNorm = 14.236 + SC1 = 6.3862e-04 + SC2 = 6.3862e-05 + SC3 = 0.011699 + W_density = 0.1000 + H_density = 0.1000 iter = 15 - elapsed = 1.0920e-03 - tTotal = 0.024006 + elapsed = 1.2650e-03 + tTotal = 0.037263 subIterW = 0 subIterH = 0 - relError = 0.4687 + relError = 0.4863 toshow = scalar structure containing the fields: - initGrNormW = 356.97 - initGrNormH = 500.18 - initNorm = 14.016 - SC1 = 0.014377 - SC2 = 2.1853e-03 - SC3 = 0.3191 - W_density = 0.1733 - H_density = 0.1200 + initGrNormW = 458.61 + initGrNormH = 637.45 + initNorm = 14.236 + SC1 = 3.6671e-04 + SC2 = 3.6671e-05 + SC3 = 8.0798e-03 + W_density = 0.1000 + H_density = 0.1000 iter = 16 - elapsed = 1.1630e-03 - tTotal = 0.025169 - subIterW = 1 + elapsed = 6.7900e-04 + tTotal = 0.037942 + subIterW = 0 subIterH = 0 - relError = 0.4750 + relError = 0.4863 toshow = scalar structure containing the fields: - initGrNormW = 356.97 - initGrNormH = 500.18 - initNorm = 14.016 - SC1 = 0.010027 - SC2 = 1.3637e-03 - SC3 = 0.085309 - W_density = 0.1567 - H_density = 0.1050 + initGrNormW = 458.61 + initGrNormH = 637.45 + initNorm = 14.236 + SC1 = 2.3955e-04 + SC2 = 2.3955e-05 + SC3 = 5.8746e-03 + W_density = 0.1000 + H_density = 0.1000 iter = 17 - elapsed = 1.2730e-03 - tTotal = 0.026442 - subIterW = 1 - subIterH = 1 - relError = 0.4791 + elapsed = 6.6200e-04 + tTotal = 0.038604 + subIterW = 0 + subIterH = 0 + relError = 0.4862 toshow = scalar structure containing the fields: - initGrNormW = 356.97 - initGrNormH = 500.18 - initNorm = 14.016 - SC1 = 9.1906e-04 - SC2 = 9.1906e-05 - SC3 = 0.018660 + initGrNormW = 458.61 + initGrNormH = 637.45 + initNorm = 14.236 + SC1 = 1.5658e-04 + SC2 = 1.5658e-05 + SC3 = 3.7030e-03 W_density = 0.1000 H_density = 0.1000 iter = 18 - elapsed = 8.5800e-04 - tTotal = 0.027300 + elapsed = 6.4400e-04 + tTotal = 0.039248 subIterW = 0 subIterH = 0 - relError = 0.4794 + relError = 0.4862 toshow = scalar structure containing the fields: - initGrNormW = 356.97 - initGrNormH = 500.18 - initNorm = 14.016 - SC1 = 5.5772e-04 - SC2 = 5.5772e-05 - SC3 = 0.011627 + initGrNormW = 458.61 + initGrNormH = 637.45 + initNorm = 14.236 + SC1 = 1.0236e-04 + SC2 = 1.0236e-05 + SC3 = 2.5104e-03 W_density = 0.1000 H_density = 0.1000 iter = 19 - elapsed = 4.7500e-04 - tTotal = 0.027775 + elapsed = 6.7800e-04 + tTotal = 0.039926 subIterW = 0 subIterH = 0 - relError = 0.4794 + relError = 0.4862 toshow = scalar structure containing the fields: - initGrNormW = 356.97 - initGrNormH = 500.18 - initNorm = 14.016 - SC1 = 3.6243e-04 - SC2 = 3.6243e-05 - SC3 = 7.8458e-03 + initGrNormW = 458.61 + initGrNormH = 637.45 + initNorm = 14.236 + SC1 = 6.6918e-05 + SC2 = 6.6918e-06 + SC3 = 1.6412e-03 W_density = 0.1000 H_density = 0.1000 iter = 20 - elapsed = 4.2900e-04 - tTotal = 0.028204 + elapsed = 6.3600e-04 + tTotal = 0.040562 subIterW = 0 subIterH = 0 - relError = 0.4793 + relError = 0.4862 toshow = scalar structure containing the fields: - initGrNormW = 356.97 - initGrNormH = 500.18 - initNorm = 14.016 - SC1 = 2.3563e-04 - SC2 = 2.3563e-05 - SC3 = 5.3048e-03 + initGrNormW = 458.61 + initGrNormH = 637.45 + initNorm = 14.236 + SC1 = 4.3751e-05 + SC2 = 4.3751e-06 + SC3 = 1.0347e-03 W_density = 0.1000 H_density = 0.1000 iter = 21 - elapsed = 4.0600e-04 - tTotal = 0.028610 + elapsed = 6.4600e-04 + tTotal = 0.041208 subIterW = 0 subIterH = 0 - relError = 0.4793 + relError = 0.4862 toshow = scalar structure containing the fields: - initGrNormW = 356.97 - initGrNormH = 500.18 - initNorm = 14.016 - SC1 = 1.5322e-04 - SC2 = 1.5322e-05 - SC3 = 3.3167e-03 + initGrNormW = 458.61 + initGrNormH = 637.45 + initNorm = 14.236 + SC1 = 2.8606e-05 + SC2 = 2.8606e-06 + SC3 = 6.5318e-04 W_density = 0.1000 H_density = 0.1000 iter = 22 - elapsed = 4.9800e-04 - tTotal = 0.029108 + elapsed = 7.7900e-04 + tTotal = 0.041987 subIterW = 0 subIterH = 0 - relError = 0.4793 + relError = 0.4862 toshow = scalar structure containing the fields: - initGrNormW = 356.97 - initGrNormH = 500.18 - initNorm = 14.016 - SC1 = 9.9638e-05 - SC2 = 9.9638e-06 - SC3 = 2.3366e-03 + initGrNormW = 458.61 + initGrNormH = 637.45 + initNorm = 14.236 + SC1 = 1.8703e-05 + SC2 = 1.8703e-06 + SC3 = 4.9541e-04 W_density = 0.1000 H_density = 0.1000 iter = 23 - elapsed = 7.7900e-04 - tTotal = 0.029887 + elapsed = 7.3500e-04 + tTotal = 0.042722 subIterW = 0 subIterH = 0 - relError = 0.4793 + relError = 0.4862 Stop: tolerance reached. final = @@ -4850,19 +4901,19 @@ scalar structure containing the fields: iterations = 23 - elapsed_total = 0.029887 - relative_error = 0.4793 + elapsed_total = 0.042722 + relative_error = 0.4862 W_density = 0.1000 H_density = 0.1000 ***** test [W,H,iter,HIS]=nmf_bpas(A,k,'verbose','type','plain','winit',rand(m,k)); - initGrNormW = 373.38 - initGrNormH = 438.79 - initNorm = 14.016 - SC1 = 1.1523 - SC2 = 576.15 - SC3 = 0.4968 + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 1.0893 + SC2 = 544.66 + SC3 = 0.4878 W_density = 1 H_density = 1 par = @@ -4873,8 +4924,8 @@ n = 20 type = plain nnls_solver = bp - alpha = 0.4979 - beta = 0.4979 + alpha = 0.5052 + beta = 0.5052 max_iter = 100 min_iter = 20 max_time = 1000 @@ -4886,468 +4937,525 @@ scalar structure containing the fields: - initGrNormW = 373.38 - initGrNormH = 438.79 - initNorm = 14.016 - SC1 = 0.013177 - SC2 = 9.5140e-03 - SC3 = 0.4830 + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 0.017865 + SC2 = 0.013506 + SC3 = 0.5816 W_density = 0.7500 - H_density = 0.6300 + H_density = 0.6700 iter = 1 - elapsed = 5.3480e-03 - tTotal = 5.3480e-03 - subIterW = 44 - subIterH = 33 - relError = 0.3263 + elapsed = 6.2510e-03 + tTotal = 6.2510e-03 + subIterW = 41 + subIterH = 26 + relError = 0.3512 toshow = scalar structure containing the fields: - initGrNormW = 373.38 - initGrNormH = 438.79 - initNorm = 14.016 - SC1 = 8.2105e-03 - SC2 = 5.6653e-03 - SC3 = 0.3027 - W_density = 0.7300 - H_density = 0.6050 + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 0.012353 + SC2 = 8.9929e-03 + SC3 = 0.4066 + W_density = 0.7467 + H_density = 0.6700 iter = 2 - elapsed = 2.9050e-03 - tTotal = 8.2530e-03 - subIterW = 19 - subIterH = 18 - relError = 0.2802 + elapsed = 5.7250e-03 + tTotal = 0.011976 + subIterW = 28 + subIterH = 20 + relError = 0.2879 toshow = scalar structure containing the fields: - initGrNormW = 373.38 - initGrNormH = 438.79 - initNorm = 14.016 - SC1 = 6.0097e-03 - SC2 = 4.1707e-03 - SC3 = 0.2277 - W_density = 0.7367 - H_density = 0.5900 + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 6.7471e-03 + SC2 = 4.6285e-03 + SC3 = 0.2504 + W_density = 0.7267 + H_density = 0.6050 iter = 3 - elapsed = 2.7070e-03 - tTotal = 0.010960 - subIterW = 14 - subIterH = 17 - relError = 0.2590 + elapsed = 4.6230e-03 + tTotal = 0.016599 + subIterW = 17 + subIterH = 18 + relError = 0.2622 toshow = scalar structure containing the fields: - initGrNormW = 373.38 - initGrNormH = 438.79 - initNorm = 14.016 - SC1 = 4.4657e-03 - SC2 = 3.0724e-03 - SC3 = 0.1963 - W_density = 0.7300 - H_density = 0.5900 + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 4.3900e-03 + SC2 = 2.9326e-03 + SC3 = 0.1736 + W_density = 0.7233 + H_density = 0.5750 iter = 4 - elapsed = 2.5780e-03 - tTotal = 0.013538 - subIterW = 4 - subIterH = 15 - relError = 0.2480 + elapsed = 3.4690e-03 + tTotal = 0.020068 + subIterW = 9 + subIterH = 8 + relError = 0.2539 toshow = scalar structure containing the fields: - initGrNormW = 373.38 - initGrNormH = 438.79 - initNorm = 14.016 - SC1 = 3.3733e-03 - SC2 = 2.2803e-03 - SC3 = 0.1663 + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 3.1995e-03 + SC2 = 2.1053e-03 + SC3 = 0.1398 W_density = 0.7267 - H_density = 0.5850 + H_density = 0.5500 iter = 5 - elapsed = 2.7300e-03 - tTotal = 0.016268 - subIterW = 8 + elapsed = 3.4630e-03 + tTotal = 0.023531 + subIterW = 9 subIterH = 9 - relError = 0.2421 + relError = 0.2503 toshow = scalar structure containing the fields: - initGrNormW = 373.38 - initGrNormH = 438.79 - initNorm = 14.016 - SC1 = 2.5945e-03 - SC2 = 1.7643e-03 - SC3 = 0.1281 - W_density = 0.7200 - H_density = 0.5900 + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 2.6091e-03 + SC2 = 1.7168e-03 + SC3 = 0.1133 + W_density = 0.7233 + H_density = 0.5500 iter = 6 - elapsed = 2.0120e-03 - tTotal = 0.018280 - subIterW = 3 - subIterH = 6 - relError = 0.2389 + elapsed = 3.2770e-03 + tTotal = 0.026808 + subIterW = 5 + subIterH = 5 + relError = 0.2482 toshow = scalar structure containing the fields: - initGrNormW = 373.38 - initGrNormH = 438.79 - initNorm = 14.016 - SC1 = 2.2212e-03 - SC2 = 1.5104e-03 - SC3 = 0.1109 - W_density = 0.7200 - H_density = 0.6100 + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 2.1206e-03 + SC2 = 1.3784e-03 + SC3 = 0.1015 + W_density = 0.7233 + H_density = 0.5250 iter = 7 - elapsed = 2.3340e-03 - tTotal = 0.020614 - subIterW = 5 + elapsed = 3.2110e-03 + tTotal = 0.030019 + subIterW = 3 subIterH = 5 - relError = 0.2368 + relError = 0.2468 toshow = scalar structure containing the fields: - initGrNormW = 373.38 - initGrNormH = 438.79 - initNorm = 14.016 - SC1 = 2.0297e-03 - SC2 = 1.3924e-03 - SC3 = 0.098257 - W_density = 0.7233 - H_density = 0.6300 + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 1.9463e-03 + SC2 = 1.2690e-03 + SC3 = 0.091730 + W_density = 0.7267 + H_density = 0.5400 iter = 8 - elapsed = 1.9920e-03 - tTotal = 0.022606 - subIterW = 3 - subIterH = 6 - relError = 0.2352 + elapsed = 3.2630e-03 + tTotal = 0.033282 + subIterW = 5 + subIterH = 5 + relError = 0.2459 toshow = scalar structure containing the fields: - initGrNormW = 373.38 - initGrNormH = 438.79 - initNorm = 14.016 - SC1 = 1.9402e-03 - SC2 = 1.3116e-03 - SC3 = 0.092583 - W_density = 0.7200 - H_density = 0.6050 + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 1.8145e-03 + SC2 = 1.1831e-03 + SC3 = 0.083610 + W_density = 0.7267 + H_density = 0.5350 iter = 9 - elapsed = 1.9680e-03 - tTotal = 0.024574 - subIterW = 3 - subIterH = 5 - relError = 0.2337 + elapsed = 3.0310e-03 + tTotal = 0.036313 + subIterW = 2 + subIterH = 0 + relError = 0.2450 toshow = scalar structure containing the fields: - initGrNormW = 373.38 - initGrNormH = 438.79 - initNorm = 14.016 - SC1 = 1.8625e-03 - SC2 = 1.2404e-03 - SC3 = 0.088680 - W_density = 0.7167 - H_density = 0.5900 + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 1.7393e-03 + SC2 = 1.1305e-03 + SC3 = 0.076260 + W_density = 0.7233 + H_density = 0.5400 iter = 10 - elapsed = 2.0760e-03 - tTotal = 0.026650 - subIterW = 3 - subIterH = 2 - relError = 0.2325 + elapsed = 3.0260e-03 + tTotal = 0.039339 + subIterW = 0 + subIterH = 3 + relError = 0.2443 toshow = scalar structure containing the fields: - initGrNormW = 373.38 - initGrNormH = 438.79 - initNorm = 14.016 - SC1 = 1.7362e-03 - SC2 = 1.1737e-03 - SC3 = 0.081030 - W_density = 0.7233 - H_density = 0.5900 + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 1.6978e-03 + SC2 = 1.1002e-03 + SC3 = 0.069970 + W_density = 0.7133 + H_density = 0.5400 iter = 11 - elapsed = 1.9680e-03 - tTotal = 0.028618 - subIterW = 4 + elapsed = 3.1560e-03 + tTotal = 0.042495 + subIterW = 3 subIterH = 2 - relError = 0.2314 + relError = 0.2436 toshow = scalar structure containing the fields: - initGrNormW = 373.38 - initGrNormH = 438.79 - initNorm = 14.016 - SC1 = 1.6751e-03 - SC2 = 1.1357e-03 - SC3 = 0.079312 - W_density = 0.7267 - H_density = 0.6000 + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 1.7272e-03 + SC2 = 1.1192e-03 + SC3 = 0.068398 + W_density = 0.7167 + H_density = 0.5400 iter = 12 - elapsed = 2.0060e-03 - tTotal = 0.030624 + elapsed = 3.1660e-03 + tTotal = 0.045661 subIterW = 3 - subIterH = 4 - relError = 0.2305 + subIterH = 2 + relError = 0.2429 toshow = scalar structure containing the fields: - initGrNormW = 373.38 - initGrNormH = 438.79 - initNorm = 14.016 - SC1 = 1.5731e-03 - SC2 = 1.0602e-03 - SC3 = 0.078042 - W_density = 0.7300 - H_density = 0.5800 + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 1.7830e-03 + SC2 = 1.1661e-03 + SC3 = 0.066376 + W_density = 0.7233 + H_density = 0.5400 iter = 13 - elapsed = 1.9910e-03 - tTotal = 0.032615 - subIterW = 3 - subIterH = 4 - relError = 0.2295 + elapsed = 2.4780e-03 + tTotal = 0.048139 + subIterW = 4 + subIterH = 0 + relError = 0.2422 toshow = scalar structure containing the fields: - initGrNormW = 373.38 - initGrNormH = 438.79 - initNorm = 14.016 - SC1 = 1.4649e-03 - SC2 = 9.9612e-04 - SC3 = 0.072467 - W_density = 0.7333 - H_density = 0.5850 + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 1.7905e-03 + SC2 = 1.1710e-03 + SC3 = 0.065656 + W_density = 0.7233 + H_density = 0.5450 iter = 14 - elapsed = 2.0170e-03 - tTotal = 0.034632 - subIterW = 3 - subIterH = 3 - relError = 0.2287 + elapsed = 3.1900e-03 + tTotal = 0.051329 + subIterW = 2 + subIterH = 4 + relError = 0.2415 toshow = scalar structure containing the fields: - initGrNormW = 373.38 - initGrNormH = 438.79 - initNorm = 14.016 - SC1 = 1.3749e-03 - SC2 = 9.3219e-04 - SC3 = 0.067508 - W_density = 0.7333 - H_density = 0.5950 + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 1.5407e-03 + SC2 = 1.0045e-03 + SC3 = 0.062870 + W_density = 0.7267 + H_density = 0.5300 iter = 15 - elapsed = 2.0800e-03 - tTotal = 0.036712 - subIterW = 6 - subIterH = 4 - relError = 0.2280 + elapsed = 3.7720e-03 + tTotal = 0.055101 + subIterW = 3 + subIterH = 5 + relError = 0.2409 toshow = scalar structure containing the fields: - initGrNormW = 373.38 - initGrNormH = 438.79 - initNorm = 14.016 - SC1 = 1.3179e-03 - SC2 = 8.9615e-04 - SC3 = 0.063452 - W_density = 0.7367 - H_density = 0.5900 + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 1.4732e-03 + SC2 = 9.6054e-04 + SC3 = 0.060159 + W_density = 0.7300 + H_density = 0.5350 iter = 16 - elapsed = 1.7830e-03 - tTotal = 0.038495 - subIterW = 0 - subIterH = 0 - relError = 0.2273 + elapsed = 3.7650e-03 + tTotal = 0.058866 + subIterW = 5 + subIterH = 4 + relError = 0.2403 toshow = scalar structure containing the fields: - initGrNormW = 373.38 - initGrNormH = 438.79 - initNorm = 14.016 - SC1 = 1.2896e-03 - SC2 = 8.7178e-04 - SC3 = 0.060902 - W_density = 0.7300 - H_density = 0.5850 + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 1.3859e-03 + SC2 = 9.1195e-04 + SC3 = 0.056860 + W_density = 0.7333 + H_density = 0.5300 iter = 17 - elapsed = 1.9640e-03 - tTotal = 0.040459 - subIterW = 2 - subIterH = 3 - relError = 0.2267 + elapsed = 3.4060e-03 + tTotal = 0.062272 + subIterW = 3 + subIterH = 0 + relError = 0.2399 toshow = scalar structure containing the fields: - initGrNormW = 373.38 - initGrNormH = 438.79 - initNorm = 14.016 - SC1 = 1.2577e-03 - SC2 = 8.3764e-04 - SC3 = 0.059475 - W_density = 0.7233 - H_density = 0.5800 + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 1.3609e-03 + SC2 = 8.9817e-04 + SC3 = 0.055587 + W_density = 0.7367 + H_density = 0.5300 iter = 18 - elapsed = 1.8910e-03 - tTotal = 0.042350 - subIterW = 2 + elapsed = 3.2940e-03 + tTotal = 0.065566 + subIterW = 3 subIterH = 2 - relError = 0.2262 + relError = 0.2395 toshow = scalar structure containing the fields: - initGrNormW = 373.38 - initGrNormH = 438.79 - initNorm = 14.016 - SC1 = 1.2246e-03 - SC2 = 8.0336e-04 - SC3 = 0.057162 - W_density = 0.7100 - H_density = 0.5750 + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 1.3985e-03 + SC2 = 9.2020e-04 + SC3 = 0.056014 + W_density = 0.7333 + H_density = 0.5350 iter = 19 - elapsed = 1.9270e-03 - tTotal = 0.044277 - subIterW = 4 + elapsed = 3.2010e-03 + tTotal = 0.068767 + subIterW = 0 subIterH = 0 - relError = 0.2257 + relError = 0.2391 toshow = scalar structure containing the fields: - initGrNormW = 373.38 - initGrNormH = 438.79 - initNorm = 14.016 - SC1 = 1.1513e-03 - SC2 = 7.5295e-04 - SC3 = 0.052654 - W_density = 0.7067 - H_density = 0.5750 + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 1.3293e-03 + SC2 = 8.6933e-04 + SC3 = 0.056141 + W_density = 0.7333 + H_density = 0.5300 iter = 20 - elapsed = 1.5680e-03 - tTotal = 0.045845 - subIterW = 0 - subIterH = 0 - relError = 0.2252 + elapsed = 3.1960e-03 + tTotal = 0.071963 + subIterW = 2 + subIterH = 3 + relError = 0.2386 toshow = scalar structure containing the fields: - initGrNormW = 373.38 - initGrNormH = 438.79 - initNorm = 14.016 - SC1 = 1.0018e-03 - SC2 = 6.5716e-04 - SC3 = 0.048999 - W_density = 0.7133 - H_density = 0.5700 + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 1.2891e-03 + SC2 = 8.3793e-04 + SC3 = 0.055052 + W_density = 0.7367 + H_density = 0.5150 iter = 21 - elapsed = 2.0660e-03 - tTotal = 0.047911 - subIterW = 1 + elapsed = 3.4720e-03 + tTotal = 0.075435 + subIterW = 0 subIterH = 3 - relError = 0.2248 + relError = 0.2382 toshow = scalar structure containing the fields: - initGrNormW = 373.38 - initGrNormH = 438.79 - initNorm = 14.016 - SC1 = 9.4876e-04 - SC2 = 6.2239e-04 - SC3 = 0.045301 - W_density = 0.7167 - H_density = 0.5650 + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 1.2582e-03 + SC2 = 8.1028e-04 + SC3 = 0.053344 + W_density = 0.7367 + H_density = 0.5050 iter = 22 - elapsed = 1.9120e-03 - tTotal = 0.049823 - subIterW = 0 - subIterH = 0 - relError = 0.2245 + elapsed = 3.3840e-03 + tTotal = 0.078819 + subIterW = 2 + subIterH = 2 + relError = 0.2379 toshow = scalar structure containing the fields: - initGrNormW = 373.38 - initGrNormH = 438.79 - initNorm = 14.016 - SC1 = 7.9536e-04 - SC2 = 5.2016e-04 - SC3 = 0.040420 - W_density = 0.7167 - H_density = 0.5550 + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 1.1344e-03 + SC2 = 7.2826e-04 + SC3 = 0.049871 + W_density = 0.7333 + H_density = 0.5000 iter = 23 - elapsed = 2.1020e-03 - tTotal = 0.051925 - subIterW = 2 + elapsed = 3.1900e-03 + tTotal = 0.082009 + subIterW = 0 subIterH = 3 - relError = 0.2243 + relError = 0.2375 toshow = scalar structure containing the fields: - initGrNormW = 373.38 - initGrNormH = 438.79 - initNorm = 14.016 - SC1 = 7.2052e-04 - SC2 = 4.6978e-04 - SC3 = 0.036673 - W_density = 0.7133 - H_density = 0.5550 + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 1.0418e-03 + SC2 = 6.7092e-04 + SC3 = 0.044909 + W_density = 0.7367 + H_density = 0.5050 iter = 24 - elapsed = 1.6390e-03 - tTotal = 0.053564 + elapsed = 3.2110e-03 + tTotal = 0.085220 subIterW = 0 subIterH = 0 - relError = 0.2241 + relError = 0.2372 + +toshow = + + scalar structure containing the fields: + + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 9.6410e-04 + SC2 = 6.1895e-04 + SC3 = 0.041523 + W_density = 0.7367 + H_density = 0.5000 + iter = 25 + elapsed = 3.1650e-03 + tTotal = 0.088385 + subIterW = 2 + subIterH = 0 + relError = 0.2370 + +toshow = + + scalar structure containing the fields: + + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 7.2702e-04 + SC2 = 4.6529e-04 + SC3 = 0.016509 + W_density = 0.7367 + H_density = 0.4950 + iter = 26 + elapsed = 2.3540e-03 + tTotal = 0.090739 + subIterW = 0 + subIterH = 0 + relError = 0.2368 + +toshow = + + scalar structure containing the fields: + + initGrNormW = 346.99 + initGrNormH = 419.82 + initNorm = 14.236 + SC1 = 5.6045e-04 + SC2 = 3.5756e-04 + SC3 = 0.013118 + W_density = 0.7367 + H_density = 0.4900 + iter = 27 + elapsed = 2.7620e-03 + tTotal = 0.093501 + subIterW = 0 + subIterH = 0 + relError = 0.2366 Stop: tolerance reached. final = scalar structure containing the fields: - iterations = 24 - elapsed_total = 0.053564 - relative_error = 0.2241 - W_density = 0.7133 - H_density = 0.5550 + iterations = 27 + elapsed_total = 0.093501 + relative_error = 0.2366 + W_density = 0.7367 + H_density = 0.4900 ***** demo m = 300; @@ -5395,38 +5503,23 @@ recovery_error_W = norm(W_org-W)/norm(W_org) recovery_error_H = norm(H_org-H)/norm(H_org) 6 tests, 6 passed, 0 known failure, 0 skipped -[inst/circulant_make_matrix.m] ->>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/circulant_make_matrix.m -***** shared v,C - v = [1 2 3]'; C = [1 3 2; 2 1 3; 3 2 1]; -***** assert (circulant_make_matrix(v), C); -1 test, 1 passed, 0 known failure, 0 skipped -[inst/funm.m] ->>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/funm.m -***** function b = fsin (a) - b = sin (a); -***** endfunction +[inst/@blksparse/ctranspose.m] +>>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/@blksparse/ctranspose.m ***** test - mtx = randn (100); - assert (funm (mtx, "fsin"), thfm (mtx, "sin"), 1e-9) + r = blksparse ([1,2],[1,2],cat(3,eye(2),[1 2; -2 1])); + rt = r'; + assert (full(rt'),full(r)); 1 test, 1 passed, 0 known failure, 0 skipped -[inst/vec_projection.m] ->>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/vec_projection.m +[inst/smwsolve.m] +>>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/smwsolve.m ***** test - assert (vec_projection ([1,0,0], [0.5,0.5,0]), 0.70711,1e-5); - assert (vec_projection ([1,2000,0], [0.5,15,0]), 1998.9, 1e-1); - assert (vec_projection ([1,-2000,0], [0.5,15,0]), -1998.9, 1e-1); - assert (vec_projection ([7,7,0], [15,0,0]), 7.000, 1e-10); - assert (vec_projection ([1,1,0], [1.05,0.94,0]), 1.4121, 1e-4); - assert (vec_projection ([1,1.1,0], [1.05,0.94,0]), 1.4788, 1e-4); -1 test, 1 passed, 0 known failure, 0 skipped -[inst/circulant_eig.m] ->>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/circulant_eig.m -***** shared v,C,vs,lambda - v = [1 2 3]'; - C = gallery("circul", v)'; - [vs lambda] = circulant_eig(v); -***** assert (vs*lambda, C*vs, 100*eps); + A = 2.1*eye (10); + u = rand (10, 2); u /= diag (norm (u, "cols")); + v = rand (10, 2); v /= diag (norm (v, "cols")); + b = rand (10, 2); + x1 = (A + u*v') \ b; + x2 = smwsolve (A, u, v, b); + assert (x1, x2, 1e-13); 1 test, 1 passed, 0 known failure, 0 skipped Checking C++ files ... Done running the unit tests. @@ -5458,12 +5551,14 @@ dpkg-buildpackage: info: binary-only upload (no source included) dpkg-genchanges: info: not including original source code in upload I: copying local configuration +I: user script /srv/workspace/pbuilder/842317/tmp/hooks/B01_cleanup starting +I: user script /srv/workspace/pbuilder/842317/tmp/hooks/B01_cleanup finished I: unmounting dev/ptmx filesystem I: unmounting dev/pts filesystem I: unmounting dev/shm filesystem I: unmounting proc filesystem I: unmounting sys filesystem I: cleaning the build env -I: removing directory /srv/workspace/pbuilder/994964 and its subdirectories -I: Current time: Wed Oct 22 17:38:47 -12 2025 -I: pbuilder-time-stamp: 1761197927 +I: removing directory /srv/workspace/pbuilder/842317 and its subdirectories +I: Current time: Thu Nov 26 02:05:04 +14 2026 +I: pbuilder-time-stamp: 1795608304