Diff of the two buildlogs: -- --- b1/build.log 2025-03-22 02:10:24.382525035 +0000 +++ b2/build.log 2025-03-22 02:13:00.103487090 +0000 @@ -1,6 +1,6 @@ I: pbuilder: network access will be disabled during build -I: Current time: Fri Mar 21 14:01:18 -12 2025 -I: pbuilder-time-stamp: 1742608878 +I: Current time: Fri Apr 24 22:33:26 +14 2026 +I: pbuilder-time-stamp: 1777019606 I: Building the build Environment I: extracting base tarball [/var/cache/pbuilder/unstable-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/3091070/tmp/hooks/D02_print_environment starting +I: user script /srv/workspace/pbuilder/1344041/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 Apr 24 08:33 /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/1344041/tmp/hooks/D01_modify_environment finished +I: user script /srv/workspace/pbuilder/1344041/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]="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=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='c821eafbbdfd48f98ce3bbc7cedd3e51' - 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='3091070' - PS1='# ' - PS2='> ' + INVOCATION_ID=c75a6a9fc3e244a38c2a858e2dc85cf0 + 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=1344041 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.kuTU7Jwd/pbuilderrc_LeWf --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.kuTU7Jwd/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://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.kuTU7Jwd/pbuilderrc_y1Io --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.kuTU7Jwd/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 ionos11-amd64 6.1.0-32-amd64 #1 SMP PREEMPT_DYNAMIC Debian 6.1.129-1 (2025-03-06) x86_64 GNU/Linux + Linux i-capture-the-hostname 6.12.12+bpo-amd64 #1 SMP PREEMPT_DYNAMIC Debian 6.12.12-1~bpo12+1 (2025-02-23) x86_64 GNU/Linux I: ls -l /bin - lrwxrwxrwx 1 root root 7 Mar 4 11:20 /bin -> usr/bin -I: user script /srv/workspace/pbuilder/3091070/tmp/hooks/D02_print_environment finished + lrwxrwxrwx 1 root root 7 Mar 4 2025 /bin -> usr/bin +I: user script /srv/workspace/pbuilder/1344041/tmp/hooks/D02_print_environment finished -> Attempting to satisfy build-dependencies -> Creating pbuilder-satisfydepends-dummy package Package: pbuilder-satisfydepends-dummy @@ -635,7 +667,7 @@ Get: 528 http://deb.debian.org/debian unstable/main amd64 gfortran amd64 4:14.2.0-1 [1436 B] Get: 529 http://deb.debian.org/debian unstable/main amd64 octave-dev amd64 9.4.0-1 [1007 kB] Get: 530 http://deb.debian.org/debian unstable/main amd64 dh-octave all 1.8.0 [22.7 kB] -Fetched 208 MB in 25s (8271 kB/s) +Fetched 208 MB in 3s (61.0 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 ... 19807 files and directories currently installed.) @@ -2808,7 +2840,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/1344041/tmp/hooks/A99_set_merged_usr starting +Not re-configuring usrmerge for unstable +I: user script /srv/workspace/pbuilder/1344041/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 @@ -2848,25 +2884,6 @@ 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.UUiTdyWfjU 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 @@ -2878,15 +2895,26 @@ x2 = smwsolve (A, u, v, b); assert (x1, x2, 1e-13); 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); +[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.tHOrU5vGNv 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 @@ -2900,12 +2928,12 @@ Stop: tolerance reached. ***** test [W,H,iter,HIS]=nmf_bpas(A,k,'verbose'); - initGrNormW = 375.95 - initGrNormH = 438.89 - initNorm = 13.825 - SC1 = 1.1558 - SC2 = 577.90 - SC3 = 0.4917 + initGrNormW = 426.79 + initGrNormH = 483.57 + initNorm = 14.419 + SC1 = 1.2899 + SC2 = 644.97 + SC3 = 0.5122 W_density = 1 H_density = 1 par = @@ -2916,8 +2944,8 @@ n = 20 type = regularized nnls_solver = bp - alpha = 0.4885 - beta = 0.4885 + alpha = 0.5093 + beta = 0.5093 max_iter = 100 min_iter = 20 max_time = 1000 @@ -2929,438 +2957,438 @@ scalar structure containing the fields: - initGrNormW = 375.95 - initGrNormH = 438.89 - initNorm = 13.825 - SC1 = 0.025209 - SC2 = 0.023747 - SC3 = 1.7037 - W_density = 0.9300 - H_density = 0.6300 + initGrNormW = 426.79 + initGrNormH = 483.57 + initNorm = 14.419 + SC1 = 0.026979 + SC2 = 0.026493 + SC3 = 2.2049 + W_density = 0.9700 + H_density = 0.7100 iter = 1 - elapsed = 0.019988 - tTotal = 0.019988 - subIterW = 7 + elapsed = 4.6650e-03 + tTotal = 4.6650e-03 + subIterW = 6 subIterH = 24 - relError = 0.4313 + relError = 0.4298 toshow = scalar structure containing the fields: - initGrNormW = 375.95 - initGrNormH = 438.89 - initNorm = 13.825 - SC1 = 0.014074 - SC2 = 0.012892 - SC3 = 0.9621 - W_density = 0.8867 - H_density = 0.5700 + initGrNormW = 426.79 + initGrNormH = 483.57 + initNorm = 14.419 + SC1 = 0.015622 + SC2 = 0.014591 + SC3 = 1.2295 + W_density = 0.9200 + H_density = 0.5250 iter = 2 - elapsed = 0.024076 - tTotal = 0.044064 - subIterW = 12 - subIterH = 22 - relError = 0.3459 + elapsed = 5.5490e-03 + tTotal = 0.010214 + subIterW = 11 + subIterH = 26 + relError = 0.3474 toshow = scalar structure containing the fields: - initGrNormW = 375.95 - initGrNormH = 438.89 - initNorm = 13.825 - SC1 = 9.5747e-03 - SC2 = 8.6172e-03 - SC3 = 0.6644 - W_density = 0.8900 - H_density = 0.6150 + initGrNormW = 426.79 + initGrNormH = 483.57 + initNorm = 14.419 + SC1 = 0.010100 + SC2 = 9.1507e-03 + SC3 = 0.7819 + W_density = 0.8933 + H_density = 0.5250 iter = 3 - elapsed = 0.020470 - tTotal = 0.064534 - subIterW = 8 + elapsed = 6.1550e-03 + tTotal = 0.016369 + subIterW = 7 subIterH = 17 - relError = 0.3108 + relError = 0.3166 toshow = scalar structure containing the fields: - initGrNormW = 375.95 - initGrNormH = 438.89 - initNorm = 13.825 - SC1 = 6.8797e-03 - SC2 = 6.0541e-03 - SC3 = 0.4721 - W_density = 0.8667 - H_density = 0.6350 + initGrNormW = 426.79 + initGrNormH = 483.57 + initNorm = 14.419 + SC1 = 7.1944e-03 + SC2 = 6.2879e-03 + SC3 = 0.5561 + W_density = 0.8700 + H_density = 0.5700 iter = 4 - elapsed = 0.020262 - tTotal = 0.084796 - subIterW = 9 - subIterH = 10 - relError = 0.2942 + elapsed = 5.3190e-03 + tTotal = 0.021688 + subIterW = 8 + subIterH = 12 + relError = 0.3024 toshow = scalar structure containing the fields: - initGrNormW = 375.95 - initGrNormH = 438.89 - initNorm = 13.825 - SC1 = 5.1031e-03 - SC2 = 4.3989e-03 - SC3 = 0.3505 - W_density = 0.8500 - H_density = 0.6800 + initGrNormW = 426.79 + initGrNormH = 483.57 + initNorm = 14.419 + SC1 = 5.2745e-03 + SC2 = 4.5361e-03 + SC3 = 0.4043 + W_density = 0.8600 + H_density = 0.6000 iter = 5 - elapsed = 0.016293 - tTotal = 0.1011 - subIterW = 4 - subIterH = 12 - relError = 0.2861 + elapsed = 4.6100e-03 + tTotal = 0.026298 + subIterW = 3 + subIterH = 11 + relError = 0.2946 toshow = scalar structure containing the fields: - initGrNormW = 375.95 - initGrNormH = 438.89 - initNorm = 13.825 - SC1 = 3.8712e-03 - SC2 = 3.3215e-03 - SC3 = 0.2654 - W_density = 0.8433 - H_density = 0.7050 + initGrNormW = 426.79 + initGrNormH = 483.57 + initNorm = 14.419 + SC1 = 4.0947e-03 + SC2 = 3.4068e-03 + SC3 = 0.3115 + W_density = 0.8367 + H_density = 0.6400 iter = 6 - elapsed = 0.014386 - tTotal = 0.1155 - subIterW = 2 - subIterH = 5 - relError = 0.2815 + elapsed = 4.5780e-03 + tTotal = 0.030876 + subIterW = 6 + subIterH = 7 + relError = 0.2895 toshow = scalar structure containing the fields: - initGrNormW = 375.95 - initGrNormH = 438.89 - initNorm = 13.825 - SC1 = 3.0466e-03 - SC2 = 2.5835e-03 - SC3 = 0.2084 - W_density = 0.8333 - H_density = 0.7250 + initGrNormW = 426.79 + initGrNormH = 483.57 + initNorm = 14.419 + SC1 = 3.2494e-03 + SC2 = 2.6450e-03 + SC3 = 0.2475 + W_density = 0.8300 + H_density = 0.6350 iter = 7 - elapsed = 0.015036 - tTotal = 0.1305 - subIterW = 5 - subIterH = 4 - relError = 0.2786 + elapsed = 3.7120e-03 + tTotal = 0.034588 + subIterW = 2 + subIterH = 3 + relError = 0.2856 toshow = scalar structure containing the fields: - initGrNormW = 375.95 - initGrNormH = 438.89 - initNorm = 13.825 - SC1 = 2.4481e-03 - SC2 = 2.0564e-03 - SC3 = 0.1681 + initGrNormW = 426.79 + initGrNormH = 483.57 + initNorm = 14.419 + SC1 = 2.5692e-03 + SC2 = 2.1067e-03 + SC3 = 0.1905 W_density = 0.8333 - H_density = 0.7350 + H_density = 0.6550 iter = 8 - elapsed = 0.016157 - tTotal = 0.1467 + elapsed = 3.8590e-03 + tTotal = 0.038447 subIterW = 4 - subIterH = 2 - relError = 0.2762 + subIterH = 4 + relError = 0.2822 toshow = scalar structure containing the fields: - initGrNormW = 375.95 - initGrNormH = 438.89 - initNorm = 13.825 - SC1 = 2.0335e-03 - SC2 = 1.6838e-03 - SC3 = 0.1377 - W_density = 0.8233 - H_density = 0.7350 + initGrNormW = 426.79 + initGrNormH = 483.57 + initNorm = 14.419 + SC1 = 2.0839e-03 + SC2 = 1.6879e-03 + SC3 = 0.1513 + W_density = 0.8267 + H_density = 0.6650 iter = 9 - elapsed = 0.018466 - tTotal = 0.1651 - subIterW = 3 - subIterH = 2 - relError = 0.2741 + elapsed = 3.8330e-03 + tTotal = 0.042280 + subIterW = 4 + subIterH = 5 + relError = 0.2792 toshow = scalar structure containing the fields: - initGrNormW = 375.95 - initGrNormH = 438.89 - initNorm = 13.825 - SC1 = 1.7291e-03 - SC2 = 1.4248e-03 - SC3 = 0.1132 + initGrNormW = 426.79 + initGrNormH = 483.57 + initNorm = 14.419 + SC1 = 1.7163e-03 + SC2 = 1.3662e-03 + SC3 = 0.1222 W_density = 0.8133 - H_density = 0.7450 + H_density = 0.6600 iter = 10 - elapsed = 0.015080 - tTotal = 0.1802 - subIterW = 3 - subIterH = 5 - relError = 0.2718 + elapsed = 3.6240e-03 + tTotal = 0.045904 + subIterW = 4 + subIterH = 0 + relError = 0.2764 toshow = scalar structure containing the fields: - initGrNormW = 375.95 - initGrNormH = 438.89 - initNorm = 13.825 - SC1 = 1.4743e-03 - SC2 = 1.2030e-03 - SC3 = 0.092192 - W_density = 0.7967 - H_density = 0.7750 + initGrNormW = 426.79 + initGrNormH = 483.57 + initNorm = 14.419 + SC1 = 1.3878e-03 + SC2 = 1.1075e-03 + SC3 = 0.097329 + W_density = 0.8067 + H_density = 0.6800 iter = 11 - elapsed = 0.015852 - tTotal = 0.1961 - subIterW = 5 - subIterH = 5 - relError = 0.2695 + elapsed = 3.7930e-03 + tTotal = 0.049697 + subIterW = 2 + subIterH = 3 + relError = 0.2742 toshow = scalar structure containing the fields: - initGrNormW = 375.95 - initGrNormH = 438.89 - initNorm = 13.825 - SC1 = 1.2578e-03 - SC2 = 1.0163e-03 - SC3 = 0.077200 - W_density = 0.7900 - H_density = 0.7650 + initGrNormW = 426.79 + initGrNormH = 483.57 + initNorm = 14.419 + SC1 = 1.1467e-03 + SC2 = 9.1738e-04 + SC3 = 0.079175 + W_density = 0.8033 + H_density = 0.7050 iter = 12 - elapsed = 0.018752 - tTotal = 0.2148 + elapsed = 3.8430e-03 + tTotal = 0.053540 subIterW = 4 - subIterH = 5 - relError = 0.2672 + subIterH = 4 + relError = 0.2724 toshow = scalar structure containing the fields: - initGrNormW = 375.95 - initGrNormH = 438.89 - initNorm = 13.825 - SC1 = 1.0506e-03 - SC2 = 8.4467e-04 - SC3 = 0.064311 - W_density = 0.7867 - H_density = 0.7850 + initGrNormW = 426.79 + initGrNormH = 483.57 + initNorm = 14.419 + SC1 = 9.8026e-04 + SC2 = 7.7636e-04 + SC3 = 0.038391 + W_density = 0.8033 + H_density = 0.7150 iter = 13 - elapsed = 0.018231 - tTotal = 0.2330 - subIterW = 5 + elapsed = 3.0210e-03 + tTotal = 0.056561 + subIterW = 0 subIterH = 4 - relError = 0.2653 + relError = 0.2711 toshow = scalar structure containing the fields: - initGrNormW = 375.95 - initGrNormH = 438.89 - initNorm = 13.825 - SC1 = 8.7042e-04 - SC2 = 6.9807e-04 - SC3 = 0.053598 - W_density = 0.7900 - H_density = 0.7850 + initGrNormW = 426.79 + initGrNormH = 483.57 + initNorm = 14.419 + SC1 = 8.5358e-04 + SC2 = 6.6921e-04 + SC3 = 0.058435 + W_density = 0.7967 + H_density = 0.7100 iter = 14 - elapsed = 0.015355 - tTotal = 0.2484 - subIterW = 3 - subIterH = 4 - relError = 0.2637 + elapsed = 4.2800e-03 + tTotal = 0.060841 + subIterW = 4 + subIterH = 0 + relError = 0.2701 toshow = scalar structure containing the fields: - initGrNormW = 375.95 - initGrNormH = 438.89 - initNorm = 13.825 - SC1 = 7.3315e-04 - SC2 = 5.8066e-04 - SC3 = 0.045385 - W_density = 0.7900 - H_density = 0.7800 + initGrNormW = 426.79 + initGrNormH = 483.57 + initNorm = 14.419 + SC1 = 7.4787e-04 + SC2 = 5.8035e-04 + SC3 = 0.051482 + W_density = 0.7967 + H_density = 0.7250 iter = 15 - elapsed = 0.015120 - tTotal = 0.2635 - subIterW = 4 + elapsed = 3.8090e-03 + tTotal = 0.064650 + subIterW = 2 subIterH = 3 - relError = 0.2625 + relError = 0.2694 toshow = scalar structure containing the fields: - initGrNormW = 375.95 - initGrNormH = 438.89 - initNorm = 13.825 - SC1 = 6.1726e-04 - SC2 = 4.8393e-04 - SC3 = 0.037274 - W_density = 0.7800 - H_density = 0.7750 + initGrNormW = 426.79 + initGrNormH = 483.57 + initNorm = 14.419 + SC1 = 6.5978e-04 + SC2 = 5.1595e-04 + SC3 = 0.044886 + W_density = 0.8067 + H_density = 0.7400 iter = 16 - elapsed = 0.014156 - tTotal = 0.2777 - subIterW = 4 - subIterH = 3 - relError = 0.2616 + elapsed = 3.7960e-03 + tTotal = 0.068446 + subIterW = 3 + subIterH = 2 + relError = 0.2688 toshow = scalar structure containing the fields: - initGrNormW = 375.95 - initGrNormH = 438.89 - initNorm = 13.825 - SC1 = 5.1283e-04 - SC2 = 3.9898e-04 - SC3 = 0.017693 - W_density = 0.7800 - H_density = 0.7600 + initGrNormW = 426.79 + initGrNormH = 483.57 + initNorm = 14.419 + SC1 = 5.9686e-04 + SC2 = 4.6436e-04 + SC3 = 0.039927 + W_density = 0.8033 + H_density = 0.7250 iter = 17 - elapsed = 0.010582 - tTotal = 0.2883 + elapsed = 3.6350e-03 + tTotal = 0.072081 subIterW = 0 - subIterH = 3 - relError = 0.2611 + subIterH = 2 + relError = 0.2683 toshow = scalar structure containing the fields: - initGrNormW = 375.95 - initGrNormH = 438.89 - initNorm = 13.825 - SC1 = 4.3282e-04 - SC2 = 3.3500e-04 - SC3 = 0.026235 - W_density = 0.7767 - H_density = 0.7550 + initGrNormW = 426.79 + initGrNormH = 483.57 + initNorm = 14.419 + SC1 = 5.4512e-04 + SC2 = 4.2410e-04 + SC3 = 0.035592 + W_density = 0.8000 + H_density = 0.7150 iter = 18 - elapsed = 0.012516 - tTotal = 0.3008 - subIterW = 0 - subIterH = 0 - relError = 0.2607 + elapsed = 3.7840e-03 + tTotal = 0.075865 + subIterW = 3 + subIterH = 2 + relError = 0.2679 toshow = scalar structure containing the fields: - initGrNormW = 375.95 - initGrNormH = 438.89 - initNorm = 13.825 - SC1 = 3.7592e-04 - SC2 = 2.9097e-04 - SC3 = 0.012678 - W_density = 0.7767 - H_density = 0.7600 + initGrNormW = 426.79 + initGrNormH = 483.57 + initNorm = 14.419 + SC1 = 5.0988e-04 + SC2 = 3.9363e-04 + SC3 = 0.033335 + W_density = 0.7933 + H_density = 0.7100 iter = 19 - elapsed = 0.010199 - tTotal = 0.3110 - subIterW = 0 + elapsed = 3.6270e-03 + tTotal = 0.079492 + subIterW = 2 subIterH = 0 - relError = 0.2604 + relError = 0.2674 toshow = scalar structure containing the fields: - initGrNormW = 375.95 - initGrNormH = 438.89 - initNorm = 13.825 - SC1 = 3.3630e-04 - SC2 = 2.6030e-04 - SC3 = 0.011238 - W_density = 0.7767 - H_density = 0.7600 + initGrNormW = 426.79 + initGrNormH = 483.57 + initNorm = 14.419 + SC1 = 4.7702e-04 + SC2 = 3.6826e-04 + SC3 = 0.031152 + W_density = 0.7900 + H_density = 0.7200 iter = 20 - elapsed = 6.9890e-03 - tTotal = 0.3180 - subIterW = 0 - subIterH = 0 - relError = 0.2601 + elapsed = 3.7850e-03 + tTotal = 0.083277 + subIterW = 3 + subIterH = 2 + relError = 0.2670 toshow = scalar structure containing the fields: - initGrNormW = 375.95 - initGrNormH = 438.89 - initNorm = 13.825 - SC1 = 3.0716e-04 - SC2 = 2.3774e-04 - SC3 = 0.018257 + initGrNormW = 426.79 + initGrNormH = 483.57 + initNorm = 14.419 + SC1 = 4.5843e-04 + SC2 = 3.5116e-04 + SC3 = 0.029622 W_density = 0.7800 - H_density = 0.7600 + H_density = 0.7350 iter = 21 - elapsed = 0.012660 - tTotal = 0.3306 - subIterW = 0 - subIterH = 2 - relError = 0.2598 + elapsed = 3.8970e-03 + tTotal = 0.087174 + subIterW = 3 + subIterH = 3 + relError = 0.2666 toshow = scalar structure containing the fields: - initGrNormW = 375.95 - initGrNormH = 438.89 - initNorm = 13.825 - SC1 = 2.8491e-04 - SC2 = 2.1995e-04 - SC3 = 0.016739 - W_density = 0.7767 - H_density = 0.7550 + initGrNormW = 426.79 + initGrNormH = 483.57 + initNorm = 14.419 + SC1 = 4.3233e-04 + SC2 = 3.3030e-04 + SC3 = 0.027794 + W_density = 0.7733 + H_density = 0.7350 iter = 22 - elapsed = 0.013114 - tTotal = 0.3437 - subIterW = 0 - subIterH = 0 - relError = 0.2596 + elapsed = 3.9010e-03 + tTotal = 0.091075 + subIterW = 2 + subIterH = 2 + relError = 0.2662 toshow = scalar structure containing the fields: - initGrNormW = 375.95 - initGrNormH = 438.89 - initNorm = 13.825 - SC1 = 2.6027e-04 - SC2 = 2.0093e-04 - SC3 = 0.015105 - W_density = 0.7767 - H_density = 0.7550 + initGrNormW = 426.79 + initGrNormH = 483.57 + initNorm = 14.419 + SC1 = 4.0992e-04 + SC2 = 3.1154e-04 + SC3 = 0.025822 + W_density = 0.7667 + H_density = 0.7450 iter = 23 - elapsed = 0.010515 - tTotal = 0.3543 - subIterW = 2 - subIterH = 0 - relError = 0.2594 + elapsed = 4.4170e-03 + tTotal = 0.095492 + subIterW = 0 + subIterH = 2 + relError = 0.2659 Stop: tolerance reached. final = @@ -3368,19 +3396,19 @@ scalar structure containing the fields: iterations = 23 - elapsed_total = 0.3543 - relative_error = 0.2594 - W_density = 0.7767 - H_density = 0.7550 + elapsed_total = 0.095492 + relative_error = 0.2659 + W_density = 0.7667 + H_density = 0.7450 ***** test [W,H,iter,HIS]=nmf_bpas(A,k,'verbose','nnlssolver','as'); - initGrNormW = 433.70 - initGrNormH = 532.78 - initNorm = 13.825 - SC1 = 1.3740 - SC2 = 686.99 - SC3 = 0.5177 + initGrNormW = 302.05 + initGrNormH = 414.73 + initNorm = 14.419 + SC1 = 1.0261 + SC2 = 513.07 + SC3 = 0.4807 W_density = 1 H_density = 1 par = @@ -3391,8 +3419,8 @@ n = 20 type = regularized nnls_solver = as - alpha = 0.4885 - beta = 0.4885 + alpha = 0.5093 + beta = 0.5093 max_iter = 100 min_iter = 20 max_time = 1000 @@ -3404,438 +3432,438 @@ scalar structure containing the fields: - initGrNormW = 433.70 - initGrNormH = 532.78 - initNorm = 13.825 - SC1 = 0.028526 - SC2 = 0.028126 - SC3 = 1.3422 - W_density = 0.9800 + initGrNormW = 302.05 + initGrNormH = 414.73 + initNorm = 14.419 + SC1 = 0.036055 + SC2 = 0.034901 + SC3 = 1.3088 + W_density = 0.9600 H_density = 0.6150 iter = 1 - elapsed = 0.013967 - tTotal = 0.013967 - subIterW = 5 - subIterH = 42 - relError = 0.4425 + elapsed = 4.2140e-03 + tTotal = 4.2140e-03 + subIterW = 7 + subIterH = 31 + relError = 0.4324 toshow = scalar structure containing the fields: - initGrNormW = 433.70 - initGrNormH = 532.78 - initNorm = 13.825 - SC1 = 0.015212 - SC2 = 0.014481 - SC3 = 0.7127 - W_density = 0.9267 - H_density = 0.5200 + initGrNormW = 302.05 + initGrNormH = 414.73 + initNorm = 14.419 + SC1 = 0.019358 + SC2 = 0.018545 + SC3 = 0.7040 + W_density = 0.9333 + H_density = 0.5500 iter = 2 - elapsed = 0.020090 - tTotal = 0.034057 - subIterW = 12 - subIterH = 40 - relError = 0.3450 + elapsed = 4.1760e-03 + tTotal = 8.3900e-03 + subIterW = 13 + subIterH = 36 + relError = 0.3436 toshow = scalar structure containing the fields: - initGrNormW = 433.70 - initGrNormH = 532.78 - initNorm = 13.825 - SC1 = 9.4019e-03 - SC2 = 8.8566e-03 - SC3 = 0.4388 - W_density = 0.9200 - H_density = 0.5550 + initGrNormW = 302.05 + initGrNormH = 414.73 + initNorm = 14.419 + SC1 = 0.012393 + SC2 = 0.011600 + SC3 = 0.4488 + W_density = 0.9167 + H_density = 0.5750 iter = 3 - elapsed = 0.015972 - tTotal = 0.050029 - subIterW = 15 - subIterH = 42 - relError = 0.3073 + elapsed = 4.5340e-03 + tTotal = 0.012924 + subIterW = 14 + subIterH = 32 + relError = 0.3082 toshow = scalar structure containing the fields: - initGrNormW = 433.70 - initGrNormH = 532.78 - initNorm = 13.825 - SC1 = 6.6270e-03 - SC2 = 6.0439e-03 - SC3 = 0.3051 - W_density = 0.8867 - H_density = 0.6050 + initGrNormW = 302.05 + initGrNormH = 414.73 + initNorm = 14.419 + SC1 = 9.0424e-03 + SC2 = 8.0296e-03 + SC3 = 0.3203 + W_density = 0.8700 + H_density = 0.6100 iter = 4 - elapsed = 0.016239 - tTotal = 0.066268 - subIterW = 22 - subIterH = 37 - relError = 0.2911 + elapsed = 5.8110e-03 + tTotal = 0.018735 + subIterW = 24 + subIterH = 41 + relError = 0.2915 toshow = scalar structure containing the fields: - initGrNormW = 433.70 - initGrNormH = 532.78 - initNorm = 13.825 - SC1 = 4.8799e-03 - SC2 = 4.3627e-03 - SC3 = 0.2224 - W_density = 0.8667 - H_density = 0.6250 + initGrNormW = 302.05 + initGrNormH = 414.73 + initNorm = 14.419 + SC1 = 6.7471e-03 + SC2 = 5.8295e-03 + SC3 = 0.2529 + W_density = 0.8567 + H_density = 0.6450 iter = 5 - elapsed = 0.017617 - tTotal = 0.083885 - subIterW = 23 - subIterH = 34 - relError = 0.2828 + elapsed = 4.4960e-03 + tTotal = 0.023231 + subIterW = 24 + subIterH = 33 + relError = 0.2824 toshow = scalar structure containing the fields: - initGrNormW = 433.70 - initGrNormH = 532.78 - initNorm = 13.825 - SC1 = 3.7450e-03 - SC2 = 3.2582e-03 - SC3 = 0.1731 + initGrNormW = 302.05 + initGrNormH = 414.73 + initNorm = 14.419 + SC1 = 5.1631e-03 + SC2 = 4.3576e-03 + SC3 = 0.1922 W_density = 0.8400 - H_density = 0.6600 + H_density = 0.6650 iter = 6 - elapsed = 0.017385 - tTotal = 0.1013 - subIterW = 26 - subIterH = 35 - relError = 0.2776 + elapsed = 4.4880e-03 + tTotal = 0.027719 + subIterW = 24 + subIterH = 30 + relError = 0.2771 toshow = scalar structure containing the fields: - initGrNormW = 433.70 - initGrNormH = 532.78 - initNorm = 13.825 - SC1 = 2.9679e-03 - SC2 = 2.4990e-03 - SC3 = 0.1352 - W_density = 0.8200 + initGrNormW = 302.05 + initGrNormH = 414.73 + initNorm = 14.419 + SC1 = 4.0309e-03 + SC2 = 3.3456e-03 + SC3 = 0.1456 + W_density = 0.8233 H_density = 0.7000 iter = 7 - elapsed = 0.016340 - tTotal = 0.1176 - subIterW = 29 - subIterH = 34 - relError = 0.2744 + elapsed = 4.3850e-03 + tTotal = 0.032104 + subIterW = 28 + subIterH = 26 + relError = 0.2736 toshow = scalar structure containing the fields: - initGrNormW = 433.70 - initGrNormH = 532.78 - initNorm = 13.825 - SC1 = 2.3574e-03 - SC2 = 1.9661e-03 - SC3 = 0.1047 - W_density = 0.8100 - H_density = 0.7150 + initGrNormW = 302.05 + initGrNormH = 414.73 + initNorm = 14.419 + SC1 = 3.2267e-03 + SC2 = 2.6137e-03 + SC3 = 0.1143 + W_density = 0.8000 + H_density = 0.7100 iter = 8 - elapsed = 0.015907 - tTotal = 0.1335 - subIterW = 30 - subIterH = 29 - relError = 0.2723 + elapsed = 4.0000e-03 + tTotal = 0.036104 + subIterW = 29 + subIterH = 27 + relError = 0.2710 toshow = scalar structure containing the fields: - initGrNormW = 433.70 - initGrNormH = 532.78 - initNorm = 13.825 - SC1 = 1.9377e-03 - SC2 = 1.5851e-03 - SC3 = 0.082667 - W_density = 0.8067 - H_density = 0.7300 + initGrNormW = 302.05 + initGrNormH = 414.73 + initNorm = 14.419 + SC1 = 2.5713e-03 + SC2 = 2.0673e-03 + SC3 = 0.091923 + W_density = 0.7900 + H_density = 0.7100 iter = 9 - elapsed = 0.014203 - tTotal = 0.1477 + elapsed = 4.4770e-03 + tTotal = 0.040581 subIterW = 30 - subIterH = 28 - relError = 0.2708 + subIterH = 29 + relError = 0.2691 toshow = scalar structure containing the fields: - initGrNormW = 433.70 - initGrNormH = 532.78 - initNorm = 13.825 - SC1 = 1.6149e-03 - SC2 = 1.2984e-03 - SC3 = 0.069858 - W_density = 0.7967 - H_density = 0.7250 + initGrNormW = 302.05 + initGrNormH = 414.73 + initNorm = 14.419 + SC1 = 2.0994e-03 + SC2 = 1.6711e-03 + SC3 = 0.074484 + W_density = 0.7767 + H_density = 0.7200 iter = 10 - elapsed = 0.016116 - tTotal = 0.1638 - subIterW = 32 - subIterH = 30 - relError = 0.2697 + elapsed = 4.7750e-03 + tTotal = 0.045356 + subIterW = 35 + subIterH = 31 + relError = 0.2676 toshow = scalar structure containing the fields: - initGrNormW = 433.70 - initGrNormH = 532.78 - initNorm = 13.825 - SC1 = 1.3411e-03 - SC2 = 1.0702e-03 - SC3 = 0.058978 - W_density = 0.7867 - H_density = 0.7300 + initGrNormW = 302.05 + initGrNormH = 414.73 + initNorm = 14.419 + SC1 = 1.7478e-03 + SC2 = 1.3843e-03 + SC3 = 0.060205 + W_density = 0.7733 + H_density = 0.7350 iter = 11 - elapsed = 0.016607 - tTotal = 0.1804 - subIterW = 30 - subIterH = 29 - relError = 0.2688 + elapsed = 4.7260e-03 + tTotal = 0.050082 + subIterW = 33 + subIterH = 30 + relError = 0.2662 toshow = scalar structure containing the fields: - initGrNormW = 433.70 - initGrNormH = 532.78 - initNorm = 13.825 - SC1 = 1.1254e-03 - SC2 = 8.8903e-04 - SC3 = 0.050699 - W_density = 0.7800 - H_density = 0.7400 + initGrNormW = 302.05 + initGrNormH = 414.73 + initNorm = 14.419 + SC1 = 1.4945e-03 + SC2 = 1.1597e-03 + SC3 = 0.051122 + W_density = 0.7633 + H_density = 0.7550 iter = 12 - elapsed = 0.015093 - tTotal = 0.1955 - subIterW = 30 - subIterH = 25 - relError = 0.2680 + elapsed = 4.0410e-03 + tTotal = 0.054123 + subIterW = 36 + subIterH = 23 + relError = 0.2651 toshow = scalar structure containing the fields: - initGrNormW = 433.70 - initGrNormH = 532.78 - initNorm = 13.825 - SC1 = 9.4867e-04 - SC2 = 7.4565e-04 - SC3 = 0.041914 - W_density = 0.7733 - H_density = 0.7300 + initGrNormW = 302.05 + initGrNormH = 414.73 + initNorm = 14.419 + SC1 = 1.2714e-03 + SC2 = 9.8405e-04 + SC3 = 0.045491 + W_density = 0.7700 + H_density = 0.7400 iter = 13 - elapsed = 0.018576 - tTotal = 0.2141 - subIterW = 30 - subIterH = 25 - relError = 0.2673 + elapsed = 4.4470e-03 + tTotal = 0.058570 + subIterW = 33 + subIterH = 24 + relError = 0.2641 toshow = scalar structure containing the fields: - initGrNormW = 433.70 - initGrNormH = 532.78 - initNorm = 13.825 - SC1 = 8.1023e-04 - SC2 = 6.3360e-04 - SC3 = 0.035872 - W_density = 0.7700 - H_density = 0.7300 + initGrNormW = 302.05 + initGrNormH = 414.73 + initNorm = 14.419 + SC1 = 1.1014e-03 + SC2 = 8.3489e-04 + SC3 = 0.036111 + W_density = 0.7567 + H_density = 0.7350 iter = 14 - elapsed = 0.017891 - tTotal = 0.2320 - subIterW = 31 - subIterH = 28 - relError = 0.2666 + elapsed = 4.5000e-03 + tTotal = 0.063070 + subIterW = 36 + subIterH = 23 + relError = 0.2634 toshow = scalar structure containing the fields: - initGrNormW = 433.70 - initGrNormH = 532.78 - initNorm = 13.825 - SC1 = 7.0423e-04 - SC2 = 5.4789e-04 - SC3 = 0.030129 - W_density = 0.7667 - H_density = 0.7350 + initGrNormW = 302.05 + initGrNormH = 414.73 + initNorm = 14.419 + SC1 = 9.3720e-04 + SC2 = 7.0478e-04 + SC3 = 0.030655 + W_density = 0.7533 + H_density = 0.7300 iter = 15 - elapsed = 0.018624 - tTotal = 0.2506 - subIterW = 32 + elapsed = 4.0970e-03 + tTotal = 0.067167 + subIterW = 36 subIterH = 24 - relError = 0.2660 + relError = 0.2629 toshow = scalar structure containing the fields: - initGrNormW = 433.70 - initGrNormH = 532.78 - initNorm = 13.825 - SC1 = 6.2537e-04 - SC2 = 4.8529e-04 - SC3 = 0.026367 - W_density = 0.7667 + initGrNormW = 302.05 + initGrNormH = 414.73 + initNorm = 14.419 + SC1 = 7.9732e-04 + SC2 = 5.9639e-04 + SC3 = 0.028127 + W_density = 0.7500 H_density = 0.7350 iter = 16 - elapsed = 0.019051 - tTotal = 0.2697 - subIterW = 32 - subIterH = 27 - relError = 0.2653 + elapsed = 4.5170e-03 + tTotal = 0.071684 + subIterW = 37 + subIterH = 23 + relError = 0.2625 toshow = scalar structure containing the fields: - initGrNormW = 433.70 - initGrNormH = 532.78 - initNorm = 13.825 - SC1 = 5.6245e-04 - SC2 = 4.3196e-04 - SC3 = 0.022202 - W_density = 0.7600 - H_density = 0.7450 + initGrNormW = 302.05 + initGrNormH = 414.73 + initNorm = 14.419 + SC1 = 6.7685e-04 + SC2 = 5.0629e-04 + SC3 = 0.021466 + W_density = 0.7500 + H_density = 0.7350 iter = 17 - elapsed = 0.019331 - tTotal = 0.2890 - subIterW = 35 - subIterH = 29 - relError = 0.2647 + elapsed = 4.1640e-03 + tTotal = 0.075848 + subIterW = 38 + subIterH = 24 + relError = 0.2623 toshow = scalar structure containing the fields: - initGrNormW = 433.70 - initGrNormH = 532.78 - initNorm = 13.825 - SC1 = 5.1106e-04 - SC2 = 3.8943e-04 - SC3 = 0.019917 - W_density = 0.7567 - H_density = 0.7500 + initGrNormW = 302.05 + initGrNormH = 414.73 + initNorm = 14.419 + SC1 = 5.8796e-04 + SC2 = 4.3979e-04 + SC3 = 0.020906 + W_density = 0.7500 + H_density = 0.7350 iter = 18 - elapsed = 0.017745 - tTotal = 0.3068 - subIterW = 32 - subIterH = 27 - relError = 0.2641 + elapsed = 4.0330e-03 + tTotal = 0.079881 + subIterW = 37 + subIterH = 22 + relError = 0.2621 toshow = scalar structure containing the fields: - initGrNormW = 433.70 - initGrNormH = 532.78 - initNorm = 13.825 - SC1 = 4.6586e-04 - SC2 = 3.5685e-04 - SC3 = 0.017764 - W_density = 0.7633 - H_density = 0.7600 + initGrNormW = 302.05 + initGrNormH = 414.73 + initNorm = 14.419 + SC1 = 5.2038e-04 + SC2 = 3.8716e-04 + SC3 = 0.016449 + W_density = 0.7433 + H_density = 0.7350 iter = 19 - elapsed = 0.014386 - tTotal = 0.3211 - subIterW = 29 - subIterH = 23 - relError = 0.2635 + elapsed = 4.0770e-03 + tTotal = 0.083958 + subIterW = 39 + subIterH = 21 + relError = 0.2619 toshow = scalar structure containing the fields: - initGrNormW = 433.70 - initGrNormH = 532.78 - initNorm = 13.825 - SC1 = 4.3206e-04 - SC2 = 3.3009e-04 - SC3 = 0.016145 - W_density = 0.7600 - H_density = 0.7650 + initGrNormW = 302.05 + initGrNormH = 414.73 + initNorm = 14.419 + SC1 = 4.5715e-04 + SC2 = 3.4012e-04 + SC3 = 0.014293 + W_density = 0.7467 + H_density = 0.7350 iter = 20 - elapsed = 0.013863 - tTotal = 0.3350 - subIterW = 30 + elapsed = 3.9620e-03 + tTotal = 0.087920 + subIterW = 34 subIterH = 23 - relError = 0.2629 + relError = 0.2618 toshow = scalar structure containing the fields: - initGrNormW = 433.70 - initGrNormH = 532.78 - initNorm = 13.825 - SC1 = 3.9999e-04 - SC2 = 3.0720e-04 - SC3 = 0.015249 - W_density = 0.7633 - H_density = 0.7700 + initGrNormW = 302.05 + initGrNormH = 414.73 + initNorm = 14.419 + SC1 = 4.0718e-04 + SC2 = 3.0132e-04 + SC3 = 0.012626 + W_density = 0.7400 + H_density = 0.7350 iter = 21 - elapsed = 0.013268 - tTotal = 0.3483 - subIterW = 30 + elapsed = 4.5100e-03 + tTotal = 0.092430 + subIterW = 35 subIterH = 23 - relError = 0.2623 + relError = 0.2617 toshow = scalar structure containing the fields: - initGrNormW = 433.70 - initGrNormH = 532.78 - initNorm = 13.825 - SC1 = 3.7145e-04 - SC2 = 2.8304e-04 - SC3 = 0.014655 - W_density = 0.7533 - H_density = 0.7700 + initGrNormW = 302.05 + initGrNormH = 414.73 + initNorm = 14.419 + SC1 = 3.5237e-04 + SC2 = 2.5793e-04 + SC3 = 0.011107 + W_density = 0.7300 + H_density = 0.7350 iter = 22 - elapsed = 0.014335 - tTotal = 0.3626 - subIterW = 33 + elapsed = 4.1100e-03 + tTotal = 0.096540 + subIterW = 37 subIterH = 22 - relError = 0.2617 + relError = 0.2616 toshow = scalar structure containing the fields: - initGrNormW = 433.70 - initGrNormH = 532.78 - initNorm = 13.825 - SC1 = 3.4212e-04 - SC2 = 2.6069e-04 - SC3 = 0.013411 - W_density = 0.7500 - H_density = 0.7700 + initGrNormW = 302.05 + initGrNormH = 414.73 + initNorm = 14.419 + SC1 = 3.0859e-04 + SC2 = 2.2527e-04 + SC3 = 9.3612e-03 + W_density = 0.7300 + H_density = 0.7300 iter = 23 - elapsed = 0.015573 - tTotal = 0.3782 - subIterW = 35 - subIterH = 22 - relError = 0.2611 + elapsed = 4.1690e-03 + tTotal = 0.1007 + subIterW = 37 + subIterH = 23 + relError = 0.2616 Stop: tolerance reached. final = @@ -3843,19 +3871,19 @@ scalar structure containing the fields: iterations = 23 - elapsed_total = 0.3782 - relative_error = 0.2611 - W_density = 0.7500 - H_density = 0.7700 + elapsed_total = 0.1007 + relative_error = 0.2616 + W_density = 0.7300 + H_density = 0.7300 ***** test [W,H,iter,HIS]=nmf_bpas(A,k,'verbose','type','sparse'); - initGrNormW = 378.08 - initGrNormH = 483.35 - initNorm = 13.825 - SC1 = 1.2273 - SC2 = 613.65 - SC3 = 0.4963 + initGrNormW = 349.45 + initGrNormH = 457.65 + initNorm = 14.419 + SC1 = 1.1516 + SC2 = 575.81 + SC3 = 0.4892 W_density = 1 H_density = 1 par = @@ -3866,8 +3894,8 @@ n = 20 type = sparse nnls_solver = bp - alpha = 0.4885 - beta = 0.4885 + alpha = 0.5093 + beta = 0.5093 max_iter = 100 min_iter = 20 max_time = 1000 @@ -3879,458 +3907,496 @@ scalar structure containing the fields: - initGrNormW = 378.08 - initGrNormH = 483.35 - initNorm = 13.825 - SC1 = 0.022508 - SC2 = 0.019987 - SC3 = 1.6119 - W_density = 0.9333 - H_density = 0.5400 + initGrNormW = 349.45 + initGrNormH = 457.65 + initNorm = 14.419 + SC1 = 0.023804 + SC2 = 0.022994 + SC3 = 1.6858 + W_density = 0.9700 + H_density = 0.5900 iter = 1 - elapsed = 0.019471 - tTotal = 0.019471 - subIterW = 7 - subIterH = 38 - relError = 0.4364 + elapsed = 6.0240e-03 + tTotal = 6.0240e-03 + subIterW = 5 + subIterH = 33 + relError = 0.4287 toshow = scalar structure containing the fields: - initGrNormW = 378.08 - initGrNormH = 483.35 - initNorm = 13.825 - SC1 = 0.012318 - SC2 = 9.1892e-03 - SC3 = 0.7946 - W_density = 0.8033 - H_density = 0.3550 + initGrNormW = 349.45 + initGrNormH = 457.65 + initNorm = 14.419 + SC1 = 0.011513 + SC2 = 9.8552e-03 + SC3 = 0.7975 + W_density = 0.9367 + H_density = 0.3700 iter = 2 - elapsed = 0.017739 - tTotal = 0.037210 - subIterW = 12 - subIterH = 23 - relError = 0.3709 + elapsed = 6.1170e-03 + tTotal = 0.012141 + subIterW = 10 + subIterH = 34 + relError = 0.3514 toshow = scalar structure containing the fields: - initGrNormW = 378.08 - initGrNormH = 483.35 - initNorm = 13.825 - SC1 = 8.0092e-03 - SC2 = 5.3501e-03 - SC3 = 0.5470 - W_density = 0.7900 - H_density = 0.2750 + initGrNormW = 349.45 + initGrNormH = 457.65 + initNorm = 14.419 + SC1 = 6.2548e-03 + SC2 = 4.8912e-03 + SC3 = 0.4578 + W_density = 0.9300 + H_density = 0.3550 iter = 3 - elapsed = 0.016935 - tTotal = 0.054145 - subIterW = 7 - subIterH = 17 - relError = 0.3540 + elapsed = 3.6600e-03 + tTotal = 0.015801 + subIterW = 5 + subIterH = 12 + relError = 0.3241 toshow = scalar structure containing the fields: - initGrNormW = 378.08 - initGrNormH = 483.35 - initNorm = 13.825 - SC1 = 6.1790e-03 - SC2 = 3.6209e-03 - SC3 = 0.4017 - W_density = 0.7167 - H_density = 0.2600 + initGrNormW = 349.45 + initGrNormH = 457.65 + initNorm = 14.419 + SC1 = 3.7768e-03 + SC2 = 2.8401e-03 + SC3 = 0.2791 + W_density = 0.9333 + H_density = 0.3300 iter = 4 - elapsed = 0.011951 - tTotal = 0.066096 - subIterW = 4 + elapsed = 3.5630e-03 + tTotal = 0.019364 + subIterW = 3 subIterH = 8 - relError = 0.3500 + relError = 0.3140 toshow = scalar structure containing the fields: - initGrNormW = 378.08 - initGrNormH = 483.35 - initNorm = 13.825 - SC1 = 5.0808e-03 - SC2 = 2.6217e-03 - SC3 = 0.3155 - W_density = 0.6400 - H_density = 0.2550 + initGrNormW = 349.45 + initGrNormH = 457.65 + initNorm = 14.419 + SC1 = 2.5262e-03 + SC2 = 1.8239e-03 + SC3 = 0.1886 + W_density = 0.9300 + H_density = 0.3200 iter = 5 - elapsed = 0.013945 - tTotal = 0.080041 + elapsed = 3.4700e-03 + tTotal = 0.022834 subIterW = 2 - subIterH = 10 - relError = 0.3476 + subIterH = 6 + relError = 0.3104 toshow = scalar structure containing the fields: - initGrNormW = 378.08 - initGrNormH = 483.35 - initNorm = 13.825 - SC1 = 3.9066e-03 - SC2 = 1.9064e-03 - SC3 = 0.2476 - W_density = 0.6267 - H_density = 0.2400 + initGrNormW = 349.45 + initGrNormH = 457.65 + initNorm = 14.419 + SC1 = 1.9703e-03 + SC2 = 1.3517e-03 + SC3 = 0.1543 + W_density = 0.9167 + H_density = 0.2900 iter = 6 - elapsed = 0.011081 - tTotal = 0.091122 + elapsed = 4.1820e-03 + tTotal = 0.027016 subIterW = 5 - subIterH = 6 - relError = 0.3470 + subIterH = 4 + relError = 0.3094 toshow = scalar structure containing the fields: - initGrNormW = 378.08 - initGrNormH = 483.35 - initNorm = 13.825 - SC1 = 2.9325e-03 - SC2 = 1.4017e-03 - SC3 = 0.1851 - W_density = 0.6233 - H_density = 0.2450 + initGrNormW = 349.45 + initGrNormH = 457.65 + initNorm = 14.419 + SC1 = 1.9486e-03 + SC2 = 1.3212e-03 + SC3 = 0.038478 + W_density = 0.9167 + H_density = 0.2800 iter = 7 - elapsed = 0.010378 - tTotal = 0.1015 - subIterW = 2 - subIterH = 5 - relError = 0.3472 + elapsed = 2.5790e-03 + tTotal = 0.029595 + subIterW = 0 + subIterH = 3 + relError = 0.3095 toshow = scalar structure containing the fields: - initGrNormW = 378.08 - initGrNormH = 483.35 - initNorm = 13.825 - SC1 = 2.3074e-03 - SC2 = 1.0752e-03 - SC3 = 0.1471 - W_density = 0.6133 - H_density = 0.2300 + initGrNormW = 349.45 + initGrNormH = 457.65 + initNorm = 14.419 + SC1 = 1.7704e-03 + SC2 = 1.1897e-03 + SC3 = 0.1355 + W_density = 0.9167 + H_density = 0.2550 iter = 8 - elapsed = 0.010102 - tTotal = 0.1116 + elapsed = 4.5480e-03 + tTotal = 0.034143 subIterW = 2 - subIterH = 3 - relError = 0.3476 + subIterH = 9 + relError = 0.3099 toshow = scalar structure containing the fields: - initGrNormW = 378.08 - initGrNormH = 483.35 - initNorm = 13.825 - SC1 = 1.8670e-03 - SC2 = 8.4389e-04 - SC3 = 0.1185 - W_density = 0.6000 - H_density = 0.2250 + initGrNormW = 349.45 + initGrNormH = 457.65 + initNorm = 14.419 + SC1 = 1.8358e-03 + SC2 = 1.1932e-03 + SC3 = 0.1433 + W_density = 0.8967 + H_density = 0.2450 iter = 9 - elapsed = 9.7840e-03 - tTotal = 0.1214 - subIterW = 4 - subIterH = 0 - relError = 0.3484 + elapsed = 3.4630e-03 + tTotal = 0.037606 + subIterW = 5 + subIterH = 6 + relError = 0.3102 toshow = scalar structure containing the fields: - initGrNormW = 378.08 - initGrNormH = 483.35 - initNorm = 13.825 - SC1 = 1.4724e-03 - SC2 = 6.6259e-04 - SC3 = 0.029559 - W_density = 0.6000 - H_density = 0.2200 + initGrNormW = 349.45 + initGrNormH = 457.65 + initNorm = 14.419 + SC1 = 1.8730e-03 + SC2 = 1.1875e-03 + SC3 = 0.1421 + W_density = 0.8767 + H_density = 0.2500 iter = 10 - elapsed = 7.0150e-03 - tTotal = 0.1284 - subIterW = 0 - subIterH = 0 - relError = 0.3492 + elapsed = 3.5770e-03 + tTotal = 0.041183 + subIterW = 5 + subIterH = 5 + relError = 0.3109 toshow = scalar structure containing the fields: - initGrNormW = 378.08 - initGrNormH = 483.35 - initNorm = 13.825 - SC1 = 1.1673e-03 - SC2 = 5.1828e-04 - SC3 = 0.072458 - W_density = 0.5900 - H_density = 0.2200 + initGrNormW = 349.45 + initGrNormH = 457.65 + initNorm = 14.419 + SC1 = 2.0316e-03 + SC2 = 1.2637e-03 + SC3 = 0.1510 + W_density = 0.8667 + H_density = 0.2400 iter = 11 - elapsed = 7.6000e-03 - tTotal = 0.1360 - subIterW = 1 - subIterH = 0 - relError = 0.3500 + elapsed = 3.4850e-03 + tTotal = 0.044668 + subIterW = 2 + subIterH = 3 + relError = 0.3119 toshow = scalar structure containing the fields: - initGrNormW = 378.08 - initGrNormH = 483.35 - initNorm = 13.825 - SC1 = 9.4193e-04 - SC2 = 4.1822e-04 - SC3 = 0.057794 - W_density = 0.5867 - H_density = 0.2200 + initGrNormW = 349.45 + initGrNormH = 457.65 + initNorm = 14.419 + SC1 = 2.1622e-03 + SC2 = 1.3319e-03 + SC3 = 0.1572 + W_density = 0.8600 + H_density = 0.2300 iter = 12 - elapsed = 7.5970e-03 - tTotal = 0.1436 - subIterW = 0 - subIterH = 0 - relError = 0.3507 + elapsed = 3.6020e-03 + tTotal = 0.048270 + subIterW = 4 + subIterH = 4 + relError = 0.3132 toshow = scalar structure containing the fields: - initGrNormW = 378.08 - initGrNormH = 483.35 - initNorm = 13.825 - SC1 = 8.1336e-04 - SC2 = 3.6113e-04 - SC3 = 0.015025 - W_density = 0.5867 - H_density = 0.2200 + initGrNormW = 349.45 + initGrNormH = 457.65 + initNorm = 14.419 + SC1 = 2.3371e-03 + SC2 = 1.4396e-03 + SC3 = 0.1603 + W_density = 0.8633 + H_density = 0.2250 iter = 13 - elapsed = 5.2430e-03 - tTotal = 0.1488 - subIterW = 0 + elapsed = 4.1000e-03 + tTotal = 0.052370 + subIterW = 5 subIterH = 0 - relError = 0.3514 + relError = 0.3147 toshow = scalar structure containing the fields: - initGrNormW = 378.08 - initGrNormH = 483.35 - initNorm = 13.825 - SC1 = 7.5278e-04 - SC2 = 3.3273e-04 - SC3 = 0.015081 - W_density = 0.5867 - H_density = 0.2200 + initGrNormW = 349.45 + initGrNormH = 457.65 + initNorm = 14.419 + SC1 = 2.2987e-03 + SC2 = 1.3792e-03 + SC3 = 0.1546 + W_density = 0.8500 + H_density = 0.2100 iter = 14 - elapsed = 4.9920e-03 - tTotal = 0.1538 - subIterW = 0 - subIterH = 0 - relError = 0.3520 + elapsed = 3.5500e-03 + tTotal = 0.055920 + subIterW = 4 + subIterH = 3 + relError = 0.3163 toshow = scalar structure containing the fields: - initGrNormW = 378.08 - initGrNormH = 483.35 - initNorm = 13.825 - SC1 = 7.2673e-04 - SC2 = 3.1976e-04 - SC3 = 0.043482 - W_density = 0.5833 - H_density = 0.2200 + initGrNormW = 349.45 + initGrNormH = 457.65 + initNorm = 14.419 + SC1 = 1.8978e-03 + SC2 = 1.1235e-03 + SC3 = 0.1319 + W_density = 0.8400 + H_density = 0.2100 iter = 15 - elapsed = 0.010042 - tTotal = 0.1639 - subIterW = 0 + elapsed = 3.4850e-03 + tTotal = 0.059405 + subIterW = 4 subIterH = 2 - relError = 0.3526 + relError = 0.3176 toshow = scalar structure containing the fields: - initGrNormW = 378.08 - initGrNormH = 483.35 - initNorm = 13.825 - SC1 = 6.9298e-04 - SC2 = 3.0214e-04 - SC3 = 0.013201 - W_density = 0.5833 + initGrNormW = 349.45 + initGrNormH = 457.65 + initNorm = 14.419 + SC1 = 1.6468e-03 + SC2 = 9.6175e-04 + SC3 = 0.1184 + W_density = 0.8300 H_density = 0.2150 iter = 16 - elapsed = 7.0700e-03 - tTotal = 0.1709 - subIterW = 0 + elapsed = 3.3610e-03 + tTotal = 0.062766 + subIterW = 3 subIterH = 0 - relError = 0.3532 + relError = 0.3190 toshow = scalar structure containing the fields: - initGrNormW = 378.08 - initGrNormH = 483.35 - initNorm = 13.825 - SC1 = 6.6974e-04 - SC2 = 2.9200e-04 - SC3 = 0.012846 - W_density = 0.5833 + initGrNormW = 349.45 + initGrNormH = 457.65 + initNorm = 14.419 + SC1 = 1.4185e-03 + SC2 = 8.2558e-04 + SC3 = 0.1036 + W_density = 0.8267 H_density = 0.2150 iter = 17 - elapsed = 5.0370e-03 - tTotal = 0.1760 + elapsed = 2.5800e-03 + tTotal = 0.065346 subIterW = 0 subIterH = 0 - relError = 0.3537 + relError = 0.3204 toshow = scalar structure containing the fields: - initGrNormW = 378.08 - initGrNormH = 483.35 - initNorm = 13.825 - SC1 = 6.5901e-04 - SC2 = 2.8601e-04 - SC3 = 0.041166 - W_density = 0.5800 - H_density = 0.2150 + initGrNormW = 349.45 + initGrNormH = 457.65 + initNorm = 14.419 + SC1 = 1.2624e-03 + SC2 = 7.2712e-04 + SC3 = 0.094563 + W_density = 0.8200 + H_density = 0.2100 iter = 18 - elapsed = 7.5720e-03 - tTotal = 0.1836 - subIterW = 0 + elapsed = 3.3260e-03 + tTotal = 0.068672 + subIterW = 2 subIterH = 0 - relError = 0.3543 + relError = 0.3217 toshow = scalar structure containing the fields: - initGrNormW = 378.08 - initGrNormH = 483.35 - initNorm = 13.825 - SC1 = 6.4479e-04 - SC2 = 2.7726e-04 - SC3 = 0.040299 - W_density = 0.5767 + initGrNormW = 349.45 + initGrNormH = 457.65 + initNorm = 14.419 + SC1 = 1.1518e-03 + SC2 = 6.6345e-04 + SC3 = 0.087000 + W_density = 0.8200 H_density = 0.2100 iter = 19 - elapsed = 9.5360e-03 - tTotal = 0.1931 - subIterW = 0 + elapsed = 2.6460e-03 + tTotal = 0.071318 + subIterW = 2 subIterH = 0 - relError = 0.3548 + relError = 0.3228 toshow = scalar structure containing the fields: - initGrNormW = 378.08 - initGrNormH = 483.35 - initNorm = 13.825 - SC1 = 6.2928e-04 - SC2 = 2.7059e-04 - SC3 = 0.012177 - W_density = 0.5767 + initGrNormW = 349.45 + initGrNormH = 457.65 + initNorm = 14.419 + SC1 = 1.1299e-03 + SC2 = 6.4628e-04 + SC3 = 0.083400 + W_density = 0.8133 H_density = 0.2100 iter = 20 - elapsed = 4.8420e-03 - tTotal = 0.1979 - subIterW = 0 + elapsed = 2.6400e-03 + tTotal = 0.073958 + subIterW = 2 subIterH = 0 - relError = 0.3553 + relError = 0.3239 toshow = scalar structure containing the fields: - initGrNormW = 378.08 - initGrNormH = 483.35 - initNorm = 13.825 - SC1 = 6.1772e-04 - SC2 = 2.6562e-04 - SC3 = 0.013024 - W_density = 0.5767 - H_density = 0.2100 + initGrNormW = 349.45 + initGrNormH = 457.65 + initNorm = 14.419 + SC1 = 1.1020e-03 + SC2 = 6.2813e-04 + SC3 = 0.022117 + W_density = 0.8133 + H_density = 0.2050 iter = 21 - elapsed = 4.9740e-03 - tTotal = 0.2029 + elapsed = 2.5570e-03 + tTotal = 0.076515 subIterW = 0 subIterH = 0 - relError = 0.3558 + relError = 0.3249 toshow = scalar structure containing the fields: - initGrNormW = 378.08 - initGrNormH = 483.35 - initNorm = 13.825 - SC1 = 6.0928e-04 - SC2 = 2.6199e-04 - SC3 = 0.011129 - W_density = 0.5767 - H_density = 0.2100 + initGrNormW = 349.45 + initGrNormH = 457.65 + initNorm = 14.419 + SC1 = 1.0116e-03 + SC2 = 5.7457e-04 + SC3 = 0.073587 + W_density = 0.8100 + H_density = 0.2050 iter = 22 - elapsed = 5.7140e-03 - tTotal = 0.2086 + elapsed = 2.6800e-03 + tTotal = 0.079195 subIterW = 0 subIterH = 0 - relError = 0.3562 + relError = 0.3255 toshow = scalar structure containing the fields: - initGrNormW = 378.08 - initGrNormH = 483.35 - initNorm = 13.825 - SC1 = 6.0373e-04 - SC2 = 2.5960e-04 - SC3 = 0.011818 - W_density = 0.5767 - H_density = 0.2100 + initGrNormW = 349.45 + initGrNormH = 457.65 + initNorm = 14.419 + SC1 = 9.5931e-04 + SC2 = 5.4297e-04 + SC3 = 0.068392 + W_density = 0.8067 + H_density = 0.2050 iter = 23 - elapsed = 5.2950e-03 - tTotal = 0.2139 + elapsed = 2.6430e-03 + tTotal = 0.081838 + subIterW = 0 + subIterH = 0 + relError = 0.3261 + +toshow = + + scalar structure containing the fields: + + initGrNormW = 349.45 + initGrNormH = 457.65 + initNorm = 14.419 + SC1 = 9.1652e-04 + SC2 = 5.1691e-04 + SC3 = 0.063525 + W_density = 0.8033 + H_density = 0.2050 + iter = 24 + elapsed = 2.6190e-03 + tTotal = 0.084457 + subIterW = 0 + subIterH = 0 + relError = 0.3267 + +toshow = + + scalar structure containing the fields: + + initGrNormW = 349.45 + initGrNormH = 457.65 + initNorm = 14.419 + SC1 = 8.8474e-04 + SC2 = 4.9899e-04 + SC3 = 0.016964 + W_density = 0.8033 + H_density = 0.2050 + iter = 25 + elapsed = 1.9940e-03 + tTotal = 0.086451 subIterW = 0 subIterH = 0 - relError = 0.3567 + relError = 0.3273 Stop: tolerance reached. final = scalar structure containing the fields: - iterations = 23 - elapsed_total = 0.2139 - relative_error = 0.3567 - W_density = 0.5767 - H_density = 0.2100 + iterations = 25 + elapsed_total = 0.086451 + relative_error = 0.3273 + W_density = 0.8033 + H_density = 0.2050 ***** test [W,H,iter,HIS]=nmf_bpas(A,k,'verbose','type','sparse','nnlssolver','bp','alpha',1.1,'beta',1.3); - initGrNormW = 346.84 - initGrNormH = 520.85 - initNorm = 13.825 - SC1 = 1.2515 - SC2 = 625.77 - SC3 = 0.4933 + initGrNormW = 369.08 + initGrNormH = 527.81 + initNorm = 14.419 + SC1 = 1.2881 + SC2 = 644.05 + SC3 = 0.4891 W_density = 1 H_density = 1 par = @@ -4354,438 +4420,438 @@ scalar structure containing the fields: - initGrNormW = 346.84 - initGrNormH = 520.85 - initNorm = 13.825 - SC1 = 0.026978 - SC2 = 0.025683 - SC3 = 0.9973 - W_density = 1 - H_density = 0.5700 + initGrNormW = 369.08 + initGrNormH = 527.81 + initNorm = 14.419 + SC1 = 0.026501 + SC2 = 0.024010 + SC3 = 2.0021 + W_density = 0.9867 + H_density = 0.5350 iter = 1 - elapsed = 0.018047 - tTotal = 0.018047 - subIterW = 0 - subIterH = 37 - relError = 0.5664 + elapsed = 5.9080e-03 + tTotal = 5.9080e-03 + subIterW = 3 + subIterH = 39 + relError = 0.5448 toshow = scalar structure containing the fields: - initGrNormW = 346.84 - initGrNormH = 520.85 - initNorm = 13.825 - SC1 = 0.015566 - SC2 = 8.8727e-03 - SC3 = 0.9350 - W_density = 0.7400 - H_density = 0.2600 + initGrNormW = 369.08 + initGrNormH = 527.81 + initNorm = 14.419 + SC1 = 0.015943 + SC2 = 7.9076e-03 + SC3 = 0.9764 + W_density = 0.6300 + H_density = 0.2200 iter = 2 - elapsed = 0.022747 - tTotal = 0.040794 - subIterW = 5 - subIterH = 37 - relError = 0.4632 + elapsed = 6.3430e-03 + tTotal = 0.012251 + subIterW = 9 + subIterH = 33 + relError = 0.4596 toshow = scalar structure containing the fields: - initGrNormW = 346.84 - initGrNormH = 520.85 - initNorm = 13.825 - SC1 = 0.019563 - SC2 = 6.6122e-03 - SC3 = 0.8895 - W_density = 0.4300 - H_density = 0.1850 + initGrNormW = 369.08 + initGrNormH = 527.81 + initNorm = 14.419 + SC1 = 0.017734 + SC2 = 6.7388e-03 + SC3 = 0.9640 + W_density = 0.5067 + H_density = 0.1800 iter = 3 - elapsed = 0.019652 - tTotal = 0.060446 - subIterW = 10 - subIterH = 16 - relError = 0.4587 + elapsed = 4.4870e-03 + tTotal = 0.016738 + subIterW = 9 + subIterH = 9 + relError = 0.4640 toshow = scalar structure containing the fields: - initGrNormW = 346.84 - initGrNormH = 520.85 - initNorm = 13.825 - SC1 = 0.017443 - SC2 = 3.6978e-03 - SC3 = 0.6659 - W_density = 0.2667 - H_density = 0.1300 + initGrNormW = 369.08 + initGrNormH = 527.81 + initNorm = 14.419 + SC1 = 0.018172 + SC2 = 4.4704e-03 + SC3 = 0.6555 + W_density = 0.3133 + H_density = 0.1450 iter = 4 - elapsed = 0.011880 - tTotal = 0.072326 - subIterW = 4 - subIterH = 2 - relError = 0.4590 + elapsed = 3.8760e-03 + tTotal = 0.020614 + subIterW = 9 + subIterH = 3 + relError = 0.4698 toshow = scalar structure containing the fields: - initGrNormW = 346.84 - initGrNormH = 520.85 - initNorm = 13.825 - SC1 = 0.013991 - SC2 = 2.7142e-03 - SC3 = 0.5480 - W_density = 0.2400 - H_density = 0.1250 + initGrNormW = 369.08 + initGrNormH = 527.81 + initNorm = 14.419 + SC1 = 0.010686 + SC2 = 1.9663e-03 + SC3 = 0.3447 + W_density = 0.2267 + H_density = 0.1200 iter = 5 - elapsed = 8.2850e-03 - tTotal = 0.080611 - subIterW = 3 - subIterH = 0 - relError = 0.4637 + elapsed = 3.5560e-03 + tTotal = 0.024170 + subIterW = 4 + subIterH = 3 + relError = 0.4767 toshow = scalar structure containing the fields: - initGrNormW = 346.84 - initGrNormH = 520.85 - initNorm = 13.825 - SC1 = 0.011242 - SC2 = 1.7313e-03 - SC3 = 0.3440 - W_density = 0.1767 - H_density = 0.1200 + initGrNormW = 369.08 + initGrNormH = 527.81 + initNorm = 14.419 + SC1 = 9.2689e-03 + SC2 = 1.2606e-03 + SC3 = 0.1809 + W_density = 0.1533 + H_density = 0.1100 iter = 6 - elapsed = 7.9250e-03 - tTotal = 0.088536 - subIterW = 0 - subIterH = 0 - relError = 0.4684 + elapsed = 2.6110e-03 + tTotal = 0.026781 + subIterW = 2 + subIterH = 1 + relError = 0.4836 toshow = scalar structure containing the fields: - initGrNormW = 346.84 - initGrNormH = 520.85 - initNorm = 13.825 - SC1 = 9.8724e-03 - SC2 = 1.4414e-03 - SC3 = 0.2225 - W_density = 0.1667 - H_density = 0.1150 + initGrNormW = 369.08 + initGrNormH = 527.81 + initNorm = 14.419 + SC1 = 8.9857e-03 + SC2 = 1.2221e-03 + SC3 = 0.050179 + W_density = 0.1567 + H_density = 0.1050 iter = 7 - elapsed = 8.6100e-03 - tTotal = 0.097146 - subIterW = 1 + elapsed = 2.4420e-03 + tTotal = 0.029223 + subIterW = 2 subIterH = 0 - relError = 0.4724 + relError = 0.4868 toshow = scalar structure containing the fields: - initGrNormW = 346.84 - initGrNormH = 520.85 - initNorm = 13.825 - SC1 = 6.3969e-03 - SC2 = 9.2115e-04 - SC3 = 0.072207 - W_density = 0.1667 - H_density = 0.1100 + initGrNormW = 369.08 + initGrNormH = 527.81 + initNorm = 14.419 + SC1 = 1.0758e-03 + SC2 = 1.0758e-04 + SC3 = 0.021005 + W_density = 0.1000 + H_density = 0.1000 iter = 8 - elapsed = 5.3030e-03 - tTotal = 0.1024 + elapsed = 1.5800e-03 + tTotal = 0.030803 subIterW = 0 subIterH = 0 - relError = 0.4764 + relError = 0.4867 toshow = scalar structure containing the fields: - initGrNormW = 346.84 - initGrNormH = 520.85 - initNorm = 13.825 - SC1 = 8.7594e-03 - SC2 = 1.2263e-03 - SC3 = 0.1170 - W_density = 0.1600 - H_density = 0.1100 + initGrNormW = 369.08 + initGrNormH = 527.81 + initNorm = 14.419 + SC1 = 7.0379e-04 + SC2 = 7.0379e-05 + SC3 = 0.015854 + W_density = 0.1000 + H_density = 0.1000 iter = 9 - elapsed = 5.6680e-03 - tTotal = 0.1081 - subIterW = 1 + elapsed = 9.0500e-04 + tTotal = 0.031708 + subIterW = 0 subIterH = 0 - relError = 0.4812 + relError = 0.4867 toshow = scalar structure containing the fields: - initGrNormW = 346.84 - initGrNormH = 520.85 - initNorm = 13.825 - SC1 = 7.8176e-04 - SC2 = 7.8176e-05 - SC3 = 0.012313 + initGrNormW = 369.08 + initGrNormH = 527.81 + initNorm = 14.419 + SC1 = 4.6248e-04 + SC2 = 4.6248e-05 + SC3 = 9.3386e-03 W_density = 0.1000 H_density = 0.1000 iter = 10 - elapsed = 5.9200e-03 - tTotal = 0.1140 + elapsed = 8.7400e-04 + tTotal = 0.032582 subIterW = 0 - subIterH = 1 - relError = 0.4825 + subIterH = 0 + relError = 0.4867 toshow = scalar structure containing the fields: - initGrNormW = 346.84 - initGrNormH = 520.85 - initNorm = 13.825 - SC1 = 2.9292e-04 - SC2 = 2.9292e-05 - SC3 = 7.1816e-03 + initGrNormW = 369.08 + initGrNormH = 527.81 + initNorm = 14.419 + SC1 = 3.0400e-04 + SC2 = 3.0400e-05 + SC3 = 5.7424e-03 W_density = 0.1000 H_density = 0.1000 iter = 11 - elapsed = 2.8560e-03 - tTotal = 0.1169 + elapsed = 8.9900e-04 + tTotal = 0.033481 subIterW = 0 subIterH = 0 - relError = 0.4825 + relError = 0.4866 toshow = scalar structure containing the fields: - initGrNormW = 346.84 - initGrNormH = 520.85 - initNorm = 13.825 - SC1 = 1.8842e-04 - SC2 = 1.8842e-05 - SC3 = 3.6592e-03 + initGrNormW = 369.08 + initGrNormH = 527.81 + initNorm = 14.419 + SC1 = 1.9986e-04 + SC2 = 1.9986e-05 + SC3 = 5.3198e-03 W_density = 0.1000 H_density = 0.1000 iter = 12 - elapsed = 2.8820e-03 - tTotal = 0.1198 + elapsed = 8.7400e-04 + tTotal = 0.034355 subIterW = 0 subIterH = 0 - relError = 0.4825 + relError = 0.4866 toshow = scalar structure containing the fields: - initGrNormW = 346.84 - initGrNormH = 520.85 - initNorm = 13.825 - SC1 = 1.2161e-04 - SC2 = 1.2161e-05 - SC3 = 2.4454e-03 + initGrNormW = 369.08 + initGrNormH = 527.81 + initNorm = 14.419 + SC1 = 1.3142e-04 + SC2 = 1.3142e-05 + SC3 = 2.8502e-03 W_density = 0.1000 H_density = 0.1000 iter = 13 - elapsed = 2.4930e-03 - tTotal = 0.1223 + elapsed = 8.5800e-04 + tTotal = 0.035213 subIterW = 0 subIterH = 0 - relError = 0.4825 + relError = 0.4866 toshow = scalar structure containing the fields: - initGrNormW = 346.84 - initGrNormH = 520.85 - initNorm = 13.825 - SC1 = 7.8492e-05 - SC2 = 7.8492e-06 - SC3 = 1.8414e-03 + initGrNormW = 369.08 + initGrNormH = 527.81 + initNorm = 14.419 + SC1 = 8.6418e-05 + SC2 = 8.6418e-06 + SC3 = 2.2002e-03 W_density = 0.1000 H_density = 0.1000 iter = 14 - elapsed = 2.7100e-03 - tTotal = 0.1250 + elapsed = 8.5300e-04 + tTotal = 0.036066 subIterW = 0 subIterH = 0 - relError = 0.4825 + relError = 0.4866 toshow = scalar structure containing the fields: - initGrNormW = 346.84 - initGrNormH = 520.85 - initNorm = 13.825 - SC1 = 5.0666e-05 - SC2 = 5.0666e-06 - SC3 = 1.1886e-03 + initGrNormW = 369.08 + initGrNormH = 527.81 + initNorm = 14.419 + SC1 = 5.6831e-05 + SC2 = 5.6831e-06 + SC3 = 1.0400e-03 W_density = 0.1000 H_density = 0.1000 iter = 15 - elapsed = 2.8760e-03 - tTotal = 0.1279 + elapsed = 8.5800e-04 + tTotal = 0.036924 subIterW = 0 subIterH = 0 - relError = 0.4825 + relError = 0.4866 toshow = scalar structure containing the fields: - initGrNormW = 346.84 - initGrNormH = 520.85 - initNorm = 13.825 - SC1 = 3.2705e-05 - SC2 = 3.2705e-06 - SC3 = 6.5764e-04 + initGrNormW = 369.08 + initGrNormH = 527.81 + initNorm = 14.419 + SC1 = 3.7374e-05 + SC2 = 3.7374e-06 + SC3 = 7.5467e-04 W_density = 0.1000 H_density = 0.1000 iter = 16 - elapsed = 2.7360e-03 - tTotal = 0.1306 + elapsed = 8.6600e-04 + tTotal = 0.037790 subIterW = 0 subIterH = 0 - relError = 0.4824 + relError = 0.4866 toshow = scalar structure containing the fields: - initGrNormW = 346.84 - initGrNormH = 520.85 - initNorm = 13.825 - SC1 = 2.1112e-05 - SC2 = 2.1112e-06 - SC3 = 4.4025e-04 + initGrNormW = 369.08 + initGrNormH = 527.81 + initNorm = 14.419 + SC1 = 2.4580e-05 + SC2 = 2.4580e-06 + SC3 = 5.3308e-04 W_density = 0.1000 H_density = 0.1000 iter = 17 - elapsed = 2.7720e-03 - tTotal = 0.1334 + elapsed = 8.5000e-04 + tTotal = 0.038640 subIterW = 0 subIterH = 0 - relError = 0.4824 + relError = 0.4866 toshow = scalar structure containing the fields: - initGrNormW = 346.84 - initGrNormH = 520.85 - initNorm = 13.825 - SC1 = 1.3629e-05 - SC2 = 1.3629e-06 - SC3 = 3.0693e-04 + initGrNormW = 369.08 + initGrNormH = 527.81 + initNorm = 14.419 + SC1 = 1.6165e-05 + SC2 = 1.6165e-06 + SC3 = 3.2641e-04 W_density = 0.1000 H_density = 0.1000 iter = 18 - elapsed = 2.7290e-03 - tTotal = 0.1361 + elapsed = 8.0900e-04 + tTotal = 0.039449 subIterW = 0 subIterH = 0 - relError = 0.4824 + relError = 0.4866 toshow = scalar structure containing the fields: - initGrNormW = 346.84 - initGrNormH = 520.85 - initNorm = 13.825 - SC1 = 8.7978e-06 - SC2 = 8.7978e-07 - SC3 = 2.2516e-04 + initGrNormW = 369.08 + initGrNormH = 527.81 + initNorm = 14.419 + SC1 = 1.0632e-05 + SC2 = 1.0632e-06 + SC3 = 2.3945e-04 W_density = 0.1000 H_density = 0.1000 iter = 19 - elapsed = 2.7570e-03 - tTotal = 0.1388 + elapsed = 8.5300e-04 + tTotal = 0.040302 subIterW = 0 subIterH = 0 - relError = 0.4824 + relError = 0.4866 toshow = scalar structure containing the fields: - initGrNormW = 346.84 - initGrNormH = 520.85 - initNorm = 13.825 - SC1 = 5.6794e-06 - SC2 = 5.6794e-07 - SC3 = 1.2299e-04 + initGrNormW = 369.08 + initGrNormH = 527.81 + initNorm = 14.419 + SC1 = 6.9922e-06 + SC2 = 6.9922e-07 + SC3 = 1.4119e-04 W_density = 0.1000 H_density = 0.1000 iter = 20 - elapsed = 2.6070e-03 - tTotal = 0.1415 + elapsed = 9.0400e-04 + tTotal = 0.041206 subIterW = 0 subIterH = 0 - relError = 0.4824 + relError = 0.4866 toshow = scalar structure containing the fields: - initGrNormW = 346.84 - initGrNormH = 520.85 - initNorm = 13.825 - SC1 = 3.6663e-06 - SC2 = 3.6663e-07 - SC3 = 8.6009e-05 + initGrNormW = 369.08 + initGrNormH = 527.81 + initNorm = 14.419 + SC1 = 4.5987e-06 + SC2 = 4.5987e-07 + SC3 = 1.0771e-04 W_density = 0.1000 H_density = 0.1000 iter = 21 - elapsed = 2.9550e-03 - tTotal = 0.1444 + elapsed = 8.8700e-04 + tTotal = 0.042093 subIterW = 0 subIterH = 0 - relError = 0.4824 + relError = 0.4866 toshow = scalar structure containing the fields: - initGrNormW = 346.84 - initGrNormH = 520.85 - initNorm = 13.825 - SC1 = 2.3668e-06 - SC2 = 2.3668e-07 - SC3 = 4.9354e-05 + initGrNormW = 369.08 + initGrNormH = 527.81 + initNorm = 14.419 + SC1 = 3.0245e-06 + SC2 = 3.0245e-07 + SC3 = 6.8118e-05 W_density = 0.1000 H_density = 0.1000 iter = 22 - elapsed = 2.8680e-03 - tTotal = 0.1473 + elapsed = 9.7800e-04 + tTotal = 0.043071 subIterW = 0 subIterH = 0 - relError = 0.4824 + relError = 0.4866 toshow = scalar structure containing the fields: - initGrNormW = 346.84 - initGrNormH = 520.85 - initNorm = 13.825 - SC1 = 1.5279e-06 - SC2 = 1.5279e-07 - SC3 = 3.5843e-05 + initGrNormW = 369.08 + initGrNormH = 527.81 + initNorm = 14.419 + SC1 = 1.9892e-06 + SC2 = 1.9892e-07 + SC3 = 5.0644e-05 W_density = 0.1000 H_density = 0.1000 iter = 23 - elapsed = 2.7660e-03 - tTotal = 0.1500 + elapsed = 9.5800e-04 + tTotal = 0.044029 subIterW = 0 subIterH = 0 - relError = 0.4824 + relError = 0.4866 Stop: tolerance reached. final = @@ -4793,19 +4859,19 @@ scalar structure containing the fields: iterations = 23 - elapsed_total = 0.1500 - relative_error = 0.4824 + elapsed_total = 0.044029 + relative_error = 0.4866 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 = 377.29 - initGrNormH = 436.44 - initNorm = 13.825 - SC1 = 1.1538 - SC2 = 576.91 - SC3 = 0.4984 + initGrNormW = 410.83 + initGrNormH = 489.11 + initNorm = 14.419 + SC1 = 1.2775 + SC2 = 638.76 + SC3 = 0.5149 W_density = 1 H_density = 1 par = @@ -4816,8 +4882,8 @@ n = 20 type = plain nnls_solver = bp - alpha = 0.4885 - beta = 0.4885 + alpha = 0.5093 + beta = 0.5093 max_iter = 100 min_iter = 20 max_time = 1000 @@ -4829,438 +4895,438 @@ scalar structure containing the fields: - initGrNormW = 377.29 - initGrNormH = 436.44 - initNorm = 13.825 - SC1 = 0.014029 - SC2 = 0.010185 - SC3 = 0.5606 - W_density = 0.7500 + initGrNormW = 410.83 + initGrNormH = 489.11 + initNorm = 14.419 + SC1 = 0.014182 + SC2 = 0.010325 + SC3 = 0.5495 + W_density = 0.7400 H_density = 0.5800 iter = 1 - elapsed = 0.027419 - tTotal = 0.027419 - subIterW = 48 - subIterH = 29 - relError = 0.3328 + elapsed = 7.6880e-03 + tTotal = 7.6880e-03 + subIterW = 34 + subIterH = 34 + relError = 0.3413 toshow = scalar structure containing the fields: - initGrNormW = 377.29 - initGrNormH = 436.44 - initNorm = 13.825 - SC1 = 8.2543e-03 - SC2 = 5.7120e-03 - SC3 = 0.3201 - W_density = 0.7233 - H_density = 0.6100 + initGrNormW = 410.83 + initGrNormH = 489.11 + initNorm = 14.419 + SC1 = 9.9753e-03 + SC2 = 6.8830e-03 + SC3 = 0.3979 + W_density = 0.7100 + H_density = 0.6150 iter = 2 - elapsed = 0.025289 - tTotal = 0.052708 - subIterW = 19 - subIterH = 20 - relError = 0.2907 + elapsed = 7.4830e-03 + tTotal = 0.015171 + subIterW = 24 + subIterH = 23 + relError = 0.2860 toshow = scalar structure containing the fields: - initGrNormW = 377.29 - initGrNormH = 436.44 - initNorm = 13.825 - SC1 = 5.1804e-03 - SC2 = 3.5330e-03 - SC3 = 0.2440 - W_density = 0.7267 - H_density = 0.6000 + initGrNormW = 410.83 + initGrNormH = 489.11 + initNorm = 14.419 + SC1 = 6.7415e-03 + SC2 = 4.5977e-03 + SC3 = 0.2818 + W_density = 0.7167 + H_density = 0.6100 iter = 3 - elapsed = 0.022313 - tTotal = 0.075021 - subIterW = 8 - subIterH = 17 - relError = 0.2743 + elapsed = 5.9780e-03 + tTotal = 0.021149 + subIterW = 14 + subIterH = 14 + relError = 0.2626 toshow = scalar structure containing the fields: - initGrNormW = 377.29 - initGrNormH = 436.44 - initNorm = 13.825 - SC1 = 4.1872e-03 - SC2 = 2.8473e-03 - SC3 = 0.2028 - W_density = 0.7233 - H_density = 0.6050 + initGrNormW = 410.83 + initGrNormH = 489.11 + initNorm = 14.419 + SC1 = 4.9119e-03 + SC2 = 3.3499e-03 + SC3 = 0.2102 + W_density = 0.7267 + H_density = 0.5950 iter = 4 - elapsed = 0.017024 - tTotal = 0.092045 + elapsed = 4.9160e-03 + tTotal = 0.026065 subIterW = 9 - subIterH = 7 - relError = 0.2665 + subIterH = 10 + relError = 0.2518 toshow = scalar structure containing the fields: - initGrNormW = 377.29 - initGrNormH = 436.44 - initNorm = 13.825 - SC1 = 3.8662e-03 - SC2 = 2.5904e-03 - SC3 = 0.1834 - W_density = 0.7167 - H_density = 0.5950 + initGrNormW = 410.83 + initGrNormH = 489.11 + initNorm = 14.419 + SC1 = 3.8437e-03 + SC2 = 2.6291e-03 + SC3 = 0.1693 + W_density = 0.7233 + H_density = 0.5850 iter = 5 - elapsed = 0.017034 - tTotal = 0.1091 - subIterW = 12 + elapsed = 4.0250e-03 + tTotal = 0.030090 + subIterW = 4 subIterH = 6 - relError = 0.2610 + relError = 0.2459 toshow = scalar structure containing the fields: - initGrNormW = 377.29 - initGrNormH = 436.44 - initNorm = 13.825 - SC1 = 3.6494e-03 - SC2 = 2.4159e-03 - SC3 = 0.1650 - W_density = 0.6967 - H_density = 0.5950 + initGrNormW = 410.83 + initGrNormH = 489.11 + initNorm = 14.419 + SC1 = 3.0627e-03 + SC2 = 2.0398e-03 + SC3 = 0.1459 + W_density = 0.7100 + H_density = 0.5800 iter = 6 - elapsed = 0.013682 - tTotal = 0.1228 - subIterW = 8 - subIterH = 4 - relError = 0.2566 + elapsed = 4.1430e-03 + tTotal = 0.034233 + subIterW = 6 + subIterH = 8 + relError = 0.2420 toshow = scalar structure containing the fields: - initGrNormW = 377.29 - initGrNormH = 436.44 - initNorm = 13.825 - SC1 = 3.1738e-03 - SC2 = 2.1709e-03 - SC3 = 0.1523 - W_density = 0.7233 - H_density = 0.5850 + initGrNormW = 410.83 + initGrNormH = 489.11 + initNorm = 14.419 + SC1 = 2.3471e-03 + SC2 = 1.5632e-03 + SC3 = 0.1178 + W_density = 0.7167 + H_density = 0.5800 iter = 7 - elapsed = 0.014310 - tTotal = 0.1371 - subIterW = 10 - subIterH = 6 - relError = 0.2530 + elapsed = 4.6820e-03 + tTotal = 0.038915 + subIterW = 9 + subIterH = 5 + relError = 0.2398 toshow = scalar structure containing the fields: - initGrNormW = 377.29 - initGrNormH = 436.44 - initNorm = 13.825 - SC1 = 3.0390e-03 - SC2 = 2.0726e-03 - SC3 = 0.1519 - W_density = 0.7233 - H_density = 0.6200 + initGrNormW = 410.83 + initGrNormH = 489.11 + initNorm = 14.419 + SC1 = 1.8392e-03 + SC2 = 1.2175e-03 + SC3 = 0.092096 + W_density = 0.7100 + H_density = 0.5850 iter = 8 - elapsed = 0.016682 - tTotal = 0.1538 - subIterW = 9 - subIterH = 10 - relError = 0.2501 + elapsed = 3.8440e-03 + tTotal = 0.042759 + subIterW = 4 + subIterH = 2 + relError = 0.2385 toshow = scalar structure containing the fields: - initGrNormW = 377.29 - initGrNormH = 436.44 - initNorm = 13.825 - SC1 = 2.7761e-03 - SC2 = 1.8711e-03 - SC3 = 0.1337 - W_density = 0.7133 - H_density = 0.6050 + initGrNormW = 410.83 + initGrNormH = 489.11 + initNorm = 14.419 + SC1 = 1.5657e-03 + SC2 = 1.0365e-03 + SC3 = 0.081747 + W_density = 0.7100 + H_density = 0.5850 iter = 9 - elapsed = 0.014105 - tTotal = 0.1679 - subIterW = 8 - subIterH = 5 - relError = 0.2477 + elapsed = 3.8120e-03 + tTotal = 0.046571 + subIterW = 3 + subIterH = 2 + relError = 0.2376 toshow = scalar structure containing the fields: - initGrNormW = 377.29 - initGrNormH = 436.44 - initNorm = 13.825 - SC1 = 2.4883e-03 - SC2 = 1.6920e-03 - SC3 = 0.1147 - W_density = 0.7200 - H_density = 0.6000 + initGrNormW = 410.83 + initGrNormH = 489.11 + initNorm = 14.419 + SC1 = 1.4339e-03 + SC2 = 9.4636e-04 + SC3 = 0.074722 + W_density = 0.7033 + H_density = 0.5950 iter = 10 - elapsed = 0.013273 - tTotal = 0.1811 - subIterW = 4 - subIterH = 5 - relError = 0.2457 + elapsed = 3.8520e-03 + tTotal = 0.050423 + subIterW = 3 + subIterH = 3 + relError = 0.2370 toshow = scalar structure containing the fields: - initGrNormW = 377.29 - initGrNormH = 436.44 - initNorm = 13.825 - SC1 = 2.3474e-03 - SC2 = 1.6009e-03 - SC3 = 0.1040 - W_density = 0.7167 - H_density = 0.6100 + initGrNormW = 410.83 + initGrNormH = 489.11 + initNorm = 14.419 + SC1 = 1.3718e-03 + SC2 = 9.0262e-04 + SC3 = 0.037152 + W_density = 0.7033 + H_density = 0.5900 iter = 11 - elapsed = 0.013626 - tTotal = 0.1948 - subIterW = 3 - subIterH = 4 - relError = 0.2440 + elapsed = 2.9850e-03 + tTotal = 0.053408 + subIterW = 0 + subIterH = 3 + relError = 0.2364 toshow = scalar structure containing the fields: - initGrNormW = 377.29 - initGrNormH = 436.44 - initNorm = 13.825 - SC1 = 2.2544e-03 - SC2 = 1.5240e-03 - SC3 = 0.098190 - W_density = 0.7100 - H_density = 0.6200 + initGrNormW = 410.83 + initGrNormH = 489.11 + initNorm = 14.419 + SC1 = 1.3540e-03 + SC2 = 8.8822e-04 + SC3 = 0.067598 + W_density = 0.7000 + H_density = 0.5900 iter = 12 - elapsed = 0.013427 - tTotal = 0.2082 - subIterW = 6 - subIterH = 3 - relError = 0.2424 + elapsed = 2.7970e-03 + tTotal = 0.056205 + subIterW = 0 + subIterH = 0 + relError = 0.2359 toshow = scalar structure containing the fields: - initGrNormW = 377.29 - initGrNormH = 436.44 - initNorm = 13.825 - SC1 = 1.8523e-03 - SC2 = 1.2448e-03 - SC3 = 0.092645 - W_density = 0.7200 - H_density = 0.6000 + initGrNormW = 410.83 + initGrNormH = 489.11 + initNorm = 14.419 + SC1 = 1.3390e-03 + SC2 = 8.7570e-04 + SC3 = 0.064657 + W_density = 0.6933 + H_density = 0.5900 iter = 13 - elapsed = 0.013863 - tTotal = 0.2220 - subIterW = 3 - subIterH = 6 - relError = 0.2410 + elapsed = 3.0780e-03 + tTotal = 0.059283 + subIterW = 4 + subIterH = 0 + relError = 0.2354 toshow = scalar structure containing the fields: - initGrNormW = 377.29 - initGrNormH = 436.44 - initNorm = 13.825 - SC1 = 1.6393e-03 - SC2 = 1.1115e-03 - SC3 = 0.083458 - W_density = 0.7267 - H_density = 0.5950 + initGrNormW = 410.83 + initGrNormH = 489.11 + initNorm = 14.419 + SC1 = 1.2874e-03 + SC2 = 8.3941e-04 + SC3 = 0.062041 + W_density = 0.6933 + H_density = 0.5850 iter = 14 - elapsed = 0.012723 - tTotal = 0.2348 - subIterW = 2 - subIterH = 0 - relError = 0.2399 + elapsed = 3.9110e-03 + tTotal = 0.063194 + subIterW = 4 + subIterH = 3 + relError = 0.2349 toshow = scalar structure containing the fields: - initGrNormW = 377.29 - initGrNormH = 436.44 - initNorm = 13.825 - SC1 = 1.5227e-03 - SC2 = 1.0233e-03 - SC3 = 0.080874 - W_density = 0.7267 - H_density = 0.5900 + initGrNormW = 410.83 + initGrNormH = 489.11 + initNorm = 14.419 + SC1 = 1.1575e-03 + SC2 = 7.4545e-04 + SC3 = 0.057019 + W_density = 0.6933 + H_density = 0.5650 iter = 15 - elapsed = 0.013427 - tTotal = 0.2482 - subIterW = 2 - subIterH = 5 - relError = 0.2389 + elapsed = 3.9770e-03 + tTotal = 0.067171 + subIterW = 6 + subIterH = 4 + relError = 0.2345 toshow = scalar structure containing the fields: - initGrNormW = 377.29 - initGrNormH = 436.44 - initNorm = 13.825 - SC1 = 1.3882e-03 - SC2 = 9.2729e-04 - SC3 = 0.074085 - W_density = 0.7200 - H_density = 0.5900 + initGrNormW = 410.83 + initGrNormH = 489.11 + initNorm = 14.419 + SC1 = 1.0534e-03 + SC2 = 6.8051e-04 + SC3 = 0.052518 + W_density = 0.7000 + H_density = 0.5650 iter = 16 - elapsed = 0.010462 - tTotal = 0.2587 - subIterW = 6 - subIterH = 0 - relError = 0.2381 + elapsed = 3.8020e-03 + tTotal = 0.070973 + subIterW = 2 + subIterH = 2 + relError = 0.2342 toshow = scalar structure containing the fields: - initGrNormW = 377.29 - initGrNormH = 436.44 - initNorm = 13.825 - SC1 = 1.2570e-03 - SC2 = 8.3718e-04 - SC3 = 0.067501 - W_density = 0.7167 - H_density = 0.5900 + initGrNormW = 410.83 + initGrNormH = 489.11 + initNorm = 14.419 + SC1 = 9.8555e-04 + SC2 = 6.3667e-04 + SC3 = 0.049679 + W_density = 0.7000 + H_density = 0.5550 iter = 17 - elapsed = 0.012617 - tTotal = 0.2713 - subIterW = 0 + elapsed = 3.8360e-03 + tTotal = 0.074809 + subIterW = 2 subIterH = 2 - relError = 0.2374 + relError = 0.2340 toshow = scalar structure containing the fields: - initGrNormW = 377.29 - initGrNormH = 436.44 - initNorm = 13.825 - SC1 = 1.1688e-03 - SC2 = 7.8312e-04 - SC3 = 0.062106 - W_density = 0.7200 - H_density = 0.5850 + initGrNormW = 410.83 + initGrNormH = 489.11 + initNorm = 14.419 + SC1 = 9.3437e-04 + SC2 = 6.0734e-04 + SC3 = 0.046923 + W_density = 0.7033 + H_density = 0.5600 iter = 18 - elapsed = 0.012946 - tTotal = 0.2842 - subIterW = 0 + elapsed = 3.7000e-03 + tTotal = 0.078509 + subIterW = 3 subIterH = 0 - relError = 0.2368 + relError = 0.2338 toshow = scalar structure containing the fields: - initGrNormW = 377.29 - initGrNormH = 436.44 - initNorm = 13.825 - SC1 = 1.0995e-03 - SC2 = 7.4106e-04 - SC3 = 0.058980 - W_density = 0.7300 - H_density = 0.5850 + initGrNormW = 410.83 + initGrNormH = 489.11 + initNorm = 14.419 + SC1 = 9.0922e-04 + SC2 = 5.9281e-04 + SC3 = 0.046018 + W_density = 0.7067 + H_density = 0.5650 iter = 19 - elapsed = 0.013063 - tTotal = 0.2973 - subIterW = 3 - subIterH = 4 - relError = 0.2363 + elapsed = 3.5720e-03 + tTotal = 0.082081 + subIterW = 0 + subIterH = 0 + relError = 0.2336 toshow = scalar structure containing the fields: - initGrNormW = 377.29 - initGrNormH = 436.44 - initNorm = 13.825 - SC1 = 1.0238e-03 - SC2 = 6.9007e-04 - SC3 = 0.056245 - W_density = 0.7333 - H_density = 0.5800 + initGrNormW = 410.83 + initGrNormH = 489.11 + initNorm = 14.419 + SC1 = 9.1115e-04 + SC2 = 5.9407e-04 + SC3 = 0.045435 + W_density = 0.7033 + H_density = 0.5650 iter = 20 - elapsed = 0.012805 - tTotal = 0.3101 - subIterW = 3 - subIterH = 3 - relError = 0.2358 + elapsed = 3.6730e-03 + tTotal = 0.085754 + subIterW = 0 + subIterH = 2 + relError = 0.2334 toshow = scalar structure containing the fields: - initGrNormW = 377.29 - initGrNormH = 436.44 - initNorm = 13.825 - SC1 = 9.6623e-04 - SC2 = 6.5317e-04 - SC3 = 0.052798 - W_density = 0.7367 - H_density = 0.5850 + initGrNormW = 410.83 + initGrNormH = 489.11 + initNorm = 14.419 + SC1 = 9.1770e-04 + SC2 = 6.0017e-04 + SC3 = 0.045799 + W_density = 0.7067 + H_density = 0.5700 iter = 21 - elapsed = 0.012879 - tTotal = 0.3230 - subIterW = 3 - subIterH = 2 - relError = 0.2354 + elapsed = 3.5470e-03 + tTotal = 0.089301 + subIterW = 0 + subIterH = 0 + relError = 0.2332 toshow = scalar structure containing the fields: - initGrNormW = 377.29 - initGrNormH = 436.44 - initNorm = 13.825 - SC1 = 9.2471e-04 - SC2 = 6.2325e-04 - SC3 = 0.025006 - W_density = 0.7367 - H_density = 0.5800 + initGrNormW = 410.83 + initGrNormH = 489.11 + initNorm = 14.419 + SC1 = 9.3026e-04 + SC2 = 6.0467e-04 + SC3 = 0.046516 + W_density = 0.7033 + H_density = 0.5700 iter = 22 - elapsed = 9.9100e-03 - tTotal = 0.3329 + elapsed = 3.7540e-03 + tTotal = 0.093055 subIterW = 0 - subIterH = 0 - relError = 0.2351 + subIterH = 2 + relError = 0.2330 toshow = scalar structure containing the fields: - initGrNormW = 377.29 - initGrNormH = 436.44 - initNorm = 13.825 - SC1 = 8.6271e-04 - SC2 = 5.7456e-04 - SC3 = 0.047487 - W_density = 0.7267 - H_density = 0.5750 + initGrNormW = 410.83 + initGrNormH = 489.11 + initNorm = 14.419 + SC1 = 9.3931e-04 + SC2 = 6.1055e-04 + SC3 = 0.047361 + W_density = 0.7067 + H_density = 0.5650 iter = 23 - elapsed = 0.013731 - tTotal = 0.3466 - subIterW = 4 - subIterH = 3 - relError = 0.2347 + elapsed = 3.6410e-03 + tTotal = 0.096696 + subIterW = 0 + subIterH = 0 + relError = 0.2327 Stop: tolerance reached. final = @@ -5268,10 +5334,10 @@ scalar structure containing the fields: iterations = 23 - elapsed_total = 0.3466 - relative_error = 0.2347 - W_density = 0.7267 - H_density = 0.5750 + elapsed_total = 0.096696 + relative_error = 0.2327 + W_density = 0.7067 + H_density = 0.5650 ***** demo m = 300; @@ -5319,24 +5385,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_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 @@ -5352,6 +5400,39 @@ [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) @@ -5361,13 +5442,6 @@ 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 @@ -5398,12 +5472,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/1344041/tmp/hooks/B01_cleanup starting +I: user script /srv/workspace/pbuilder/1344041/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/3091070 and its subdirectories -I: Current time: Fri Mar 21 14:10:23 -12 2025 -I: pbuilder-time-stamp: 1742609423 +I: removing directory /srv/workspace/pbuilder/1344041 and its subdirectories +I: Current time: Fri Apr 24 22:35:59 +14 2026 +I: pbuilder-time-stamp: 1777019759