Diff of the two buildlogs: -- --- b1/build.log 2025-02-26 14:37:22.429495248 +0000 +++ b2/build.log 2025-02-26 14:47:05.139212458 +0000 @@ -1,6 +1,6 @@ I: pbuilder: network access will be disabled during build -I: Current time: Tue Mar 31 08:55:40 -12 2026 -I: pbuilder-time-stamp: 1774990540 +I: Current time: Thu Feb 27 04:37:25 +14 2025 +I: pbuilder-time-stamp: 1740580645 I: Building the build Environment I: extracting base tarball [/var/cache/pbuilder/trixie-reproducible-base.tgz] I: copying local configuration @@ -22,52 +22,84 @@ dpkg-source: info: unpacking octave-linear-algebra_2.2.3-4.debian.tar.xz I: Not using root during the build. I: Installing the build-deps -I: user script /srv/workspace/pbuilder/956707/tmp/hooks/D02_print_environment starting +I: user script /srv/workspace/pbuilder/2992669/tmp/hooks/D01_modify_environment starting +debug: Running on ionos11-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 Feb 26 14:37 /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/2992669/tmp/hooks/D01_modify_environment finished +I: user script /srv/workspace/pbuilder/2992669/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=42 ' - DISTRIBUTION='trixie' - 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]="2" [2]="37" [3]="1" [4]="release" [5]="x86_64-pc-linux-gnu") + BASH_VERSION='5.2.37(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=20 ' + DIRSTACK=() + DISTRIBUTION=trixie + EUID=0 + FUNCNAME=([0]="Echo" [1]="main") + GROUPS=() + HOME=/root + HOSTNAME=i-capture-the-hostname + HOSTTYPE=x86_64 + HOST_ARCH=amd64 IFS=' ' - INVOCATION_ID='65bfe91082f041038a61b02742a235c1' - 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='956707' - PS1='# ' - PS2='> ' + INVOCATION_ID=6186f4e186d44bb6aa5fefa1aead0345 + 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=2992669 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.XlgZPRpl/pbuilderrc_CCO1 --distribution trixie --hookdir /etc/pbuilder/first-build-hooks --debbuildopts -b --basetgz /var/cache/pbuilder/trixie-reproducible-base.tgz --buildresult /srv/reproducible-results/rbuild-debian/r-b-build.XlgZPRpl/b1 --logfile b1/build.log octave-linear-algebra_2.2.3-4.dsc' - SUDO_GID='111' - SUDO_UID='106' - SUDO_USER='jenkins' - TERM='unknown' - TZ='/usr/share/zoneinfo/Etc/GMT+12' - USER='root' - _='/usr/bin/systemd-run' - http_proxy='http://213.165.73.152:3128' + 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.XlgZPRpl/pbuilderrc_OaCI --distribution trixie --hookdir /etc/pbuilder/rebuild-hooks --debbuildopts -b --basetgz /var/cache/pbuilder/trixie-reproducible-base.tgz --buildresult /srv/reproducible-results/rbuild-debian/r-b-build.XlgZPRpl/b2 --logfile b2/build.log octave-linear-algebra_2.2.3-4.dsc' + SUDO_GID=111 + SUDO_UID=106 + SUDO_USER=jenkins + TERM=unknown + TZ=/usr/share/zoneinfo/Etc/GMT-14 + UID=0 + USER=root + _='I: set' + http_proxy=http://46.16.76.132:3128 I: uname -a - Linux ionos15-amd64 6.12.9+bpo-amd64 #1 SMP PREEMPT_DYNAMIC Debian 6.12.9-1~bpo12+1 (2025-01-19) x86_64 GNU/Linux + Linux i-capture-the-hostname 6.1.0-31-amd64 #1 SMP PREEMPT_DYNAMIC Debian 6.1.128-1 (2025-02-07) x86_64 GNU/Linux I: ls -l /bin - lrwxrwxrwx 1 root root 7 Nov 22 2024 /bin -> usr/bin -I: user script /srv/workspace/pbuilder/956707/tmp/hooks/D02_print_environment finished + lrwxrwxrwx 1 root root 7 Nov 22 14:40 /bin -> usr/bin +I: user script /srv/workspace/pbuilder/2992669/tmp/hooks/D02_print_environment finished -> Attempting to satisfy build-dependencies -> Creating pbuilder-satisfydepends-dummy package Package: pbuilder-satisfydepends-dummy @@ -640,7 +672,7 @@ Get: 533 http://deb.debian.org/debian trixie/main amd64 gfortran amd64 4:14.2.0-1 [1436 B] Get: 534 http://deb.debian.org/debian trixie/main amd64 octave-dev amd64 9.4.0-1 [1007 kB] Get: 535 http://deb.debian.org/debian trixie/main amd64 dh-octave all 1.8.0 [22.7 kB] -Fetched 208 MB in 22s (9649 kB/s) +Fetched 208 MB in 19s (11.2 MB/s) Preconfiguring packages ... Selecting previously unselected package netbase. (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 ... 19800 files and directories currently installed.) @@ -2833,7 +2865,11 @@ Building tag database... -> Finished parsing the build-deps 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-4_source.changes +I: user script /srv/workspace/pbuilder/2992669/tmp/hooks/A99_set_merged_usr starting +Not re-configuring usrmerge for trixie +I: user script /srv/workspace/pbuilder/2992669/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-4_source.changes dpkg-buildpackage: info: source package octave-linear-algebra dpkg-buildpackage: info: source version 2.2.3-4 dpkg-buildpackage: info: source distribution unstable @@ -2873,6 +2909,25 @@ Checking package... Run the unit tests... Checking m files ... +[inst/circulant_matrix_vector_product.m] +>>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/circulant_matrix_vector_product.m +***** shared v,x + v = [1 2 3]'; x = [2 5 6]'; +***** assert (circulant_matrix_vector_product(v, x), circulant_make_matrix(v)*x, eps); +warning: `circulant_make_matrix (V)' has been deprecated in favor of `gallery ("circul", V)''. This function will be removed from future versions of the `linear-algebra' package +warning: called from + circulant_make_matrix at line 43 column 5 + __test__ at line 2 column 1 + test at line 682 column 11 + /tmp/tmp.Y63sZ4l5rj at line 14 column 31 + +1 test, 1 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/smwsolve.m] >>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/smwsolve.m ***** test @@ -2884,26 +2939,15 @@ x2 = smwsolve (A, u, v, b); assert (x1, x2, 1e-13); 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); -1 test, 1 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); -warning: `circulant_make_matrix (V)' has been deprecated in favor of `gallery ("circul", V)''. This function will be removed from future versions of the `linear-algebra' package -warning: called from - circulant_make_matrix at line 43 column 5 - __test__ at line 2 column 1 - test at line 682 column 11 - /tmp/tmp.TRDMcvUGdS at line 30 column 31 - +[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/nmf_bpas.m] >>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/nmf_bpas.m @@ -2917,12 +2961,12 @@ Stop: tolerance reached. ***** test [W,H,iter,HIS]=nmf_bpas(A,k,'verbose'); - initGrNormW = 441.51 - initGrNormH = 544.06 - initNorm = 14.240 - SC1 = 1.4013 - SC2 = 700.66 - SC3 = 0.5264 + initGrNormW = 445.22 + initGrNormH = 482.13 + initNorm = 14.367 + SC1 = 1.3125 + SC2 = 656.26 + SC3 = 0.5123 W_density = 1 H_density = 1 par = @@ -2933,8 +2977,8 @@ n = 20 type = regularized nnls_solver = bp - alpha = 0.5074 - beta = 0.5074 + alpha = 0.5111 + beta = 0.5111 max_iter = 100 min_iter = 20 max_time = 1000 @@ -2946,438 +2990,438 @@ scalar structure containing the fields: - initGrNormW = 441.51 - initGrNormH = 544.06 - initNorm = 14.240 - SC1 = 0.027997 - SC2 = 0.027885 - SC3 = 2.4623 - W_density = 0.9933 - H_density = 0.7150 + initGrNormW = 445.22 + initGrNormH = 482.13 + initNorm = 14.367 + SC1 = 0.024190 + SC2 = 0.023368 + SC3 = 1.8555 + W_density = 0.9467 + H_density = 0.6500 iter = 1 - elapsed = 0.012328 - tTotal = 0.012328 - subIterW = 1 - subIterH = 28 - relError = 0.4332 + elapsed = 0.024463 + tTotal = 0.024463 + subIterW = 8 + subIterH = 27 + relError = 0.4143 toshow = scalar structure containing the fields: - initGrNormW = 441.51 - initGrNormH = 544.06 - initNorm = 14.240 - SC1 = 0.015331 - SC2 = 0.014779 - SC3 = 1.2929 - W_density = 0.9467 - H_density = 0.5300 + initGrNormW = 445.22 + initGrNormH = 482.13 + initNorm = 14.367 + SC1 = 0.013962 + SC2 = 0.013068 + SC3 = 1.0214 + W_density = 0.9133 + H_density = 0.5550 iter = 2 - elapsed = 0.010088 - tTotal = 0.022416 - subIterW = 10 + elapsed = 0.021039 + tTotal = 0.045502 + subIterW = 9 subIterH = 24 - relError = 0.3352 + relError = 0.3269 toshow = scalar structure containing the fields: - initGrNormW = 441.51 - initGrNormH = 544.06 - initNorm = 14.240 - SC1 = 9.9382e-03 - SC2 = 9.3816e-03 - SC3 = 0.8444 - W_density = 0.9400 - H_density = 0.5400 + initGrNormW = 445.22 + initGrNormH = 482.13 + initNorm = 14.367 + SC1 = 9.0021e-03 + SC2 = 8.2099e-03 + SC3 = 0.6754 + W_density = 0.8867 + H_density = 0.5800 iter = 3 - elapsed = 9.4120e-03 - tTotal = 0.031828 - subIterW = 2 - subIterH = 14 - relError = 0.3018 + elapsed = 0.019955 + tTotal = 0.065457 + subIterW = 9 + subIterH = 12 + relError = 0.2962 toshow = scalar structure containing the fields: - initGrNormW = 441.51 - initGrNormH = 544.06 - initNorm = 14.240 - SC1 = 7.0134e-03 - SC2 = 6.4243e-03 - SC3 = 0.5892 - W_density = 0.9067 - H_density = 0.6250 + initGrNormW = 445.22 + initGrNormH = 482.13 + initNorm = 14.367 + SC1 = 6.2349e-03 + SC2 = 5.6738e-03 + SC3 = 0.4725 + W_density = 0.8800 + H_density = 0.6300 iter = 4 - elapsed = 9.8710e-03 - tTotal = 0.041699 - subIterW = 9 + elapsed = 0.020425 + tTotal = 0.085882 + subIterW = 6 subIterH = 14 - relError = 0.2851 + relError = 0.2809 toshow = scalar structure containing the fields: - initGrNormW = 441.51 - initGrNormH = 544.06 - initNorm = 14.240 - SC1 = 5.1181e-03 - SC2 = 4.6268e-03 - SC3 = 0.4244 - W_density = 0.8900 - H_density = 0.6400 + initGrNormW = 445.22 + initGrNormH = 482.13 + initNorm = 14.367 + SC1 = 4.6303e-03 + SC2 = 4.1025e-03 + SC3 = 0.3582 + W_density = 0.8633 + H_density = 0.6550 iter = 5 - elapsed = 9.1940e-03 - tTotal = 0.050893 - subIterW = 5 - subIterH = 5 - relError = 0.2760 + elapsed = 0.020272 + tTotal = 0.1062 + subIterW = 7 + subIterH = 8 + relError = 0.2722 toshow = scalar structure containing the fields: - initGrNormW = 441.51 - initGrNormH = 544.06 - initNorm = 14.240 - SC1 = 3.9435e-03 - SC2 = 3.4861e-03 - SC3 = 0.3238 - W_density = 0.8700 - H_density = 0.6750 + initGrNormW = 445.22 + initGrNormH = 482.13 + initNorm = 14.367 + SC1 = 3.5370e-03 + SC2 = 3.0772e-03 + SC3 = 0.2749 + W_density = 0.8500 + H_density = 0.6800 iter = 6 - elapsed = 9.7190e-03 - tTotal = 0.060612 + elapsed = 0.018946 + tTotal = 0.1251 subIterW = 6 - subIterH = 8 - relError = 0.2705 + subIterH = 7 + relError = 0.2670 toshow = scalar structure containing the fields: - initGrNormW = 441.51 - initGrNormH = 544.06 - initNorm = 14.240 - SC1 = 3.0986e-03 - SC2 = 2.7144e-03 - SC3 = 0.2519 - W_density = 0.8633 - H_density = 0.7150 + initGrNormW = 445.22 + initGrNormH = 482.13 + initNorm = 14.367 + SC1 = 2.7830e-03 + SC2 = 2.3822e-03 + SC3 = 0.2152 + W_density = 0.8367 + H_density = 0.6750 iter = 7 - elapsed = 8.1680e-03 - tTotal = 0.068780 + elapsed = 0.015825 + tTotal = 0.1409 subIterW = 4 - subIterH = 6 - relError = 0.2666 + subIterH = 5 + relError = 0.2640 toshow = scalar structure containing the fields: - initGrNormW = 441.51 - initGrNormH = 544.06 - initNorm = 14.240 - SC1 = 2.5447e-03 - SC2 = 2.1579e-03 - SC3 = 0.2042 - W_density = 0.8367 - H_density = 0.7200 + initGrNormW = 445.22 + initGrNormH = 482.13 + initNorm = 14.367 + SC1 = 2.2333e-03 + SC2 = 1.8939e-03 + SC3 = 0.1695 + W_density = 0.8267 + H_density = 0.7000 iter = 8 - elapsed = 8.2370e-03 - tTotal = 0.077017 - subIterW = 7 - subIterH = 3 - relError = 0.2636 + elapsed = 0.016390 + tTotal = 0.1573 + subIterW = 3 + subIterH = 7 + relError = 0.2621 toshow = scalar structure containing the fields: - initGrNormW = 441.51 - initGrNormH = 544.06 - initNorm = 14.240 - SC1 = 2.0757e-03 - SC2 = 1.7519e-03 - SC3 = 0.092489 - W_density = 0.8367 - H_density = 0.7350 + initGrNormW = 445.22 + initGrNormH = 482.13 + initNorm = 14.367 + SC1 = 1.8361e-03 + SC2 = 1.5387e-03 + SC3 = 0.1380 + W_density = 0.8233 + H_density = 0.7100 iter = 9 - elapsed = 7.5310e-03 - tTotal = 0.084548 + elapsed = 0.015375 + tTotal = 0.1727 subIterW = 0 subIterH = 3 - relError = 0.2610 + relError = 0.2606 toshow = scalar structure containing the fields: - initGrNormW = 441.51 - initGrNormH = 544.06 - initNorm = 14.240 - SC1 = 1.7613e-03 - SC2 = 1.4478e-03 - SC3 = 0.1367 - W_density = 0.8167 - H_density = 0.7500 + initGrNormW = 445.22 + initGrNormH = 482.13 + initNorm = 14.367 + SC1 = 1.5344e-03 + SC2 = 1.2582e-03 + SC3 = 0.1142 + W_density = 0.8100 + H_density = 0.7300 iter = 10 - elapsed = 8.3310e-03 - tTotal = 0.092879 - subIterW = 8 - subIterH = 5 - relError = 0.2587 + elapsed = 0.015941 + tTotal = 0.1886 + subIterW = 3 + subIterH = 4 + relError = 0.2595 toshow = scalar structure containing the fields: - initGrNormW = 441.51 - initGrNormH = 544.06 - initNorm = 14.240 - SC1 = 1.4747e-03 - SC2 = 1.1916e-03 - SC3 = 0.1104 - W_density = 0.7933 + initGrNormW = 445.22 + initGrNormH = 482.13 + initNorm = 14.367 + SC1 = 1.2861e-03 + SC2 = 1.0392e-03 + SC3 = 0.093843 + W_density = 0.8000 H_density = 0.7400 iter = 11 - elapsed = 8.1040e-03 - tTotal = 0.1010 - subIterW = 8 + elapsed = 0.016576 + tTotal = 0.2052 + subIterW = 5 subIterH = 2 - relError = 0.2566 + relError = 0.2586 toshow = scalar structure containing the fields: - initGrNormW = 441.51 - initGrNormH = 544.06 - initNorm = 14.240 - SC1 = 1.2284e-03 - SC2 = 9.8271e-04 - SC3 = 0.055963 - W_density = 0.7933 - H_density = 0.7350 + initGrNormW = 445.22 + initGrNormH = 482.13 + initNorm = 14.367 + SC1 = 1.0940e-03 + SC2 = 8.6423e-04 + SC3 = 0.079376 + W_density = 0.7900 + H_density = 0.7400 iter = 12 - elapsed = 6.4320e-03 - tTotal = 0.1074 - subIterW = 0 - subIterH = 6 - relError = 0.2549 + elapsed = 0.016730 + tTotal = 0.2219 + subIterW = 5 + subIterH = 4 + relError = 0.2578 toshow = scalar structure containing the fields: - initGrNormW = 441.51 - initGrNormH = 544.06 - initNorm = 14.240 - SC1 = 1.0356e-03 - SC2 = 8.2019e-04 - SC3 = 0.076257 + initGrNormW = 445.22 + initGrNormH = 482.13 + initNorm = 14.367 + SC1 = 9.2534e-04 + SC2 = 7.2547e-04 + SC3 = 0.065806 W_density = 0.7833 - H_density = 0.7350 + H_density = 0.7450 iter = 13 - elapsed = 6.5820e-03 - tTotal = 0.1140 - subIterW = 5 + elapsed = 0.015322 + tTotal = 0.2373 + subIterW = 2 subIterH = 0 - relError = 0.2534 + relError = 0.2572 toshow = scalar structure containing the fields: - initGrNormW = 441.51 - initGrNormH = 544.06 - initNorm = 14.240 - SC1 = 8.7511e-04 - SC2 = 6.9309e-04 - SC3 = 0.064155 - W_density = 0.7833 - H_density = 0.7550 + initGrNormW = 445.22 + initGrNormH = 482.13 + initNorm = 14.367 + SC1 = 7.8598e-04 + SC2 = 6.0835e-04 + SC3 = 0.054932 + W_density = 0.7700 + H_density = 0.7600 iter = 14 - elapsed = 8.0620e-03 - tTotal = 0.1221 - subIterW = 2 + elapsed = 0.016807 + tTotal = 0.2541 + subIterW = 4 subIterH = 3 - relError = 0.2522 + relError = 0.2568 toshow = scalar structure containing the fields: - initGrNormW = 441.51 - initGrNormH = 544.06 - initNorm = 14.240 - SC1 = 7.4921e-04 - SC2 = 5.9188e-04 - SC3 = 0.054445 - W_density = 0.7800 - H_density = 0.7700 + initGrNormW = 445.22 + initGrNormH = 482.13 + initNorm = 14.367 + SC1 = 6.6310e-04 + SC2 = 5.1324e-04 + SC3 = 0.027792 + W_density = 0.7700 + H_density = 0.7550 iter = 15 - elapsed = 8.0260e-03 - tTotal = 0.1301 - subIterW = 3 - subIterH = 4 - relError = 0.2512 + elapsed = 0.012360 + tTotal = 0.2664 + subIterW = 0 + subIterH = 0 + relError = 0.2565 toshow = scalar structure containing the fields: - initGrNormW = 441.51 - initGrNormH = 544.06 - initNorm = 14.240 - SC1 = 6.4149e-04 - SC2 = 5.0293e-04 - SC3 = 0.045922 - W_density = 0.7733 - H_density = 0.7750 + initGrNormW = 445.22 + initGrNormH = 482.13 + initNorm = 14.367 + SC1 = 5.6303e-04 + SC2 = 4.3354e-04 + SC3 = 0.038565 + W_density = 0.7667 + H_density = 0.7550 iter = 16 - elapsed = 7.9320e-03 - tTotal = 0.1380 - subIterW = 5 - subIterH = 0 - relError = 0.2504 + elapsed = 0.014721 + tTotal = 0.2811 + subIterW = 0 + subIterH = 2 + relError = 0.2563 toshow = scalar structure containing the fields: - initGrNormW = 441.51 - initGrNormH = 544.06 - initNorm = 14.240 - SC1 = 5.4634e-04 - SC2 = 4.2942e-04 - SC3 = 0.038433 - W_density = 0.7767 - H_density = 0.7750 + initGrNormW = 445.22 + initGrNormH = 482.13 + initNorm = 14.367 + SC1 = 4.8181e-04 + SC2 = 3.7099e-04 + SC3 = 0.019395 + W_density = 0.7667 + H_density = 0.7550 iter = 17 - elapsed = 6.0010e-03 - tTotal = 0.1440 + elapsed = 8.4170e-03 + tTotal = 0.2896 subIterW = 0 subIterH = 0 - relError = 0.2498 + relError = 0.2561 toshow = scalar structure containing the fields: - initGrNormW = 441.51 - initGrNormH = 544.06 - initNorm = 14.240 - SC1 = 4.7360e-04 - SC2 = 3.7415e-04 - SC3 = 0.032716 - W_density = 0.7733 - H_density = 0.7750 + initGrNormW = 445.22 + initGrNormH = 482.13 + initNorm = 14.367 + SC1 = 4.1797e-04 + SC2 = 3.2184e-04 + SC3 = 0.017037 + W_density = 0.7667 + H_density = 0.7550 iter = 18 - elapsed = 6.0060e-03 - tTotal = 0.1500 + elapsed = 9.1260e-03 + tTotal = 0.2987 subIterW = 0 subIterH = 0 - relError = 0.2492 + relError = 0.2559 toshow = scalar structure containing the fields: - initGrNormW = 441.51 - initGrNormH = 544.06 - initNorm = 14.240 - SC1 = 4.1486e-04 - SC2 = 3.2857e-04 - SC3 = 0.028553 - W_density = 0.7767 - H_density = 0.7850 + initGrNormW = 445.22 + initGrNormH = 482.13 + initNorm = 14.367 + SC1 = 3.6977e-04 + SC2 = 2.8324e-04 + SC3 = 0.024617 + W_density = 0.7633 + H_density = 0.7650 iter = 19 - elapsed = 7.6990e-03 - tTotal = 0.1577 + elapsed = 0.015566 + tTotal = 0.3143 subIterW = 0 subIterH = 2 - relError = 0.2487 + relError = 0.2558 toshow = scalar structure containing the fields: - initGrNormW = 441.51 - initGrNormH = 544.06 - initNorm = 14.240 - SC1 = 3.6944e-04 - SC2 = 2.9186e-04 - SC3 = 0.025152 - W_density = 0.7733 - H_density = 0.7900 + initGrNormW = 445.22 + initGrNormH = 482.13 + initNorm = 14.367 + SC1 = 3.2921e-04 + SC2 = 2.5086e-04 + SC3 = 0.021741 + W_density = 0.7600 + H_density = 0.7600 iter = 20 - elapsed = 7.4530e-03 - tTotal = 0.1652 + elapsed = 0.014550 + tTotal = 0.3288 subIterW = 0 subIterH = 0 - relError = 0.2483 + relError = 0.2556 toshow = scalar structure containing the fields: - initGrNormW = 441.51 - initGrNormH = 544.06 - initNorm = 14.240 - SC1 = 3.3095e-04 - SC2 = 2.6212e-04 - SC3 = 0.022212 - W_density = 0.7767 - H_density = 0.8100 + initGrNormW = 445.22 + initGrNormH = 482.13 + initNorm = 14.367 + SC1 = 2.9694e-04 + SC2 = 2.2568e-04 + SC3 = 0.011404 + W_density = 0.7600 + H_density = 0.7600 iter = 21 - elapsed = 7.8100e-03 - tTotal = 0.1730 + elapsed = 0.012616 + tTotal = 0.3414 subIterW = 0 - subIterH = 3 - relError = 0.2480 + subIterH = 2 + relError = 0.2555 toshow = scalar structure containing the fields: - initGrNormW = 441.51 - initGrNormH = 544.06 - initNorm = 14.240 - SC1 = 2.9961e-04 - SC2 = 2.3729e-04 - SC3 = 0.019847 - W_density = 0.7767 - H_density = 0.8100 + initGrNormW = 445.22 + initGrNormH = 482.13 + initNorm = 14.367 + SC1 = 2.7093e-04 + SC2 = 2.0591e-04 + SC3 = 0.017398 + W_density = 0.7567 + H_density = 0.7600 iter = 22 - elapsed = 6.7830e-03 - tTotal = 0.1798 - subIterW = 4 + elapsed = 0.013006 + tTotal = 0.3544 + subIterW = 3 subIterH = 0 - relError = 0.2477 + relError = 0.2554 toshow = scalar structure containing the fields: - initGrNormW = 441.51 - initGrNormH = 544.06 - initNorm = 14.240 - SC1 = 2.6898e-04 - SC2 = 2.1357e-04 - SC3 = 0.017507 - W_density = 0.7767 - H_density = 0.8100 + initGrNormW = 445.22 + initGrNormH = 482.13 + initNorm = 14.367 + SC1 = 2.5124e-04 + SC2 = 1.9094e-04 + SC3 = 9.0359e-03 + W_density = 0.7567 + H_density = 0.7600 iter = 23 - elapsed = 0.013030 - tTotal = 0.1928 - subIterW = 2 + elapsed = 8.8940e-03 + tTotal = 0.3633 + subIterW = 0 subIterH = 0 - relError = 0.2474 + relError = 0.2553 Stop: tolerance reached. final = @@ -3385,19 +3429,19 @@ scalar structure containing the fields: iterations = 23 - elapsed_total = 0.1928 - relative_error = 0.2474 - W_density = 0.7767 - H_density = 0.8100 + elapsed_total = 0.3633 + relative_error = 0.2553 + W_density = 0.7567 + H_density = 0.7600 ***** test [W,H,iter,HIS]=nmf_bpas(A,k,'verbose','nnlssolver','as'); - initGrNormW = 353.39 - initGrNormH = 440.78 - initNorm = 14.240 - SC1 = 1.1299 - SC2 = 564.95 - SC3 = 0.4928 + initGrNormW = 403.56 + initGrNormH = 521.61 + initNorm = 14.367 + SC1 = 1.3190 + SC2 = 659.50 + SC3 = 0.5200 W_density = 1 H_density = 1 par = @@ -3408,8 +3452,8 @@ n = 20 type = regularized nnls_solver = as - alpha = 0.5074 - beta = 0.5074 + alpha = 0.5111 + beta = 0.5111 max_iter = 100 min_iter = 20 max_time = 1000 @@ -3421,438 +3465,438 @@ scalar structure containing the fields: - initGrNormW = 353.39 - initGrNormH = 440.78 - initNorm = 14.240 - SC1 = 0.030417 - SC2 = 0.029322 - SC3 = 1.1950 - W_density = 0.9500 + initGrNormW = 403.56 + initGrNormH = 521.61 + initNorm = 14.367 + SC1 = 0.028376 + SC2 = 0.027752 + SC3 = 1.2022 + W_density = 0.9667 H_density = 0.6750 iter = 1 - elapsed = 0.010833 - tTotal = 0.010833 - subIterW = 10 - subIterH = 29 - relError = 0.4268 + elapsed = 0.015598 + tTotal = 0.015598 + subIterW = 5 + subIterH = 26 + relError = 0.4352 toshow = scalar structure containing the fields: - initGrNormW = 353.39 - initGrNormH = 440.78 - initNorm = 14.240 - SC1 = 0.017468 - SC2 = 0.016245 - SC3 = 0.6375 - W_density = 0.9233 - H_density = 0.5600 + initGrNormW = 403.56 + initGrNormH = 521.61 + initNorm = 14.367 + SC1 = 0.015592 + SC2 = 0.014688 + SC3 = 0.6507 + W_density = 0.9200 + H_density = 0.5450 iter = 2 - elapsed = 9.2050e-03 - tTotal = 0.020038 - subIterW = 13 - subIterH = 36 - relError = 0.3376 + elapsed = 0.016674 + tTotal = 0.032272 + subIterW = 12 + subIterH = 35 + relError = 0.3386 toshow = scalar structure containing the fields: - initGrNormW = 353.39 - initGrNormH = 440.78 - initNorm = 14.240 - SC1 = 0.010663 - SC2 = 9.6605e-03 - SC3 = 0.4137 - W_density = 0.8967 - H_density = 0.5600 + initGrNormW = 403.56 + initGrNormH = 521.61 + initNorm = 14.367 + SC1 = 9.9606e-03 + SC2 = 9.1638e-03 + SC3 = 0.4407 + W_density = 0.8900 + H_density = 0.5950 iter = 3 - elapsed = 0.011910 - tTotal = 0.031948 + elapsed = 0.018337 + tTotal = 0.050609 subIterW = 18 - subIterH = 40 - relError = 0.3044 + subIterH = 34 + relError = 0.3005 toshow = scalar structure containing the fields: - initGrNormW = 353.39 - initGrNormH = 440.78 - initNorm = 14.240 - SC1 = 7.4864e-03 - SC2 = 6.6779e-03 - SC3 = 0.2853 - W_density = 0.8900 - H_density = 0.5850 + initGrNormW = 403.56 + initGrNormH = 521.61 + initNorm = 14.367 + SC1 = 6.8576e-03 + SC2 = 6.1582e-03 + SC3 = 0.3029 + W_density = 0.8633 + H_density = 0.6200 iter = 4 - elapsed = 0.011908 - tTotal = 0.043856 - subIterW = 19 - subIterH = 38 - relError = 0.2865 + elapsed = 0.017481 + tTotal = 0.068090 + subIterW = 20 + subIterH = 31 + relError = 0.2822 toshow = scalar structure containing the fields: - initGrNormW = 353.39 - initGrNormH = 440.78 - initNorm = 14.240 - SC1 = 5.6027e-03 - SC2 = 4.9304e-03 - SC3 = 0.2150 - W_density = 0.8867 - H_density = 0.6100 + initGrNormW = 403.56 + initGrNormH = 521.61 + initNorm = 14.367 + SC1 = 5.0343e-03 + SC2 = 4.4504e-03 + SC3 = 0.2337 + W_density = 0.8567 + H_density = 0.6250 iter = 5 - elapsed = 9.9670e-03 - tTotal = 0.053823 - subIterW = 17 - subIterH = 35 - relError = 0.2749 + elapsed = 0.018812 + tTotal = 0.086902 + subIterW = 21 + subIterH = 36 + relError = 0.2731 toshow = scalar structure containing the fields: - initGrNormW = 353.39 - initGrNormH = 440.78 - initNorm = 14.240 - SC1 = 4.4357e-03 - SC2 = 3.7615e-03 - SC3 = 0.1631 - W_density = 0.8533 - H_density = 0.6450 + initGrNormW = 403.56 + initGrNormH = 521.61 + initNorm = 14.367 + SC1 = 3.8720e-03 + SC2 = 3.3609e-03 + SC3 = 0.1673 + W_density = 0.8467 + H_density = 0.6550 iter = 6 - elapsed = 0.010855 - tTotal = 0.064678 - subIterW = 24 - subIterH = 32 - relError = 0.2672 + elapsed = 0.018193 + tTotal = 0.1051 + subIterW = 23 + subIterH = 31 + relError = 0.2678 toshow = scalar structure containing the fields: - initGrNormW = 353.39 - initGrNormH = 440.78 - initNorm = 14.240 - SC1 = 3.5028e-03 - SC2 = 2.9073e-03 - SC3 = 0.1354 - W_density = 0.8300 - H_density = 0.6750 + initGrNormW = 403.56 + initGrNormH = 521.61 + initNorm = 14.367 + SC1 = 3.1003e-03 + SC2 = 2.5980e-03 + SC3 = 0.1417 + W_density = 0.8233 + H_density = 0.6800 iter = 7 - elapsed = 8.8960e-03 - tTotal = 0.073574 - subIterW = 27 - subIterH = 28 - relError = 0.2622 + elapsed = 0.020740 + tTotal = 0.1258 + subIterW = 28 + subIterH = 29 + relError = 0.2643 toshow = scalar structure containing the fields: - initGrNormW = 353.39 - initGrNormH = 440.78 - initNorm = 14.240 - SC1 = 2.7976e-03 - SC2 = 2.2828e-03 - SC3 = 0.1018 - W_density = 0.8200 - H_density = 0.6850 + initGrNormW = 403.56 + initGrNormH = 521.61 + initNorm = 14.367 + SC1 = 2.5141e-03 + SC2 = 2.0515e-03 + SC3 = 0.1099 + W_density = 0.8133 + H_density = 0.7050 iter = 8 - elapsed = 9.3730e-03 - tTotal = 0.082947 - subIterW = 27 - subIterH = 33 - relError = 0.2587 + elapsed = 0.019092 + tTotal = 0.1449 + subIterW = 31 + subIterH = 25 + relError = 0.2621 toshow = scalar structure containing the fields: - initGrNormW = 353.39 - initGrNormH = 440.78 - initNorm = 14.240 - SC1 = 2.2826e-03 - SC2 = 1.8489e-03 - SC3 = 0.084514 - W_density = 0.8100 - H_density = 0.6850 + initGrNormW = 403.56 + initGrNormH = 521.61 + initNorm = 14.367 + SC1 = 2.0401e-03 + SC2 = 1.6565e-03 + SC3 = 0.090528 + W_density = 0.8067 + H_density = 0.7150 iter = 9 - elapsed = 9.6330e-03 - tTotal = 0.092580 - subIterW = 28 - subIterH = 31 - relError = 0.2559 + elapsed = 0.017982 + tTotal = 0.1629 + subIterW = 30 + subIterH = 27 + relError = 0.2606 toshow = scalar structure containing the fields: - initGrNormW = 353.39 - initGrNormH = 440.78 - initNorm = 14.240 - SC1 = 1.9283e-03 - SC2 = 1.5349e-03 - SC3 = 0.069489 + initGrNormW = 403.56 + initGrNormH = 521.61 + initNorm = 14.367 + SC1 = 1.6824e-03 + SC2 = 1.3527e-03 + SC3 = 0.074424 W_density = 0.8033 - H_density = 0.6950 + H_density = 0.7250 iter = 10 - elapsed = 8.8740e-03 - tTotal = 0.1015 - subIterW = 31 - subIterH = 27 - relError = 0.2536 + elapsed = 0.016510 + tTotal = 0.1794 + subIterW = 30 + subIterH = 28 + relError = 0.2595 toshow = scalar structure containing the fields: - initGrNormW = 353.39 - initGrNormH = 440.78 - initNorm = 14.240 - SC1 = 1.6389e-03 - SC2 = 1.2881e-03 - SC3 = 0.059951 - W_density = 0.7900 - H_density = 0.7150 + initGrNormW = 403.56 + initGrNormH = 521.61 + initNorm = 14.367 + SC1 = 1.3909e-03 + SC2 = 1.1100e-03 + SC3 = 0.060076 + W_density = 0.7933 + H_density = 0.7300 iter = 11 - elapsed = 9.6070e-03 - tTotal = 0.1111 - subIterW = 34 - subIterH = 26 - relError = 0.2514 + elapsed = 0.016860 + tTotal = 0.1963 + subIterW = 32 + subIterH = 27 + relError = 0.2588 toshow = scalar structure containing the fields: - initGrNormW = 353.39 - initGrNormH = 440.78 - initNorm = 14.240 - SC1 = 1.3833e-03 - SC2 = 1.0762e-03 - SC3 = 0.051268 - W_density = 0.7733 - H_density = 0.7150 + initGrNormW = 403.56 + initGrNormH = 521.61 + initNorm = 14.367 + SC1 = 1.1641e-03 + SC2 = 9.1728e-04 + SC3 = 0.049549 + W_density = 0.7833 + H_density = 0.7300 iter = 12 - elapsed = 9.6400e-03 - tTotal = 0.1207 - subIterW = 35 - subIterH = 27 - relError = 0.2498 + elapsed = 0.018191 + tTotal = 0.2145 + subIterW = 33 + subIterH = 28 + relError = 0.2582 toshow = scalar structure containing the fields: - initGrNormW = 353.39 - initGrNormH = 440.78 - initNorm = 14.240 - SC1 = 1.1765e-03 - SC2 = 9.0120e-04 - SC3 = 0.041531 - W_density = 0.7667 - H_density = 0.7150 + initGrNormW = 403.56 + initGrNormH = 521.61 + initNorm = 14.367 + SC1 = 9.7789e-04 + SC2 = 7.6667e-04 + SC3 = 0.039289 + W_density = 0.7833 + H_density = 0.7300 iter = 13 - elapsed = 9.6160e-03 - tTotal = 0.1303 - subIterW = 34 - subIterH = 25 - relError = 0.2486 + elapsed = 0.017121 + tTotal = 0.2316 + subIterW = 33 + subIterH = 28 + relError = 0.2577 toshow = scalar structure containing the fields: - initGrNormW = 353.39 - initGrNormH = 440.78 - initNorm = 14.240 - SC1 = 9.9215e-04 - SC2 = 7.5602e-04 - SC3 = 0.033741 - W_density = 0.7633 - H_density = 0.7150 + initGrNormW = 403.56 + initGrNormH = 521.61 + initNorm = 14.367 + SC1 = 8.4013e-04 + SC2 = 6.5026e-04 + SC3 = 0.035191 + W_density = 0.7767 + H_density = 0.7300 iter = 14 - elapsed = 9.3770e-03 - tTotal = 0.1397 - subIterW = 34 - subIterH = 23 - relError = 0.2478 + elapsed = 0.018429 + tTotal = 0.2500 + subIterW = 33 + subIterH = 29 + relError = 0.2572 toshow = scalar structure containing the fields: - initGrNormW = 353.39 - initGrNormH = 440.78 - initNorm = 14.240 - SC1 = 8.3650e-04 - SC2 = 6.3742e-04 - SC3 = 0.029253 - W_density = 0.7633 - H_density = 0.7200 + initGrNormW = 403.56 + initGrNormH = 521.61 + initNorm = 14.367 + SC1 = 7.3224e-04 + SC2 = 5.5650e-04 + SC3 = 0.029912 + W_density = 0.7600 + H_density = 0.7300 iter = 15 - elapsed = 8.6830e-03 - tTotal = 0.1484 - subIterW = 33 - subIterH = 24 - relError = 0.2472 + elapsed = 0.022188 + tTotal = 0.2722 + subIterW = 39 + subIterH = 29 + relError = 0.2568 toshow = scalar structure containing the fields: - initGrNormW = 353.39 - initGrNormH = 440.78 - initNorm = 14.240 - SC1 = 7.0302e-04 - SC2 = 5.3430e-04 - SC3 = 0.024351 - W_density = 0.7600 - H_density = 0.7200 + initGrNormW = 403.56 + initGrNormH = 521.61 + initNorm = 14.367 + SC1 = 6.3645e-04 + SC2 = 4.8243e-04 + SC3 = 0.026046 + W_density = 0.7567 + H_density = 0.7350 iter = 16 - elapsed = 9.0130e-03 - tTotal = 0.1574 - subIterW = 35 - subIterH = 25 - relError = 0.2467 + elapsed = 0.019054 + tTotal = 0.2913 + subIterW = 39 + subIterH = 28 + relError = 0.2564 toshow = scalar structure containing the fields: - initGrNormW = 353.39 - initGrNormH = 440.78 - initNorm = 14.240 - SC1 = 5.8119e-04 - SC2 = 4.4170e-04 - SC3 = 0.020602 + initGrNormW = 403.56 + initGrNormH = 521.61 + initNorm = 14.367 + SC1 = 5.5944e-04 + SC2 = 4.2406e-04 + SC3 = 0.021312 W_density = 0.7567 - H_density = 0.7250 + H_density = 0.7400 iter = 17 - elapsed = 9.9140e-03 - tTotal = 0.1673 - subIterW = 36 - subIterH = 25 - relError = 0.2463 + elapsed = 0.018083 + tTotal = 0.3093 + subIterW = 39 + subIterH = 26 + relError = 0.2561 toshow = scalar structure containing the fields: - initGrNormW = 353.39 - initGrNormH = 440.78 - initNorm = 14.240 - SC1 = 4.8776e-04 - SC2 = 3.7167e-04 - SC3 = 0.016301 - W_density = 0.7600 - H_density = 0.7450 + initGrNormW = 403.56 + initGrNormH = 521.61 + initNorm = 14.367 + SC1 = 5.0446e-04 + SC2 = 3.7733e-04 + SC3 = 0.018989 + W_density = 0.7433 + H_density = 0.7400 iter = 18 - elapsed = 8.6480e-03 - tTotal = 0.1760 - subIterW = 33 - subIterH = 25 - relError = 0.2460 + elapsed = 0.018463 + tTotal = 0.3278 + subIterW = 42 + subIterH = 26 + relError = 0.2557 toshow = scalar structure containing the fields: - initGrNormW = 353.39 - initGrNormH = 440.78 - initNorm = 14.240 - SC1 = 4.1418e-04 - SC2 = 3.1478e-04 - SC3 = 0.014560 - W_density = 0.7567 - H_density = 0.7450 + initGrNormW = 403.56 + initGrNormH = 521.61 + initNorm = 14.367 + SC1 = 4.3932e-04 + SC2 = 3.2949e-04 + SC3 = 0.017190 + W_density = 0.7433 + H_density = 0.7500 iter = 19 - elapsed = 8.8690e-03 - tTotal = 0.1848 - subIterW = 34 - subIterH = 24 - relError = 0.2458 + elapsed = 0.018281 + tTotal = 0.3461 + subIterW = 40 + subIterH = 26 + relError = 0.2555 toshow = scalar structure containing the fields: - initGrNormW = 353.39 - initGrNormH = 440.78 - initNorm = 14.240 - SC1 = 3.5138e-04 - SC2 = 2.6705e-04 - SC3 = 0.011430 - W_density = 0.7600 - H_density = 0.7500 + initGrNormW = 403.56 + initGrNormH = 521.61 + initNorm = 14.367 + SC1 = 3.9207e-04 + SC2 = 2.9405e-04 + SC3 = 0.014729 + W_density = 0.7433 + H_density = 0.7450 iter = 20 - elapsed = 6.4600e-03 - tTotal = 0.1913 - subIterW = 34 - subIterH = 25 - relError = 0.2457 + elapsed = 0.019308 + tTotal = 0.3654 + subIterW = 39 + subIterH = 27 + relError = 0.2552 toshow = scalar structure containing the fields: - initGrNormW = 353.39 - initGrNormH = 440.78 - initNorm = 14.240 - SC1 = 3.0330e-04 - SC2 = 2.3111e-04 - SC3 = 0.010187 - W_density = 0.7600 - H_density = 0.7500 + initGrNormW = 403.56 + initGrNormH = 521.61 + initNorm = 14.367 + SC1 = 3.5367e-04 + SC2 = 2.6596e-04 + SC3 = 0.012698 + W_density = 0.7433 + H_density = 0.7400 iter = 21 - elapsed = 6.8520e-03 - tTotal = 0.1981 - subIterW = 32 - subIterH = 26 - relError = 0.2455 + elapsed = 0.020510 + tTotal = 0.3859 + subIterW = 39 + subIterH = 31 + relError = 0.2548 toshow = scalar structure containing the fields: - initGrNormW = 353.39 - initGrNormH = 440.78 - initNorm = 14.240 - SC1 = 2.6723e-04 - SC2 = 2.0309e-04 - SC3 = 9.1690e-03 - W_density = 0.7533 - H_density = 0.7500 + initGrNormW = 403.56 + initGrNormH = 521.61 + initNorm = 14.367 + SC1 = 3.2453e-04 + SC2 = 2.4275e-04 + SC3 = 0.012290 + W_density = 0.7367 + H_density = 0.7350 iter = 22 - elapsed = 7.7730e-03 - tTotal = 0.2059 - subIterW = 34 - subIterH = 26 - relError = 0.2454 + elapsed = 0.021762 + tTotal = 0.4077 + subIterW = 42 + subIterH = 31 + relError = 0.2545 toshow = scalar structure containing the fields: - initGrNormW = 353.39 - initGrNormH = 440.78 - initNorm = 14.240 - SC1 = 2.3642e-04 - SC2 = 1.7968e-04 - SC3 = 7.4565e-03 - W_density = 0.7533 - H_density = 0.7550 + initGrNormW = 403.56 + initGrNormH = 521.61 + initNorm = 14.367 + SC1 = 2.9943e-04 + SC2 = 2.2277e-04 + SC3 = 0.010242 + W_density = 0.7333 + H_density = 0.7350 iter = 23 - elapsed = 9.0130e-03 - tTotal = 0.2149 - subIterW = 37 - subIterH = 25 - relError = 0.2453 + elapsed = 0.021733 + tTotal = 0.4294 + subIterW = 43 + subIterH = 32 + relError = 0.2542 Stop: tolerance reached. final = @@ -3860,19 +3904,19 @@ scalar structure containing the fields: iterations = 23 - elapsed_total = 0.2149 - relative_error = 0.2453 - W_density = 0.7533 - H_density = 0.7550 + elapsed_total = 0.4294 + relative_error = 0.2542 + W_density = 0.7333 + H_density = 0.7350 ***** test [W,H,iter,HIS]=nmf_bpas(A,k,'verbose','type','sparse'); - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 1.0998 - SC2 = 549.89 - SC3 = 0.4822 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 1.1643 + SC2 = 582.16 + SC3 = 0.4945 W_density = 1 H_density = 1 par = @@ -3883,8 +3927,8 @@ n = 20 type = sparse nnls_solver = bp - alpha = 0.5074 - beta = 0.5074 + alpha = 0.5111 + beta = 0.5111 max_iter = 100 min_iter = 20 max_time = 1000 @@ -3896,572 +3940,553 @@ scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 0.022359 - SC2 = 0.020347 - SC3 = 1.4972 - W_density = 0.9600 - H_density = 0.6200 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 0.023195 + SC2 = 0.021710 + SC3 = 1.5969 + W_density = 0.9667 + H_density = 0.5850 iter = 1 - elapsed = 0.010539 - tTotal = 0.010539 - subIterW = 7 - subIterH = 28 - relError = 0.4200 + elapsed = 0.024016 + tTotal = 0.024016 + subIterW = 5 + subIterH = 33 + relError = 0.4166 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 0.011743 - SC2 = 8.9245e-03 - SC3 = 0.7712 - W_density = 0.8467 - H_density = 0.3700 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 0.013214 + SC2 = 0.010281 + SC3 = 0.8879 + W_density = 0.8867 + H_density = 0.3450 iter = 2 - elapsed = 0.012618 - tTotal = 0.023157 - subIterW = 6 - subIterH = 32 - relError = 0.3455 + elapsed = 0.028042 + tTotal = 0.052058 + subIterW = 14 + subIterH = 37 + relError = 0.3504 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 7.1799e-03 - SC2 = 5.0547e-03 - SC3 = 0.4836 - W_density = 0.8433 - H_density = 0.3200 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 8.2647e-03 + SC2 = 5.4382e-03 + SC3 = 0.6110 + W_density = 0.8033 + H_density = 0.3000 iter = 3 - elapsed = 8.5040e-03 - tTotal = 0.031661 - subIterW = 6 - subIterH = 11 - relError = 0.3277 + elapsed = 0.020713 + tTotal = 0.072771 + subIterW = 10 + subIterH = 15 + relError = 0.3314 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 4.8836e-03 - SC2 = 3.2134e-03 - SC3 = 0.3466 - W_density = 0.8400 - H_density = 0.3000 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 5.5241e-03 + SC2 = 3.4581e-03 + SC3 = 0.4030 + W_density = 0.7833 + H_density = 0.2950 iter = 4 - elapsed = 6.6690e-03 - tTotal = 0.038330 - subIterW = 3 - subIterH = 8 - relError = 0.3218 + elapsed = 0.011534 + tTotal = 0.084305 + subIterW = 5 + subIterH = 3 + relError = 0.3269 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 3.5893e-03 - SC2 = 2.2469e-03 - SC3 = 0.2549 - W_density = 0.8200 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 3.9863e-03 + SC2 = 2.4078e-03 + SC3 = 0.2949 + W_density = 0.7767 H_density = 0.2900 iter = 5 - elapsed = 8.1280e-03 - tTotal = 0.046458 - subIterW = 6 - subIterH = 6 - relError = 0.3186 + elapsed = 0.012906 + tTotal = 0.097211 + subIterW = 2 + subIterH = 5 + relError = 0.3254 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 2.6761e-03 - SC2 = 1.6699e-03 - SC3 = 0.1817 - W_density = 0.8267 - H_density = 0.2750 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 3.0008e-03 + SC2 = 1.8005e-03 + SC3 = 0.2110 + W_density = 0.7800 + H_density = 0.2850 iter = 6 - elapsed = 6.6920e-03 - tTotal = 0.053150 - subIterW = 2 - subIterH = 4 - relError = 0.3174 + elapsed = 0.012399 + tTotal = 0.1096 + subIterW = 3 + subIterH = 0 + relError = 0.3246 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 2.2927e-03 - SC2 = 1.3894e-03 - SC3 = 0.1563 - W_density = 0.8100 - H_density = 0.2600 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 2.1975e-03 + SC2 = 1.2745e-03 + SC3 = 0.1559 + W_density = 0.7700 + H_density = 0.2800 iter = 7 - elapsed = 6.8840e-03 - tTotal = 0.060034 - subIterW = 6 - subIterH = 5 - relError = 0.3167 + elapsed = 0.012967 + tTotal = 0.1226 + subIterW = 4 + subIterH = 3 + relError = 0.3232 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 2.1718e-03 - SC2 = 1.2987e-03 - SC3 = 0.1472 - W_density = 0.8000 - H_density = 0.2600 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 1.5938e-03 + SC2 = 9.1485e-04 + SC3 = 0.035162 + W_density = 0.7700 + H_density = 0.2550 iter = 8 - elapsed = 6.4620e-03 - tTotal = 0.066496 - subIterW = 2 - subIterH = 2 - relError = 0.3169 + elapsed = 9.6140e-03 + tTotal = 0.1322 + subIterW = 0 + subIterH = 5 + relError = 0.3226 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 2.3164e-03 - SC2 = 1.3852e-03 - SC3 = 0.1513 - W_density = 0.7967 - H_density = 0.2700 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 1.4012e-03 + SC2 = 7.9587e-04 + SC3 = 0.098475 + W_density = 0.7633 + H_density = 0.2650 iter = 9 - elapsed = 6.6410e-03 - tTotal = 0.073137 - subIterW = 3 + elapsed = 0.012518 + tTotal = 0.1447 + subIterW = 2 subIterH = 4 - relError = 0.3174 + relError = 0.3229 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 2.4540e-03 - SC2 = 1.4381e-03 - SC3 = 0.1755 - W_density = 0.7900 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 1.3833e-03 + SC2 = 7.8297e-04 + SC3 = 0.092874 + W_density = 0.7633 H_density = 0.2600 iter = 10 - elapsed = 6.7990e-03 - tTotal = 0.079936 + elapsed = 0.011992 + tTotal = 0.1567 subIterW = 2 - subIterH = 6 - relError = 0.3177 + subIterH = 0 + relError = 0.3233 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 2.5823e-03 - SC2 = 1.4978e-03 - SC3 = 0.1676 - W_density = 0.7767 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 1.3283e-03 + SC2 = 7.4121e-04 + SC3 = 0.089682 + W_density = 0.7533 H_density = 0.2550 iter = 11 - elapsed = 6.4640e-03 - tTotal = 0.086400 + elapsed = 0.012771 + tTotal = 0.1695 subIterW = 4 subIterH = 0 - relError = 0.3182 + relError = 0.3238 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 2.5961e-03 - SC2 = 1.4798e-03 - SC3 = 0.1721 - W_density = 0.7700 - H_density = 0.2450 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 1.3528e-03 + SC2 = 7.4943e-04 + SC3 = 0.089351 + W_density = 0.7433 + H_density = 0.2600 iter = 12 - elapsed = 6.7420e-03 - tTotal = 0.093142 - subIterW = 2 - subIterH = 3 - relError = 0.3193 + elapsed = 0.013444 + tTotal = 0.1829 + subIterW = 3 + subIterH = 0 + relError = 0.3243 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 2.4744e-03 - SC2 = 1.3956e-03 - SC3 = 0.1621 - W_density = 0.7667 - H_density = 0.2400 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 1.3294e-03 + SC2 = 7.3384e-04 + SC3 = 0.027297 + W_density = 0.7433 + H_density = 0.2550 iter = 13 - elapsed = 6.8200e-03 - tTotal = 0.099962 - subIterW = 3 - subIterH = 3 - relError = 0.3207 + elapsed = 6.5650e-03 + tTotal = 0.1895 + subIterW = 0 + subIterH = 0 + relError = 0.3245 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 2.3961e-03 - SC2 = 1.3322e-03 - SC3 = 0.1566 - W_density = 0.7600 - H_density = 0.2450 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 1.3463e-03 + SC2 = 7.4317e-04 + SC3 = 0.090442 + W_density = 0.7400 + H_density = 0.2550 iter = 14 - elapsed = 7.6980e-03 - tTotal = 0.1077 - subIterW = 2 - subIterH = 3 - relError = 0.3223 + elapsed = 8.3130e-03 + tTotal = 0.1978 + subIterW = 3 + subIterH = 0 + relError = 0.3246 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 2.3578e-03 - SC2 = 1.2921e-03 - SC3 = 0.1516 - W_density = 0.7533 - H_density = 0.2350 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 1.3400e-03 + SC2 = 7.2898e-04 + SC3 = 0.093182 + W_density = 0.7367 + H_density = 0.2500 iter = 15 - elapsed = 6.5910e-03 - tTotal = 0.1143 - subIterW = 2 - subIterH = 4 - relError = 0.3241 + elapsed = 9.9440e-03 + tTotal = 0.2077 + subIterW = 3 + subIterH = 5 + relError = 0.3250 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 2.3540e-03 - SC2 = 1.2712e-03 - SC3 = 0.1460 - W_density = 0.7367 - H_density = 0.2400 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 1.3779e-03 + SC2 = 7.4132e-04 + SC3 = 0.091480 + W_density = 0.7300 + H_density = 0.2500 iter = 16 - elapsed = 6.7710e-03 - tTotal = 0.1210 - subIterW = 7 - subIterH = 0 - relError = 0.3260 + elapsed = 0.010916 + tTotal = 0.2187 + subIterW = 2 + subIterH = 2 + relError = 0.3257 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 2.3660e-03 - SC2 = 1.2729e-03 - SC3 = 0.1456 - W_density = 0.7333 - H_density = 0.2350 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 1.4245e-03 + SC2 = 7.6070e-04 + SC3 = 0.093081 + W_density = 0.7300 + H_density = 0.2400 iter = 17 - elapsed = 6.2940e-03 - tTotal = 0.1273 - subIterW = 0 - subIterH = 0 - relError = 0.3284 + elapsed = 0.010371 + tTotal = 0.2290 + subIterW = 2 + subIterH = 2 + relError = 0.3265 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 2.3929e-03 - SC2 = 1.2634e-03 - SC3 = 0.1500 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 1.4334e-03 + SC2 = 7.6257e-04 + SC3 = 0.096669 W_density = 0.7267 - H_density = 0.2150 + H_density = 0.2350 iter = 18 - elapsed = 8.0720e-03 - tTotal = 0.1354 - subIterW = 2 - subIterH = 3 - relError = 0.3308 + elapsed = 8.7710e-03 + tTotal = 0.2378 + subIterW = 0 + subIterH = 0 + relError = 0.3275 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 2.1237e-03 - SC2 = 1.0958e-03 - SC3 = 0.1336 - W_density = 0.7133 - H_density = 0.2100 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 1.4487e-03 + SC2 = 7.6780e-04 + SC3 = 0.030501 + W_density = 0.7267 + H_density = 0.2250 iter = 19 - elapsed = 6.7040e-03 - tTotal = 0.1421 - subIterW = 4 - subIterH = 0 - relError = 0.3329 + elapsed = 7.4590e-03 + tTotal = 0.2453 + subIterW = 0 + subIterH = 2 + relError = 0.3285 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 1.8267e-03 - SC2 = 9.2796e-04 - SC3 = 0.1138 - W_density = 0.7033 - H_density = 0.2100 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 1.4820e-03 + SC2 = 7.8844e-04 + SC3 = 0.1032 + W_density = 0.7267 + H_density = 0.2250 iter = 20 - elapsed = 5.3930e-03 - tTotal = 0.1475 - subIterW = 3 + elapsed = 7.4450e-03 + tTotal = 0.2527 + subIterW = 2 subIterH = 0 - relError = 0.3350 + relError = 0.3298 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 1.5770e-03 - SC2 = 7.9166e-04 - SC3 = 0.095795 - W_density = 0.6933 - H_density = 0.2100 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 1.5321e-03 + SC2 = 8.0589e-04 + SC3 = 0.1043 + W_density = 0.7233 + H_density = 0.2150 iter = 21 - elapsed = 5.3970e-03 - tTotal = 0.1529 - subIterW = 3 - subIterH = 0 - relError = 0.3370 + elapsed = 9.6520e-03 + tTotal = 0.2624 + subIterW = 0 + subIterH = 2 + relError = 0.3311 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 1.3752e-03 - SC2 = 6.8485e-04 - SC3 = 0.084116 - W_density = 0.6900 - H_density = 0.2050 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 1.4099e-03 + SC2 = 7.3313e-04 + SC3 = 0.095995 + W_density = 0.7233 + H_density = 0.2150 iter = 22 - elapsed = 6.5980e-03 - tTotal = 0.1595 - subIterW = 0 - subIterH = 0 - relError = 0.3387 + elapsed = 0.011720 + tTotal = 0.2741 + subIterW = 2 + subIterH = 4 + relError = 0.3321 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 1.2304e-03 - SC2 = 6.0782e-04 - SC3 = 0.074794 - W_density = 0.6867 - H_density = 0.2050 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 1.2325e-03 + SC2 = 6.3595e-04 + SC3 = 0.083210 + W_density = 0.7200 + H_density = 0.2100 iter = 23 - elapsed = 5.2150e-03 - tTotal = 0.1647 - subIterW = 0 + elapsed = 9.7520e-03 + tTotal = 0.2838 + subIterW = 3 subIterH = 0 - relError = 0.3403 + relError = 0.3328 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 1.1309e-03 - SC2 = 5.5640e-04 - SC3 = 0.066871 - W_density = 0.6833 - H_density = 0.2050 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 1.1183e-03 + SC2 = 5.7926e-04 + SC3 = 0.072782 + W_density = 0.7233 + H_density = 0.2100 iter = 24 - elapsed = 5.2600e-03 - tTotal = 0.1700 + elapsed = 7.5200e-03 + tTotal = 0.2913 subIterW = 0 subIterH = 0 - relError = 0.3417 + relError = 0.3335 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 1.0849e-03 - SC2 = 5.3161e-04 - SC3 = 0.059308 - W_density = 0.6800 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 1.0044e-03 + SC2 = 5.1829e-04 + SC3 = 0.016758 + W_density = 0.7233 H_density = 0.2050 iter = 25 - elapsed = 5.2390e-03 - tTotal = 0.1752 + elapsed = 7.4490e-03 + tTotal = 0.2988 subIterW = 0 subIterH = 0 - relError = 0.3432 + relError = 0.3342 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 1.0098e-03 - SC2 = 4.9685e-04 - SC3 = 0.053392 - W_density = 0.6867 - H_density = 0.2000 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 9.2677e-04 + SC2 = 4.7636e-04 + SC3 = 0.058460 + W_density = 0.7200 + H_density = 0.2050 iter = 26 - elapsed = 6.9270e-03 - tTotal = 0.1821 - subIterW = 2 + elapsed = 9.7930e-03 + tTotal = 0.3086 + subIterW = 0 subIterH = 0 - relError = 0.3445 + relError = 0.3349 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 9.6607e-04 - SC2 = 4.7338e-04 - SC3 = 0.011162 - W_density = 0.6867 - H_density = 0.1950 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 8.3539e-04 + SC2 = 4.2605e-04 + SC3 = 0.052049 + W_density = 0.7167 + H_density = 0.2000 iter = 27 - elapsed = 5.1360e-03 - tTotal = 0.1873 + elapsed = 0.013370 + tTotal = 0.3220 subIterW = 0 subIterH = 0 - relError = 0.3458 + relError = 0.3354 toshow = scalar structure containing the fields: - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 6.9990e-04 - SC2 = 3.4015e-04 - SC3 = 8.8157e-03 - W_density = 0.6867 - H_density = 0.1850 + initGrNormW = 355.75 + initGrNormH = 460.81 + initNorm = 14.367 + SC1 = 7.1389e-04 + SC2 = 3.6408e-04 + SC3 = 0.011189 + W_density = 0.7167 + H_density = 0.2000 iter = 28 - elapsed = 5.5540e-03 - tTotal = 0.1928 - subIterW = 0 - subIterH = 2 - relError = 0.3466 - -toshow = - - scalar structure containing the fields: - - initGrNormW = 348.10 - initGrNormH = 425.68 - initNorm = 14.240 - SC1 = 4.9655e-04 - SC2 = 2.4132e-04 - SC3 = 7.4628e-03 - W_density = 0.6867 - H_density = 0.1850 - iter = 29 - elapsed = 3.9020e-03 - tTotal = 0.1967 + elapsed = 7.4190e-03 + tTotal = 0.3294 subIterW = 0 subIterH = 0 - relError = 0.3472 + relError = 0.3359 Stop: tolerance reached. final = scalar structure containing the fields: - iterations = 29 - elapsed_total = 0.1967 - relative_error = 0.3472 - W_density = 0.6867 - H_density = 0.1850 + iterations = 28 + elapsed_total = 0.3294 + relative_error = 0.3359 + W_density = 0.7167 + H_density = 0.2000 ***** test [W,H,iter,HIS]=nmf_bpas(A,k,'verbose','type','sparse','nnlssolver','bp','alpha',1.1,'beta',1.3); - initGrNormW = 402.54 - initGrNormH = 572.48 - initNorm = 14.240 - SC1 = 1.3997 - SC2 = 699.83 - SC3 = 0.5099 + initGrNormW = 395.57 + initGrNormH = 552.16 + initNorm = 14.367 + SC1 = 1.3585 + SC2 = 679.23 + SC3 = 0.5051 W_density = 1 H_density = 1 par = @@ -4485,438 +4510,438 @@ scalar structure containing the fields: - initGrNormW = 402.54 - initGrNormH = 572.48 - initNorm = 14.240 - SC1 = 0.026545 - SC2 = 0.023944 - SC3 = 2.2406 - W_density = 0.9933 - H_density = 0.5650 + initGrNormW = 395.57 + initGrNormH = 552.16 + initNorm = 14.367 + SC1 = 0.025094 + SC2 = 0.023388 + SC3 = 1.9995 + W_density = 0.9867 + H_density = 0.5350 iter = 1 - elapsed = 0.010719 - tTotal = 0.010719 + elapsed = 0.023016 + tTotal = 0.023016 subIterW = 2 - subIterH = 35 - relError = 0.5488 + subIterH = 36 + relError = 0.5391 toshow = scalar structure containing the fields: - initGrNormW = 402.54 - initGrNormH = 572.48 - initNorm = 14.240 - SC1 = 0.013306 - SC2 = 5.5088e-03 - SC3 = 0.7689 - W_density = 0.4900 - H_density = 0.2200 + initGrNormW = 395.57 + initGrNormH = 552.16 + initNorm = 14.367 + SC1 = 0.015213 + SC2 = 5.9028e-03 + SC3 = 0.8761 + W_density = 0.4700 + H_density = 0.2050 iter = 2 - elapsed = 0.014318 - tTotal = 0.025037 - subIterW = 2 - subIterH = 36 - relError = 0.4329 + elapsed = 0.023727 + tTotal = 0.046743 + subIterW = 3 + subIterH = 29 + relError = 0.4401 toshow = scalar structure containing the fields: - initGrNormW = 402.54 - initGrNormH = 572.48 - initNorm = 14.240 - SC1 = 0.012472 - SC2 = 4.5647e-03 - SC3 = 0.6509 - W_density = 0.4567 + initGrNormW = 395.57 + initGrNormH = 552.16 + initNorm = 14.367 + SC1 = 0.016665 + SC2 = 5.2329e-03 + SC3 = 0.7051 + W_density = 0.3767 H_density = 0.2100 iter = 3 - elapsed = 7.2700e-03 - tTotal = 0.032307 + elapsed = 0.015724 + tTotal = 0.062467 subIterW = 4 - subIterH = 4 - relError = 0.4304 + subIterH = 6 + relError = 0.4406 toshow = scalar structure containing the fields: - initGrNormW = 402.54 - initGrNormH = 572.48 - initNorm = 14.240 - SC1 = 0.013113 - SC2 = 3.8291e-03 - SC3 = 0.6139 - W_density = 0.3600 - H_density = 0.1800 + initGrNormW = 395.57 + initGrNormH = 552.16 + initNorm = 14.367 + SC1 = 0.017387 + SC2 = 4.6250e-03 + SC3 = 0.6093 + W_density = 0.3367 + H_density = 0.1550 iter = 4 - elapsed = 5.7010e-03 - tTotal = 0.038008 - subIterW = 3 + elapsed = 0.015043 + tTotal = 0.077510 + subIterW = 5 subIterH = 3 - relError = 0.4378 + relError = 0.4503 toshow = scalar structure containing the fields: - initGrNormW = 402.54 - initGrNormH = 572.48 - initNorm = 14.240 - SC1 = 0.014473 - SC2 = 3.8787e-03 - SC3 = 0.7046 - W_density = 0.3400 - H_density = 0.1600 + initGrNormW = 395.57 + initGrNormH = 552.16 + initNorm = 14.367 + SC1 = 0.015495 + SC2 = 3.1609e-03 + SC3 = 0.4769 + W_density = 0.2533 + H_density = 0.1300 iter = 5 - elapsed = 5.1310e-03 - tTotal = 0.043139 - subIterW = 2 - subIterH = 2 - relError = 0.4482 + elapsed = 0.015304 + tTotal = 0.092814 + subIterW = 4 + subIterH = 5 + relError = 0.4566 toshow = scalar structure containing the fields: - initGrNormW = 402.54 - initGrNormH = 572.48 - initNorm = 14.240 - SC1 = 0.015758 - SC2 = 3.2147e-03 - SC3 = 0.5435 - W_density = 0.2500 - H_density = 0.1350 + initGrNormW = 395.57 + initGrNormH = 552.16 + initNorm = 14.367 + SC1 = 0.012068 + SC2 = 2.2446e-03 + SC3 = 0.3369 + W_density = 0.2300 + H_density = 0.1200 iter = 6 - elapsed = 5.1000e-03 - tTotal = 0.048239 - subIterW = 4 + elapsed = 0.010187 + tTotal = 0.1030 + subIterW = 2 subIterH = 1 - relError = 0.4589 + relError = 0.4639 toshow = scalar structure containing the fields: - initGrNormW = 402.54 - initGrNormH = 572.48 - initNorm = 14.240 - SC1 = 0.011639 - SC2 = 2.0949e-03 - SC3 = 0.2963 - W_density = 0.2200 - H_density = 0.1200 + initGrNormW = 395.57 + initGrNormH = 552.16 + initNorm = 14.367 + SC1 = 9.1072e-03 + SC2 = 1.2568e-03 + SC3 = 0.1166 + W_density = 0.1600 + H_density = 0.1050 iter = 7 - elapsed = 5.0620e-03 - tTotal = 0.053301 - subIterW = 3 + elapsed = 0.010950 + tTotal = 0.1140 + subIterW = 2 subIterH = 1 - relError = 0.4656 + relError = 0.4689 toshow = scalar structure containing the fields: - initGrNormW = 402.54 - initGrNormH = 572.48 - initNorm = 14.240 - SC1 = 9.7178e-03 - SC2 = 1.3216e-03 - SC3 = 0.093522 - W_density = 0.1533 - H_density = 0.1100 + initGrNormW = 395.57 + initGrNormH = 552.16 + initNorm = 14.367 + SC1 = 1.3875e-03 + SC2 = 1.3875e-04 + SC3 = 0.026919 + W_density = 0.1000 + H_density = 0.1000 iter = 8 - elapsed = 4.9550e-03 - tTotal = 0.058256 - subIterW = 2 - subIterH = 1 - relError = 0.4704 + elapsed = 5.9770e-03 + tTotal = 0.1199 + subIterW = 0 + subIterH = 0 + relError = 0.4698 toshow = scalar structure containing the fields: - initGrNormW = 402.54 - initGrNormH = 572.48 - initNorm = 14.240 - SC1 = 8.2596e-04 - SC2 = 8.2596e-05 - SC3 = 0.016657 + initGrNormW = 395.57 + initGrNormH = 552.16 + initNorm = 14.367 + SC1 = 8.3518e-04 + SC2 = 8.3518e-05 + SC3 = 0.018470 W_density = 0.1000 H_density = 0.1000 iter = 9 - elapsed = 3.1770e-03 - tTotal = 0.061433 + elapsed = 2.7080e-03 + tTotal = 0.1226 subIterW = 0 - subIterH = 1 - relError = 0.4711 + subIterH = 0 + relError = 0.4697 toshow = scalar structure containing the fields: - initGrNormW = 402.54 - initGrNormH = 572.48 - initNorm = 14.240 - SC1 = 4.8925e-04 - SC2 = 4.8925e-05 - SC3 = 0.010650 + initGrNormW = 395.57 + initGrNormH = 552.16 + initNorm = 14.367 + SC1 = 5.5030e-04 + SC2 = 5.5030e-05 + SC3 = 0.011743 W_density = 0.1000 H_density = 0.1000 iter = 10 - elapsed = 1.5180e-03 - tTotal = 0.062951 + elapsed = 2.9580e-03 + tTotal = 0.1256 subIterW = 0 subIterH = 0 - relError = 0.4710 + relError = 0.4697 toshow = scalar structure containing the fields: - initGrNormW = 402.54 - initGrNormH = 572.48 - initNorm = 14.240 - SC1 = 3.2110e-04 - SC2 = 3.2110e-05 - SC3 = 9.3103e-03 + initGrNormW = 395.57 + initGrNormH = 552.16 + initNorm = 14.367 + SC1 = 3.6283e-04 + SC2 = 3.6283e-05 + SC3 = 7.7424e-03 W_density = 0.1000 H_density = 0.1000 iter = 11 - elapsed = 1.4910e-03 - tTotal = 0.064442 + elapsed = 2.9820e-03 + tTotal = 0.1286 subIterW = 0 subIterH = 0 - relError = 0.4710 + relError = 0.4696 toshow = scalar structure containing the fields: - initGrNormW = 402.54 - initGrNormH = 572.48 - initNorm = 14.240 - SC1 = 2.1083e-04 - SC2 = 2.1083e-05 - SC3 = 4.9372e-03 + initGrNormW = 395.57 + initGrNormH = 552.16 + initNorm = 14.367 + SC1 = 2.3928e-04 + SC2 = 2.3928e-05 + SC3 = 6.2161e-03 W_density = 0.1000 H_density = 0.1000 iter = 12 - elapsed = 1.4990e-03 - tTotal = 0.065941 + elapsed = 3.2270e-03 + tTotal = 0.1318 subIterW = 0 subIterH = 0 - relError = 0.4710 + relError = 0.4696 toshow = scalar structure containing the fields: - initGrNormW = 402.54 - initGrNormH = 572.48 - initNorm = 14.240 - SC1 = 1.3845e-04 - SC2 = 1.3845e-05 - SC3 = 3.1221e-03 + initGrNormW = 395.57 + initGrNormH = 552.16 + initNorm = 14.367 + SC1 = 1.5783e-04 + SC2 = 1.5783e-05 + SC3 = 4.1001e-03 W_density = 0.1000 H_density = 0.1000 iter = 13 - elapsed = 1.5250e-03 - tTotal = 0.067466 + elapsed = 2.7690e-03 + tTotal = 0.1346 subIterW = 0 subIterH = 0 - relError = 0.4710 + relError = 0.4696 toshow = scalar structure containing the fields: - initGrNormW = 402.54 - initGrNormH = 572.48 - initNorm = 14.240 - SC1 = 9.0927e-05 - SC2 = 9.0927e-06 - SC3 = 1.7300e-03 + initGrNormW = 395.57 + initGrNormH = 552.16 + initNorm = 14.367 + SC1 = 1.0411e-04 + SC2 = 1.0411e-05 + SC3 = 2.2217e-03 W_density = 0.1000 H_density = 0.1000 iter = 14 - elapsed = 1.5440e-03 - tTotal = 0.069010 + elapsed = 2.2290e-03 + tTotal = 0.1368 subIterW = 0 subIterH = 0 - relError = 0.4710 + relError = 0.4696 toshow = scalar structure containing the fields: - initGrNormW = 402.54 - initGrNormH = 572.48 - initNorm = 14.240 - SC1 = 5.9720e-05 - SC2 = 5.9720e-06 - SC3 = 1.4544e-03 + initGrNormW = 395.57 + initGrNormH = 552.16 + initNorm = 14.367 + SC1 = 6.8683e-05 + SC2 = 6.8683e-06 + SC3 = 1.4656e-03 W_density = 0.1000 H_density = 0.1000 iter = 15 - elapsed = 1.5160e-03 - tTotal = 0.070526 + elapsed = 2.7280e-03 + tTotal = 0.1395 subIterW = 0 subIterH = 0 - relError = 0.4710 + relError = 0.4696 toshow = scalar structure containing the fields: - initGrNormW = 402.54 - initGrNormH = 572.48 - initNorm = 14.240 - SC1 = 3.9225e-05 - SC2 = 3.9225e-06 - SC3 = 8.5293e-04 + initGrNormW = 395.57 + initGrNormH = 552.16 + initNorm = 14.367 + SC1 = 4.5312e-05 + SC2 = 4.5312e-06 + SC3 = 9.6693e-04 W_density = 0.1000 H_density = 0.1000 iter = 16 - elapsed = 1.4900e-03 - tTotal = 0.072016 + elapsed = 2.8320e-03 + tTotal = 0.1424 subIterW = 0 subIterH = 0 - relError = 0.4710 + relError = 0.4696 toshow = scalar structure containing the fields: - initGrNormW = 402.54 - initGrNormH = 572.48 - initNorm = 14.240 - SC1 = 2.5764e-05 - SC2 = 2.5764e-06 - SC3 = 6.0333e-04 + initGrNormW = 395.57 + initGrNormH = 552.16 + initNorm = 14.367 + SC1 = 2.9895e-05 + SC2 = 2.9895e-06 + SC3 = 5.9540e-04 W_density = 0.1000 H_density = 0.1000 iter = 17 - elapsed = 1.4970e-03 - tTotal = 0.073513 + elapsed = 2.7120e-03 + tTotal = 0.1451 subIterW = 0 subIterH = 0 - relError = 0.4710 + relError = 0.4696 toshow = scalar structure containing the fields: - initGrNormW = 402.54 - initGrNormH = 572.48 - initNorm = 14.240 - SC1 = 1.6923e-05 - SC2 = 1.6923e-06 - SC3 = 4.2932e-04 + initGrNormW = 395.57 + initGrNormH = 552.16 + initNorm = 14.367 + SC1 = 1.9724e-05 + SC2 = 1.9724e-06 + SC3 = 4.7139e-04 W_density = 0.1000 H_density = 0.1000 iter = 18 - elapsed = 1.5210e-03 - tTotal = 0.075034 + elapsed = 3.0090e-03 + tTotal = 0.1481 subIterW = 0 subIterH = 0 - relError = 0.4710 + relError = 0.4696 toshow = scalar structure containing the fields: - initGrNormW = 402.54 - initGrNormH = 572.48 - initNorm = 14.240 - SC1 = 1.1116e-05 - SC2 = 1.1116e-06 - SC3 = 2.3338e-04 + initGrNormW = 395.57 + initGrNormH = 552.16 + initNorm = 14.367 + SC1 = 1.3013e-05 + SC2 = 1.3013e-06 + SC3 = 2.7769e-04 W_density = 0.1000 H_density = 0.1000 iter = 19 - elapsed = 1.5350e-03 - tTotal = 0.076569 + elapsed = 3.1910e-03 + tTotal = 0.1513 subIterW = 0 subIterH = 0 - relError = 0.4710 + relError = 0.4696 toshow = scalar structure containing the fields: - initGrNormW = 402.54 - initGrNormH = 572.48 - initNorm = 14.240 - SC1 = 7.3015e-06 - SC2 = 7.3015e-07 - SC3 = 1.6465e-04 + initGrNormW = 395.57 + initGrNormH = 552.16 + initNorm = 14.367 + SC1 = 8.5857e-06 + SC2 = 8.5857e-07 + SC3 = 2.0520e-04 W_density = 0.1000 H_density = 0.1000 iter = 20 - elapsed = 1.5100e-03 - tTotal = 0.078079 + elapsed = 2.7420e-03 + tTotal = 0.1540 subIterW = 0 subIterH = 0 - relError = 0.4710 + relError = 0.4696 toshow = scalar structure containing the fields: - initGrNormW = 402.54 - initGrNormH = 572.48 - initNorm = 14.240 - SC1 = 4.7961e-06 - SC2 = 4.7961e-07 - SC3 = 1.2167e-04 + initGrNormW = 395.57 + initGrNormH = 552.16 + initNorm = 14.367 + SC1 = 5.6647e-06 + SC2 = 5.6647e-07 + SC3 = 1.3539e-04 W_density = 0.1000 H_density = 0.1000 iter = 21 - elapsed = 1.4950e-03 - tTotal = 0.079574 + elapsed = 2.7730e-03 + tTotal = 0.1568 subIterW = 0 subIterH = 0 - relError = 0.4710 + relError = 0.4696 toshow = scalar structure containing the fields: - initGrNormW = 402.54 - initGrNormH = 572.48 - initNorm = 14.240 - SC1 = 3.1504e-06 - SC2 = 3.1504e-07 - SC3 = 6.8504e-05 + initGrNormW = 395.57 + initGrNormH = 552.16 + initNorm = 14.367 + SC1 = 3.7375e-06 + SC2 = 3.7375e-07 + SC3 = 8.2709e-05 W_density = 0.1000 H_density = 0.1000 iter = 22 - elapsed = 1.7040e-03 - tTotal = 0.081278 + elapsed = 3.1900e-03 + tTotal = 0.1600 subIterW = 0 subIterH = 0 - relError = 0.4710 + relError = 0.4696 toshow = scalar structure containing the fields: - initGrNormW = 402.54 - initGrNormH = 572.48 - initNorm = 14.240 - SC1 = 2.0694e-06 - SC2 = 2.0694e-07 - SC3 = 4.8459e-05 + initGrNormW = 395.57 + initGrNormH = 552.16 + initNorm = 14.367 + SC1 = 2.4659e-06 + SC2 = 2.4659e-07 + SC3 = 5.6669e-05 W_density = 0.1000 H_density = 0.1000 iter = 23 - elapsed = 1.7040e-03 - tTotal = 0.082982 + elapsed = 2.8920e-03 + tTotal = 0.1629 subIterW = 0 subIterH = 0 - relError = 0.4710 + relError = 0.4696 Stop: tolerance reached. final = @@ -4924,19 +4949,19 @@ scalar structure containing the fields: iterations = 23 - elapsed_total = 0.082982 - relative_error = 0.4710 + elapsed_total = 0.1629 + relative_error = 0.4696 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 = 361.80 - initGrNormH = 474.39 - initNorm = 14.240 - SC1 = 1.1932 - SC2 = 596.61 - SC3 = 0.5035 + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 1.0260 + SC2 = 512.98 + SC3 = 0.4831 W_density = 1 H_density = 1 par = @@ -4947,8 +4972,8 @@ n = 20 type = plain nnls_solver = bp - alpha = 0.5074 - beta = 0.5074 + alpha = 0.5111 + beta = 0.5111 max_iter = 100 min_iter = 20 max_time = 1000 @@ -4960,467 +4985,562 @@ scalar structure containing the fields: - initGrNormW = 361.80 - initGrNormH = 474.39 - initNorm = 14.240 - SC1 = 0.019259 - SC2 = 0.014329 - SC3 = 0.7037 - W_density = 0.7533 - H_density = 0.5750 + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 0.015980 + SC2 = 0.011857 + SC3 = 0.5315 + W_density = 0.7567 + H_density = 0.5700 iter = 1 - elapsed = 0.017206 - tTotal = 0.017206 + elapsed = 0.033259 + tTotal = 0.033259 subIterW = 38 - subIterH = 37 - relError = 0.3198 + subIterH = 30 + relError = 0.3029 toshow = scalar structure containing the fields: - initGrNormW = 361.80 - initGrNormH = 474.39 - initNorm = 14.240 - SC1 = 0.010521 - SC2 = 7.4281e-03 - SC3 = 0.3829 - W_density = 0.7400 - H_density = 0.6200 + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 7.7801e-03 + SC2 = 5.4150e-03 + SC3 = 0.3016 + W_density = 0.7267 + H_density = 0.6050 iter = 2 - elapsed = 0.015785 - tTotal = 0.032991 - subIterW = 14 - subIterH = 29 - relError = 0.2680 + elapsed = 0.033224 + tTotal = 0.066483 + subIterW = 21 + subIterH = 22 + relError = 0.2596 toshow = scalar structure containing the fields: - initGrNormW = 361.80 - initGrNormH = 474.39 - initNorm = 14.240 - SC1 = 6.3730e-03 - SC2 = 4.4484e-03 - SC3 = 0.2897 - W_density = 0.7400 - H_density = 0.5950 + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 4.6233e-03 + SC2 = 3.2270e-03 + SC3 = 0.1880 + W_density = 0.7333 + H_density = 0.6150 iter = 3 - elapsed = 0.012263 - tTotal = 0.045254 - subIterW = 10 - subIterH = 13 - relError = 0.2518 + elapsed = 0.023083 + tTotal = 0.089566 + subIterW = 14 + subIterH = 10 + relError = 0.2500 toshow = scalar structure containing the fields: - initGrNormW = 361.80 - initGrNormH = 474.39 - initNorm = 14.240 - SC1 = 4.6779e-03 - SC2 = 3.2464e-03 - SC3 = 0.2158 - W_density = 0.7333 - H_density = 0.6100 + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 3.7327e-03 + SC2 = 2.6278e-03 + SC3 = 0.1500 + W_density = 0.7500 + H_density = 0.6050 iter = 4 - elapsed = 9.7760e-03 - tTotal = 0.055030 - subIterW = 8 - subIterH = 9 - relError = 0.2434 + elapsed = 0.020778 + tTotal = 0.1103 + subIterW = 5 + subIterH = 12 + relError = 0.2451 toshow = scalar structure containing the fields: - initGrNormW = 361.80 - initGrNormH = 474.39 - initNorm = 14.240 - SC1 = 3.8260e-03 - SC2 = 2.6400e-03 - SC3 = 0.1793 - W_density = 0.7300 - H_density = 0.6200 + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 3.3646e-03 + SC2 = 2.3081e-03 + SC3 = 0.1342 + W_density = 0.7333 + H_density = 0.5950 iter = 5 - elapsed = 8.3750e-03 - tTotal = 0.063405 - subIterW = 9 + elapsed = 0.016982 + tTotal = 0.1273 + subIterW = 6 subIterH = 6 - relError = 0.2378 + relError = 0.2409 toshow = scalar structure containing the fields: - initGrNormW = 361.80 - initGrNormH = 474.39 - initNorm = 14.240 - SC1 = 3.2258e-03 - SC2 = 2.2065e-03 - SC3 = 0.1505 - W_density = 0.7267 - H_density = 0.6150 + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 2.4109e-03 + SC2 = 1.6394e-03 + SC3 = 0.1111 + W_density = 0.7333 + H_density = 0.5900 iter = 6 - elapsed = 8.1920e-03 - tTotal = 0.071597 - subIterW = 3 - subIterH = 6 - relError = 0.2341 + elapsed = 0.020441 + tTotal = 0.1478 + subIterW = 6 + subIterH = 10 + relError = 0.2380 toshow = scalar structure containing the fields: - initGrNormW = 361.80 - initGrNormH = 474.39 - initNorm = 14.240 - SC1 = 2.8234e-03 - SC2 = 1.9425e-03 - SC3 = 0.1235 - W_density = 0.7267 - H_density = 0.6200 + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 1.9355e-03 + SC2 = 1.3239e-03 + SC3 = 0.092084 + W_density = 0.7367 + H_density = 0.5950 iter = 7 - elapsed = 8.1110e-03 - tTotal = 0.079708 - subIterW = 3 - subIterH = 4 - relError = 0.2317 + elapsed = 0.016992 + tTotal = 0.1648 + subIterW = 4 + subIterH = 5 + relError = 0.2366 toshow = scalar structure containing the fields: - initGrNormW = 361.80 - initGrNormH = 474.39 - initNorm = 14.240 - SC1 = 2.6023e-03 - SC2 = 1.8060e-03 - SC3 = 0.1205 - W_density = 0.7333 - H_density = 0.6150 + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 1.8206e-03 + SC2 = 1.2489e-03 + SC3 = 0.084545 + W_density = 0.7367 + H_density = 0.5950 iter = 8 - elapsed = 9.2980e-03 - tTotal = 0.089006 - subIterW = 5 - subIterH = 3 - relError = 0.2299 + elapsed = 0.016167 + tTotal = 0.1809 + subIterW = 7 + subIterH = 1 + relError = 0.2355 toshow = scalar structure containing the fields: - initGrNormW = 361.80 - initGrNormH = 474.39 - initNorm = 14.240 - SC1 = 2.6260e-03 - SC2 = 1.8277e-03 - SC3 = 0.1232 - W_density = 0.7333 - H_density = 0.6400 + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 1.7527e-03 + SC2 = 1.2164e-03 + SC3 = 0.081333 + W_density = 0.7433 + H_density = 0.6100 iter = 9 - elapsed = 7.8450e-03 - tTotal = 0.096851 - subIterW = 2 - subIterH = 8 - relError = 0.2282 + elapsed = 0.016808 + tTotal = 0.1977 + subIterW = 4 + subIterH = 7 + relError = 0.2345 toshow = scalar structure containing the fields: - initGrNormW = 361.80 - initGrNormH = 474.39 - initNorm = 14.240 - SC1 = 2.7609e-03 - SC2 = 1.9271e-03 - SC3 = 0.1304 - W_density = 0.7400 - H_density = 0.6250 + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 1.7704e-03 + SC2 = 1.2251e-03 + SC3 = 0.083148 + W_density = 0.7433 + H_density = 0.6100 iter = 10 - elapsed = 7.4920e-03 - tTotal = 0.1043 - subIterW = 2 - subIterH = 3 - relError = 0.2265 + elapsed = 0.016231 + tTotal = 0.2140 + subIterW = 4 + subIterH = 4 + relError = 0.2336 toshow = scalar structure containing the fields: - initGrNormW = 361.80 - initGrNormH = 474.39 - initNorm = 14.240 - SC1 = 2.7539e-03 - SC2 = 1.9112e-03 - SC3 = 0.1322 - W_density = 0.7333 - H_density = 0.6250 + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 1.7437e-03 + SC2 = 1.2101e-03 + SC3 = 0.080936 + W_density = 0.7500 + H_density = 0.6050 iter = 11 - elapsed = 7.5450e-03 - tTotal = 0.1119 - subIterW = 2 - subIterH = 4 - relError = 0.2248 + elapsed = 0.015415 + tTotal = 0.2294 + subIterW = 3 + subIterH = 3 + relError = 0.2327 toshow = scalar structure containing the fields: - initGrNormW = 361.80 - initGrNormH = 474.39 - initNorm = 14.240 - SC1 = 2.6509e-03 - SC2 = 1.8185e-03 - SC3 = 0.1244 - W_density = 0.7233 - H_density = 0.6250 + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 1.7528e-03 + SC2 = 1.2199e-03 + SC3 = 0.081888 + W_density = 0.7500 + H_density = 0.6050 iter = 12 - elapsed = 7.8040e-03 - tTotal = 0.1197 - subIterW = 8 + elapsed = 0.015218 + tTotal = 0.2446 + subIterW = 2 subIterH = 2 - relError = 0.2233 + relError = 0.2319 toshow = scalar structure containing the fields: - initGrNormW = 361.80 - initGrNormH = 474.39 - initNorm = 14.240 - SC1 = 2.4696e-03 - SC2 = 1.6991e-03 - SC3 = 0.1067 - W_density = 0.7300 - H_density = 0.6250 + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 1.7747e-03 + SC2 = 1.2245e-03 + SC3 = 0.081306 + W_density = 0.7367 + H_density = 0.6150 iter = 13 - elapsed = 7.6380e-03 - tTotal = 0.1273 - subIterW = 4 - subIterH = 3 - relError = 0.2221 + elapsed = 0.016012 + tTotal = 0.2606 + subIterW = 6 + subIterH = 2 + relError = 0.2310 toshow = scalar structure containing the fields: - initGrNormW = 361.80 - initGrNormH = 474.39 - initNorm = 14.240 - SC1 = 2.3963e-03 - SC2 = 1.6438e-03 - SC3 = 0.095627 + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 1.7922e-03 + SC2 = 1.2187e-03 + SC3 = 0.083866 W_density = 0.7300 - H_density = 0.6150 + H_density = 0.6050 iter = 14 - elapsed = 7.5810e-03 - tTotal = 0.1349 - subIterW = 4 + elapsed = 0.015439 + tTotal = 0.2760 + subIterW = 2 subIterH = 2 - relError = 0.2211 + relError = 0.2302 toshow = scalar structure containing the fields: - initGrNormW = 361.80 - initGrNormH = 474.39 - initNorm = 14.240 - SC1 = 2.3884e-03 - SC2 = 1.6241e-03 - SC3 = 0.089570 - W_density = 0.7233 - H_density = 0.6100 + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 1.6951e-03 + SC2 = 1.1357e-03 + SC3 = 0.082165 + W_density = 0.7200 + H_density = 0.5900 iter = 15 - elapsed = 8.7760e-03 - tTotal = 0.1437 - subIterW = 6 - subIterH = 0 - relError = 0.2204 + elapsed = 0.015849 + tTotal = 0.2919 + subIterW = 3 + subIterH = 3 + relError = 0.2295 toshow = scalar structure containing the fields: - initGrNormW = 361.80 - initGrNormH = 474.39 - initNorm = 14.240 - SC1 = 2.2775e-03 - SC2 = 1.5396e-03 - SC3 = 0.079714 - W_density = 0.7200 - H_density = 0.6000 + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 1.6291e-03 + SC2 = 1.0948e-03 + SC3 = 0.078362 + W_density = 0.7167 + H_density = 0.5950 iter = 16 - elapsed = 7.8060e-03 - tTotal = 0.1515 - subIterW = 5 - subIterH = 3 - relError = 0.2198 + elapsed = 0.014174 + tTotal = 0.3061 + subIterW = 3 + subIterH = 0 + relError = 0.2288 toshow = scalar structure containing the fields: - initGrNormW = 361.80 - initGrNormH = 474.39 - initNorm = 14.240 - SC1 = 1.8482e-03 - SC2 = 1.2346e-03 - SC3 = 0.070878 + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 1.4928e-03 + SC2 = 9.9421e-04 + SC3 = 0.073049 W_density = 0.7133 - H_density = 0.6000 + H_density = 0.5950 iter = 17 - elapsed = 7.5680e-03 - tTotal = 0.1591 - subIterW = 4 - subIterH = 2 - relError = 0.2193 + elapsed = 0.015562 + tTotal = 0.3216 + subIterW = 0 + subIterH = 4 + relError = 0.2283 toshow = scalar structure containing the fields: - initGrNormW = 361.80 - initGrNormH = 474.39 - initNorm = 14.240 - SC1 = 1.6880e-03 - SC2 = 1.1242e-03 - SC3 = 0.062962 - W_density = 0.7100 - H_density = 0.6000 + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 1.3920e-03 + SC2 = 9.2709e-04 + SC3 = 0.068342 + W_density = 0.7133 + H_density = 0.5900 iter = 18 - elapsed = 5.6940e-03 - tTotal = 0.1648 - subIterW = 0 + elapsed = 0.015548 + tTotal = 0.3372 + subIterW = 4 subIterH = 0 - relError = 0.2188 + relError = 0.2278 toshow = scalar structure containing the fields: - initGrNormW = 361.80 - initGrNormH = 474.39 - initNorm = 14.240 - SC1 = 1.5752e-03 - SC2 = 1.0428e-03 - SC3 = 0.057417 - W_density = 0.7067 - H_density = 0.5900 + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 1.3205e-03 + SC2 = 8.8207e-04 + SC3 = 0.063885 + W_density = 0.7167 + H_density = 0.5950 iter = 19 - elapsed = 9.1290e-03 - tTotal = 0.1739 + elapsed = 0.015190 + tTotal = 0.3524 subIterW = 3 - subIterH = 2 - relError = 0.2185 + subIterH = 0 + relError = 0.2273 toshow = scalar structure containing the fields: - initGrNormW = 361.80 - initGrNormH = 474.39 - initNorm = 14.240 - SC1 = 1.4523e-03 - SC2 = 9.6722e-04 - SC3 = 0.051273 - W_density = 0.7100 + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 1.2935e-03 + SC2 = 8.6664e-04 + SC3 = 0.060929 + W_density = 0.7167 H_density = 0.5950 iter = 20 - elapsed = 8.9990e-03 - tTotal = 0.1829 - subIterW = 4 + elapsed = 0.015456 + tTotal = 0.3678 + subIterW = 7 subIterH = 0 - relError = 0.2182 + relError = 0.2269 toshow = scalar structure containing the fields: - initGrNormW = 361.80 - initGrNormH = 474.39 - initNorm = 14.240 - SC1 = 1.1158e-03 - SC2 = 7.3644e-04 - SC3 = 0.047629 - W_density = 0.7067 - H_density = 0.5850 + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 1.3022e-03 + SC2 = 8.8026e-04 + SC3 = 0.060800 + W_density = 0.7300 + H_density = 0.5950 iter = 21 - elapsed = 7.4620e-03 - tTotal = 0.1903 - subIterW = 0 + elapsed = 0.014557 + tTotal = 0.3824 + subIterW = 4 subIterH = 2 - relError = 0.2179 + relError = 0.2266 toshow = scalar structure containing the fields: - initGrNormW = 361.80 - initGrNormH = 474.39 - initNorm = 14.240 - SC1 = 8.3485e-04 - SC2 = 5.4098e-04 - SC3 = 0.043267 - W_density = 0.6967 - H_density = 0.5750 + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 1.3041e-03 + SC2 = 8.8418e-04 + SC3 = 0.060727 + W_density = 0.7400 + H_density = 0.5850 iter = 22 - elapsed = 7.9540e-03 - tTotal = 0.1983 + elapsed = 0.015429 + tTotal = 0.3978 + subIterW = 4 + subIterH = 2 + relError = 0.2262 + +toshow = + + scalar structure containing the fields: + + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 1.2118e-03 + SC2 = 8.1190e-04 + SC3 = 0.056818 + W_density = 0.7300 + H_density = 0.5800 + iter = 23 + elapsed = 0.014793 + tTotal = 0.4126 subIterW = 3 - subIterH = 4 - relError = 0.2178 + subIterH = 3 + relError = 0.2258 toshow = scalar structure containing the fields: - initGrNormW = 361.80 - initGrNormH = 474.39 - initNorm = 14.240 - SC1 = 7.2473e-04 - SC2 = 4.6673e-04 - SC3 = 0.039414 - W_density = 0.6933 + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 1.1141e-03 + SC2 = 7.4866e-04 + SC3 = 0.052450 + W_density = 0.7367 H_density = 0.5700 - iter = 23 - elapsed = 7.1350e-03 - tTotal = 0.2054 + iter = 24 + elapsed = 0.018357 + tTotal = 0.4310 + subIterW = 2 + subIterH = 2 + relError = 0.2255 + +toshow = + + scalar structure containing the fields: + + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 1.0654e-03 + SC2 = 7.1805e-04 + SC3 = 0.048874 + W_density = 0.7400 + H_density = 0.5750 + iter = 25 + elapsed = 0.015974 + tTotal = 0.4469 subIterW = 0 subIterH = 0 - relError = 0.2176 + relError = 0.2252 toshow = scalar structure containing the fields: - initGrNormW = 361.80 - initGrNormH = 474.39 - initNorm = 14.240 - SC1 = 6.8528e-04 - SC2 = 4.4132e-04 - SC3 = 0.019097 - W_density = 0.6933 + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 1.0538e-03 + SC2 = 7.0816e-04 + SC3 = 0.046180 + W_density = 0.7367 + H_density = 0.5750 + iter = 26 + elapsed = 0.011109 + tTotal = 0.4580 + subIterW = 0 + subIterH = 0 + relError = 0.2249 + +toshow = + + scalar structure containing the fields: + + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 9.9101e-04 + SC2 = 6.5803e-04 + SC3 = 0.044801 + W_density = 0.7300 + H_density = 0.5650 + iter = 27 + elapsed = 0.014896 + tTotal = 0.4729 + subIterW = 2 + subIterH = 2 + relError = 0.2247 + +toshow = + + scalar structure containing the fields: + + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 9.4053e-04 + SC2 = 6.2639e-04 + SC3 = 0.021007 + W_density = 0.7300 H_density = 0.5700 - iter = 24 - elapsed = 4.6480e-03 - tTotal = 0.2101 + iter = 28 + elapsed = 0.011668 + tTotal = 0.4846 subIterW = 0 subIterH = 0 - relError = 0.2175 + relError = 0.2245 + +toshow = + + scalar structure containing the fields: + + initGrNormW = 326.84 + initGrNormH = 395.37 + initNorm = 14.367 + SC1 = 8.9858e-04 + SC2 = 5.9846e-04 + SC3 = 0.020735 + W_density = 0.7300 + H_density = 0.5700 + iter = 29 + elapsed = 8.5780e-03 + tTotal = 0.4932 + subIterW = 0 + subIterH = 0 + relError = 0.2242 Stop: tolerance reached. final = scalar structure containing the fields: - iterations = 24 - elapsed_total = 0.2101 - relative_error = 0.2175 - W_density = 0.6933 + iterations = 29 + elapsed_total = 0.4932 + relative_error = 0.2242 + W_density = 0.7300 H_density = 0.5700 ***** demo @@ -5469,6 +5589,24 @@ 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_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); +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/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 ***** test @@ -5484,39 +5622,6 @@ [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/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/@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)); -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 -***** shared v,x - v = [1 2 3]'; x = [2 5 6]'; -***** assert (circulant_matrix_vector_product(v, x), circulant_make_matrix(v)*x, eps); -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/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/funm.m] >>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/funm.m ***** function b = fsin (a) @@ -5526,6 +5631,13 @@ mtx = randn (100); assert (funm (mtx, "fsin"), thfm (mtx, "sin"), 1e-9) 1 test, 1 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)); +1 test, 1 passed, 0 known failure, 0 skipped Checking C++ files ... Done running the unit tests. Summary: 18 tests, 18 passed, 0 known failures, 0 skipped @@ -5556,12 +5668,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/2992669/tmp/hooks/B01_cleanup starting +I: user script /srv/workspace/pbuilder/2992669/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/956707 and its subdirectories -I: Current time: Tue Mar 31 09:00:20 -12 2026 -I: pbuilder-time-stamp: 1774990820 +I: removing directory /srv/workspace/pbuilder/2992669 and its subdirectories +I: Current time: Thu Feb 27 04:47:04 +14 2025 +I: pbuilder-time-stamp: 1740581224