Diff of the two buildlogs: -- --- b1/build.log 2024-05-18 21:37:36.632775053 +0000 +++ b2/build.log 2024-05-18 21:42:37.879949095 +0000 @@ -1,6 +1,6 @@ I: pbuilder: network access will be disabled during build -I: Current time: Sat May 18 09:28:45 -12 2024 -I: pbuilder-time-stamp: 1716067725 +I: Current time: Sat Jun 21 18:00:39 +14 2025 +I: pbuilder-time-stamp: 1750478439 I: Building the build Environment I: extracting base tarball [/var/cache/pbuilder/unstable-reproducible-base.tgz] I: copying local configuration @@ -26,52 +26,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/1023116/tmp/hooks/D02_print_environment starting +I: user script /srv/workspace/pbuilder/1800908/tmp/hooks/D01_modify_environment starting +debug: Running on ionos15-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 Jun 21 04:00 /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/1800908/tmp/hooks/D01_modify_environment finished +I: user script /srv/workspace/pbuilder/1800908/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=20 ' - DISTRIBUTION='unstable' - 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]="21" [3]="1" [4]="release" [5]="x86_64-pc-linux-gnu") + BASH_VERSION='5.2.21(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=42 ' + DIRSTACK=() + DISTRIBUTION=unstable + EUID=0 + FUNCNAME=([0]="Echo" [1]="main") + GROUPS=() + HOME=/root + HOSTNAME=i-capture-the-hostname + HOSTTYPE=x86_64 + HOST_ARCH=amd64 IFS=' ' - INVOCATION_ID='e90d8204c40c4b5697dd0db4ae854b93' - 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='1023116' - PS1='# ' - PS2='> ' + INVOCATION_ID=312f255ab3f542698fb4e58762d51659 + 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=1800908 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.Xl5cZTdO/pbuilderrc_5n9Z --distribution unstable --hookdir /etc/pbuilder/first-build-hooks --debbuildopts -b --basetgz /var/cache/pbuilder/unstable-reproducible-base.tgz --buildresult /srv/reproducible-results/rbuild-debian/r-b-build.Xl5cZTdO/b1 --logfile b1/build.log octave-linear-algebra_2.2.3-4.dsc' - SUDO_GID='110' - SUDO_UID='105' - SUDO_USER='jenkins' - TERM='unknown' - TZ='/usr/share/zoneinfo/Etc/GMT+12' - USER='root' - _='/usr/bin/systemd-run' - http_proxy='http://46.16.76.132: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.Xl5cZTdO/pbuilderrc_OYQF --distribution unstable --hookdir /etc/pbuilder/rebuild-hooks --debbuildopts -b --basetgz /var/cache/pbuilder/unstable-reproducible-base.tgz --buildresult /srv/reproducible-results/rbuild-debian/r-b-build.Xl5cZTdO/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://213.165.73.152:3128 I: uname -a - Linux ionos1-amd64 6.1.0-21-amd64 #1 SMP PREEMPT_DYNAMIC Debian 6.1.90-1 (2024-05-03) x86_64 GNU/Linux + Linux i-capture-the-hostname 6.6.13+bpo-amd64 #1 SMP PREEMPT_DYNAMIC Debian 6.6.13-1~bpo12+1 (2024-02-15) x86_64 GNU/Linux I: ls -l /bin - lrwxrwxrwx 1 root root 7 May 17 07:43 /bin -> usr/bin -I: user script /srv/workspace/pbuilder/1023116/tmp/hooks/D02_print_environment finished + lrwxrwxrwx 1 root root 7 Jun 20 14:05 /bin -> usr/bin +I: user script /srv/workspace/pbuilder/1800908/tmp/hooks/D02_print_environment finished -> Attempting to satisfy build-dependencies -> Creating pbuilder-satisfydepends-dummy package Package: pbuilder-satisfydepends-dummy @@ -651,7 +683,7 @@ Get: 540 http://deb.debian.org/debian unstable/main amd64 gfortran amd64 4:13.2.0-7 [1440 B] Get: 541 http://deb.debian.org/debian unstable/main amd64 octave-dev amd64 9.1.0-3 [991 kB] Get: 542 http://deb.debian.org/debian unstable/main amd64 dh-octave all 1.8.0 [22.7 kB] -Fetched 234 MB in 22s (10.7 MB/s) +Fetched 234 MB in 5s (44.7 MB/s) debconf: delaying package configuration, since apt-utils is not installed 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 ... 19718 files and directories currently installed.) @@ -2876,7 +2908,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/1800908/tmp/hooks/A99_set_merged_usr starting +Not re-configuring usrmerge for unstable +I: user script /srv/workspace/pbuilder/1800908/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 @@ -2914,6 +2950,38 @@ Checking package... Run the unit tests... Checking m files ... +[inst/smwsolve.m] +>>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/smwsolve.m +***** 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); +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.UtAnmInGl0 at line 30 column 31 + +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 ***** shared m, n, k, A @@ -2926,12 +2994,12 @@ Stop: tolerance reached. ***** test [W,H,iter,HIS]=nmf_bpas(A,k,'verbose'); - initGrNormW = 352.95 - initGrNormH = 426.19 - initNorm = 14.263 - SC1 = 1.1067 - SC2 = 553.37 - SC3 = 0.4895 + initGrNormW = 365.90 + initGrNormH = 418.01 + initNorm = 13.836 + SC1 = 1.1111 + SC2 = 555.53 + SC3 = 0.4851 W_density = 1 H_density = 1 par = @@ -2942,8 +3010,8 @@ n = 20 type = regularized nnls_solver = bp - alpha = 0.5082 - beta = 0.5082 + alpha = 0.4893 + beta = 0.4893 max_iter = 100 min_iter = 20 max_time = 1000 @@ -2955,438 +3023,438 @@ scalar structure containing the fields: - initGrNormW = 352.95 - initGrNormH = 426.19 - initNorm = 14.263 - SC1 = 0.030812 - SC2 = 0.030010 - SC3 = 2.1956 - W_density = 0.9633 - H_density = 0.7150 + initGrNormW = 365.90 + initGrNormH = 418.01 + initNorm = 13.836 + SC1 = 0.027244 + SC2 = 0.026263 + SC3 = 1.8675 + W_density = 0.9500 + H_density = 0.6800 iter = 1 - elapsed = 0.027773 - tTotal = 0.027773 - subIterW = 6 - subIterH = 24 - relError = 0.4282 + elapsed = 0.017143 + tTotal = 0.017143 + subIterW = 7 + subIterH = 31 + relError = 0.4368 toshow = scalar structure containing the fields: - initGrNormW = 352.95 - initGrNormH = 426.19 - initNorm = 14.263 - SC1 = 0.016821 - SC2 = 0.015946 - SC3 = 1.1856 - W_density = 0.9233 - H_density = 0.5250 + initGrNormW = 365.90 + initGrNormH = 418.01 + initNorm = 13.836 + SC1 = 0.015332 + SC2 = 0.014197 + SC3 = 1.0380 + W_density = 0.8900 + H_density = 0.5600 iter = 2 - elapsed = 0.029295 - tTotal = 0.057068 - subIterW = 11 - subIterH = 25 - relError = 0.3362 + elapsed = 0.016777 + tTotal = 0.033920 + subIterW = 12 + subIterH = 21 + relError = 0.3488 toshow = scalar structure containing the fields: - initGrNormW = 352.95 - initGrNormH = 426.19 - initNorm = 14.263 - SC1 = 0.010995 - SC2 = 0.010292 - SC3 = 0.7657 - W_density = 0.9033 - H_density = 0.5650 + initGrNormW = 365.90 + initGrNormH = 418.01 + initNorm = 13.836 + SC1 = 0.010329 + SC2 = 9.3577e-03 + SC3 = 0.6919 + W_density = 0.8767 + H_density = 0.6050 iter = 3 - elapsed = 0.026477 - tTotal = 0.083545 - subIterW = 4 - subIterH = 15 - relError = 0.3092 + elapsed = 0.012673 + tTotal = 0.046593 + subIterW = 6 + subIterH = 14 + relError = 0.3195 toshow = scalar structure containing the fields: - initGrNormW = 352.95 - initGrNormH = 426.19 - initNorm = 14.263 - SC1 = 7.8249e-03 - SC2 = 7.1832e-03 - SC3 = 0.5339 - W_density = 0.8767 - H_density = 0.6400 + initGrNormW = 365.90 + initGrNormH = 418.01 + initNorm = 13.836 + SC1 = 7.7721e-03 + SC2 = 6.6995e-03 + SC3 = 0.5159 + W_density = 0.8400 + H_density = 0.6200 iter = 4 - elapsed = 0.021861 - tTotal = 0.1054 - subIterW = 7 - subIterH = 9 - relError = 0.2931 + elapsed = 0.012570 + tTotal = 0.059163 + subIterW = 9 + subIterH = 8 + relError = 0.3014 toshow = scalar structure containing the fields: - initGrNormW = 352.95 - initGrNormH = 426.19 - initNorm = 14.263 - SC1 = 5.7505e-03 - SC2 = 5.1525e-03 - SC3 = 0.3920 - W_density = 0.8567 - H_density = 0.6650 + initGrNormW = 365.90 + initGrNormH = 418.01 + initNorm = 13.836 + SC1 = 5.8393e-03 + SC2 = 4.9167e-03 + SC3 = 0.3966 + W_density = 0.8333 + H_density = 0.6400 iter = 5 - elapsed = 0.023690 - tTotal = 0.1291 + elapsed = 0.010954 + tTotal = 0.070117 subIterW = 8 - subIterH = 8 - relError = 0.2832 + subIterH = 7 + relError = 0.2900 toshow = scalar structure containing the fields: - initGrNormW = 352.95 - initGrNormH = 426.19 - initNorm = 14.263 - SC1 = 4.3641e-03 - SC2 = 3.8404e-03 - SC3 = 0.2977 - W_density = 0.8433 - H_density = 0.6900 + initGrNormW = 365.90 + initGrNormH = 418.01 + initNorm = 13.836 + SC1 = 4.4618e-03 + SC2 = 3.6765e-03 + SC3 = 0.3004 + W_density = 0.8167 + H_density = 0.6550 iter = 6 - elapsed = 0.018196 - tTotal = 0.1473 - subIterW = 4 - subIterH = 6 - relError = 0.2778 + elapsed = 0.012736 + tTotal = 0.082853 + subIterW = 7 + subIterH = 4 + relError = 0.2836 toshow = scalar structure containing the fields: - initGrNormW = 352.95 - initGrNormH = 426.19 - initNorm = 14.263 - SC1 = 3.4792e-03 - SC2 = 2.9713e-03 - SC3 = 0.2363 - W_density = 0.8200 - H_density = 0.7200 + initGrNormW = 365.90 + initGrNormH = 418.01 + initNorm = 13.836 + SC1 = 3.4705e-03 + SC2 = 2.8181e-03 + SC3 = 0.2296 + W_density = 0.8033 + H_density = 0.6800 iter = 7 - elapsed = 0.019288 - tTotal = 0.1666 + elapsed = 0.010959 + tTotal = 0.093812 subIterW = 6 - subIterH = 6 - relError = 0.2747 + subIterH = 7 + relError = 0.2797 toshow = scalar structure containing the fields: - initGrNormW = 352.95 - initGrNormH = 426.19 - initNorm = 14.263 - SC1 = 2.7819e-03 - SC2 = 2.3591e-03 - SC3 = 0.1884 - W_density = 0.8200 - H_density = 0.7250 + initGrNormW = 365.90 + initGrNormH = 418.01 + initNorm = 13.836 + SC1 = 2.8036e-03 + SC2 = 2.2149e-03 + SC3 = 0.1862 + W_density = 0.7900 + H_density = 0.6900 iter = 8 - elapsed = 0.019615 - tTotal = 0.1862 - subIterW = 2 - subIterH = 5 - relError = 0.2729 + elapsed = 0.010413 + tTotal = 0.1042 + subIterW = 6 + subIterH = 0 + relError = 0.2769 toshow = scalar structure containing the fields: - initGrNormW = 352.95 - initGrNormH = 426.19 - initNorm = 14.263 - SC1 = 2.2970e-03 - SC2 = 1.9203e-03 - SC3 = 0.1524 - W_density = 0.8067 - H_density = 0.7600 + initGrNormW = 365.90 + initGrNormH = 418.01 + initNorm = 13.836 + SC1 = 2.2356e-03 + SC2 = 1.7527e-03 + SC3 = 0.1481 + W_density = 0.7867 + H_density = 0.6950 iter = 9 - elapsed = 0.019168 - tTotal = 0.2054 - subIterW = 4 - subIterH = 7 - relError = 0.2717 + elapsed = 0.010427 + tTotal = 0.1147 + subIterW = 5 + subIterH = 3 + relError = 0.2748 toshow = scalar structure containing the fields: - initGrNormW = 352.95 - initGrNormH = 426.19 - initNorm = 14.263 - SC1 = 1.9341e-03 - SC2 = 1.5899e-03 - SC3 = 0.1280 - W_density = 0.8033 - H_density = 0.7700 + initGrNormW = 365.90 + initGrNormH = 418.01 + initNorm = 13.836 + SC1 = 1.8008e-03 + SC2 = 1.4047e-03 + SC3 = 0.1172 + W_density = 0.7800 + H_density = 0.6950 iter = 10 - elapsed = 0.016932 - tTotal = 0.2223 - subIterW = 5 - subIterH = 3 - relError = 0.2708 + elapsed = 7.9060e-03 + tTotal = 0.1226 + subIterW = 1 + subIterH = 0 + relError = 0.2730 toshow = scalar structure containing the fields: - initGrNormW = 352.95 - initGrNormH = 426.19 - initNorm = 14.263 - SC1 = 1.6496e-03 - SC2 = 1.3395e-03 - SC3 = 0.1081 - W_density = 0.8000 - H_density = 0.7700 + initGrNormW = 365.90 + initGrNormH = 418.01 + initNorm = 13.836 + SC1 = 1.4713e-03 + SC2 = 1.1270e-03 + SC3 = 0.094935 + W_density = 0.7700 + H_density = 0.7000 iter = 11 - elapsed = 0.015034 - tTotal = 0.2373 - subIterW = 0 + elapsed = 0.010828 + tTotal = 0.1334 + subIterW = 7 subIterH = 3 - relError = 0.2700 + relError = 0.2716 toshow = scalar structure containing the fields: - initGrNormW = 352.95 - initGrNormH = 426.19 - initNorm = 14.263 - SC1 = 1.4189e-03 - SC2 = 1.1380e-03 - SC3 = 0.091113 - W_density = 0.7900 - H_density = 0.7750 + initGrNormW = 365.90 + initGrNormH = 418.01 + initNorm = 13.836 + SC1 = 1.1978e-03 + SC2 = 9.0795e-04 + SC3 = 0.075972 + W_density = 0.7600 + H_density = 0.7050 iter = 12 - elapsed = 0.015409 - tTotal = 0.2527 + elapsed = 9.9100e-03 + tTotal = 0.1433 subIterW = 3 - subIterH = 3 - relError = 0.2694 + subIterH = 0 + relError = 0.2704 toshow = scalar structure containing the fields: - initGrNormW = 352.95 - initGrNormH = 426.19 - initNorm = 14.263 - SC1 = 1.2285e-03 - SC2 = 9.7541e-04 - SC3 = 0.077640 - W_density = 0.7867 - H_density = 0.7750 + initGrNormW = 365.90 + initGrNormH = 418.01 + initNorm = 13.836 + SC1 = 9.8170e-04 + SC2 = 7.4020e-04 + SC3 = 0.062182 + W_density = 0.7633 + H_density = 0.7200 iter = 13 - elapsed = 0.017153 - tTotal = 0.2699 - subIterW = 0 - subIterH = 2 - relError = 0.2688 + elapsed = 9.9450e-03 + tTotal = 0.1532 + subIterW = 3 + subIterH = 3 + relError = 0.2695 toshow = scalar structure containing the fields: - initGrNormW = 352.95 - initGrNormH = 426.19 - initNorm = 14.263 - SC1 = 1.0624e-03 - SC2 = 8.3931e-04 - SC3 = 0.065293 - W_density = 0.7800 - H_density = 0.7750 + initGrNormW = 365.90 + initGrNormH = 418.01 + initNorm = 13.836 + SC1 = 8.0844e-04 + SC2 = 6.1441e-04 + SC3 = 0.049566 + W_density = 0.7633 + H_density = 0.7250 iter = 14 - elapsed = 0.013651 - tTotal = 0.2835 + elapsed = 8.7740e-03 + tTotal = 0.1620 subIterW = 2 subIterH = 0 - relError = 0.2682 + relError = 0.2686 toshow = scalar structure containing the fields: - initGrNormW = 352.95 - initGrNormH = 426.19 - initNorm = 14.263 - SC1 = 9.1086e-04 - SC2 = 7.2140e-04 - SC3 = 0.033116 - W_density = 0.7800 - H_density = 0.7700 + initGrNormW = 365.90 + initGrNormH = 418.01 + initNorm = 13.836 + SC1 = 6.7694e-04 + SC2 = 5.1582e-04 + SC3 = 0.040676 + W_density = 0.7600 + H_density = 0.7250 iter = 15 - elapsed = 0.012994 - tTotal = 0.2965 + elapsed = 9.9680e-03 + tTotal = 0.1720 subIterW = 0 - subIterH = 3 - relError = 0.2676 + subIterH = 2 + relError = 0.2679 toshow = scalar structure containing the fields: - initGrNormW = 352.95 - initGrNormH = 426.19 - initNorm = 14.263 - SC1 = 7.9136e-04 - SC2 = 6.2359e-04 - SC3 = 0.047797 - W_density = 0.7800 - H_density = 0.7650 + initGrNormW = 365.90 + initGrNormH = 418.01 + initNorm = 13.836 + SC1 = 5.9225e-04 + SC2 = 4.5129e-04 + SC3 = 0.020962 + W_density = 0.7600 + H_density = 0.7350 iter = 16 - elapsed = 0.015785 - tTotal = 0.3123 - subIterW = 2 - subIterH = 3 - relError = 0.2672 + elapsed = 8.4630e-03 + tTotal = 0.1804 + subIterW = 0 + subIterH = 5 + relError = 0.2673 toshow = scalar structure containing the fields: - initGrNormW = 352.95 - initGrNormH = 426.19 - initNorm = 14.263 - SC1 = 6.7980e-04 - SC2 = 5.3432e-04 - SC3 = 0.040555 - W_density = 0.7733 - H_density = 0.7650 + initGrNormW = 365.90 + initGrNormH = 418.01 + initNorm = 13.836 + SC1 = 5.4257e-04 + SC2 = 4.1344e-04 + SC3 = 0.031762 + W_density = 0.7567 + H_density = 0.7400 iter = 17 - elapsed = 0.014506 - tTotal = 0.3268 - subIterW = 2 + elapsed = 9.1350e-03 + tTotal = 0.1896 + subIterW = 0 subIterH = 0 - relError = 0.2668 + relError = 0.2667 toshow = scalar structure containing the fields: - initGrNormW = 352.95 - initGrNormH = 426.19 - initNorm = 14.263 - SC1 = 5.8390e-04 - SC2 = 4.5894e-04 - SC3 = 0.019595 - W_density = 0.7733 - H_density = 0.7600 + initGrNormW = 365.90 + initGrNormH = 418.01 + initNorm = 13.836 + SC1 = 5.1330e-04 + SC2 = 3.8908e-04 + SC3 = 0.029804 + W_density = 0.7567 + H_density = 0.7400 iter = 18 - elapsed = 0.013299 - tTotal = 0.3401 - subIterW = 0 - subIterH = 0 - relError = 0.2666 + elapsed = 0.010096 + tTotal = 0.1997 + subIterW = 2 + subIterH = 2 + relError = 0.2662 toshow = scalar structure containing the fields: - initGrNormW = 352.95 - initGrNormH = 426.19 - initNorm = 14.263 - SC1 = 5.0909e-04 - SC2 = 3.9811e-04 - SC3 = 0.030174 - W_density = 0.7700 - H_density = 0.7650 + initGrNormW = 365.90 + initGrNormH = 418.01 + initNorm = 13.836 + SC1 = 4.9860e-04 + SC2 = 3.7794e-04 + SC3 = 0.028529 + W_density = 0.7533 + H_density = 0.7550 iter = 19 - elapsed = 0.018430 - tTotal = 0.3586 + elapsed = 0.010003 + tTotal = 0.2097 subIterW = 0 - subIterH = 0 - relError = 0.2663 + subIterH = 3 + relError = 0.2656 toshow = scalar structure containing the fields: - initGrNormW = 352.95 - initGrNormH = 426.19 - initNorm = 14.263 - SC1 = 4.4877e-04 - SC2 = 3.4914e-04 - SC3 = 0.015937 - W_density = 0.7700 - H_density = 0.7600 + initGrNormW = 365.90 + initGrNormH = 418.01 + initNorm = 13.836 + SC1 = 4.9283e-04 + SC2 = 3.7061e-04 + SC3 = 0.028133 + W_density = 0.7500 + H_density = 0.7500 iter = 20 - elapsed = 0.013845 - tTotal = 0.3724 + elapsed = 0.010002 + tTotal = 0.2197 subIterW = 0 - subIterH = 0 - relError = 0.2662 + subIterH = 3 + relError = 0.2650 toshow = scalar structure containing the fields: - initGrNormW = 352.95 - initGrNormH = 426.19 - initNorm = 14.263 - SC1 = 4.0116e-04 - SC2 = 3.1050e-04 - SC3 = 0.024045 - W_density = 0.7700 - H_density = 0.7700 + initGrNormW = 365.90 + initGrNormH = 418.01 + initNorm = 13.836 + SC1 = 4.8264e-04 + SC2 = 3.6391e-04 + SC3 = 0.027794 + W_density = 0.7567 + H_density = 0.7450 iter = 21 - elapsed = 0.020951 - tTotal = 0.3934 - subIterW = 2 - subIterH = 2 - relError = 0.2660 + elapsed = 9.1920e-03 + tTotal = 0.2289 + subIterW = 0 + subIterH = 0 + relError = 0.2644 toshow = scalar structure containing the fields: - initGrNormW = 352.95 - initGrNormH = 426.19 - initNorm = 14.263 - SC1 = 3.6124e-04 - SC2 = 2.7888e-04 - SC3 = 0.021470 - W_density = 0.7667 - H_density = 0.7750 + initGrNormW = 365.90 + initGrNormH = 418.01 + initNorm = 13.836 + SC1 = 4.6764e-04 + SC2 = 3.5634e-04 + SC3 = 0.026773 + W_density = 0.7600 + H_density = 0.7500 iter = 22 - elapsed = 0.015754 - tTotal = 0.4091 - subIterW = 0 - subIterH = 2 - relError = 0.2659 + elapsed = 0.010665 + tTotal = 0.2395 + subIterW = 3 + subIterH = 3 + relError = 0.2637 toshow = scalar structure containing the fields: - initGrNormW = 352.95 - initGrNormH = 426.19 - initNorm = 14.263 - SC1 = 3.2848e-04 - SC2 = 2.5227e-04 - SC3 = 0.019200 - W_density = 0.7600 - H_density = 0.7700 + initGrNormW = 365.90 + initGrNormH = 418.01 + initNorm = 13.836 + SC1 = 4.5021e-04 + SC2 = 3.4666e-04 + SC3 = 0.026065 + W_density = 0.7733 + H_density = 0.7550 iter = 23 - elapsed = 0.014860 - tTotal = 0.4240 - subIterW = 2 - subIterH = 0 - relError = 0.2658 + elapsed = 0.010450 + tTotal = 0.2500 + subIterW = 6 + subIterH = 1 + relError = 0.2631 Stop: tolerance reached. final = @@ -3394,19 +3462,19 @@ scalar structure containing the fields: iterations = 23 - elapsed_total = 0.4240 - relative_error = 0.2658 - W_density = 0.7600 - H_density = 0.7700 + elapsed_total = 0.2500 + relative_error = 0.2631 + W_density = 0.7733 + H_density = 0.7550 ***** test [W,H,iter,HIS]=nmf_bpas(A,k,'verbose','nnlssolver','as'); - initGrNormW = 379.23 - initGrNormH = 485.73 - initNorm = 14.263 - SC1 = 1.2325 - SC2 = 616.24 - SC3 = 0.5092 + initGrNormW = 366.15 + initGrNormH = 442.83 + initNorm = 13.836 + SC1 = 1.1492 + SC2 = 574.60 + SC3 = 0.4939 W_density = 1 H_density = 1 par = @@ -3417,8 +3485,8 @@ n = 20 type = regularized nnls_solver = as - alpha = 0.5082 - beta = 0.5082 + alpha = 0.4893 + beta = 0.4893 max_iter = 100 min_iter = 20 max_time = 1000 @@ -3430,438 +3498,438 @@ scalar structure containing the fields: - initGrNormW = 379.23 - initGrNormH = 485.73 - initNorm = 14.263 - SC1 = 0.029465 - SC2 = 0.028934 - SC3 = 1.2733 - W_density = 0.9700 - H_density = 0.6500 + initGrNormW = 366.15 + initGrNormH = 442.83 + initNorm = 13.836 + SC1 = 0.027740 + SC2 = 0.026520 + SC3 = 1.0592 + W_density = 0.9600 + H_density = 0.7100 iter = 1 - elapsed = 0.015737 - tTotal = 0.015737 + elapsed = 0.015574 + tTotal = 0.015574 subIterW = 7 - subIterH = 32 - relError = 0.4305 + subIterH = 29 + relError = 0.4497 toshow = scalar structure containing the fields: - initGrNormW = 379.23 - initGrNormH = 485.73 - initNorm = 14.263 - SC1 = 0.015745 - SC2 = 0.015084 - SC3 = 0.6530 - W_density = 0.9333 - H_density = 0.5250 + initGrNormW = 366.15 + initGrNormH = 442.83 + initNorm = 13.836 + SC1 = 0.015081 + SC2 = 0.013935 + SC3 = 0.5656 + W_density = 0.9167 + H_density = 0.5700 iter = 2 - elapsed = 0.017227 - tTotal = 0.032964 - subIterW = 12 - subIterH = 42 - relError = 0.3362 + elapsed = 0.012608 + tTotal = 0.028182 + subIterW = 13 + subIterH = 34 + relError = 0.3549 toshow = scalar structure containing the fields: - initGrNormW = 379.23 - initGrNormH = 485.73 - initNorm = 14.263 - SC1 = 0.010220 - SC2 = 9.6272e-03 - SC3 = 0.4555 - W_density = 0.9100 - H_density = 0.5800 + initGrNormW = 366.15 + initGrNormH = 442.83 + initNorm = 13.836 + SC1 = 9.9482e-03 + SC2 = 9.0529e-03 + SC3 = 0.3867 + W_density = 0.8967 + H_density = 0.6100 iter = 3 - elapsed = 0.018860 - tTotal = 0.051824 + elapsed = 0.011888 + tTotal = 0.040070 subIterW = 14 - subIterH = 42 - relError = 0.3058 + subIterH = 35 + relError = 0.3242 toshow = scalar structure containing the fields: - initGrNormW = 379.23 - initGrNormH = 485.73 - initNorm = 14.263 - SC1 = 7.3129e-03 - SC2 = 6.6255e-03 - SC3 = 0.3124 - W_density = 0.8667 + initGrNormW = 366.15 + initGrNormH = 442.83 + initNorm = 13.836 + SC1 = 7.1875e-03 + SC2 = 6.3250e-03 + SC3 = 0.2779 + W_density = 0.8700 H_density = 0.6400 iter = 4 - elapsed = 0.016791 - tTotal = 0.068615 - subIterW = 23 - subIterH = 35 - relError = 0.2912 + elapsed = 0.012141 + tTotal = 0.052211 + subIterW = 22 + subIterH = 32 + relError = 0.3061 toshow = scalar structure containing the fields: - initGrNormW = 379.23 - initGrNormH = 485.73 - initNorm = 14.263 - SC1 = 5.3425e-03 - SC2 = 4.7976e-03 - SC3 = 0.2236 - W_density = 0.8600 - H_density = 0.6500 + initGrNormW = 366.15 + initGrNormH = 442.83 + initNorm = 13.836 + SC1 = 5.3211e-03 + SC2 = 4.5442e-03 + SC3 = 0.2066 + W_density = 0.8333 + H_density = 0.6750 iter = 5 - elapsed = 0.018634 - tTotal = 0.087249 - subIterW = 22 - subIterH = 35 - relError = 0.2833 + elapsed = 0.012140 + tTotal = 0.064351 + subIterW = 26 + subIterH = 28 + relError = 0.2936 toshow = scalar structure containing the fields: - initGrNormW = 379.23 - initGrNormH = 485.73 - initNorm = 14.263 - SC1 = 4.0711e-03 - SC2 = 3.6151e-03 - SC3 = 0.1769 - W_density = 0.8500 - H_density = 0.6650 + initGrNormW = 366.15 + initGrNormH = 442.83 + initNorm = 13.836 + SC1 = 3.9380e-03 + SC2 = 3.3552e-03 + SC3 = 0.1526 + W_density = 0.8367 + H_density = 0.6800 iter = 6 - elapsed = 0.020803 - tTotal = 0.1081 - subIterW = 22 - subIterH = 34 - relError = 0.2783 + elapsed = 0.011657 + tTotal = 0.076008 + subIterW = 24 + subIterH = 28 + relError = 0.2845 toshow = scalar structure containing the fields: - initGrNormW = 379.23 - initGrNormH = 485.73 - initNorm = 14.263 - SC1 = 3.2387e-03 - SC2 = 2.8112e-03 - SC3 = 0.1302 - W_density = 0.8367 - H_density = 0.7100 + initGrNormW = 366.15 + initGrNormH = 442.83 + initNorm = 13.836 + SC1 = 3.0452e-03 + SC2 = 2.5519e-03 + SC3 = 0.1167 + W_density = 0.8333 + H_density = 0.7150 iter = 7 - elapsed = 0.018207 - tTotal = 0.1263 - subIterW = 25 - subIterH = 26 - relError = 0.2750 + elapsed = 0.010507 + tTotal = 0.086515 + subIterW = 27 + subIterH = 24 + relError = 0.2782 toshow = scalar structure containing the fields: - initGrNormW = 379.23 - initGrNormH = 485.73 - initNorm = 14.263 - SC1 = 2.6119e-03 - SC2 = 2.2201e-03 - SC3 = 0.1083 - W_density = 0.8267 - H_density = 0.7300 + initGrNormW = 366.15 + initGrNormH = 442.83 + initNorm = 13.836 + SC1 = 2.4400e-03 + SC2 = 2.0106e-03 + SC3 = 0.091593 + W_density = 0.8133 + H_density = 0.7200 iter = 8 - elapsed = 0.016388 - tTotal = 0.1426 - subIterW = 26 - subIterH = 28 - relError = 0.2731 + elapsed = 0.012275 + tTotal = 0.098790 + subIterW = 30 + subIterH = 30 + relError = 0.2740 toshow = scalar structure containing the fields: - initGrNormW = 379.23 - initGrNormH = 485.73 - initNorm = 14.263 - SC1 = 2.1142e-03 - SC2 = 1.7844e-03 - SC3 = 0.089997 - W_density = 0.8200 - H_density = 0.7550 + initGrNormW = 366.15 + initGrNormH = 442.83 + initNorm = 13.836 + SC1 = 1.9960e-03 + SC2 = 1.6247e-03 + SC3 = 0.073457 + W_density = 0.8067 + H_density = 0.7500 iter = 9 - elapsed = 0.015861 - tTotal = 0.1585 - subIterW = 26 + elapsed = 0.010242 + tTotal = 0.1090 + subIterW = 32 subIterH = 25 - relError = 0.2719 + relError = 0.2712 toshow = scalar structure containing the fields: - initGrNormW = 379.23 - initGrNormH = 485.73 - initNorm = 14.263 - SC1 = 1.7568e-03 - SC2 = 1.4581e-03 - SC3 = 0.078311 - W_density = 0.8100 - H_density = 0.7600 + initGrNormW = 366.15 + initGrNormH = 442.83 + initNorm = 13.836 + SC1 = 1.6407e-03 + SC2 = 1.3257e-03 + SC3 = 0.064830 + W_density = 0.8033 + H_density = 0.7700 iter = 10 - elapsed = 0.015976 - tTotal = 0.1745 - subIterW = 24 - subIterH = 21 - relError = 0.2710 + elapsed = 0.012689 + tTotal = 0.1217 + subIterW = 30 + subIterH = 23 + relError = 0.2691 toshow = scalar structure containing the fields: - initGrNormW = 379.23 - initGrNormH = 485.73 - initNorm = 14.263 - SC1 = 1.4933e-03 - SC2 = 1.2095e-03 - SC3 = 0.063728 - W_density = 0.7933 - H_density = 0.7700 + initGrNormW = 366.15 + initGrNormH = 442.83 + initNorm = 13.836 + SC1 = 1.3505e-03 + SC2 = 1.0858e-03 + SC3 = 0.050574 + W_density = 0.8000 + H_density = 0.7750 iter = 11 - elapsed = 0.023253 - tTotal = 0.1977 - subIterW = 29 - subIterH = 20 - relError = 0.2703 + elapsed = 0.011159 + tTotal = 0.1329 + subIterW = 31 + subIterH = 21 + relError = 0.2677 toshow = scalar structure containing the fields: - initGrNormW = 379.23 - initGrNormH = 485.73 - initNorm = 14.263 - SC1 = 1.2526e-03 - SC2 = 1.0071e-03 - SC3 = 0.052486 + initGrNormW = 366.15 + initGrNormH = 442.83 + initNorm = 13.836 + SC1 = 1.1199e-03 + SC2 = 8.8245e-04 + SC3 = 0.041579 W_density = 0.7867 - H_density = 0.7650 + H_density = 0.7550 iter = 12 - elapsed = 0.019405 - tTotal = 0.2171 - subIterW = 32 - subIterH = 22 - relError = 0.2699 + elapsed = 9.6370e-03 + tTotal = 0.1425 + subIterW = 31 + subIterH = 25 + relError = 0.2669 toshow = scalar structure containing the fields: - initGrNormW = 379.23 - initGrNormH = 485.73 - initNorm = 14.263 - SC1 = 1.0663e-03 - SC2 = 8.4879e-04 - SC3 = 0.043760 + initGrNormW = 366.15 + initGrNormH = 442.83 + initNorm = 13.836 + SC1 = 9.2165e-04 + SC2 = 7.2626e-04 + SC3 = 0.034716 W_density = 0.7867 - H_density = 0.7750 + H_density = 0.7550 iter = 13 - elapsed = 0.015900 - tTotal = 0.2330 - subIterW = 28 - subIterH = 21 - relError = 0.2694 + elapsed = 9.9320e-03 + tTotal = 0.1524 + subIterW = 31 + subIterH = 26 + relError = 0.2662 toshow = scalar structure containing the fields: - initGrNormW = 379.23 - initGrNormH = 485.73 - initNorm = 14.263 - SC1 = 9.2083e-04 - SC2 = 7.3114e-04 - SC3 = 0.037082 - W_density = 0.7867 - H_density = 0.7750 + initGrNormW = 366.15 + initGrNormH = 442.83 + initNorm = 13.836 + SC1 = 7.7446e-04 + SC2 = 6.0563e-04 + SC3 = 0.027931 + W_density = 0.7800 + H_density = 0.7450 iter = 14 - elapsed = 0.013343 - tTotal = 0.2464 - subIterW = 28 - subIterH = 21 - relError = 0.2690 + elapsed = 0.010161 + tTotal = 0.1626 + subIterW = 32 + subIterH = 28 + relError = 0.2657 toshow = scalar structure containing the fields: - initGrNormW = 379.23 - initGrNormH = 485.73 - initNorm = 14.263 - SC1 = 8.0790e-04 - SC2 = 6.4471e-04 - SC3 = 0.032071 - W_density = 0.7900 - H_density = 0.7700 + initGrNormW = 366.15 + initGrNormH = 442.83 + initNorm = 13.836 + SC1 = 6.5593e-04 + SC2 = 5.0769e-04 + SC3 = 0.023134 + W_density = 0.7733 + H_density = 0.7500 iter = 15 - elapsed = 0.013616 - tTotal = 0.2600 - subIterW = 27 - subIterH = 22 - relError = 0.2685 + elapsed = 0.010312 + tTotal = 0.1729 + subIterW = 33 + subIterH = 29 + relError = 0.2654 toshow = scalar structure containing the fields: - initGrNormW = 379.23 - initGrNormH = 485.73 - initNorm = 14.263 - SC1 = 7.3325e-04 - SC2 = 5.8220e-04 - SC3 = 0.029966 - W_density = 0.7867 - H_density = 0.7800 + initGrNormW = 366.15 + initGrNormH = 442.83 + initNorm = 13.836 + SC1 = 5.5936e-04 + SC2 = 4.3182e-04 + SC3 = 0.020079 + W_density = 0.7767 + H_density = 0.7400 iter = 16 - elapsed = 0.015088 - tTotal = 0.2751 - subIterW = 31 - subIterH = 21 - relError = 0.2679 + elapsed = 0.011651 + tTotal = 0.1846 + subIterW = 32 + subIterH = 27 + relError = 0.2652 toshow = scalar structure containing the fields: - initGrNormW = 379.23 - initGrNormH = 485.73 - initNorm = 14.263 - SC1 = 6.7877e-04 - SC2 = 5.3623e-04 - SC3 = 0.025356 + initGrNormW = 366.15 + initGrNormH = 442.83 + initNorm = 13.836 + SC1 = 4.8527e-04 + SC2 = 3.7463e-04 + SC3 = 0.017508 W_density = 0.7800 - H_density = 0.7850 + H_density = 0.7450 iter = 17 - elapsed = 0.017379 - tTotal = 0.2925 - subIterW = 32 - subIterH = 18 - relError = 0.2673 + elapsed = 0.010556 + tTotal = 0.1951 + subIterW = 30 + subIterH = 25 + relError = 0.2650 toshow = scalar structure containing the fields: - initGrNormW = 379.23 - initGrNormH = 485.73 - initNorm = 14.263 - SC1 = 6.3643e-04 - SC2 = 5.0150e-04 - SC3 = 0.023844 + initGrNormW = 366.15 + initGrNormH = 442.83 + initNorm = 13.836 + SC1 = 4.3020e-04 + SC2 = 3.3211e-04 + SC3 = 0.015079 W_density = 0.7767 - H_density = 0.7850 + H_density = 0.7400 iter = 18 - elapsed = 0.019249 - tTotal = 0.3117 - subIterW = 32 - subIterH = 21 - relError = 0.2665 + elapsed = 0.011459 + tTotal = 0.2066 + subIterW = 31 + subIterH = 27 + relError = 0.2648 toshow = scalar structure containing the fields: - initGrNormW = 379.23 - initGrNormH = 485.73 - initNorm = 14.263 - SC1 = 5.9001e-04 - SC2 = 4.6492e-04 - SC3 = 0.021192 + initGrNormW = 366.15 + initGrNormH = 442.83 + initNorm = 13.836 + SC1 = 3.7986e-04 + SC2 = 2.9174e-04 + SC3 = 0.013227 W_density = 0.7767 - H_density = 0.7850 + H_density = 0.7400 iter = 19 - elapsed = 0.024221 - tTotal = 0.3359 - subIterW = 32 - subIterH = 22 - relError = 0.2656 + elapsed = 0.011501 + tTotal = 0.2181 + subIterW = 30 + subIterH = 26 + relError = 0.2647 toshow = scalar structure containing the fields: - initGrNormW = 379.23 - initGrNormH = 485.73 - initNorm = 14.263 - SC1 = 5.5031e-04 - SC2 = 4.3365e-04 - SC3 = 0.019474 + initGrNormW = 366.15 + initGrNormH = 442.83 + initNorm = 13.836 + SC1 = 3.3788e-04 + SC2 = 2.5814e-04 + SC3 = 0.011997 W_density = 0.7767 - H_density = 0.7800 + H_density = 0.7350 iter = 20 - elapsed = 0.022199 - tTotal = 0.3581 - subIterW = 33 - subIterH = 22 - relError = 0.2648 + elapsed = 0.011775 + tTotal = 0.2299 + subIterW = 30 + subIterH = 27 + relError = 0.2646 toshow = scalar structure containing the fields: - initGrNormW = 379.23 - initGrNormH = 485.73 - initNorm = 14.263 - SC1 = 5.0957e-04 - SC2 = 3.9644e-04 - SC3 = 0.017237 + initGrNormW = 366.15 + initGrNormH = 442.83 + initNorm = 13.836 + SC1 = 2.9988e-04 + SC2 = 2.2791e-04 + SC3 = 0.011026 W_density = 0.7733 - H_density = 0.7650 + H_density = 0.7300 iter = 21 - elapsed = 0.022290 - tTotal = 0.3804 - subIterW = 33 - subIterH = 24 - relError = 0.2640 + elapsed = 0.012018 + tTotal = 0.2419 + subIterW = 30 + subIterH = 28 + relError = 0.2646 toshow = scalar structure containing the fields: - initGrNormW = 379.23 - initGrNormH = 485.73 - initNorm = 14.263 - SC1 = 4.4773e-04 - SC2 = 3.4744e-04 - SC3 = 0.015447 - W_density = 0.7700 - H_density = 0.7650 + initGrNormW = 366.15 + initGrNormH = 442.83 + initNorm = 13.836 + SC1 = 2.6858e-04 + SC2 = 2.0412e-04 + SC3 = 9.7223e-03 + W_density = 0.7733 + H_density = 0.7300 iter = 22 - elapsed = 0.016448 - tTotal = 0.3969 - subIterW = 34 - subIterH = 22 - relError = 0.2634 + elapsed = 0.012513 + tTotal = 0.2544 + subIterW = 30 + subIterH = 29 + relError = 0.2645 toshow = scalar structure containing the fields: - initGrNormW = 379.23 - initGrNormH = 485.73 - initNorm = 14.263 - SC1 = 3.9888e-04 - SC2 = 3.0793e-04 - SC3 = 0.014590 - W_density = 0.7633 - H_density = 0.7650 + initGrNormW = 366.15 + initGrNormH = 442.83 + initNorm = 13.836 + SC1 = 2.4523e-04 + SC2 = 1.8637e-04 + SC3 = 8.1450e-03 + W_density = 0.7767 + H_density = 0.7300 iter = 23 - elapsed = 0.022077 - tTotal = 0.4190 - subIterW = 35 - subIterH = 22 - relError = 0.2629 + elapsed = 0.012338 + tTotal = 0.2667 + subIterW = 29 + subIterH = 29 + relError = 0.2645 Stop: tolerance reached. final = @@ -3869,19 +3937,19 @@ scalar structure containing the fields: iterations = 23 - elapsed_total = 0.4190 - relative_error = 0.2629 - W_density = 0.7633 - H_density = 0.7650 + elapsed_total = 0.2667 + relative_error = 0.2645 + W_density = 0.7767 + H_density = 0.7300 ***** test [W,H,iter,HIS]=nmf_bpas(A,k,'verbose','type','sparse'); - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 1.3251 - SC2 = 662.55 - SC3 = 0.5113 + initGrNormW = 385.10 + initGrNormH = 490.96 + initNorm = 13.836 + SC1 = 1.2480 + SC2 = 623.98 + SC3 = 0.5035 W_density = 1 H_density = 1 par = @@ -3892,8 +3960,8 @@ n = 20 type = sparse nnls_solver = bp - alpha = 0.5082 - beta = 0.5082 + alpha = 0.4893 + beta = 0.4893 max_iter = 100 min_iter = 20 max_time = 1000 @@ -3905,534 +3973,496 @@ scalar structure containing the fields: - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 0.023223 - SC2 = 0.022387 - SC3 = 1.8020 - W_density = 0.9767 - H_density = 0.5700 + initGrNormW = 385.10 + initGrNormH = 490.96 + initNorm = 13.836 + SC1 = 0.023087 + SC2 = 0.020917 + SC3 = 1.5672 + W_density = 0.9367 + H_density = 0.5500 iter = 1 - elapsed = 0.023476 - tTotal = 0.023476 - subIterW = 5 - subIterH = 32 - relError = 0.4434 + elapsed = 0.011876 + tTotal = 0.011876 + subIterW = 8 + subIterH = 34 + relError = 0.4512 toshow = scalar structure containing the fields: - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 0.011153 - SC2 = 9.2120e-03 - SC3 = 0.8541 - W_density = 0.9100 + initGrNormW = 385.10 + initGrNormH = 490.96 + initNorm = 13.836 + SC1 = 0.011261 + SC2 = 8.7838e-03 + SC3 = 0.7999 + W_density = 0.8900 H_density = 0.3600 iter = 2 - elapsed = 0.032391 - tTotal = 0.055867 + elapsed = 0.011394 + tTotal = 0.023270 subIterW = 11 - subIterH = 29 - relError = 0.3587 + subIterH = 25 + relError = 0.3725 toshow = scalar structure containing the fields: - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 6.7361e-03 - SC2 = 4.8230e-03 - SC3 = 0.4972 - W_density = 0.8333 - H_density = 0.3600 + initGrNormW = 385.10 + initGrNormH = 490.96 + initNorm = 13.836 + SC1 = 7.2748e-03 + SC2 = 5.2088e-03 + SC3 = 0.5367 + W_density = 0.8833 + H_density = 0.3400 iter = 3 - elapsed = 0.020299 - tTotal = 0.076166 - subIterW = 8 - subIterH = 15 - relError = 0.3348 + elapsed = 9.6530e-03 + tTotal = 0.032923 + subIterW = 9 + subIterH = 10 + relError = 0.3466 toshow = scalar structure containing the fields: - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 4.1853e-03 - SC2 = 2.8377e-03 - SC3 = 0.3157 - W_density = 0.8233 - H_density = 0.3450 + initGrNormW = 385.10 + initGrNormH = 490.96 + initNorm = 13.836 + SC1 = 4.9554e-03 + SC2 = 3.3003e-03 + SC3 = 0.3943 + W_density = 0.8667 + H_density = 0.3150 iter = 4 - elapsed = 0.015130 - tTotal = 0.091296 - subIterW = 5 - subIterH = 7 - relError = 0.3224 + elapsed = 0.010255 + tTotal = 0.043178 + subIterW = 6 + subIterH = 10 + relError = 0.3359 toshow = scalar structure containing the fields: - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 2.6733e-03 - SC2 = 1.7430e-03 - SC3 = 0.2116 - W_density = 0.8267 - H_density = 0.3150 + initGrNormW = 385.10 + initGrNormH = 490.96 + initNorm = 13.836 + SC1 = 3.5250e-03 + SC2 = 2.2419e-03 + SC3 = 0.2898 + W_density = 0.8533 + H_density = 0.3000 iter = 5 - elapsed = 0.018867 - tTotal = 0.1102 - subIterW = 2 - subIterH = 9 - relError = 0.3169 + elapsed = 7.4300e-03 + tTotal = 0.050608 + subIterW = 7 + subIterH = 5 + relError = 0.3330 toshow = scalar structure containing the fields: - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 1.8906e-03 - SC2 = 1.2100e-03 - SC3 = 0.1445 - W_density = 0.8267 - H_density = 0.3100 + initGrNormW = 385.10 + initGrNormH = 490.96 + initNorm = 13.836 + SC1 = 2.7193e-03 + SC2 = 1.7186e-03 + SC3 = 0.2124 + W_density = 0.8533 + H_density = 0.2850 iter = 6 - elapsed = 0.017298 - tTotal = 0.1275 - subIterW = 2 - subIterH = 4 - relError = 0.3133 + elapsed = 8.0670e-03 + tTotal = 0.058675 + subIterW = 6 + subIterH = 3 + relError = 0.3328 toshow = scalar structure containing the fields: - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 1.5279e-03 - SC2 = 9.6561e-04 - SC3 = 0.1156 - W_density = 0.8333 - H_density = 0.2900 + initGrNormW = 385.10 + initGrNormH = 490.96 + initNorm = 13.836 + SC1 = 2.0798e-03 + SC2 = 1.2770e-03 + SC3 = 0.1647 + W_density = 0.8467 + H_density = 0.2550 iter = 7 - elapsed = 0.018261 - tTotal = 0.1457 - subIterW = 3 + elapsed = 6.4510e-03 + tTotal = 0.065126 + subIterW = 4 subIterH = 5 - relError = 0.3109 + relError = 0.3336 toshow = scalar structure containing the fields: - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 1.2775e-03 - SC2 = 7.8947e-04 - SC3 = 0.096619 - W_density = 0.8267 - H_density = 0.2750 + initGrNormW = 385.10 + initGrNormH = 490.96 + initNorm = 13.836 + SC1 = 1.6465e-03 + SC2 = 9.9121e-04 + SC3 = 0.1304 + W_density = 0.8333 + H_density = 0.2500 iter = 8 - elapsed = 0.018769 - tTotal = 0.1645 + elapsed = 6.3120e-03 + tTotal = 0.071438 subIterW = 4 subIterH = 2 - relError = 0.3094 + relError = 0.3348 toshow = scalar structure containing the fields: - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 1.0828e-03 - SC2 = 6.6487e-04 - SC3 = 0.023763 - W_density = 0.8267 - H_density = 0.2650 + initGrNormW = 385.10 + initGrNormH = 490.96 + initNorm = 13.836 + SC1 = 1.5203e-03 + SC2 = 9.0306e-04 + SC3 = 0.1228 + W_density = 0.8200 + H_density = 0.2450 iter = 9 - elapsed = 0.010965 - tTotal = 0.1755 - subIterW = 0 - subIterH = 2 - relError = 0.3088 + elapsed = 6.1630e-03 + tTotal = 0.077601 + subIterW = 4 + subIterH = 0 + relError = 0.3365 toshow = scalar structure containing the fields: - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 9.6597e-04 - SC2 = 5.9310e-04 - SC3 = 0.020921 - W_density = 0.8267 - H_density = 0.2700 + initGrNormW = 385.10 + initGrNormH = 490.96 + initNorm = 13.836 + SC1 = 1.6080e-03 + SC2 = 9.3584e-04 + SC3 = 0.1320 + W_density = 0.8067 + H_density = 0.2350 iter = 10 - elapsed = 0.010753 - tTotal = 0.1862 - subIterW = 0 - subIterH = 3 - relError = 0.3086 + elapsed = 7.6690e-03 + tTotal = 0.085270 + subIterW = 4 + subIterH = 2 + relError = 0.3389 toshow = scalar structure containing the fields: - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 9.0677e-04 - SC2 = 5.5857e-04 - SC3 = 0.069353 - W_density = 0.8267 - H_density = 0.2750 + initGrNormW = 385.10 + initGrNormH = 490.96 + initNorm = 13.836 + SC1 = 1.7757e-03 + SC2 = 1.0370e-03 + SC3 = 0.1417 + W_density = 0.8167 + H_density = 0.2250 iter = 11 - elapsed = 0.015828 - tTotal = 0.2020 - subIterW = 2 - subIterH = 3 - relError = 0.3087 + elapsed = 8.7610e-03 + tTotal = 0.094031 + subIterW = 5 + subIterH = 2 + relError = 0.3416 toshow = scalar structure containing the fields: - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 9.8335e-04 - SC2 = 6.0771e-04 - SC3 = 0.021501 - W_density = 0.8267 - H_density = 0.2800 + initGrNormW = 385.10 + initGrNormH = 490.96 + initNorm = 13.836 + SC1 = 1.9229e-03 + SC2 = 1.0383e-03 + SC3 = 0.1375 + W_density = 0.7467 + H_density = 0.2200 iter = 12 - elapsed = 9.1020e-03 - tTotal = 0.2111 - subIterW = 0 + elapsed = 7.9200e-03 + tTotal = 0.1020 + subIterW = 4 subIterH = 0 - relError = 0.3089 + relError = 0.3438 toshow = scalar structure containing the fields: - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 1.1234e-03 - SC2 = 6.8527e-04 - SC3 = 0.087836 - W_density = 0.8233 - H_density = 0.2650 + initGrNormW = 385.10 + initGrNormH = 490.96 + initNorm = 13.836 + SC1 = 1.9980e-03 + SC2 = 1.0789e-03 + SC3 = 0.1346 + W_density = 0.7533 + H_density = 0.2100 iter = 13 - elapsed = 0.017760 - tTotal = 0.2289 - subIterW = 0 - subIterH = 3 - relError = 0.3094 + elapsed = 9.2120e-03 + tTotal = 0.1112 + subIterW = 4 + subIterH = 2 + relError = 0.3453 toshow = scalar structure containing the fields: - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 1.3607e-03 - SC2 = 8.2188e-04 - SC3 = 0.1075 - W_density = 0.8200 - H_density = 0.2600 + initGrNormW = 385.10 + initGrNormH = 490.96 + initNorm = 13.836 + SC1 = 1.9549e-03 + SC2 = 1.0400e-03 + SC3 = 0.1330 + W_density = 0.7467 + H_density = 0.2050 iter = 14 - elapsed = 0.015764 - tTotal = 0.2447 - subIterW = 0 - subIterH = 3 - relError = 0.3101 + elapsed = 6.9750e-03 + tTotal = 0.1181 + subIterW = 2 + subIterH = 0 + relError = 0.3464 toshow = scalar structure containing the fields: - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 1.5976e-03 - SC2 = 9.6492e-04 - SC3 = 0.1186 - W_density = 0.8133 - H_density = 0.2700 + initGrNormW = 385.10 + initGrNormH = 490.96 + initNorm = 13.836 + SC1 = 1.8365e-03 + SC2 = 9.6600e-04 + SC3 = 0.1326 + W_density = 0.7433 + H_density = 0.2000 iter = 15 - elapsed = 0.015817 - tTotal = 0.2605 - subIterW = 3 - subIterH = 2 - relError = 0.3111 + elapsed = 6.6620e-03 + tTotal = 0.1248 + subIterW = 0 + subIterH = 3 + relError = 0.3474 toshow = scalar structure containing the fields: - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 1.8603e-03 - SC2 = 1.1236e-03 - SC3 = 0.1401 - W_density = 0.8167 - H_density = 0.2700 + initGrNormW = 385.10 + initGrNormH = 490.96 + initNorm = 13.836 + SC1 = 1.8094e-03 + SC2 = 9.3729e-04 + SC3 = 0.1255 + W_density = 0.7300 + H_density = 0.2000 iter = 16 - elapsed = 0.016286 - tTotal = 0.2768 - subIterW = 0 - subIterH = 2 - relError = 0.3129 + elapsed = 5.3440e-03 + tTotal = 0.1301 + subIterW = 3 + subIterH = 0 + relError = 0.3481 toshow = scalar structure containing the fields: - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 2.2035e-03 - SC2 = 1.3089e-03 - SC3 = 0.1544 - W_density = 0.8033 - H_density = 0.2700 + initGrNormW = 385.10 + initGrNormH = 490.96 + initNorm = 13.836 + SC1 = 1.4551e-03 + SC2 = 7.4502e-04 + SC3 = 0.1043 + W_density = 0.7200 + H_density = 0.1950 iter = 17 - elapsed = 0.019672 - tTotal = 0.2964 - subIterW = 4 - subIterH = 2 - relError = 0.3154 + elapsed = 7.0390e-03 + tTotal = 0.1372 + subIterW = 3 + subIterH = 0 + relError = 0.3488 toshow = scalar structure containing the fields: - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 2.3839e-03 - SC2 = 1.3826e-03 - SC3 = 0.1730 - W_density = 0.8000 - H_density = 0.2450 + initGrNormW = 385.10 + initGrNormH = 490.96 + initNorm = 13.836 + SC1 = 1.2529e-03 + SC2 = 6.4398e-04 + SC3 = 0.093264 + W_density = 0.7233 + H_density = 0.2000 iter = 18 - elapsed = 0.014614 - tTotal = 0.3111 - subIterW = 3 - subIterH = 4 - relError = 0.3178 + elapsed = 6.8790e-03 + tTotal = 0.1441 + subIterW = 0 + subIterH = 0 + relError = 0.3495 toshow = scalar structure containing the fields: - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 2.3156e-03 - SC2 = 1.3106e-03 - SC3 = 0.1718 - W_density = 0.7800 - H_density = 0.2450 + initGrNormW = 385.10 + initGrNormH = 490.96 + initNorm = 13.836 + SC1 = 1.1465e-03 + SC2 = 5.8472e-04 + SC3 = 0.079778 + W_density = 0.7167 + H_density = 0.2000 iter = 19 - elapsed = 0.013549 - tTotal = 0.3246 - subIterW = 4 - subIterH = 2 - relError = 0.3203 + elapsed = 5.5330e-03 + tTotal = 0.1496 + subIterW = 2 + subIterH = 0 + relError = 0.3503 toshow = scalar structure containing the fields: - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 2.3524e-03 - SC2 = 1.2985e-03 - SC3 = 0.1621 - W_density = 0.7567 - H_density = 0.2400 + initGrNormW = 385.10 + initGrNormH = 490.96 + initNorm = 13.836 + SC1 = 9.8845e-04 + SC2 = 4.9620e-04 + SC3 = 0.070255 + W_density = 0.7100 + H_density = 0.1900 iter = 20 - elapsed = 0.014275 - tTotal = 0.3389 - subIterW = 7 - subIterH = 0 - relError = 0.3234 + elapsed = 7.1990e-03 + tTotal = 0.1568 + subIterW = 2 + subIterH = 2 + relError = 0.3510 toshow = scalar structure containing the fields: - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 2.2031e-03 - SC2 = 1.1985e-03 - SC3 = 0.1479 - W_density = 0.7467 - H_density = 0.2350 + initGrNormW = 385.10 + initGrNormH = 490.96 + initNorm = 13.836 + SC1 = 9.8118e-04 + SC2 = 4.8667e-04 + SC3 = 0.065210 + W_density = 0.7000 + H_density = 0.1900 iter = 21 - elapsed = 0.012519 - tTotal = 0.3514 + elapsed = 7.4030e-03 + tTotal = 0.1642 subIterW = 3 - subIterH = 0 - relError = 0.3269 + subIterH = 2 + relError = 0.3516 toshow = scalar structure containing the fields: - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 1.8442e-03 - SC2 = 9.7004e-04 - SC3 = 0.1248 - W_density = 0.7233 - H_density = 0.2300 + initGrNormW = 385.10 + initGrNormH = 490.96 + initNorm = 13.836 + SC1 = 1.0457e-03 + SC2 = 5.2284e-04 + SC3 = 0.060314 + W_density = 0.7067 + H_density = 0.1900 iter = 22 - elapsed = 0.011663 - tTotal = 0.3631 - subIterW = 6 + elapsed = 6.2840e-03 + tTotal = 0.1705 + subIterW = 2 subIterH = 0 - relError = 0.3303 + relError = 0.3521 toshow = scalar structure containing the fields: - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 1.4107e-03 - SC2 = 7.1663e-04 - SC3 = 0.1030 - W_density = 0.7000 - H_density = 0.2200 + initGrNormW = 385.10 + initGrNormH = 490.96 + initNorm = 13.836 + SC1 = 8.9661e-04 + SC2 = 4.4651e-04 + SC3 = 0.013938 + W_density = 0.7067 + H_density = 0.1850 iter = 23 - elapsed = 0.013951 - tTotal = 0.3770 - subIterW = 8 - subIterH = 2 - relError = 0.3330 - -toshow = - - scalar structure containing the fields: - - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 1.1058e-03 - SC2 = 5.5510e-04 - SC3 = 0.084950 - W_density = 0.6933 - H_density = 0.2150 - iter = 24 - elapsed = 0.011290 - tTotal = 0.3883 - subIterW = 2 + elapsed = 5.9350e-03 + tTotal = 0.1764 + subIterW = 0 subIterH = 0 - relError = 0.3352 + relError = 0.3522 toshow = scalar structure containing the fields: - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 9.1528e-04 - SC2 = 4.5764e-04 - SC3 = 0.072112 - W_density = 0.6900 - H_density = 0.2150 - iter = 25 - elapsed = 9.1150e-03 - tTotal = 0.3974 + initGrNormW = 385.10 + initGrNormH = 490.96 + initNorm = 13.836 + SC1 = 7.7026e-04 + SC2 = 3.8051e-04 + SC3 = 0.050780 + W_density = 0.7000 + H_density = 0.1850 + iter = 24 + elapsed = 8.1850e-03 + tTotal = 0.1846 subIterW = 3 - subIterH = 0 - relError = 0.3369 + subIterH = 2 + relError = 0.3521 toshow = scalar structure containing the fields: - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 8.0567e-04 - SC2 = 4.0122e-04 - SC3 = 0.063758 - W_density = 0.6867 - H_density = 0.2150 - iter = 26 - elapsed = 8.0090e-03 - tTotal = 0.4054 + initGrNormW = 385.10 + initGrNormH = 490.96 + initNorm = 13.836 + SC1 = 6.5544e-04 + SC2 = 3.2117e-04 + SC3 = 0.037267 + W_density = 0.6967 + H_density = 0.1800 + iter = 25 + elapsed = 6.4930e-03 + tTotal = 0.1911 subIterW = 0 subIterH = 0 - relError = 0.3384 - -toshow = - - scalar structure containing the fields: - - initGrNormW = 401.81 - initGrNormH = 526.81 - initNorm = 14.263 - SC1 = 7.8132e-04 - SC2 = 3.9222e-04 - SC3 = 0.059675 - W_density = 0.6933 - H_density = 0.2150 - iter = 27 - elapsed = 0.010931 - tTotal = 0.4164 - subIterW = 2 - subIterH = 0 - relError = 0.3398 + relError = 0.3518 Stop: tolerance reached. final = scalar structure containing the fields: - iterations = 27 - elapsed_total = 0.4164 - relative_error = 0.3398 - W_density = 0.6933 - H_density = 0.2150 + iterations = 25 + elapsed_total = 0.1911 + relative_error = 0.3518 + W_density = 0.6967 + H_density = 0.1800 ***** test [W,H,iter,HIS]=nmf_bpas(A,k,'verbose','type','sparse','nnlssolver','bp','alpha',1.1,'beta',1.3); - initGrNormW = 386.35 - initGrNormH = 616.44 - initNorm = 14.263 - SC1 = 1.4550 - SC2 = 727.51 - SC3 = 0.5231 + initGrNormW = 383.96 + initGrNormH = 564.73 + initNorm = 13.836 + SC1 = 1.3658 + SC2 = 682.89 + SC3 = 0.5015 W_density = 1 H_density = 1 par = @@ -4456,438 +4486,438 @@ scalar structure containing the fields: - initGrNormW = 386.35 - initGrNormH = 616.44 - initNorm = 14.263 - SC1 = 0.028736 - SC2 = 0.026322 - SC3 = 2.5231 - W_density = 0.9733 - H_density = 0.5500 + initGrNormW = 383.96 + initGrNormH = 564.73 + initNorm = 13.836 + SC1 = 0.026485 + SC2 = 0.024843 + SC3 = 2.2176 + W_density = 0.9900 + H_density = 0.5600 iter = 1 - elapsed = 0.021717 - tTotal = 0.021717 + elapsed = 0.011172 + tTotal = 0.011172 subIterW = 2 - subIterH = 38 - relError = 0.5684 + subIterH = 36 + relError = 0.5751 toshow = scalar structure containing the fields: - initGrNormW = 386.35 - initGrNormH = 616.44 - initNorm = 14.263 - SC1 = 0.016113 - SC2 = 7.5088e-03 - SC3 = 0.8738 - W_density = 0.6067 - H_density = 0.2150 + initGrNormW = 383.96 + initGrNormH = 564.73 + initNorm = 13.836 + SC1 = 0.013658 + SC2 = 7.1569e-03 + SC3 = 0.8672 + W_density = 0.6600 + H_density = 0.2350 iter = 2 - elapsed = 0.025134 - tTotal = 0.046851 - subIterW = 13 - subIterH = 29 - relError = 0.4538 + elapsed = 0.012026 + tTotal = 0.023198 + subIterW = 4 + subIterH = 36 + relError = 0.4623 toshow = scalar structure containing the fields: - initGrNormW = 386.35 - initGrNormH = 616.44 - initNorm = 14.263 - SC1 = 0.011856 - SC2 = 3.7229e-03 - SC3 = 0.5708 - W_density = 0.4000 - H_density = 0.1700 + initGrNormW = 383.96 + initGrNormH = 564.73 + initNorm = 13.836 + SC1 = 0.013979 + SC2 = 5.1724e-03 + SC3 = 0.7567 + W_density = 0.4700 + H_density = 0.1850 iter = 3 - elapsed = 0.021489 - tTotal = 0.068340 - subIterW = 14 - subIterH = 7 - relError = 0.4532 + elapsed = 9.2500e-03 + tTotal = 0.032448 + subIterW = 4 + subIterH = 10 + relError = 0.4587 toshow = scalar structure containing the fields: - initGrNormW = 386.35 - initGrNormH = 616.44 - initNorm = 14.263 - SC1 = 0.011357 - SC2 = 2.4985e-03 - SC3 = 0.4807 - W_density = 0.2633 - H_density = 0.1500 + initGrNormW = 383.96 + initGrNormH = 564.73 + initNorm = 13.836 + SC1 = 0.012619 + SC2 = 3.9625e-03 + SC3 = 0.5131 + W_density = 0.4033 + H_density = 0.1700 iter = 4 - elapsed = 0.015831 - tTotal = 0.084171 - subIterW = 3 - subIterH = 2 - relError = 0.4580 + elapsed = 6.6440e-03 + tTotal = 0.039092 + subIterW = 6 + subIterH = 3 + relError = 0.4653 toshow = scalar structure containing the fields: - initGrNormW = 386.35 - initGrNormH = 616.44 - initNorm = 14.263 - SC1 = 0.014592 - SC2 = 3.0935e-03 - SC3 = 0.4621 - W_density = 0.2500 - H_density = 0.1550 + initGrNormW = 383.96 + initGrNormH = 564.73 + initNorm = 13.836 + SC1 = 0.012095 + SC2 = 2.6125e-03 + SC3 = 0.4375 + W_density = 0.2600 + H_density = 0.1500 iter = 5 - elapsed = 0.010166 - tTotal = 0.094337 - subIterW = 2 - subIterH = 0 - relError = 0.4632 + elapsed = 6.7610e-03 + tTotal = 0.045853 + subIterW = 3 + subIterH = 5 + relError = 0.4689 toshow = scalar structure containing the fields: - initGrNormW = 386.35 - initGrNormH = 616.44 - initNorm = 14.263 - SC1 = 0.018984 - SC2 = 3.1513e-03 - SC3 = 0.4182 - W_density = 0.1867 + initGrNormW = 383.96 + initGrNormH = 564.73 + initNorm = 13.836 + SC1 = 0.014492 + SC2 = 2.4637e-03 + SC3 = 0.4617 + W_density = 0.1933 H_density = 0.1350 iter = 6 - elapsed = 0.010595 - tTotal = 0.1049 - subIterW = 2 - subIterH = 1 - relError = 0.4665 + elapsed = 5.6310e-03 + tTotal = 0.051484 + subIterW = 0 + subIterH = 2 + relError = 0.4714 toshow = scalar structure containing the fields: - initGrNormW = 386.35 - initGrNormH = 616.44 - initNorm = 14.263 - SC1 = 0.012132 - SC2 = 1.7956e-03 - SC3 = 0.2347 - W_density = 0.1733 - H_density = 0.1100 + initGrNormW = 383.96 + initGrNormH = 564.73 + initNorm = 13.836 + SC1 = 0.014145 + SC2 = 2.3480e-03 + SC3 = 0.4728 + W_density = 0.1900 + H_density = 0.1300 iter = 7 - elapsed = 8.9800e-03 - tTotal = 0.1139 - subIterW = 1 - subIterH = 1 - relError = 0.4687 + elapsed = 4.6260e-03 + tTotal = 0.056110 + subIterW = 0 + subIterH = 0 + relError = 0.4750 toshow = scalar structure containing the fields: - initGrNormW = 386.35 - initGrNormH = 616.44 - initNorm = 14.263 - SC1 = 7.2035e-03 - SC2 = 1.0229e-03 - SC3 = 0.1589 - W_density = 0.1667 - H_density = 0.1050 + initGrNormW = 383.96 + initGrNormH = 564.73 + initNorm = 13.836 + SC1 = 0.014859 + SC2 = 2.2586e-03 + SC3 = 0.3310 + W_density = 0.1700 + H_density = 0.1250 iter = 8 - elapsed = 9.2990e-03 - tTotal = 0.1232 + elapsed = 5.1060e-03 + tTotal = 0.061216 subIterW = 1 subIterH = 0 - relError = 0.4724 + relError = 0.4811 toshow = scalar structure containing the fields: - initGrNormW = 386.35 - initGrNormH = 616.44 - initNorm = 14.263 - SC1 = 8.0445e-03 - SC2 = 1.0941e-03 - SC3 = 0.063811 - W_density = 0.1567 + initGrNormW = 383.96 + initGrNormH = 564.73 + initNorm = 13.836 + SC1 = 9.6141e-03 + SC2 = 1.3460e-03 + SC3 = 0.1334 + W_density = 0.1633 H_density = 0.1050 iter = 9 - elapsed = 6.7430e-03 - tTotal = 0.1300 - subIterW = 1 - subIterH = 0 - relError = 0.4761 + elapsed = 5.0450e-03 + tTotal = 0.066261 + subIterW = 2 + subIterH = 1 + relError = 0.4861 toshow = scalar structure containing the fields: - initGrNormW = 386.35 - initGrNormH = 616.44 - initNorm = 14.263 - SC1 = 2.8805e-04 - SC2 = 2.8805e-05 - SC3 = 4.7772e-03 + initGrNormW = 383.96 + initGrNormH = 564.73 + initNorm = 13.836 + SC1 = 7.5333e-04 + SC2 = 7.5333e-05 + SC3 = 0.013358 W_density = 0.1000 H_density = 0.1000 iter = 10 - elapsed = 5.5970e-03 - tTotal = 0.1356 + elapsed = 3.4090e-03 + tTotal = 0.069670 subIterW = 0 subIterH = 0 - relError = 0.4764 + relError = 0.4866 toshow = scalar structure containing the fields: - initGrNormW = 386.35 - initGrNormH = 616.44 - initNorm = 14.263 - SC1 = 1.2506e-04 - SC2 = 1.2506e-05 - SC3 = 2.5737e-03 + initGrNormW = 383.96 + initGrNormH = 564.73 + initNorm = 13.836 + SC1 = 3.6252e-04 + SC2 = 3.6252e-05 + SC3 = 0.010980 W_density = 0.1000 H_density = 0.1000 iter = 11 - elapsed = 2.8300e-03 - tTotal = 0.1384 + elapsed = 1.7280e-03 + tTotal = 0.071398 subIterW = 0 subIterH = 0 - relError = 0.4764 + relError = 0.4866 toshow = scalar structure containing the fields: - initGrNormW = 386.35 - initGrNormH = 616.44 - initNorm = 14.263 - SC1 = 8.1929e-05 - SC2 = 8.1929e-06 - SC3 = 1.8100e-03 + initGrNormW = 383.96 + initGrNormH = 564.73 + initNorm = 13.836 + SC1 = 2.3344e-04 + SC2 = 2.3344e-05 + SC3 = 5.2433e-03 W_density = 0.1000 H_density = 0.1000 iter = 12 - elapsed = 2.6070e-03 - tTotal = 0.1410 + elapsed = 1.7470e-03 + tTotal = 0.073145 subIterW = 0 subIterH = 0 - relError = 0.4763 + relError = 0.4866 toshow = scalar structure containing the fields: - initGrNormW = 386.35 - initGrNormH = 616.44 - initNorm = 14.263 - SC1 = 5.3768e-05 - SC2 = 5.3768e-06 - SC3 = 1.2319e-03 + initGrNormW = 383.96 + initGrNormH = 564.73 + initNorm = 13.836 + SC1 = 1.5052e-04 + SC2 = 1.5052e-05 + SC3 = 3.3808e-03 W_density = 0.1000 H_density = 0.1000 iter = 13 - elapsed = 3.1210e-03 - tTotal = 0.1441 + elapsed = 1.6920e-03 + tTotal = 0.074837 subIterW = 0 subIterH = 0 - relError = 0.4763 + relError = 0.4866 toshow = scalar structure containing the fields: - initGrNormW = 386.35 - initGrNormH = 616.44 - initNorm = 14.263 - SC1 = 3.5289e-05 - SC2 = 3.5289e-06 - SC3 = 6.6149e-04 + initGrNormW = 383.96 + initGrNormH = 564.73 + initNorm = 13.836 + SC1 = 9.7062e-05 + SC2 = 9.7062e-06 + SC3 = 2.4526e-03 W_density = 0.1000 H_density = 0.1000 iter = 14 - elapsed = 2.2600e-03 - tTotal = 0.1464 + elapsed = 1.5950e-03 + tTotal = 0.076432 subIterW = 0 subIterH = 0 - relError = 0.4763 + relError = 0.4866 toshow = scalar structure containing the fields: - initGrNormW = 386.35 - initGrNormH = 616.44 - initNorm = 14.263 - SC1 = 2.3161e-05 - SC2 = 2.3161e-06 - SC3 = 5.7308e-04 + initGrNormW = 383.96 + initGrNormH = 564.73 + initNorm = 13.836 + SC1 = 6.2595e-05 + SC2 = 6.2595e-06 + SC3 = 1.5184e-03 W_density = 0.1000 H_density = 0.1000 iter = 15 - elapsed = 3.5760e-03 - tTotal = 0.1499 + elapsed = 1.6200e-03 + tTotal = 0.078052 subIterW = 0 subIterH = 0 - relError = 0.4763 + relError = 0.4866 toshow = scalar structure containing the fields: - initGrNormW = 386.35 - initGrNormH = 616.44 - initNorm = 14.263 - SC1 = 1.5201e-05 - SC2 = 1.5201e-06 - SC3 = 4.0884e-04 + initGrNormW = 383.96 + initGrNormH = 564.73 + initNorm = 13.836 + SC1 = 4.0369e-05 + SC2 = 4.0369e-06 + SC3 = 1.0201e-03 W_density = 0.1000 H_density = 0.1000 iter = 16 - elapsed = 2.7970e-03 - tTotal = 0.1527 + elapsed = 2.0580e-03 + tTotal = 0.080110 subIterW = 0 subIterH = 0 - relError = 0.4763 + relError = 0.4866 toshow = scalar structure containing the fields: - initGrNormW = 386.35 - initGrNormH = 616.44 - initNorm = 14.263 - SC1 = 9.9772e-06 - SC2 = 9.9772e-07 - SC3 = 2.3738e-04 + initGrNormW = 383.96 + initGrNormH = 564.73 + initNorm = 13.836 + SC1 = 2.6036e-05 + SC2 = 2.6036e-06 + SC3 = 6.3158e-04 W_density = 0.1000 H_density = 0.1000 iter = 17 - elapsed = 3.0400e-03 - tTotal = 0.1558 + elapsed = 1.7660e-03 + tTotal = 0.081876 subIterW = 0 subIterH = 0 - relError = 0.4763 + relError = 0.4866 toshow = scalar structure containing the fields: - initGrNormW = 386.35 - initGrNormH = 616.44 - initNorm = 14.263 - SC1 = 6.5485e-06 - SC2 = 6.5485e-07 - SC3 = 1.6879e-04 + initGrNormW = 383.96 + initGrNormH = 564.73 + initNorm = 13.836 + SC1 = 1.6792e-05 + SC2 = 1.6792e-06 + SC3 = 3.9167e-04 W_density = 0.1000 H_density = 0.1000 iter = 18 - elapsed = 2.6890e-03 - tTotal = 0.1585 + elapsed = 1.7700e-03 + tTotal = 0.083646 subIterW = 0 subIterH = 0 - relError = 0.4763 + relError = 0.4866 toshow = scalar structure containing the fields: - initGrNormW = 386.35 - initGrNormH = 616.44 - initNorm = 14.263 - SC1 = 4.2981e-06 - SC2 = 4.2981e-07 - SC3 = 1.0226e-04 + initGrNormW = 383.96 + initGrNormH = 564.73 + initNorm = 13.836 + SC1 = 1.0830e-05 + SC2 = 1.0830e-06 + SC3 = 2.5262e-04 W_density = 0.1000 H_density = 0.1000 iter = 19 - elapsed = 3.3430e-03 - tTotal = 0.1618 + elapsed = 1.6760e-03 + tTotal = 0.085322 subIterW = 0 subIterH = 0 - relError = 0.4763 + relError = 0.4866 toshow = scalar structure containing the fields: - initGrNormW = 386.35 - initGrNormH = 616.44 - initNorm = 14.263 - SC1 = 2.8211e-06 - SC2 = 2.8211e-07 - SC3 = 7.2712e-05 + initGrNormW = 383.96 + initGrNormH = 564.73 + initNorm = 13.836 + SC1 = 6.9853e-06 + SC2 = 6.9853e-07 + SC3 = 1.4608e-04 W_density = 0.1000 H_density = 0.1000 iter = 20 - elapsed = 3.5430e-03 - tTotal = 0.1654 + elapsed = 1.6740e-03 + tTotal = 0.086996 subIterW = 0 subIterH = 0 - relError = 0.4763 + relError = 0.4866 toshow = scalar structure containing the fields: - initGrNormW = 386.35 - initGrNormH = 616.44 - initNorm = 14.263 - SC1 = 1.8516e-06 - SC2 = 1.8516e-07 - SC3 = 4.5816e-05 + initGrNormW = 383.96 + initGrNormH = 564.73 + initNorm = 13.836 + SC1 = 4.5053e-06 + SC2 = 4.5053e-07 + SC3 = 1.0509e-04 W_density = 0.1000 H_density = 0.1000 iter = 21 - elapsed = 3.2290e-03 - tTotal = 0.1686 + elapsed = 1.7950e-03 + tTotal = 0.088791 subIterW = 0 subIterH = 0 - relError = 0.4763 + relError = 0.4866 toshow = scalar structure containing the fields: - initGrNormW = 386.35 - initGrNormH = 616.44 - initNorm = 14.263 - SC1 = 1.2153e-06 - SC2 = 1.2153e-07 - SC3 = 2.6850e-05 + initGrNormW = 383.96 + initGrNormH = 564.73 + initNorm = 13.836 + SC1 = 2.9058e-06 + SC2 = 2.9058e-07 + SC3 = 7.0489e-05 W_density = 0.1000 H_density = 0.1000 iter = 22 - elapsed = 3.4660e-03 - tTotal = 0.1721 + elapsed = 1.9110e-03 + tTotal = 0.090702 subIterW = 0 subIterH = 0 - relError = 0.4763 + relError = 0.4865 toshow = scalar structure containing the fields: - initGrNormW = 386.35 - initGrNormH = 616.44 - initNorm = 14.263 - SC1 = 7.9769e-07 - SC2 = 7.9769e-08 - SC3 = 1.6448e-05 + initGrNormW = 383.96 + initGrNormH = 564.73 + initNorm = 13.836 + SC1 = 1.8742e-06 + SC2 = 1.8742e-07 + SC3 = 4.2096e-05 W_density = 0.1000 H_density = 0.1000 iter = 23 - elapsed = 3.5190e-03 - tTotal = 0.1756 + elapsed = 1.8950e-03 + tTotal = 0.092597 subIterW = 0 subIterH = 0 - relError = 0.4763 + relError = 0.4865 Stop: tolerance reached. final = @@ -4895,19 +4925,19 @@ scalar structure containing the fields: iterations = 23 - elapsed_total = 0.1756 - relative_error = 0.4763 + elapsed_total = 0.092597 + relative_error = 0.4865 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.43 - initGrNormH = 433.24 - initNorm = 14.263 - SC1 = 1.1284 - SC2 = 564.21 - SC3 = 0.4974 + initGrNormW = 378.98 + initGrNormH = 458.15 + initNorm = 13.836 + SC1 = 1.1892 + SC2 = 594.58 + SC3 = 0.5034 W_density = 1 H_density = 1 par = @@ -4918,8 +4948,8 @@ n = 20 type = plain nnls_solver = bp - alpha = 0.5082 - beta = 0.5082 + alpha = 0.4893 + beta = 0.4893 max_iter = 100 min_iter = 20 max_time = 1000 @@ -4931,438 +4961,438 @@ scalar structure containing the fields: - initGrNormW = 361.43 - initGrNormH = 433.24 - initNorm = 14.263 - SC1 = 0.016091 - SC2 = 0.011682 - SC3 = 0.5270 + initGrNormW = 378.98 + initGrNormH = 458.15 + initNorm = 13.836 + SC1 = 0.018134 + SC2 = 0.013238 + SC3 = 0.5986 W_density = 0.7167 - H_density = 0.6300 + H_density = 0.6550 iter = 1 - elapsed = 0.036738 - tTotal = 0.036738 - subIterW = 52 - subIterH = 32 - relError = 0.3251 + elapsed = 0.020257 + tTotal = 0.020257 + subIterW = 42 + subIterH = 31 + relError = 0.3382 toshow = scalar structure containing the fields: - initGrNormW = 361.43 - initGrNormH = 433.24 - initNorm = 14.263 - SC1 = 9.0591e-03 - SC2 = 6.5226e-03 - SC3 = 0.2980 - W_density = 0.7167 - H_density = 0.6750 + initGrNormW = 378.98 + initGrNormH = 458.15 + initNorm = 13.836 + SC1 = 0.010653 + SC2 = 7.7342e-03 + SC3 = 0.3621 + W_density = 0.7400 + H_density = 0.6550 iter = 2 - elapsed = 0.030018 - tTotal = 0.066756 - subIterW = 23 - subIterH = 14 - relError = 0.2805 + elapsed = 0.013497 + tTotal = 0.033754 + subIterW = 21 + subIterH = 18 + relError = 0.2903 toshow = scalar structure containing the fields: - initGrNormW = 361.43 - initGrNormH = 433.24 - initNorm = 14.263 - SC1 = 5.4603e-03 - SC2 = 3.8768e-03 - SC3 = 0.2033 + initGrNormW = 378.98 + initGrNormH = 458.15 + initNorm = 13.836 + SC1 = 7.3539e-03 + SC2 = 5.2507e-03 + SC3 = 0.2541 W_density = 0.7267 H_density = 0.6650 iter = 3 - elapsed = 0.021007 - tTotal = 0.087763 - subIterW = 11 - subIterH = 14 - relError = 0.2630 + elapsed = 0.014724 + tTotal = 0.048478 + subIterW = 18 + subIterH = 13 + relError = 0.2644 toshow = scalar structure containing the fields: - initGrNormW = 361.43 - initGrNormH = 433.24 - initNorm = 14.263 - SC1 = 3.7526e-03 - SC2 = 2.6193e-03 - SC3 = 0.1649 - W_density = 0.7200 - H_density = 0.6450 + initGrNormW = 378.98 + initGrNormH = 458.15 + initNorm = 13.836 + SC1 = 4.6139e-03 + SC2 = 3.2851e-03 + SC3 = 0.1675 + W_density = 0.7300 + H_density = 0.6700 iter = 4 - elapsed = 0.022724 - tTotal = 0.1105 - subIterW = 11 - subIterH = 10 - relError = 0.2558 + elapsed = 0.011543 + tTotal = 0.060021 + subIterW = 10 + subIterH = 11 + relError = 0.2546 toshow = scalar structure containing the fields: - initGrNormW = 361.43 - initGrNormH = 433.24 - initNorm = 14.263 - SC1 = 3.0150e-03 - SC2 = 2.1105e-03 - SC3 = 0.1335 - W_density = 0.7233 - H_density = 0.6650 + initGrNormW = 378.98 + initGrNormH = 458.15 + initNorm = 13.836 + SC1 = 3.6290e-03 + SC2 = 2.5766e-03 + SC3 = 0.1375 + W_density = 0.7333 + H_density = 0.6600 iter = 5 - elapsed = 0.017478 - tTotal = 0.1280 - subIterW = 7 - subIterH = 8 - relError = 0.2518 + elapsed = 0.010780 + tTotal = 0.070801 + subIterW = 11 + subIterH = 4 + relError = 0.2502 toshow = scalar structure containing the fields: - initGrNormW = 361.43 - initGrNormH = 433.24 - initNorm = 14.263 - SC1 = 2.6780e-03 - SC2 = 1.8532e-03 - SC3 = 0.1183 - W_density = 0.7133 - H_density = 0.6600 + initGrNormW = 378.98 + initGrNormH = 458.15 + initNorm = 13.836 + SC1 = 3.1183e-03 + SC2 = 2.2078e-03 + SC3 = 0.1267 + W_density = 0.7333 + H_density = 0.6650 iter = 6 - elapsed = 0.015269 - tTotal = 0.1432 - subIterW = 3 - subIterH = 3 - relError = 0.2489 + elapsed = 7.9210e-03 + tTotal = 0.078722 + subIterW = 6 + subIterH = 5 + relError = 0.2474 toshow = scalar structure containing the fields: - initGrNormW = 361.43 - initGrNormH = 433.24 - initNorm = 14.263 - SC1 = 2.5600e-03 - SC2 = 1.7767e-03 - SC3 = 0.1155 + initGrNormW = 378.98 + initGrNormH = 458.15 + initNorm = 13.836 + SC1 = 2.6390e-03 + SC2 = 1.8525e-03 + SC3 = 0.1084 W_density = 0.7300 - H_density = 0.6300 + H_density = 0.6550 iter = 7 - elapsed = 0.020244 - tTotal = 0.1635 - subIterW = 10 - subIterH = 5 - relError = 0.2464 + elapsed = 8.6210e-03 + tTotal = 0.087343 + subIterW = 5 + subIterH = 4 + relError = 0.2453 toshow = scalar structure containing the fields: - initGrNormW = 361.43 - initGrNormH = 433.24 - initNorm = 14.263 - SC1 = 2.5656e-03 - SC2 = 1.7959e-03 - SC3 = 0.1140 - W_density = 0.7400 - H_density = 0.6250 + initGrNormW = 378.98 + initGrNormH = 458.15 + initNorm = 13.836 + SC1 = 2.3875e-03 + SC2 = 1.6426e-03 + SC3 = 0.096667 + W_density = 0.7167 + H_density = 0.6450 iter = 8 - elapsed = 0.016681 - tTotal = 0.1802 - subIterW = 3 - subIterH = 3 - relError = 0.2441 + elapsed = 8.9720e-03 + tTotal = 0.096315 + subIterW = 5 + subIterH = 7 + relError = 0.2438 toshow = scalar structure containing the fields: - initGrNormW = 361.43 - initGrNormH = 433.24 - initNorm = 14.263 - SC1 = 2.3634e-03 - SC2 = 1.6497e-03 - SC3 = 0.099923 - W_density = 0.7433 - H_density = 0.6200 + initGrNormW = 378.98 + initGrNormH = 458.15 + initNorm = 13.836 + SC1 = 2.3085e-03 + SC2 = 1.5744e-03 + SC3 = 0.092775 + W_density = 0.7133 + H_density = 0.6350 iter = 9 - elapsed = 0.023447 - tTotal = 0.2036 - subIterW = 9 - subIterH = 7 - relError = 0.2419 + elapsed = 8.9550e-03 + tTotal = 0.1053 + subIterW = 5 + subIterH = 2 + relError = 0.2423 toshow = scalar structure containing the fields: - initGrNormW = 361.43 - initGrNormH = 433.24 - initNorm = 14.263 - SC1 = 2.1062e-03 - SC2 = 1.4533e-03 - SC3 = 0.090869 - W_density = 0.7433 - H_density = 0.6050 + initGrNormW = 378.98 + initGrNormH = 458.15 + initNorm = 13.836 + SC1 = 1.8683e-03 + SC2 = 1.2629e-03 + SC3 = 0.082135 + W_density = 0.7133 + H_density = 0.6150 iter = 10 - elapsed = 0.015660 - tTotal = 0.2193 + elapsed = 9.2030e-03 + tTotal = 0.1145 subIterW = 2 - subIterH = 3 - relError = 0.2400 + subIterH = 4 + relError = 0.2411 toshow = scalar structure containing the fields: - initGrNormW = 361.43 - initGrNormH = 433.24 - initNorm = 14.263 - SC1 = 1.8466e-03 - SC2 = 1.2852e-03 - SC3 = 0.082031 - W_density = 0.7567 - H_density = 0.5950 + initGrNormW = 378.98 + initGrNormH = 458.15 + initNorm = 13.836 + SC1 = 1.5891e-03 + SC2 = 1.0679e-03 + SC3 = 0.076138 + W_density = 0.7133 + H_density = 0.6050 iter = 11 - elapsed = 0.013929 - tTotal = 0.2332 - subIterW = 5 - subIterH = 3 - relError = 0.2384 + elapsed = 8.2660e-03 + tTotal = 0.1227 + subIterW = 0 + subIterH = 4 + relError = 0.2402 toshow = scalar structure containing the fields: - initGrNormW = 361.43 - initGrNormH = 433.24 - initNorm = 14.263 - SC1 = 1.5918e-03 - SC2 = 1.0983e-03 - SC3 = 0.073032 - W_density = 0.7567 - H_density = 0.5900 + initGrNormW = 378.98 + initGrNormH = 458.15 + initNorm = 13.836 + SC1 = 1.4316e-03 + SC2 = 9.6490e-04 + SC3 = 0.071293 + W_density = 0.7200 + H_density = 0.6050 iter = 12 - elapsed = 0.012826 - tTotal = 0.2460 - subIterW = 5 - subIterH = 4 - relError = 0.2371 + elapsed = 9.5410e-03 + tTotal = 0.1323 + subIterW = 4 + subIterH = 3 + relError = 0.2396 toshow = scalar structure containing the fields: - initGrNormW = 361.43 - initGrNormH = 433.24 - initNorm = 14.263 - SC1 = 1.3345e-03 - SC2 = 9.1813e-04 - SC3 = 0.062592 - W_density = 0.7600 - H_density = 0.5800 + initGrNormW = 378.98 + initGrNormH = 458.15 + initNorm = 13.836 + SC1 = 1.3242e-03 + SC2 = 8.8984e-04 + SC3 = 0.066783 + W_density = 0.7167 + H_density = 0.6000 iter = 13 - elapsed = 0.011731 - tTotal = 0.2578 - subIterW = 2 - subIterH = 2 - relError = 0.2362 + elapsed = 8.5820e-03 + tTotal = 0.1409 + subIterW = 3 + subIterH = 0 + relError = 0.2390 toshow = scalar structure containing the fields: - initGrNormW = 361.43 - initGrNormH = 433.24 - initNorm = 14.263 - SC1 = 1.1473e-03 - SC2 = 7.8703e-04 - SC3 = 0.053095 - W_density = 0.7600 - H_density = 0.5700 + initGrNormW = 378.98 + initGrNormH = 458.15 + initNorm = 13.836 + SC1 = 1.2017e-03 + SC2 = 8.0514e-04 + SC3 = 0.061984 + W_density = 0.7200 + H_density = 0.5900 iter = 14 - elapsed = 0.012658 - tTotal = 0.2704 + elapsed = 9.3530e-03 + tTotal = 0.1502 subIterW = 4 - subIterH = 4 - relError = 0.2356 + subIterH = 2 + relError = 0.2386 toshow = scalar structure containing the fields: - initGrNormW = 361.43 - initGrNormH = 433.24 - initNorm = 14.263 - SC1 = 1.0009e-03 - SC2 = 6.8863e-04 - SC3 = 0.045090 - W_density = 0.7633 - H_density = 0.5750 + initGrNormW = 378.98 + initGrNormH = 458.15 + initNorm = 13.836 + SC1 = 1.1127e-03 + SC2 = 7.3660e-04 + SC3 = 0.058805 + W_density = 0.7167 + H_density = 0.5800 iter = 15 - elapsed = 0.011780 - tTotal = 0.2822 + elapsed = 6.6260e-03 + tTotal = 0.1568 subIterW = 0 subIterH = 3 - relError = 0.2352 + relError = 0.2382 toshow = scalar structure containing the fields: - initGrNormW = 361.43 - initGrNormH = 433.24 - initNorm = 14.263 - SC1 = 8.8793e-04 - SC2 = 6.1267e-04 - SC3 = 0.041337 - W_density = 0.7733 - H_density = 0.5650 + initGrNormW = 378.98 + initGrNormH = 458.15 + initNorm = 13.836 + SC1 = 1.0356e-03 + SC2 = 6.8144e-04 + SC3 = 0.054248 + W_density = 0.7100 + H_density = 0.5800 iter = 16 - elapsed = 0.014352 - tTotal = 0.2965 - subIterW = 3 - subIterH = 2 - relError = 0.2348 + elapsed = 5.3310e-03 + tTotal = 0.1622 + subIterW = 2 + subIterH = 0 + relError = 0.2379 toshow = scalar structure containing the fields: - initGrNormW = 361.43 - initGrNormH = 433.24 - initNorm = 14.263 - SC1 = 7.9257e-04 - SC2 = 5.4370e-04 - SC3 = 0.039218 - W_density = 0.7767 - H_density = 0.5400 + initGrNormW = 378.98 + initGrNormH = 458.15 + initNorm = 13.836 + SC1 = 9.4530e-04 + SC2 = 6.2957e-04 + SC3 = 0.051080 + W_density = 0.7233 + H_density = 0.5800 iter = 17 - elapsed = 0.011660 - tTotal = 0.3082 - subIterW = 0 - subIterH = 5 - relError = 0.2346 + elapsed = 6.3130e-03 + tTotal = 0.1685 + subIterW = 4 + subIterH = 0 + relError = 0.2376 toshow = scalar structure containing the fields: - initGrNormW = 361.43 - initGrNormH = 433.24 - initNorm = 14.263 - SC1 = 7.3429e-04 - SC2 = 4.9932e-04 - SC3 = 0.037396 - W_density = 0.7733 - H_density = 0.5400 + initGrNormW = 378.98 + initGrNormH = 458.15 + initNorm = 13.836 + SC1 = 8.9198e-04 + SC2 = 5.9049e-04 + SC3 = 0.048985 + W_density = 0.7200 + H_density = 0.5750 iter = 18 - elapsed = 0.011641 - tTotal = 0.3198 + elapsed = 6.8030e-03 + tTotal = 0.1753 subIterW = 0 - subIterH = 3 - relError = 0.2343 + subIterH = 0 + relError = 0.2374 toshow = scalar structure containing the fields: - initGrNormW = 361.43 - initGrNormH = 433.24 - initNorm = 14.263 - SC1 = 6.6541e-04 - SC2 = 4.5248e-04 - SC3 = 0.033441 - W_density = 0.7733 - H_density = 0.5400 + initGrNormW = 378.98 + initGrNormH = 458.15 + initNorm = 13.836 + SC1 = 8.4884e-04 + SC2 = 5.6193e-04 + SC3 = 0.046652 + W_density = 0.7167 + H_density = 0.5750 iter = 19 - elapsed = 0.011850 - tTotal = 0.3317 - subIterW = 2 - subIterH = 2 - relError = 0.2341 + elapsed = 5.2810e-03 + tTotal = 0.1806 + subIterW = 3 + subIterH = 0 + relError = 0.2372 toshow = scalar structure containing the fields: - initGrNormW = 361.43 - initGrNormH = 433.24 - initNorm = 14.263 - SC1 = 6.2073e-04 - SC2 = 4.2458e-04 - SC3 = 0.030592 - W_density = 0.7767 - H_density = 0.5400 + initGrNormW = 378.98 + initGrNormH = 458.15 + initNorm = 13.836 + SC1 = 7.8243e-04 + SC2 = 5.1797e-04 + SC3 = 0.043021 + W_density = 0.7133 + H_density = 0.5800 iter = 20 - elapsed = 9.4930e-03 - tTotal = 0.3412 + elapsed = 6.5570e-03 + tTotal = 0.1871 subIterW = 3 subIterH = 0 - relError = 0.2340 + relError = 0.2370 toshow = scalar structure containing the fields: - initGrNormW = 361.43 - initGrNormH = 433.24 - initNorm = 14.263 - SC1 = 5.9630e-04 - SC2 = 4.0787e-04 - SC3 = 0.012788 - W_density = 0.7767 - H_density = 0.5450 + initGrNormW = 378.98 + initGrNormH = 458.15 + initNorm = 13.836 + SC1 = 7.4388e-04 + SC2 = 4.9096e-04 + SC3 = 0.041174 + W_density = 0.7100 + H_density = 0.5850 iter = 21 - elapsed = 9.0030e-03 - tTotal = 0.3502 - subIterW = 0 + elapsed = 6.8080e-03 + tTotal = 0.1939 + subIterW = 3 subIterH = 0 - relError = 0.2338 + relError = 0.2369 toshow = scalar structure containing the fields: - initGrNormW = 361.43 - initGrNormH = 433.24 - initNorm = 14.263 - SC1 = 5.6425e-04 - SC2 = 3.8369e-04 - SC3 = 0.028364 - W_density = 0.7767 - H_density = 0.5350 + initGrNormW = 378.98 + initGrNormH = 458.15 + initNorm = 13.836 + SC1 = 6.7854e-04 + SC2 = 4.4784e-04 + SC3 = 0.020994 + W_density = 0.7100 + H_density = 0.5850 iter = 22 - elapsed = 0.012466 - tTotal = 0.3627 - subIterW = 2 - subIterH = 2 - relError = 0.2337 + elapsed = 4.0100e-03 + tTotal = 0.1979 + subIterW = 0 + subIterH = 0 + relError = 0.2367 toshow = scalar structure containing the fields: - initGrNormW = 361.43 - initGrNormH = 433.24 - initNorm = 14.263 - SC1 = 5.0183e-04 - SC2 = 3.4024e-04 - SC3 = 0.012613 - W_density = 0.7767 - H_density = 0.5300 + initGrNormW = 378.98 + initGrNormH = 458.15 + initNorm = 13.836 + SC1 = 6.6294e-04 + SC2 = 4.4019e-04 + SC3 = 0.037412 + W_density = 0.7167 + H_density = 0.5850 iter = 23 - elapsed = 0.010607 - tTotal = 0.3733 - subIterW = 0 + elapsed = 5.5190e-03 + tTotal = 0.2035 + subIterW = 4 subIterH = 0 - relError = 0.2336 + relError = 0.2366 Stop: tolerance reached. final = @@ -5370,10 +5400,10 @@ scalar structure containing the fields: iterations = 23 - elapsed_total = 0.3733 - relative_error = 0.2336 - W_density = 0.7767 - H_density = 0.5300 + elapsed_total = 0.2035 + relative_error = 0.2366 + W_density = 0.7167 + H_density = 0.5850 ***** demo m = 300; @@ -5421,56 +5451,6 @@ 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_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/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/smwsolve.m] ->>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/smwsolve.m -***** 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); -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 -***** test - mtx = randn (100); - assert (funm (mtx, "fsin"), thfm (mtx, "sin"), 1e-9) -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); -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.5NpRqk2mWn at line 54 column 31 - -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/cod.m] >>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/cod.m ***** test @@ -5486,19 +5466,28 @@ [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/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); +[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_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); +[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 @@ -5510,6 +5499,15 @@ 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) + b = sin (a); +***** endfunction +***** test + mtx = randn (100); + assert (funm (mtx, "fsin"), thfm (mtx, "sin"), 1e-9) +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 @@ -5540,12 +5538,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/1800908/tmp/hooks/B01_cleanup starting +I: user script /srv/workspace/pbuilder/1800908/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/1023116 and its subdirectories -I: Current time: Sat May 18 09:37:35 -12 2024 -I: pbuilder-time-stamp: 1716068255 +I: removing directory /srv/workspace/pbuilder/1800908 and its subdirectories +I: Current time: Sat Jun 21 18:05:37 +14 2025 +I: pbuilder-time-stamp: 1750478737