Diff of the two buildlogs: -- --- b1/build.log 2024-12-11 17:53:46.543642847 +0000 +++ b2/build.log 2024-12-11 17:59:46.619231549 +0000 @@ -1,6 +1,6 @@ I: pbuilder: network access will be disabled during build -I: Current time: Wed Dec 11 05:47:46 -12 2024 -I: pbuilder-time-stamp: 1733939266 +I: Current time: Thu Dec 12 07:53:55 +14 2024 +I: pbuilder-time-stamp: 1733939635 I: Building the build Environment I: extracting base tarball [/var/cache/pbuilder/trixie-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/22840/tmp/hooks/D02_print_environment starting +I: user script /srv/workspace/pbuilder/12084/tmp/hooks/D01_modify_environment starting +debug: Running on virt32c. +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 Dec 11 17:54 /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/12084/tmp/hooks/D01_modify_environment finished +I: user script /srv/workspace/pbuilder/12084/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='armhf' - DEBIAN_FRONTEND='noninteractive' - DEB_BUILD_OPTIONS='buildinfo=+all reproducible=+all parallel=5 ' - DISTRIBUTION='trixie' - HOME='/root' - HOST_ARCH='armhf' + 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]="arm-unknown-linux-gnueabihf") + 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=armhf + DEBIAN_FRONTEND=noninteractive + DEB_BUILD_OPTIONS='buildinfo=+all reproducible=+all parallel=4 ' + DIRSTACK=() + DISTRIBUTION=trixie + EUID=0 + FUNCNAME=([0]="Echo" [1]="main") + GROUPS=() + HOME=/root + HOSTNAME=i-capture-the-hostname + HOSTTYPE=arm + HOST_ARCH=armhf IFS=' ' - INVOCATION_ID='9753bbe545474550ab943e3d74a2bad3' - 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='22840' - PS1='# ' - PS2='> ' + INVOCATION_ID=7d7d1b7ebe484a7c9d7f1a75c014aa51 + LANG=C + LANGUAGE=it_CH:it + LC_ALL=C + MACHTYPE=arm-unknown-linux-gnueabihf + MAIL=/var/mail/root + OPTERR=1 + OPTIND=1 + OSTYPE=linux-gnueabihf + 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=12084 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.NLLmrxx4/pbuilderrc_P21u --distribution trixie --hookdir /etc/pbuilder/first-build-hooks --debbuildopts -b --basetgz /var/cache/pbuilder/trixie-reproducible-base.tgz --buildresult /srv/reproducible-results/rbuild-debian/r-b-build.NLLmrxx4/b1 --logfile b1/build.log octave-linear-algebra_2.2.3-4.dsc' - SUDO_GID='114' - SUDO_UID='109' - SUDO_USER='jenkins' - TERM='unknown' - TZ='/usr/share/zoneinfo/Etc/GMT+12' - USER='root' - _='/usr/bin/systemd-run' - http_proxy='http://10.0.0.15:3142/' + 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.NLLmrxx4/pbuilderrc_Qm9E --distribution trixie --hookdir /etc/pbuilder/rebuild-hooks --debbuildopts -b --basetgz /var/cache/pbuilder/trixie-reproducible-base.tgz --buildresult /srv/reproducible-results/rbuild-debian/r-b-build.NLLmrxx4/b2 --logfile b2/build.log octave-linear-algebra_2.2.3-4.dsc' + SUDO_GID=113 + SUDO_UID=107 + SUDO_USER=jenkins + TERM=unknown + TZ=/usr/share/zoneinfo/Etc/GMT-14 + UID=0 + USER=root + _='I: set' + http_proxy=http://10.0.0.15:3142/ I: uname -a - Linux ff64a 6.1.0-28-arm64 #1 SMP Debian 6.1.119-1 (2024-11-22) aarch64 GNU/Linux + Linux i-capture-the-hostname 6.1.0-28-armmp-lpae #1 SMP Debian 6.1.119-1 (2024-11-22) armv7l GNU/Linux I: ls -l /bin lrwxrwxrwx 1 root root 7 Nov 22 14:40 /bin -> usr/bin -I: user script /srv/workspace/pbuilder/22840/tmp/hooks/D02_print_environment finished +I: user script /srv/workspace/pbuilder/12084/tmp/hooks/D02_print_environment finished -> Attempting to satisfy build-dependencies -> Creating pbuilder-satisfydepends-dummy package Package: pbuilder-satisfydepends-dummy @@ -630,7 +662,7 @@ Get: 519 http://deb.debian.org/debian trixie/main armhf gfortran armhf 4:14.2.0-1 [1428 B] Get: 520 http://deb.debian.org/debian trixie/main armhf octave-dev armhf 9.2.0-3+b1 [995 kB] Get: 521 http://deb.debian.org/debian trixie/main armhf dh-octave all 1.8.0 [22.7 kB] -Fetched 175 MB in 5s (34.4 MB/s) +Fetched 175 MB in 4s (42.8 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 ... 19680 files and directories currently installed.) @@ -2767,7 +2799,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/12084/tmp/hooks/A99_set_merged_usr starting +Not re-configuring usrmerge for trixie +I: user script /srv/workspace/pbuilder/12084/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 @@ -2805,55 +2841,43 @@ Checking package... Run the unit tests... Checking m files ... -[inst/cartprod.m] ->>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/cartprod.m -***** assert(cartprod(1:2,0:1),[1 0; 2 0; 1 1; 2 1]) -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); +[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.EhE5KowJE7 at line 22 column 31 + /tmp/tmp.Wd8zUnTx0L at line 14 column 31 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/smwsolve.m] ->>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/smwsolve.m +[inst/funm.m] +>>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/funm.m +***** function b = fsin (a) + b = sin (a); +***** endfunction ***** test - A = 2.1*eye (10); - u = rand (10, 2); u /= diag (norm (u, "cols")); - v = rand (10, 2); v /= diag (norm (v, "cols")); - b = rand (10, 2); - x1 = (A + u*v') \ b; - x2 = smwsolve (A, u, v, b); - assert (x1, x2, 1e-13); -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); + 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 +[inst/cod.m] +>>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/cod.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 + a = rand (5, 10); + [q, r, z] = cod (a); + assert (norm (q*r*z' - a) / norm (a) < 1e-10); +***** test + a = rand (5, 10) + i * rand (5, 10); + [q, r, z] = cod (a); + assert (norm (q*r*z' - a) / norm (a) < 1e-10); +***** test + a = rand (5, 10); + [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/nmf_bpas.m] >>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/nmf_bpas.m ***** shared m, n, k, A @@ -2866,12 +2890,12 @@ Stop: tolerance reached. ***** test [W,H,iter,HIS]=nmf_bpas(A,k,'verbose'); - initGrNormW = 392.16 - initGrNormH = 478.38 - initNorm = 14.465 - SC1 = 1.2371 - SC2 = 618.57 - SC3 = 0.5056 + initGrNormW = 395.66 + initGrNormH = 464.40 + initNorm = 14.111 + SC1 = 1.2202 + SC2 = 610.10 + SC3 = 0.5024 W_density = 1 H_density = 1 par = @@ -2882,8 +2906,8 @@ n = 20 type = regularized nnls_solver = bp - alpha = 0.5165 - beta = 0.5165 + alpha = 0.4933 + beta = 0.4933 max_iter = 100 min_iter = 20 max_time = 1000 @@ -2895,438 +2919,438 @@ scalar structure containing the fields: - initGrNormW = 392.16 - initGrNormH = 478.38 - initNorm = 14.465 - SC1 = 0.030706 - SC2 = 0.030030 - SC3 = 2.4511 - W_density = 0.9733 - H_density = 0.6600 + initGrNormW = 395.66 + initGrNormH = 464.40 + initNorm = 14.111 + SC1 = 0.028103 + SC2 = 0.027035 + SC3 = 2.1188 + W_density = 0.9533 + H_density = 0.6300 iter = 1 - elapsed = 0.042582 - tTotal = 0.042582 - subIterW = 6 - subIterH = 31 - relError = 0.4266 + elapsed = 0.035888 + tTotal = 0.035888 + subIterW = 5 + subIterH = 35 + relError = 0.4524 toshow = scalar structure containing the fields: - initGrNormW = 392.16 - initGrNormH = 478.38 - initNorm = 14.465 - SC1 = 0.016865 - SC2 = 0.016258 - SC3 = 1.3089 - W_density = 0.9467 - H_density = 0.5600 + initGrNormW = 395.66 + initGrNormH = 464.40 + initNorm = 14.111 + SC1 = 0.016373 + SC2 = 0.015260 + SC3 = 1.1781 + W_density = 0.9233 + H_density = 0.5450 iter = 2 - elapsed = 0.052858 - tTotal = 0.095440 - subIterW = 8 - subIterH = 23 - relError = 0.3346 + elapsed = 0.027246 + tTotal = 0.063134 + subIterW = 11 + subIterH = 19 + relError = 0.3642 toshow = scalar structure containing the fields: - initGrNormW = 392.16 - initGrNormH = 478.38 - initNorm = 14.465 - SC1 = 0.010798 - SC2 = 0.010214 - SC3 = 0.8254 - W_density = 0.9167 - H_density = 0.6150 + initGrNormW = 395.66 + initGrNormH = 464.40 + initNorm = 14.111 + SC1 = 0.010149 + SC2 = 9.2156e-03 + SC3 = 0.7400 + W_density = 0.8900 + H_density = 0.5950 iter = 3 - elapsed = 0.045884 - tTotal = 0.1413 - subIterW = 8 - subIterH = 20 - relError = 0.2985 + elapsed = 0.025016 + tTotal = 0.088150 + subIterW = 10 + subIterH = 21 + relError = 0.3208 toshow = scalar structure containing the fields: - initGrNormW = 392.16 - initGrNormH = 478.38 - initNorm = 14.465 - SC1 = 7.5324e-03 - SC2 = 6.9900e-03 - SC3 = 0.5712 - W_density = 0.8967 - H_density = 0.6600 + initGrNormW = 395.66 + initGrNormH = 464.40 + initNorm = 14.111 + SC1 = 7.0967e-03 + SC2 = 6.2309e-03 + SC3 = 0.5186 + W_density = 0.8567 + H_density = 0.6350 iter = 4 - elapsed = 0.046997 - tTotal = 0.1883 - subIterW = 8 - subIterH = 12 - relError = 0.2842 + elapsed = 0.023494 + tTotal = 0.1116 + subIterW = 10 + subIterH = 14 + relError = 0.3003 toshow = scalar structure containing the fields: - initGrNormW = 392.16 - initGrNormH = 478.38 - initNorm = 14.465 - SC1 = 5.4907e-03 - SC2 = 5.0404e-03 - SC3 = 0.4131 - W_density = 0.8867 - H_density = 0.6800 + initGrNormW = 395.66 + initGrNormH = 464.40 + initNorm = 14.111 + SC1 = 5.1472e-03 + SC2 = 4.4266e-03 + SC3 = 0.3684 + W_density = 0.8300 + H_density = 0.6650 iter = 5 - elapsed = 0.035126 - tTotal = 0.2234 - subIterW = 5 - subIterH = 6 - relError = 0.2773 + elapsed = 0.025215 + tTotal = 0.1369 + subIterW = 8 + subIterH = 8 + relError = 0.2901 toshow = scalar structure containing the fields: - initGrNormW = 392.16 - initGrNormH = 478.38 - initNorm = 14.465 - SC1 = 4.1596e-03 - SC2 = 3.7770e-03 - SC3 = 0.3100 - W_density = 0.8767 - H_density = 0.6950 + initGrNormW = 395.66 + initGrNormH = 464.40 + initNorm = 14.111 + SC1 = 3.8598e-03 + SC2 = 3.2808e-03 + SC3 = 0.2716 + W_density = 0.8200 + H_density = 0.6550 iter = 6 - elapsed = 0.033449 - tTotal = 0.2569 - subIterW = 3 - subIterH = 3 - relError = 0.2735 + elapsed = 0.018662 + tTotal = 0.1555 + subIterW = 2 + subIterH = 7 + relError = 0.2841 toshow = scalar structure containing the fields: - initGrNormW = 392.16 - initGrNormH = 478.38 - initNorm = 14.465 - SC1 = 3.2847e-03 - SC2 = 2.9168e-03 - SC3 = 0.2441 - W_density = 0.8633 - H_density = 0.7100 + initGrNormW = 395.66 + initGrNormH = 464.40 + initNorm = 14.111 + SC1 = 3.0033e-03 + SC2 = 2.5168e-03 + SC3 = 0.2089 + W_density = 0.8133 + H_density = 0.6600 iter = 7 - elapsed = 0.034629 - tTotal = 0.2915 - subIterW = 4 - subIterH = 4 - relError = 0.2710 + elapsed = 0.018366 + tTotal = 0.1739 + subIterW = 5 + subIterH = 2 + relError = 0.2801 toshow = scalar structure containing the fields: - initGrNormW = 392.16 - initGrNormH = 478.38 - initNorm = 14.465 - SC1 = 2.6381e-03 - SC2 = 2.2952e-03 - SC3 = 0.1940 - W_density = 0.8467 - H_density = 0.7200 + initGrNormW = 395.66 + initGrNormH = 464.40 + initNorm = 14.111 + SC1 = 2.3976e-03 + SC2 = 1.9660e-03 + SC3 = 0.1659 + W_density = 0.8033 + H_density = 0.6700 iter = 8 - elapsed = 0.035938 - tTotal = 0.3275 - subIterW = 4 + elapsed = 0.019994 + tTotal = 0.1939 + subIterW = 7 subIterH = 6 - relError = 0.2693 + relError = 0.2773 toshow = scalar structure containing the fields: - initGrNormW = 392.16 - initGrNormH = 478.38 - initNorm = 14.465 - SC1 = 2.1463e-03 - SC2 = 1.8458e-03 - SC3 = 0.1591 - W_density = 0.8500 - H_density = 0.7350 + initGrNormW = 395.66 + initGrNormH = 464.40 + initNorm = 14.111 + SC1 = 1.9513e-03 + SC2 = 1.5610e-03 + SC3 = 0.1356 + W_density = 0.7967 + H_density = 0.6800 iter = 9 - elapsed = 0.033836 - tTotal = 0.3613 - subIterW = 0 - subIterH = 3 - relError = 0.2680 + elapsed = 0.018735 + tTotal = 0.2126 + subIterW = 2 + subIterH = 4 + relError = 0.2753 toshow = scalar structure containing the fields: - initGrNormW = 392.16 - initGrNormH = 478.38 - initNorm = 14.465 - SC1 = 1.7806e-03 - SC2 = 1.5135e-03 - SC3 = 0.1291 - W_density = 0.8367 - H_density = 0.7400 + initGrNormW = 395.66 + initGrNormH = 464.40 + initNorm = 14.111 + SC1 = 1.6205e-03 + SC2 = 1.2640e-03 + SC3 = 0.1109 + W_density = 0.7800 + H_density = 0.6850 iter = 10 - elapsed = 0.033980 - tTotal = 0.3953 - subIterW = 4 - subIterH = 0 - relError = 0.2670 + elapsed = 0.021508 + tTotal = 0.2341 + subIterW = 6 + subIterH = 4 + relError = 0.2737 toshow = scalar structure containing the fields: - initGrNormW = 392.16 - initGrNormH = 478.38 - initNorm = 14.465 - SC1 = 1.5134e-03 - SC2 = 1.2531e-03 - SC3 = 0.1109 - W_density = 0.8300 - H_density = 0.7400 + initGrNormW = 395.66 + initGrNormH = 464.40 + initNorm = 14.111 + SC1 = 1.3396e-03 + SC2 = 1.0341e-03 + SC3 = 0.090512 + W_density = 0.7733 + H_density = 0.6950 iter = 11 - elapsed = 0.035569 - tTotal = 0.4308 - subIterW = 4 - subIterH = 2 - relError = 0.2663 + elapsed = 0.022180 + tTotal = 0.2563 + subIterW = 6 + subIterH = 4 + relError = 0.2726 toshow = scalar structure containing the fields: - initGrNormW = 392.16 - initGrNormH = 478.38 - initNorm = 14.465 - SC1 = 1.2800e-03 - SC2 = 1.0496e-03 - SC3 = 0.092898 - W_density = 0.8233 - H_density = 0.7400 + initGrNormW = 395.66 + initGrNormH = 464.40 + initNorm = 14.111 + SC1 = 1.1275e-03 + SC2 = 8.6144e-04 + SC3 = 0.074172 + W_density = 0.7633 + H_density = 0.7150 iter = 12 - elapsed = 0.036211 - tTotal = 0.4671 - subIterW = 6 - subIterH = 2 - relError = 0.2656 + elapsed = 0.021607 + tTotal = 0.2779 + subIterW = 3 + subIterH = 4 + relError = 0.2716 toshow = scalar structure containing the fields: - initGrNormW = 392.16 - initGrNormH = 478.38 - initNorm = 14.465 - SC1 = 1.0978e-03 - SC2 = 8.8923e-04 - SC3 = 0.045773 - W_density = 0.8233 - H_density = 0.7400 + initGrNormW = 395.66 + initGrNormH = 464.40 + initNorm = 14.111 + SC1 = 9.4708e-04 + SC2 = 7.0841e-04 + SC3 = 0.061266 + W_density = 0.7433 + H_density = 0.7100 iter = 13 - elapsed = 0.027839 - tTotal = 0.4949 - subIterW = 0 - subIterH = 4 - relError = 0.2650 + elapsed = 0.020880 + tTotal = 0.2988 + subIterW = 5 + subIterH = 0 + relError = 0.2709 toshow = scalar structure containing the fields: - initGrNormW = 392.16 - initGrNormH = 478.38 - initNorm = 14.465 - SC1 = 9.4731e-04 - SC2 = 7.5785e-04 - SC3 = 0.067411 - W_density = 0.8067 - H_density = 0.7400 + initGrNormW = 395.66 + initGrNormH = 464.40 + initNorm = 14.111 + SC1 = 7.8917e-04 + SC2 = 5.8714e-04 + SC3 = 0.050291 + W_density = 0.7400 + H_density = 0.7100 iter = 14 - elapsed = 0.035695 - tTotal = 0.5306 - subIterW = 4 - subIterH = 2 - relError = 0.2646 + elapsed = 0.019986 + tTotal = 0.3188 + subIterW = 0 + subIterH = 4 + relError = 0.2704 toshow = scalar structure containing the fields: - initGrNormW = 392.16 - initGrNormH = 478.38 - initNorm = 14.465 - SC1 = 8.2643e-04 - SC2 = 6.5453e-04 - SC3 = 0.057809 - W_density = 0.7967 - H_density = 0.7400 + initGrNormW = 395.66 + initGrNormH = 464.40 + initNorm = 14.111 + SC1 = 6.6828e-04 + SC2 = 4.9586e-04 + SC3 = 0.026011 + W_density = 0.7400 + H_density = 0.7150 iter = 15 - elapsed = 0.027489 - tTotal = 0.5581 - subIterW = 2 + elapsed = 0.014965 + tTotal = 0.3337 + subIterW = 0 subIterH = 0 - relError = 0.2641 + relError = 0.2700 toshow = scalar structure containing the fields: - initGrNormW = 392.16 - initGrNormH = 478.38 - initNorm = 14.465 - SC1 = 7.2835e-04 - SC2 = 5.7394e-04 - SC3 = 0.050071 - W_density = 0.7933 - H_density = 0.7600 + initGrNormW = 395.66 + initGrNormH = 464.40 + initNorm = 14.111 + SC1 = 5.7617e-04 + SC2 = 4.2522e-04 + SC3 = 0.035734 + W_density = 0.7367 + H_density = 0.7150 iter = 16 - elapsed = 0.043056 - tTotal = 0.6011 - subIterW = 2 - subIterH = 6 - relError = 0.2636 + elapsed = 0.014736 + tTotal = 0.3485 + subIterW = 0 + subIterH = 0 + relError = 0.2697 toshow = scalar structure containing the fields: - initGrNormW = 392.16 - initGrNormH = 478.38 - initNorm = 14.465 - SC1 = 6.4357e-04 - SC2 = 5.0713e-04 - SC3 = 0.043181 - W_density = 0.7933 - H_density = 0.7600 + initGrNormW = 395.66 + initGrNormH = 464.40 + initNorm = 14.111 + SC1 = 5.0134e-04 + SC2 = 3.7099e-04 + SC3 = 0.030628 + W_density = 0.7367 + H_density = 0.7200 iter = 17 - elapsed = 0.035743 - tTotal = 0.6369 - subIterW = 4 - subIterH = 2 - relError = 0.2633 + elapsed = 0.018884 + tTotal = 0.3674 + subIterW = 2 + subIterH = 0 + relError = 0.2694 toshow = scalar structure containing the fields: - initGrNormW = 392.16 - initGrNormH = 478.38 - initNorm = 14.465 - SC1 = 5.7910e-04 - SC2 = 4.5402e-04 - SC3 = 0.038112 - W_density = 0.7900 - H_density = 0.7550 + initGrNormW = 395.66 + initGrNormH = 464.40 + initNorm = 14.111 + SC1 = 4.4451e-04 + SC2 = 3.2894e-04 + SC3 = 0.015361 + W_density = 0.7367 + H_density = 0.7150 iter = 18 - elapsed = 0.033919 - tTotal = 0.6708 - subIterW = 3 + elapsed = 0.014950 + tTotal = 0.3823 + subIterW = 0 subIterH = 0 - relError = 0.2629 + relError = 0.2691 toshow = scalar structure containing the fields: - initGrNormW = 392.16 - initGrNormH = 478.38 - initNorm = 14.465 - SC1 = 5.3012e-04 - SC2 = 4.1350e-04 - SC3 = 0.034803 - W_density = 0.7900 - H_density = 0.7550 + initGrNormW = 395.66 + initGrNormH = 464.40 + initNorm = 14.111 + SC1 = 3.9788e-04 + SC2 = 2.9603e-04 + SC3 = 0.023754 + W_density = 0.7400 + H_density = 0.7150 iter = 19 - elapsed = 0.026838 - tTotal = 0.6976 - subIterW = 2 - subIterH = 0 - relError = 0.2625 + elapsed = 0.021830 + tTotal = 0.4041 + subIterW = 3 + subIterH = 6 + relError = 0.2689 toshow = scalar structure containing the fields: - initGrNormW = 392.16 - initGrNormH = 478.38 - initNorm = 14.465 - SC1 = 4.8924e-04 - SC2 = 3.8258e-04 - SC3 = 0.031524 - W_density = 0.7933 - H_density = 0.7500 + initGrNormW = 395.66 + initGrNormH = 464.40 + initNorm = 14.111 + SC1 = 3.6683e-04 + SC2 = 2.7146e-04 + SC3 = 0.021912 + W_density = 0.7433 + H_density = 0.7150 iter = 20 - elapsed = 0.032450 - tTotal = 0.7301 + elapsed = 0.015608 + tTotal = 0.4197 subIterW = 0 - subIterH = 3 - relError = 0.2622 + subIterH = 0 + relError = 0.2687 toshow = scalar structure containing the fields: - initGrNormW = 392.16 - initGrNormH = 478.38 - initNorm = 14.465 - SC1 = 4.5498e-04 - SC2 = 3.5670e-04 - SC3 = 0.028829 - W_density = 0.8000 - H_density = 0.7500 + initGrNormW = 395.66 + initGrNormH = 464.40 + initNorm = 14.111 + SC1 = 3.4277e-04 + SC2 = 2.5365e-04 + SC3 = 0.020351 + W_density = 0.7467 + H_density = 0.7200 iter = 21 - elapsed = 0.033595 - tTotal = 0.7637 - subIterW = 2 - subIterH = 2 - relError = 0.2618 + elapsed = 0.018560 + tTotal = 0.4383 + subIterW = 0 + subIterH = 0 + relError = 0.2685 toshow = scalar structure containing the fields: - initGrNormW = 392.16 - initGrNormH = 478.38 - initNorm = 14.465 - SC1 = 4.1262e-04 - SC2 = 3.2019e-04 - SC3 = 0.026308 - W_density = 0.7900 - H_density = 0.7450 + initGrNormW = 395.66 + initGrNormH = 464.40 + initNorm = 14.111 + SC1 = 3.1287e-04 + SC2 = 2.3152e-04 + SC3 = 0.019287 + W_density = 0.7500 + H_density = 0.7150 iter = 22 - elapsed = 0.032667 - tTotal = 0.7963 + elapsed = 0.023829 + tTotal = 0.4621 subIterW = 3 - subIterH = 0 - relError = 0.2615 + subIterH = 2 + relError = 0.2684 toshow = scalar structure containing the fields: - initGrNormW = 392.16 - initGrNormH = 478.38 - initNorm = 14.465 - SC1 = 3.7108e-04 - SC2 = 2.8796e-04 - SC3 = 0.023438 - W_density = 0.7933 - H_density = 0.7400 + initGrNormW = 395.66 + initGrNormH = 464.40 + initNorm = 14.111 + SC1 = 2.9360e-04 + SC2 = 2.1609e-04 + SC3 = 0.018179 + W_density = 0.7433 + H_density = 0.7200 iter = 23 - elapsed = 0.033038 - tTotal = 0.8294 - subIterW = 0 + elapsed = 0.019705 + tTotal = 0.4818 + subIterW = 2 subIterH = 0 - relError = 0.2612 + relError = 0.2683 Stop: tolerance reached. final = @@ -3334,19 +3358,19 @@ scalar structure containing the fields: iterations = 23 - elapsed_total = 0.8294 - relative_error = 0.2612 - W_density = 0.7933 - H_density = 0.7400 + elapsed_total = 0.4818 + relative_error = 0.2683 + W_density = 0.7433 + H_density = 0.7200 ***** test [W,H,iter,HIS]=nmf_bpas(A,k,'verbose','nnlssolver','as'); - initGrNormW = 367.52 - initGrNormH = 439.89 - initNorm = 14.465 - SC1 = 1.1464 - SC2 = 573.21 - SC3 = 0.4948 + initGrNormW = 398.15 + initGrNormH = 457.57 + initNorm = 14.111 + SC1 = 1.2131 + SC2 = 606.55 + SC3 = 0.5017 W_density = 1 H_density = 1 par = @@ -3357,8 +3381,8 @@ n = 20 type = regularized nnls_solver = as - alpha = 0.5165 - beta = 0.5165 + alpha = 0.4933 + beta = 0.4933 max_iter = 100 min_iter = 20 max_time = 1000 @@ -3370,438 +3394,438 @@ scalar structure containing the fields: - initGrNormW = 367.52 - initGrNormH = 439.89 - initNorm = 14.465 - SC1 = 0.029299 - SC2 = 0.027834 - SC3 = 1.1242 + initGrNormW = 398.15 + initGrNormH = 457.57 + initNorm = 14.111 + SC1 = 0.028654 + SC2 = 0.027794 + SC3 = 1.1921 W_density = 0.9533 - H_density = 0.6650 + H_density = 0.6900 iter = 1 - elapsed = 0.038332 - tTotal = 0.038332 + elapsed = 0.024672 + tTotal = 0.024672 subIterW = 9 - subIterH = 29 - relError = 0.4294 + subIterH = 31 + relError = 0.4368 toshow = scalar structure containing the fields: - initGrNormW = 367.52 - initGrNormH = 439.89 - initNorm = 14.465 - SC1 = 0.014770 - SC2 = 0.013854 - SC3 = 0.5900 - W_density = 0.9200 - H_density = 0.5900 + initGrNormW = 398.15 + initGrNormH = 457.57 + initNorm = 14.111 + SC1 = 0.015416 + SC2 = 0.014461 + SC3 = 0.6271 + W_density = 0.9067 + H_density = 0.4900 iter = 2 - elapsed = 0.045572 - tTotal = 0.083904 - subIterW = 13 - subIterH = 39 - relError = 0.3399 + elapsed = 0.024653 + tTotal = 0.049325 + subIterW = 16 + subIterH = 46 + relError = 0.3383 toshow = scalar structure containing the fields: - initGrNormW = 367.52 - initGrNormH = 439.89 - initNorm = 14.465 - SC1 = 9.9067e-03 - SC2 = 9.0350e-03 - SC3 = 0.3950 - W_density = 0.8833 - H_density = 0.6350 + initGrNormW = 398.15 + initGrNormH = 457.57 + initNorm = 14.111 + SC1 = 9.9558e-03 + SC2 = 9.1195e-03 + SC3 = 0.4028 + W_density = 0.8767 + H_density = 0.5500 iter = 3 - elapsed = 0.045029 - tTotal = 0.1289 + elapsed = 0.025130 + tTotal = 0.074455 subIterW = 20 - subIterH = 32 - relError = 0.3143 + subIterH = 44 + relError = 0.3121 toshow = scalar structure containing the fields: - initGrNormW = 367.52 - initGrNormH = 439.89 - initNorm = 14.465 - SC1 = 7.2188e-03 - SC2 = 6.4680e-03 - SC3 = 0.2800 - W_density = 0.8700 - H_density = 0.6500 + initGrNormW = 398.15 + initGrNormH = 457.57 + initNorm = 14.111 + SC1 = 7.0841e-03 + SC2 = 6.2907e-03 + SC3 = 0.2906 + W_density = 0.8500 + H_density = 0.5900 iter = 4 - elapsed = 0.045327 - tTotal = 0.1743 - subIterW = 22 - subIterH = 31 - relError = 0.3010 + elapsed = 0.027238 + tTotal = 0.1017 + subIterW = 26 + subIterH = 38 + relError = 0.3000 toshow = scalar structure containing the fields: - initGrNormW = 367.52 - initGrNormH = 439.89 - initNorm = 14.465 - SC1 = 5.5033e-03 - SC2 = 4.8319e-03 - SC3 = 0.2106 - W_density = 0.8567 - H_density = 0.6600 + initGrNormW = 398.15 + initGrNormH = 457.57 + initNorm = 14.111 + SC1 = 5.2491e-03 + SC2 = 4.5457e-03 + SC3 = 0.2104 + W_density = 0.8333 + H_density = 0.6100 iter = 5 - elapsed = 0.042133 - tTotal = 0.2164 - subIterW = 24 - subIterH = 31 - relError = 0.2918 + elapsed = 0.024294 + tTotal = 0.1260 + subIterW = 28 + subIterH = 37 + relError = 0.2931 toshow = scalar structure containing the fields: - initGrNormW = 367.52 - initGrNormH = 439.89 - initNorm = 14.465 - SC1 = 4.3729e-03 - SC2 = 3.7169e-03 - SC3 = 0.1661 - W_density = 0.8533 - H_density = 0.6650 + initGrNormW = 398.15 + initGrNormH = 457.57 + initNorm = 14.111 + SC1 = 4.0886e-03 + SC2 = 3.4345e-03 + SC3 = 0.1608 + W_density = 0.8200 + H_density = 0.6400 iter = 6 - elapsed = 0.036638 - tTotal = 0.2530 - subIterW = 24 - subIterH = 30 - relError = 0.2852 + elapsed = 0.022482 + tTotal = 0.1485 + subIterW = 28 + subIterH = 35 + relError = 0.2882 toshow = scalar structure containing the fields: - initGrNormW = 367.52 - initGrNormH = 439.89 - initNorm = 14.465 - SC1 = 3.4767e-03 - SC2 = 2.9065e-03 - SC3 = 0.1314 - W_density = 0.8367 - H_density = 0.6750 + initGrNormW = 398.15 + initGrNormH = 457.57 + initNorm = 14.111 + SC1 = 3.2881e-03 + SC2 = 2.6831e-03 + SC3 = 0.1270 + W_density = 0.8133 + H_density = 0.6450 iter = 7 - elapsed = 0.042291 - tTotal = 0.2953 - subIterW = 25 - subIterH = 31 - relError = 0.2803 + elapsed = 0.018690 + tTotal = 0.1672 + subIterW = 26 + subIterH = 30 + relError = 0.2843 toshow = scalar structure containing the fields: - initGrNormW = 367.52 - initGrNormH = 439.89 - initNorm = 14.465 - SC1 = 2.7887e-03 - SC2 = 2.3034e-03 - SC3 = 0.1062 - W_density = 0.8300 - H_density = 0.6900 + initGrNormW = 398.15 + initGrNormH = 457.57 + initNorm = 14.111 + SC1 = 2.6748e-03 + SC2 = 2.1399e-03 + SC3 = 0.1056 + W_density = 0.8000 + H_density = 0.6500 iter = 8 - elapsed = 0.042802 - tTotal = 0.3381 + elapsed = 0.018930 + tTotal = 0.1861 subIterW = 27 - subIterH = 31 - relError = 0.2766 + subIterH = 29 + relError = 0.2811 toshow = scalar structure containing the fields: - initGrNormW = 367.52 - initGrNormH = 439.89 - initNorm = 14.465 - SC1 = 2.2756e-03 - SC2 = 1.8660e-03 - SC3 = 0.086512 - W_density = 0.8233 - H_density = 0.7150 + initGrNormW = 398.15 + initGrNormH = 457.57 + initNorm = 14.111 + SC1 = 2.1880e-03 + SC2 = 1.7373e-03 + SC3 = 0.084082 + W_density = 0.7933 + H_density = 0.6750 iter = 9 - elapsed = 0.039718 - tTotal = 0.3778 - subIterW = 29 - subIterH = 24 - relError = 0.2739 + elapsed = 0.021264 + tTotal = 0.2074 + subIterW = 28 + subIterH = 31 + relError = 0.2785 toshow = scalar structure containing the fields: - initGrNormW = 367.52 - initGrNormH = 439.89 - initNorm = 14.465 - SC1 = 1.8516e-03 - SC2 = 1.4998e-03 - SC3 = 0.067887 - W_density = 0.8067 - H_density = 0.7300 + initGrNormW = 398.15 + initGrNormH = 457.57 + initNorm = 14.111 + SC1 = 1.7893e-03 + SC2 = 1.4064e-03 + SC3 = 0.071645 + W_density = 0.7833 + H_density = 0.6850 iter = 10 - elapsed = 0.041817 - tTotal = 0.4197 - subIterW = 32 - subIterH = 25 - relError = 0.2718 + elapsed = 0.024672 + tTotal = 0.2320 + subIterW = 31 + subIterH = 31 + relError = 0.2766 toshow = scalar structure containing the fields: - initGrNormW = 367.52 - initGrNormH = 439.89 - initNorm = 14.465 - SC1 = 1.4921e-03 - SC2 = 1.1996e-03 - SC3 = 0.055572 - W_density = 0.8000 - H_density = 0.7300 + initGrNormW = 398.15 + initGrNormH = 457.57 + initNorm = 14.111 + SC1 = 1.4754e-03 + SC2 = 1.1596e-03 + SC3 = 0.058338 + W_density = 0.7833 + H_density = 0.7050 iter = 11 - elapsed = 0.041935 - tTotal = 0.4616 - subIterW = 33 - subIterH = 27 - relError = 0.2703 + elapsed = 0.023071 + tTotal = 0.2551 + subIterW = 29 + subIterH = 26 + relError = 0.2750 toshow = scalar structure containing the fields: - initGrNormW = 367.52 - initGrNormH = 439.89 - initNorm = 14.465 - SC1 = 1.2220e-03 - SC2 = 9.7518e-04 - SC3 = 0.045347 - W_density = 0.7967 - H_density = 0.7350 + initGrNormW = 398.15 + initGrNormH = 457.57 + initNorm = 14.111 + SC1 = 1.2652e-03 + SC2 = 9.6405e-04 + SC3 = 0.047906 + W_density = 0.7700 + H_density = 0.7050 iter = 12 - elapsed = 0.041034 - tTotal = 0.5026 - subIterW = 33 + elapsed = 0.019481 + tTotal = 0.2746 + subIterW = 32 subIterH = 24 - relError = 0.2692 + relError = 0.2738 toshow = scalar structure containing the fields: - initGrNormW = 367.52 - initGrNormH = 439.89 - initNorm = 14.465 - SC1 = 1.0150e-03 - SC2 = 8.0590e-04 - SC3 = 0.036136 - W_density = 0.7933 - H_density = 0.7550 + initGrNormW = 398.15 + initGrNormH = 457.57 + initNorm = 14.111 + SC1 = 1.0676e-03 + SC2 = 8.1351e-04 + SC3 = 0.040089 + W_density = 0.7700 + H_density = 0.7050 iter = 13 - elapsed = 0.040745 - tTotal = 0.5434 - subIterW = 33 - subIterH = 24 - relError = 0.2683 + elapsed = 0.020068 + tTotal = 0.2946 + subIterW = 28 + subIterH = 26 + relError = 0.2728 toshow = scalar structure containing the fields: - initGrNormW = 367.52 - initGrNormH = 439.89 - initNorm = 14.465 - SC1 = 8.5751e-04 - SC2 = 6.7400e-04 - SC3 = 0.031446 - W_density = 0.7833 - H_density = 0.7500 + initGrNormW = 398.15 + initGrNormH = 457.57 + initNorm = 14.111 + SC1 = 9.1690e-04 + SC2 = 6.9868e-04 + SC3 = 0.032202 + W_density = 0.7667 + H_density = 0.7100 iter = 14 - elapsed = 0.041614 - tTotal = 0.5850 - subIterW = 33 - subIterH = 26 - relError = 0.2676 + elapsed = 0.020773 + tTotal = 0.3154 + subIterW = 30 + subIterH = 27 + relError = 0.2718 toshow = scalar structure containing the fields: - initGrNormW = 367.52 - initGrNormH = 439.89 - initNorm = 14.465 - SC1 = 7.2869e-04 - SC2 = 5.7566e-04 - SC3 = 0.025072 - W_density = 0.7867 - H_density = 0.7650 + initGrNormW = 398.15 + initGrNormH = 457.57 + initNorm = 14.111 + SC1 = 8.0139e-04 + SC2 = 5.9944e-04 + SC3 = 0.028602 + W_density = 0.7633 + H_density = 0.6950 iter = 15 - elapsed = 0.038608 - tTotal = 0.6236 - subIterW = 34 - subIterH = 20 - relError = 0.2670 + elapsed = 0.020705 + tTotal = 0.3361 + subIterW = 30 + subIterH = 28 + relError = 0.2709 toshow = scalar structure containing the fields: - initGrNormW = 367.52 - initGrNormH = 439.89 - initNorm = 14.465 - SC1 = 6.3149e-04 - SC2 = 4.9130e-04 - SC3 = 0.021671 - W_density = 0.7800 - H_density = 0.7550 + initGrNormW = 398.15 + initGrNormH = 457.57 + initNorm = 14.111 + SC1 = 6.9174e-04 + SC2 = 5.1465e-04 + SC3 = 0.025679 + W_density = 0.7533 + H_density = 0.7050 iter = 16 - elapsed = 0.035742 - tTotal = 0.6593 - subIterW = 36 - subIterH = 24 - relError = 0.2666 + elapsed = 0.021206 + tTotal = 0.3573 + subIterW = 31 + subIterH = 26 + relError = 0.2703 toshow = scalar structure containing the fields: - initGrNormW = 367.52 - initGrNormH = 439.89 - initNorm = 14.465 - SC1 = 5.5003e-04 - SC2 = 4.2462e-04 - SC3 = 0.019449 - W_density = 0.7767 - H_density = 0.7500 + initGrNormW = 398.15 + initGrNormH = 457.57 + initNorm = 14.111 + SC1 = 6.0069e-04 + SC2 = 4.4932e-04 + SC3 = 0.021815 + W_density = 0.7567 + H_density = 0.7050 iter = 17 - elapsed = 0.041730 - tTotal = 0.7011 - subIterW = 35 - subIterH = 24 - relError = 0.2663 + elapsed = 0.021090 + tTotal = 0.3784 + subIterW = 31 + subIterH = 25 + relError = 0.2698 toshow = scalar structure containing the fields: - initGrNormW = 367.52 - initGrNormH = 439.89 - initNorm = 14.465 - SC1 = 4.8106e-04 - SC2 = 3.7234e-04 - SC3 = 0.015967 - W_density = 0.7767 - H_density = 0.7500 + initGrNormW = 398.15 + initGrNormH = 457.57 + initNorm = 14.111 + SC1 = 5.3681e-04 + SC2 = 3.9939e-04 + SC3 = 0.019778 + W_density = 0.7467 + H_density = 0.7200 iter = 18 - elapsed = 0.036922 - tTotal = 0.7380 + elapsed = 0.022155 + tTotal = 0.4006 subIterW = 34 - subIterH = 24 - relError = 0.2660 + subIterH = 25 + relError = 0.2693 toshow = scalar structure containing the fields: - initGrNormW = 367.52 - initGrNormH = 439.89 - initNorm = 14.465 - SC1 = 4.2474e-04 - SC2 = 3.2875e-04 - SC3 = 0.014456 - W_density = 0.7800 - H_density = 0.7600 + initGrNormW = 398.15 + initGrNormH = 457.57 + initNorm = 14.111 + SC1 = 4.8655e-04 + SC2 = 3.6296e-04 + SC3 = 0.017653 + W_density = 0.7467 + H_density = 0.7250 iter = 19 - elapsed = 0.036542 - tTotal = 0.7745 - subIterW = 34 - subIterH = 24 - relError = 0.2658 + elapsed = 0.023311 + tTotal = 0.4239 + subIterW = 35 + subIterH = 25 + relError = 0.2689 toshow = scalar structure containing the fields: - initGrNormW = 367.52 - initGrNormH = 439.89 - initNorm = 14.465 - SC1 = 3.7915e-04 - SC2 = 2.9270e-04 - SC3 = 0.013475 - W_density = 0.7767 - H_density = 0.7550 + initGrNormW = 398.15 + initGrNormH = 457.57 + initNorm = 14.111 + SC1 = 4.5089e-04 + SC2 = 3.3456e-04 + SC3 = 0.015516 + W_density = 0.7367 + H_density = 0.7250 iter = 20 - elapsed = 0.037042 - tTotal = 0.8116 - subIterW = 34 - subIterH = 24 - relError = 0.2656 + elapsed = 0.028959 + tTotal = 0.4528 + subIterW = 37 + subIterH = 27 + relError = 0.2685 toshow = scalar structure containing the fields: - initGrNormW = 367.52 - initGrNormH = 439.89 - initNorm = 14.465 - SC1 = 3.3708e-04 - SC2 = 2.6023e-04 - SC3 = 0.011808 - W_density = 0.7800 - H_density = 0.7500 + initGrNormW = 398.15 + initGrNormH = 457.57 + initNorm = 14.111 + SC1 = 4.2156e-04 + SC2 = 3.1111e-04 + SC3 = 0.014843 + W_density = 0.7300 + H_density = 0.7300 iter = 21 - elapsed = 0.036357 - tTotal = 0.8479 - subIterW = 34 - subIterH = 24 - relError = 0.2654 + elapsed = 0.026957 + tTotal = 0.4798 + subIterW = 40 + subIterH = 27 + relError = 0.2682 toshow = scalar structure containing the fields: - initGrNormW = 367.52 - initGrNormH = 439.89 - initNorm = 14.465 - SC1 = 3.0388e-04 - SC2 = 2.3338e-04 - SC3 = 0.010457 - W_density = 0.7733 - H_density = 0.7500 + initGrNormW = 398.15 + initGrNormH = 457.57 + initNorm = 14.111 + SC1 = 3.9747e-04 + SC2 = 2.9412e-04 + SC3 = 0.013377 + W_density = 0.7300 + H_density = 0.7350 iter = 22 - elapsed = 0.038536 - tTotal = 0.8865 + elapsed = 0.029599 + tTotal = 0.5094 subIterW = 38 - subIterH = 25 - relError = 0.2653 + subIterH = 30 + relError = 0.2679 toshow = scalar structure containing the fields: - initGrNormW = 367.52 - initGrNormH = 439.89 - initNorm = 14.465 - SC1 = 2.7571e-04 - SC2 = 2.1120e-04 - SC3 = 8.8845e-03 - W_density = 0.7700 - H_density = 0.7550 + initGrNormW = 398.15 + initGrNormH = 457.57 + initNorm = 14.111 + SC1 = 3.8308e-04 + SC2 = 2.8271e-04 + SC3 = 0.012441 + W_density = 0.7300 + H_density = 0.7400 iter = 23 - elapsed = 0.038706 - tTotal = 0.9252 - subIterW = 37 - subIterH = 23 - relError = 0.2652 + elapsed = 0.025868 + tTotal = 0.5353 + subIterW = 40 + subIterH = 29 + relError = 0.2677 Stop: tolerance reached. final = @@ -3809,19 +3833,19 @@ scalar structure containing the fields: iterations = 23 - elapsed_total = 0.9252 - relative_error = 0.2652 - W_density = 0.7700 - H_density = 0.7550 + elapsed_total = 0.5353 + relative_error = 0.2677 + W_density = 0.7300 + H_density = 0.7400 ***** test [W,H,iter,HIS]=nmf_bpas(A,k,'verbose','type','sparse'); - initGrNormW = 394.55 - initGrNormH = 471.39 - initNorm = 14.465 - SC1 = 1.2294 - SC2 = 614.72 - SC3 = 0.4971 + initGrNormW = 332.35 + initGrNormH = 432.98 + initNorm = 14.111 + SC1 = 1.0917 + SC2 = 545.83 + SC3 = 0.4817 W_density = 1 H_density = 1 par = @@ -3832,8 +3856,8 @@ n = 20 type = sparse nnls_solver = bp - alpha = 0.5165 - beta = 0.5165 + alpha = 0.4933 + beta = 0.4933 max_iter = 100 min_iter = 20 max_time = 1000 @@ -3845,438 +3869,438 @@ scalar structure containing the fields: - initGrNormW = 394.55 - initGrNormH = 471.39 - initNorm = 14.465 - SC1 = 0.023119 - SC2 = 0.021085 - SC3 = 1.6541 - W_density = 0.9800 - H_density = 0.6050 + initGrNormW = 332.35 + initGrNormH = 432.98 + initNorm = 14.111 + SC1 = 0.026453 + SC2 = 0.023226 + SC3 = 1.6844 + W_density = 0.9267 + H_density = 0.5550 iter = 1 - elapsed = 0.053207 - tTotal = 0.053207 - subIterW = 3 - subIterH = 30 - relError = 0.4340 + elapsed = 0.029989 + tTotal = 0.029989 + subIterW = 12 + subIterH = 37 + relError = 0.4403 toshow = scalar structure containing the fields: - initGrNormW = 394.55 - initGrNormH = 471.39 - initNorm = 14.465 - SC1 = 0.011850 - SC2 = 8.9584e-03 - SC3 = 0.7951 - W_density = 0.8800 - H_density = 0.3950 + initGrNormW = 332.35 + initGrNormH = 432.98 + initNorm = 14.111 + SC1 = 0.010735 + SC2 = 8.6526e-03 + SC3 = 0.7030 + W_density = 0.9067 + H_density = 0.3500 iter = 2 - elapsed = 0.069017 - tTotal = 0.1222 - subIterW = 18 - subIterH = 21 - relError = 0.3609 + elapsed = 0.029436 + tTotal = 0.059425 + subIterW = 10 + subIterH = 30 + relError = 0.3709 toshow = scalar structure containing the fields: - initGrNormW = 394.55 - initGrNormH = 471.39 - initNorm = 14.465 - SC1 = 6.8422e-03 - SC2 = 4.7074e-03 - SC3 = 0.4820 - W_density = 0.8533 - H_density = 0.3200 + initGrNormW = 332.35 + initGrNormH = 432.98 + initNorm = 14.111 + SC1 = 6.6446e-03 + SC2 = 4.6911e-03 + SC3 = 0.4271 + W_density = 0.8400 + H_density = 0.3000 iter = 3 - elapsed = 0.055602 - tTotal = 0.1778 - subIterW = 7 - subIterH = 19 - relError = 0.3390 + elapsed = 0.024170 + tTotal = 0.083595 + subIterW = 3 + subIterH = 18 + relError = 0.3543 toshow = scalar structure containing the fields: - initGrNormW = 394.55 - initGrNormH = 471.39 - initNorm = 14.465 - SC1 = 5.6500e-03 - SC2 = 3.3674e-03 - SC3 = 0.4131 - W_density = 0.7767 - H_density = 0.2800 + initGrNormW = 332.35 + initGrNormH = 432.98 + initNorm = 14.111 + SC1 = 4.5470e-03 + SC2 = 2.9737e-03 + SC3 = 0.3186 + W_density = 0.8400 + H_density = 0.2900 iter = 4 - elapsed = 0.044640 - tTotal = 0.2225 - subIterW = 9 - subIterH = 11 - relError = 0.3344 + elapsed = 0.017596 + tTotal = 0.1012 + subIterW = 3 + subIterH = 5 + relError = 0.3444 toshow = scalar structure containing the fields: - initGrNormW = 394.55 - initGrNormH = 471.39 - initNorm = 14.465 - SC1 = 4.7367e-03 - SC2 = 2.7378e-03 - SC3 = 0.3176 - W_density = 0.7700 - H_density = 0.2650 + initGrNormW = 332.35 + initGrNormH = 432.98 + initNorm = 14.111 + SC1 = 3.1411e-03 + SC2 = 1.9726e-03 + SC3 = 0.2219 + W_density = 0.8300 + H_density = 0.2800 iter = 5 - elapsed = 0.032528 - tTotal = 0.2550 + elapsed = 0.023033 + tTotal = 0.1242 subIterW = 5 - subIterH = 3 - relError = 0.3337 + subIterH = 8 + relError = 0.3382 toshow = scalar structure containing the fields: - initGrNormW = 394.55 - initGrNormH = 471.39 - initNorm = 14.465 - SC1 = 4.0273e-03 - SC2 = 2.2553e-03 - SC3 = 0.2773 - W_density = 0.7600 - H_density = 0.2500 + initGrNormW = 332.35 + initGrNormH = 432.98 + initNorm = 14.111 + SC1 = 2.3860e-03 + SC2 = 1.4793e-03 + SC3 = 0.1604 + W_density = 0.8200 + H_density = 0.2900 iter = 6 - elapsed = 0.032719 - tTotal = 0.2877 - subIterW = 6 + elapsed = 0.019018 + tTotal = 0.1432 + subIterW = 5 subIterH = 2 - relError = 0.3357 + relError = 0.3350 toshow = scalar structure containing the fields: - initGrNormW = 394.55 - initGrNormH = 471.39 - initNorm = 14.465 - SC1 = 3.4771e-03 - SC2 = 1.8429e-03 - SC3 = 0.2556 - W_density = 0.7367 - H_density = 0.2200 + initGrNormW = 332.35 + initGrNormH = 432.98 + initNorm = 14.111 + SC1 = 2.0208e-03 + SC2 = 1.2327e-03 + SC3 = 0.1315 + W_density = 0.8167 + H_density = 0.2750 iter = 7 - elapsed = 0.034120 - tTotal = 0.3218 - subIterW = 6 - subIterH = 6 - relError = 0.3387 + elapsed = 0.020401 + tTotal = 0.1636 + subIterW = 4 + subIterH = 0 + relError = 0.3338 toshow = scalar structure containing the fields: - initGrNormW = 394.55 - initGrNormH = 471.39 - initNorm = 14.465 - SC1 = 2.8668e-03 - SC2 = 1.5022e-03 - SC3 = 0.2052 - W_density = 0.7300 - H_density = 0.2150 + initGrNormW = 332.35 + initGrNormH = 432.98 + initNorm = 14.111 + SC1 = 1.9094e-03 + SC2 = 1.1418e-03 + SC3 = 0.1235 + W_density = 0.8067 + H_density = 0.2600 iter = 8 - elapsed = 0.029989 - tTotal = 0.3518 - subIterW = 2 - subIterH = 0 - relError = 0.3423 + elapsed = 0.019007 + tTotal = 0.1826 + subIterW = 3 + subIterH = 2 + relError = 0.3340 toshow = scalar structure containing the fields: - initGrNormW = 394.55 - initGrNormH = 471.39 - initNorm = 14.465 - SC1 = 2.3562e-03 - SC2 = 1.2064e-03 - SC3 = 0.1667 - W_density = 0.7200 - H_density = 0.2000 + initGrNormW = 332.35 + initGrNormH = 432.98 + initNorm = 14.111 + SC1 = 1.7954e-03 + SC2 = 1.0629e-03 + SC3 = 0.1073 + W_density = 0.8000 + H_density = 0.2600 iter = 9 - elapsed = 0.031666 - tTotal = 0.3835 - subIterW = 3 - subIterH = 2 - relError = 0.3459 + elapsed = 0.018509 + tTotal = 0.2012 + subIterW = 4 + subIterH = 4 + relError = 0.3346 toshow = scalar structure containing the fields: - initGrNormW = 394.55 - initGrNormH = 471.39 - initNorm = 14.465 - SC1 = 1.8616e-03 - SC2 = 9.3827e-04 - SC3 = 0.1217 - W_density = 0.7100 - H_density = 0.1950 + initGrNormW = 332.35 + initGrNormH = 432.98 + initNorm = 14.111 + SC1 = 1.5558e-03 + SC2 = 9.0549e-04 + SC3 = 0.098296 + W_density = 0.7867 + H_density = 0.2550 iter = 10 - elapsed = 0.029934 - tTotal = 0.4134 - subIterW = 3 - subIterH = 0 - relError = 0.3486 + elapsed = 0.018871 + tTotal = 0.2200 + subIterW = 4 + subIterH = 3 + relError = 0.3351 toshow = scalar structure containing the fields: - initGrNormW = 394.55 - initGrNormH = 471.39 - initNorm = 14.465 - SC1 = 1.3978e-03 - SC2 = 6.9609e-04 - SC3 = 0.089997 - W_density = 0.7067 - H_density = 0.1850 + initGrNormW = 332.35 + initGrNormH = 432.98 + initNorm = 14.111 + SC1 = 1.4099e-03 + SC2 = 8.3467e-04 + SC3 = 0.096156 + W_density = 0.8033 + H_density = 0.2600 iter = 11 - elapsed = 0.030061 - tTotal = 0.4435 - subIterW = 0 - subIterH = 2 - relError = 0.3507 + elapsed = 0.017382 + tTotal = 0.2374 + subIterW = 4 + subIterH = 0 + relError = 0.3354 toshow = scalar structure containing the fields: - initGrNormW = 394.55 - initGrNormH = 471.39 - initNorm = 14.465 - SC1 = 1.0446e-03 - SC2 = 5.1186e-04 - SC3 = 0.069419 - W_density = 0.6967 - H_density = 0.1800 + initGrNormW = 332.35 + initGrNormH = 432.98 + initNorm = 14.111 + SC1 = 1.5112e-03 + SC2 = 8.8861e-04 + SC3 = 0.029483 + W_density = 0.8033 + H_density = 0.2500 iter = 12 - elapsed = 0.030265 - tTotal = 0.4737 - subIterW = 3 - subIterH = 0 - relError = 0.3521 + elapsed = 0.012673 + tTotal = 0.2501 + subIterW = 0 + subIterH = 2 + relError = 0.3356 toshow = scalar structure containing the fields: - initGrNormW = 394.55 - initGrNormH = 471.39 - initNorm = 14.465 - SC1 = 7.7057e-04 - SC2 = 3.7141e-04 - SC3 = 0.053109 - W_density = 0.6867 - H_density = 0.1750 + initGrNormW = 332.35 + initGrNormH = 432.98 + initNorm = 14.111 + SC1 = 1.5966e-03 + SC2 = 9.2924e-04 + SC3 = 0.1107 + W_density = 0.7967 + H_density = 0.2500 iter = 13 - elapsed = 0.030559 - tTotal = 0.5043 - subIterW = 3 - subIterH = 0 - relError = 0.3531 + elapsed = 0.017444 + tTotal = 0.2675 + subIterW = 2 + subIterH = 2 + relError = 0.3357 toshow = scalar structure containing the fields: - initGrNormW = 394.55 - initGrNormH = 471.39 - initNorm = 14.465 - SC1 = 5.4108e-04 - SC2 = 2.5972e-04 - SC3 = 0.040830 - W_density = 0.6867 - H_density = 0.1700 + initGrNormW = 332.35 + initGrNormH = 432.98 + initNorm = 14.111 + SC1 = 1.7083e-03 + SC2 = 9.9767e-04 + SC3 = 0.1037 + W_density = 0.7967 + H_density = 0.2550 iter = 14 - elapsed = 0.029508 - tTotal = 0.5338 + elapsed = 0.017151 + tTotal = 0.2847 subIterW = 2 subIterH = 0 - relError = 0.3537 + relError = 0.3362 toshow = scalar structure containing the fields: - initGrNormW = 394.55 - initGrNormH = 471.39 - initNorm = 14.465 - SC1 = 3.8460e-04 - SC2 = 1.8384e-04 - SC3 = 0.030829 - W_density = 0.6800 - H_density = 0.1700 + initGrNormW = 332.35 + initGrNormH = 432.98 + initNorm = 14.111 + SC1 = 1.4736e-03 + SC2 = 8.3994e-04 + SC3 = 0.099329 + W_density = 0.7833 + H_density = 0.2400 iter = 15 - elapsed = 0.023303 - tTotal = 0.5571 - subIterW = 2 - subIterH = 0 - relError = 0.3540 + elapsed = 0.019262 + tTotal = 0.3039 + subIterW = 4 + subIterH = 3 + relError = 0.3368 toshow = scalar structure containing the fields: - initGrNormW = 394.55 - initGrNormH = 471.39 - initNorm = 14.465 - SC1 = 3.0043e-04 - SC2 = 1.4361e-04 - SC3 = 0.025551 - W_density = 0.6800 - H_density = 0.1750 + initGrNormW = 332.35 + initGrNormH = 432.98 + initNorm = 14.111 + SC1 = 1.3630e-03 + SC2 = 7.7421e-04 + SC3 = 0.093087 + W_density = 0.7800 + H_density = 0.2450 iter = 16 - elapsed = 0.029655 - tTotal = 0.5868 - subIterW = 2 + elapsed = 0.016729 + tTotal = 0.3207 + subIterW = 0 subIterH = 0 - relError = 0.3540 + relError = 0.3378 toshow = scalar structure containing the fields: - initGrNormW = 394.55 - initGrNormH = 471.39 - initNorm = 14.465 - SC1 = 2.7752e-04 - SC2 = 1.3265e-04 - SC3 = 5.3873e-03 - W_density = 0.6800 - H_density = 0.1750 + initGrNormW = 332.35 + initGrNormH = 432.98 + initNorm = 14.111 + SC1 = 1.2378e-03 + SC2 = 6.9812e-04 + SC3 = 0.083255 + W_density = 0.7700 + H_density = 0.2500 iter = 17 - elapsed = 0.015808 - tTotal = 0.6026 - subIterW = 0 + elapsed = 0.017499 + tTotal = 0.3382 + subIterW = 3 subIterH = 0 - relError = 0.3540 + relError = 0.3392 toshow = scalar structure containing the fields: - initGrNormW = 394.55 - initGrNormH = 471.39 - initNorm = 14.465 - SC1 = 3.1270e-04 - SC2 = 1.4947e-04 - SC3 = 4.9376e-03 - W_density = 0.6800 - H_density = 0.1750 + initGrNormW = 332.35 + initGrNormH = 432.98 + initNorm = 14.111 + SC1 = 1.0919e-03 + SC2 = 6.0709e-04 + SC3 = 0.078078 + W_density = 0.7600 + H_density = 0.2500 iter = 18 - elapsed = 0.015838 - tTotal = 0.6184 - subIterW = 0 - subIterH = 0 - relError = 0.3538 + elapsed = 0.018891 + tTotal = 0.3571 + subIterW = 3 + subIterH = 2 + relError = 0.3407 toshow = scalar structure containing the fields: - initGrNormW = 394.55 - initGrNormH = 471.39 - initNorm = 14.465 - SC1 = 3.9553e-04 - SC2 = 1.8986e-04 - SC3 = 0.025426 - W_density = 0.6833 - H_density = 0.1750 + initGrNormW = 332.35 + initGrNormH = 432.98 + initNorm = 14.111 + SC1 = 1.0444e-03 + SC2 = 5.7648e-04 + SC3 = 0.074342 + W_density = 0.7533 + H_density = 0.2500 iter = 19 - elapsed = 0.022377 - tTotal = 0.6408 - subIterW = 0 + elapsed = 0.014718 + tTotal = 0.3718 + subIterW = 2 subIterH = 0 - relError = 0.3535 + relError = 0.3422 toshow = scalar structure containing the fields: - initGrNormW = 394.55 - initGrNormH = 471.39 - initNorm = 14.465 - SC1 = 5.2290e-04 - SC2 = 2.5099e-04 - SC3 = 6.9513e-03 - W_density = 0.6833 - H_density = 0.1750 + initGrNormW = 332.35 + initGrNormH = 432.98 + initNorm = 14.111 + SC1 = 9.5992e-04 + SC2 = 5.2603e-04 + SC3 = 0.068525 + W_density = 0.7500 + H_density = 0.2450 iter = 20 - elapsed = 0.015755 - tTotal = 0.6566 + elapsed = 0.017524 + tTotal = 0.3893 subIterW = 0 subIterH = 0 - relError = 0.3531 + relError = 0.3436 toshow = scalar structure containing the fields: - initGrNormW = 394.55 - initGrNormH = 471.39 - initNorm = 14.465 - SC1 = 6.6949e-04 - SC2 = 3.2269e-04 - SC3 = 8.3658e-03 - W_density = 0.6833 - H_density = 0.1750 + initGrNormW = 332.35 + initGrNormH = 432.98 + initNorm = 14.111 + SC1 = 9.1604e-04 + SC2 = 4.9466e-04 + SC3 = 0.062619 + W_density = 0.7367 + H_density = 0.2450 iter = 21 - elapsed = 0.015973 - tTotal = 0.6725 - subIterW = 0 + elapsed = 0.014614 + tTotal = 0.4039 + subIterW = 4 subIterH = 0 - relError = 0.3525 + relError = 0.3449 toshow = scalar structure containing the fields: - initGrNormW = 394.55 - initGrNormH = 471.39 - initNorm = 14.465 - SC1 = 7.2761e-04 - SC2 = 3.5071e-04 - SC3 = 0.040092 - W_density = 0.6833 - H_density = 0.1800 + initGrNormW = 332.35 + initGrNormH = 432.98 + initNorm = 14.111 + SC1 = 8.5766e-04 + SC2 = 4.5628e-04 + SC3 = 0.058664 + W_density = 0.7300 + H_density = 0.2350 iter = 22 - elapsed = 0.030908 - tTotal = 0.7034 + elapsed = 0.019746 + tTotal = 0.4237 subIterW = 2 - subIterH = 0 - relError = 0.3519 + subIterH = 1 + relError = 0.3461 toshow = scalar structure containing the fields: - initGrNormW = 394.55 - initGrNormH = 471.39 - initNorm = 14.465 - SC1 = 6.1921e-04 - SC2 = 2.9474e-04 - SC3 = 0.032907 - W_density = 0.6733 - H_density = 0.1800 + initGrNormW = 332.35 + initGrNormH = 432.98 + initNorm = 14.111 + SC1 = 7.7651e-04 + SC2 = 4.0844e-04 + SC3 = 0.053063 + W_density = 0.7267 + H_density = 0.2250 iter = 23 - elapsed = 0.025475 - tTotal = 0.7289 - subIterW = 4 - subIterH = 0 - relError = 0.3515 + elapsed = 0.020755 + tTotal = 0.4444 + subIterW = 3 + subIterH = 2 + relError = 0.3469 Stop: tolerance reached. final = @@ -4284,19 +4308,19 @@ scalar structure containing the fields: iterations = 23 - elapsed_total = 0.7289 - relative_error = 0.3515 - W_density = 0.6733 - H_density = 0.1800 + elapsed_total = 0.4444 + relative_error = 0.3469 + W_density = 0.7267 + H_density = 0.2250 ***** test [W,H,iter,HIS]=nmf_bpas(A,k,'verbose','type','sparse','nnlssolver','bp','alpha',1.1,'beta',1.3); - initGrNormW = 409.83 - initGrNormH = 577.59 - initNorm = 14.465 - SC1 = 1.4164 - SC2 = 708.22 - SC3 = 0.5108 + initGrNormW = 322.41 + initGrNormH = 499.99 + initNorm = 14.111 + SC1 = 1.1899 + SC2 = 594.93 + SC3 = 0.4839 W_density = 1 H_density = 1 par = @@ -4320,438 +4344,438 @@ scalar structure containing the fields: - initGrNormW = 409.83 - initGrNormH = 577.59 - initNorm = 14.465 - SC1 = 0.025483 - SC2 = 0.023240 - SC3 = 2.0219 - W_density = 0.9900 - H_density = 0.5650 + initGrNormW = 322.41 + initGrNormH = 499.99 + initNorm = 14.111 + SC1 = 0.030288 + SC2 = 0.026654 + SC3 = 2.1469 + W_density = 0.9833 + H_density = 0.5200 iter = 1 - elapsed = 0.053049 - tTotal = 0.053049 - subIterW = 1 - subIterH = 37 - relError = 0.5442 + elapsed = 0.033331 + tTotal = 0.033331 + subIterW = 3 + subIterH = 36 + relError = 0.5676 toshow = scalar structure containing the fields: - initGrNormW = 409.83 - initGrNormH = 577.59 - initNorm = 14.465 - SC1 = 0.015280 - SC2 = 6.3869e-03 - SC3 = 0.6841 - W_density = 0.5133 - H_density = 0.2550 + initGrNormW = 322.41 + initGrNormH = 499.99 + initNorm = 14.111 + SC1 = 0.020394 + SC2 = 9.0143e-03 + SC3 = 1.1629 + W_density = 0.5800 + H_density = 0.1950 iter = 2 - elapsed = 0.058895 - tTotal = 0.1119 - subIterW = 6 - subIterH = 33 - relError = 0.4528 + elapsed = 0.037866 + tTotal = 0.071197 + subIterW = 12 + subIterH = 29 + relError = 0.4796 toshow = scalar structure containing the fields: - initGrNormW = 409.83 - initGrNormH = 577.59 - initNorm = 14.465 - SC1 = 0.014800 - SC2 = 4.4105e-03 - SC3 = 0.6277 - W_density = 0.3567 - H_density = 0.1850 + initGrNormW = 322.41 + initGrNormH = 499.99 + initNorm = 14.111 + SC1 = 0.017322 + SC2 = 5.4046e-03 + SC3 = 0.7131 + W_density = 0.4000 + H_density = 0.1800 iter = 3 - elapsed = 0.039735 - tTotal = 0.1517 - subIterW = 6 - subIterH = 8 - relError = 0.4465 + elapsed = 0.025758 + tTotal = 0.096955 + subIterW = 13 + subIterH = 7 + relError = 0.4810 toshow = scalar structure containing the fields: - initGrNormW = 409.83 - initGrNormH = 577.59 - initNorm = 14.465 - SC1 = 0.012511 - SC2 = 2.9526e-03 - SC3 = 0.4914 - W_density = 0.2733 - H_density = 0.1600 + initGrNormW = 322.41 + initGrNormH = 499.99 + initNorm = 14.111 + SC1 = 0.016239 + SC2 = 3.1503e-03 + SC3 = 0.4235 + W_density = 0.2333 + H_density = 0.1350 iter = 4 - elapsed = 0.025187 - tTotal = 0.1769 - subIterW = 2 - subIterH = 3 - relError = 0.4426 + elapsed = 0.020478 + tTotal = 0.1174 + subIterW = 5 + subIterH = 2 + relError = 0.4912 toshow = scalar structure containing the fields: - initGrNormW = 409.83 - initGrNormH = 577.59 - initNorm = 14.465 - SC1 = 0.010442 - SC2 = 2.2763e-03 - SC3 = 0.4534 - W_density = 0.2667 - H_density = 0.1450 + initGrNormW = 322.41 + initGrNormH = 499.99 + initNorm = 14.111 + SC1 = 0.012418 + SC2 = 2.1360e-03 + SC3 = 0.2696 + W_density = 0.2067 + H_density = 0.1200 iter = 5 - elapsed = 0.023458 - tTotal = 0.2003 - subIterW = 1 - subIterH = 2 - relError = 0.4441 + elapsed = 0.014128 + tTotal = 0.1316 + subIterW = 2 + subIterH = 1 + relError = 0.4997 toshow = scalar structure containing the fields: - initGrNormW = 409.83 - initGrNormH = 577.59 - initNorm = 14.465 - SC1 = 0.011120 - SC2 = 2.4241e-03 - SC3 = 0.4937 - W_density = 0.2700 - H_density = 0.1400 + initGrNormW = 322.41 + initGrNormH = 499.99 + initNorm = 14.111 + SC1 = 9.7131e-03 + SC2 = 1.3016e-03 + SC3 = 0.1908 + W_density = 0.1533 + H_density = 0.1050 iter = 6 - elapsed = 0.020030 - tTotal = 0.2204 + elapsed = 0.013994 + tTotal = 0.1456 subIterW = 0 - subIterH = 0 - relError = 0.4490 + subIterH = 1 + relError = 0.5044 toshow = scalar structure containing the fields: - initGrNormW = 409.83 - initGrNormH = 577.59 - initNorm = 14.465 - SC1 = 0.012153 - SC2 = 2.5278e-03 - SC3 = 0.4979 - W_density = 0.2533 - H_density = 0.1400 + initGrNormW = 322.41 + initGrNormH = 499.99 + initNorm = 14.111 + SC1 = 4.5541e-03 + SC2 = 4.5541e-04 + SC3 = 0.088537 + W_density = 0.1000 + H_density = 0.1000 iter = 7 - elapsed = 0.015054 - tTotal = 0.2354 - subIterW = 2 + elapsed = 8.6580e-03 + tTotal = 0.1542 + subIterW = 0 subIterH = 0 - relError = 0.4548 + relError = 0.5050 toshow = scalar structure containing the fields: - initGrNormW = 409.83 - initGrNormH = 577.59 - initNorm = 14.465 - SC1 = 0.012686 - SC2 = 2.4610e-03 - SC3 = 0.3718 - W_density = 0.2333 - H_density = 0.1350 + initGrNormW = 322.41 + initGrNormH = 499.99 + initNorm = 14.111 + SC1 = 2.9093e-03 + SC2 = 2.9093e-04 + SC3 = 0.054900 + W_density = 0.1000 + H_density = 0.1000 iter = 8 - elapsed = 0.021601 - tTotal = 0.2570 - subIterW = 3 + elapsed = 4.3320e-03 + tTotal = 0.1585 + subIterW = 0 subIterH = 0 - relError = 0.4620 + relError = 0.5048 toshow = scalar structure containing the fields: - initGrNormW = 409.83 - initGrNormH = 577.59 - initNorm = 14.465 - SC1 = 0.010374 - SC2 = 1.5768e-03 - SC3 = 0.2694 - W_density = 0.1767 - H_density = 0.1150 + initGrNormW = 322.41 + initGrNormH = 499.99 + initNorm = 14.111 + SC1 = 1.8771e-03 + SC2 = 1.8771e-04 + SC3 = 0.042796 + W_density = 0.1000 + H_density = 0.1000 iter = 9 - elapsed = 0.023076 - tTotal = 0.2801 - subIterW = 2 - subIterH = 1 - relError = 0.4673 + elapsed = 4.3250e-03 + tTotal = 0.1629 + subIterW = 0 + subIterH = 0 + relError = 0.5047 toshow = scalar structure containing the fields: - initGrNormW = 409.83 - initGrNormH = 577.59 - initNorm = 14.465 - SC1 = 9.2356e-03 - SC2 = 1.2930e-03 - SC3 = 0.1317 - W_density = 0.1600 - H_density = 0.1100 + initGrNormW = 322.41 + initGrNormH = 499.99 + initNorm = 14.111 + SC1 = 1.2126e-03 + SC2 = 1.2126e-04 + SC3 = 0.026540 + W_density = 0.1000 + H_density = 0.1000 iter = 10 - elapsed = 0.021513 - tTotal = 0.3016 - subIterW = 1 + elapsed = 4.0260e-03 + tTotal = 0.1669 + subIterW = 0 subIterH = 0 - relError = 0.4728 + relError = 0.5046 toshow = scalar structure containing the fields: - initGrNormW = 409.83 - initGrNormH = 577.59 - initNorm = 14.465 - SC1 = 1.5264e-03 - SC2 = 1.5264e-04 - SC3 = 0.026527 + initGrNormW = 322.41 + initGrNormH = 499.99 + initNorm = 14.111 + SC1 = 7.8397e-04 + SC2 = 7.8397e-05 + SC3 = 0.017874 W_density = 0.1000 H_density = 0.1000 iter = 11 - elapsed = 0.015163 - tTotal = 0.3168 + elapsed = 3.8380e-03 + tTotal = 0.1707 subIterW = 0 - subIterH = 1 - relError = 0.4743 + subIterH = 0 + relError = 0.5046 toshow = scalar structure containing the fields: - initGrNormW = 409.83 - initGrNormH = 577.59 - initNorm = 14.465 - SC1 = 8.6652e-04 - SC2 = 8.6652e-05 - SC3 = 0.022266 + initGrNormW = 322.41 + initGrNormH = 499.99 + initNorm = 14.111 + SC1 = 5.0710e-04 + SC2 = 5.0710e-05 + SC3 = 9.2490e-03 W_density = 0.1000 H_density = 0.1000 iter = 12 - elapsed = 6.9020e-03 - tTotal = 0.3237 + elapsed = 4.0700e-03 + tTotal = 0.1748 subIterW = 0 subIterH = 0 - relError = 0.4743 + relError = 0.5046 toshow = scalar structure containing the fields: - initGrNormW = 409.83 - initGrNormH = 577.59 - initNorm = 14.465 - SC1 = 5.7175e-04 - SC2 = 5.7175e-05 - SC3 = 0.013057 + initGrNormW = 322.41 + initGrNormH = 499.99 + initNorm = 14.111 + SC1 = 3.2812e-04 + SC2 = 3.2812e-05 + SC3 = 6.4120e-03 W_density = 0.1000 H_density = 0.1000 iter = 13 - elapsed = 6.6430e-03 - tTotal = 0.3303 + elapsed = 4.0710e-03 + tTotal = 0.1789 subIterW = 0 subIterH = 0 - relError = 0.4742 + relError = 0.5045 toshow = scalar structure containing the fields: - initGrNormW = 409.83 - initGrNormH = 577.59 - initNorm = 14.465 - SC1 = 3.7766e-04 - SC2 = 3.7766e-05 - SC3 = 8.6238e-03 + initGrNormW = 322.41 + initGrNormH = 499.99 + initNorm = 14.111 + SC1 = 2.1235e-04 + SC2 = 2.1235e-05 + SC3 = 4.3035e-03 W_density = 0.1000 H_density = 0.1000 iter = 14 - elapsed = 6.5530e-03 - tTotal = 0.3369 + elapsed = 4.2310e-03 + tTotal = 0.1831 subIterW = 0 subIterH = 0 - relError = 0.4742 + relError = 0.5045 toshow = scalar structure containing the fields: - initGrNormW = 409.83 - initGrNormH = 577.59 - initNorm = 14.465 - SC1 = 2.4952e-04 - SC2 = 2.4952e-05 - SC3 = 5.6977e-03 + initGrNormW = 322.41 + initGrNormH = 499.99 + initNorm = 14.111 + SC1 = 1.3745e-04 + SC2 = 1.3745e-05 + SC3 = 2.8926e-03 W_density = 0.1000 H_density = 0.1000 iter = 15 - elapsed = 6.6470e-03 - tTotal = 0.3435 + elapsed = 3.7930e-03 + tTotal = 0.1869 subIterW = 0 subIterH = 0 - relError = 0.4742 + relError = 0.5045 toshow = scalar structure containing the fields: - initGrNormW = 409.83 - initGrNormH = 577.59 - initNorm = 14.465 - SC1 = 1.6488e-04 - SC2 = 1.6488e-05 - SC3 = 4.2357e-03 + initGrNormW = 322.41 + initGrNormH = 499.99 + initNorm = 14.111 + SC1 = 8.8974e-05 + SC2 = 8.8974e-06 + SC3 = 1.8725e-03 W_density = 0.1000 H_density = 0.1000 iter = 16 - elapsed = 6.7100e-03 - tTotal = 0.3502 + elapsed = 3.8200e-03 + tTotal = 0.1907 subIterW = 0 subIterH = 0 - relError = 0.4742 + relError = 0.5045 toshow = scalar structure containing the fields: - initGrNormW = 409.83 - initGrNormH = 577.59 - initNorm = 14.465 - SC1 = 1.0897e-04 - SC2 = 1.0897e-05 - SC3 = 2.5840e-03 + initGrNormW = 322.41 + initGrNormH = 499.99 + initNorm = 14.111 + SC1 = 5.7599e-05 + SC2 = 5.7599e-06 + SC3 = 1.1256e-03 W_density = 0.1000 H_density = 0.1000 iter = 17 - elapsed = 6.7950e-03 - tTotal = 0.3570 + elapsed = 3.5040e-03 + tTotal = 0.1942 subIterW = 0 subIterH = 0 - relError = 0.4742 + relError = 0.5045 toshow = scalar structure containing the fields: - initGrNormW = 409.83 - initGrNormH = 577.59 - initNorm = 14.465 - SC1 = 7.2021e-05 - SC2 = 7.2021e-06 - SC3 = 1.6446e-03 + initGrNormW = 322.41 + initGrNormH = 499.99 + initNorm = 14.111 + SC1 = 3.7289e-05 + SC2 = 3.7289e-06 + SC3 = 6.3760e-04 W_density = 0.1000 H_density = 0.1000 iter = 18 - elapsed = 6.5290e-03 - tTotal = 0.3635 + elapsed = 3.8530e-03 + tTotal = 0.1981 subIterW = 0 subIterH = 0 - relError = 0.4742 + relError = 0.5045 toshow = scalar structure containing the fields: - initGrNormW = 409.83 - initGrNormH = 577.59 - initNorm = 14.465 - SC1 = 4.7603e-05 - SC2 = 4.7603e-06 - SC3 = 1.1288e-03 + initGrNormW = 322.41 + initGrNormH = 499.99 + initNorm = 14.111 + SC1 = 2.4141e-05 + SC2 = 2.4141e-06 + SC3 = 4.2610e-04 W_density = 0.1000 H_density = 0.1000 iter = 19 - elapsed = 6.7010e-03 - tTotal = 0.3702 + elapsed = 3.8250e-03 + tTotal = 0.2019 subIterW = 0 subIterH = 0 - relError = 0.4742 + relError = 0.5045 toshow = scalar structure containing the fields: - initGrNormW = 409.83 - initGrNormH = 577.59 - initNorm = 14.465 - SC1 = 3.1465e-05 - SC2 = 3.1465e-06 - SC3 = 8.0830e-04 + initGrNormW = 322.41 + initGrNormH = 499.99 + initNorm = 14.111 + SC1 = 1.5629e-05 + SC2 = 1.5629e-06 + SC3 = 3.0542e-04 W_density = 0.1000 H_density = 0.1000 iter = 20 - elapsed = 6.5350e-03 - tTotal = 0.3768 + elapsed = 3.7000e-03 + tTotal = 0.2056 subIterW = 0 subIterH = 0 - relError = 0.4742 + relError = 0.5045 toshow = scalar structure containing the fields: - initGrNormW = 409.83 - initGrNormH = 577.59 - initNorm = 14.465 - SC1 = 2.0798e-05 - SC2 = 2.0798e-06 - SC3 = 5.1291e-04 + initGrNormW = 322.41 + initGrNormH = 499.99 + initNorm = 14.111 + SC1 = 1.0118e-05 + SC2 = 1.0118e-06 + SC3 = 2.4072e-04 W_density = 0.1000 H_density = 0.1000 iter = 21 - elapsed = 6.6890e-03 - tTotal = 0.3835 + elapsed = 3.5340e-03 + tTotal = 0.2091 subIterW = 0 subIterH = 0 - relError = 0.4742 + relError = 0.5045 toshow = scalar structure containing the fields: - initGrNormW = 409.83 - initGrNormH = 577.59 - initNorm = 14.465 - SC1 = 1.3748e-05 - SC2 = 1.3748e-06 - SC3 = 2.9227e-04 + initGrNormW = 322.41 + initGrNormH = 499.99 + initNorm = 14.111 + SC1 = 6.5509e-06 + SC2 = 6.5509e-07 + SC3 = 1.2802e-04 W_density = 0.1000 H_density = 0.1000 iter = 22 - elapsed = 7.3830e-03 - tTotal = 0.3908 + elapsed = 4.2620e-03 + tTotal = 0.2134 subIterW = 0 subIterH = 0 - relError = 0.4742 + relError = 0.5045 toshow = scalar structure containing the fields: - initGrNormW = 409.83 - initGrNormH = 577.59 - initNorm = 14.465 - SC1 = 9.0874e-06 - SC2 = 9.0874e-07 - SC3 = 2.0751e-04 + initGrNormW = 322.41 + initGrNormH = 499.99 + initNorm = 14.111 + SC1 = 4.2412e-06 + SC2 = 4.2412e-07 + SC3 = 8.5952e-05 W_density = 0.1000 H_density = 0.1000 iter = 23 - elapsed = 7.3300e-03 - tTotal = 0.3982 + elapsed = 4.1710e-03 + tTotal = 0.2176 subIterW = 0 subIterH = 0 - relError = 0.4742 + relError = 0.5045 Stop: tolerance reached. final = @@ -4759,19 +4783,19 @@ scalar structure containing the fields: iterations = 23 - elapsed_total = 0.3982 - relative_error = 0.4742 + elapsed_total = 0.2176 + relative_error = 0.5045 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 = 370.94 - initGrNormH = 483.38 - initNorm = 14.465 - SC1 = 1.2186 - SC2 = 609.31 - SC3 = 0.5112 + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 1.1282 + SC2 = 564.12 + SC3 = 0.4920 W_density = 1 H_density = 1 par = @@ -4782,8 +4806,8 @@ n = 20 type = plain nnls_solver = bp - alpha = 0.5165 - beta = 0.5165 + alpha = 0.4933 + beta = 0.4933 max_iter = 100 min_iter = 20 max_time = 1000 @@ -4795,487 +4819,658 @@ scalar structure containing the fields: - initGrNormW = 370.94 - initGrNormH = 483.38 - initNorm = 14.465 - SC1 = 0.013056 - SC2 = 9.8183e-03 - SC3 = 0.5458 - W_density = 0.7833 - H_density = 0.5750 + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 0.020419 + SC2 = 0.015069 + SC3 = 0.6129 + W_density = 0.7533 + H_density = 0.5850 iter = 1 - elapsed = 0.068593 - tTotal = 0.068593 - subIterW = 33 + elapsed = 0.043446 + tTotal = 0.043446 + subIterW = 41 subIterH = 36 - relError = 0.3266 + relError = 0.3490 toshow = scalar structure containing the fields: - initGrNormW = 370.94 - initGrNormH = 483.38 - initNorm = 14.465 - SC1 = 8.5387e-03 - SC2 = 6.1820e-03 - SC3 = 0.3654 - W_density = 0.7667 - H_density = 0.6300 + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 0.010034 + SC2 = 6.9238e-03 + SC3 = 0.3801 + W_density = 0.7167 + H_density = 0.6100 iter = 2 - elapsed = 0.055562 - tTotal = 0.1242 - subIterW = 16 - subIterH = 22 - relError = 0.2846 + elapsed = 0.043975 + tTotal = 0.087421 + subIterW = 21 + subIterH = 21 + relError = 0.2999 toshow = scalar structure containing the fields: - initGrNormW = 370.94 - initGrNormH = 483.38 - initNorm = 14.465 - SC1 = 5.4790e-03 - SC2 = 3.9011e-03 - SC3 = 0.2371 - W_density = 0.7767 - H_density = 0.6050 + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 6.4904e-03 + SC2 = 4.4264e-03 + SC3 = 0.2528 + W_density = 0.7133 + H_density = 0.6350 iter = 3 - elapsed = 0.058024 - tTotal = 0.1822 - subIterW = 10 + elapsed = 0.027741 + tTotal = 0.1152 + subIterW = 16 subIterH = 11 - relError = 0.2655 + relError = 0.2850 toshow = scalar structure containing the fields: - initGrNormW = 370.94 - initGrNormH = 483.38 - initNorm = 14.465 - SC1 = 4.1689e-03 - SC2 = 2.9099e-03 - SC3 = 0.1848 - W_density = 0.7633 - H_density = 0.5950 + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 4.9932e-03 + SC2 = 3.3354e-03 + SC3 = 0.2092 + W_density = 0.7033 + H_density = 0.6000 iter = 4 - elapsed = 0.038839 - tTotal = 0.2210 - subIterW = 12 + elapsed = 0.032651 + tTotal = 0.1478 + subIterW = 10 subIterH = 8 - relError = 0.2565 + relError = 0.2771 toshow = scalar structure containing the fields: - initGrNormW = 370.94 - initGrNormH = 483.38 - initNorm = 14.465 - SC1 = 3.1008e-03 - SC2 = 2.1333e-03 - SC3 = 0.1455 - W_density = 0.7600 - H_density = 0.5800 + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 3.8791e-03 + SC2 = 2.5990e-03 + SC3 = 0.1619 + W_density = 0.7000 + H_density = 0.6050 iter = 5 - elapsed = 0.045453 - tTotal = 0.2665 - subIterW = 8 - subIterH = 8 - relError = 0.2520 + elapsed = 0.023427 + tTotal = 0.1712 + subIterW = 3 + subIterH = 6 + relError = 0.2718 toshow = scalar structure containing the fields: - initGrNormW = 370.94 - initGrNormH = 483.38 - initNorm = 14.465 - SC1 = 2.7181e-03 - SC2 = 1.8701e-03 - SC3 = 0.1327 - W_density = 0.7633 - H_density = 0.5750 + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 2.9954e-03 + SC2 = 2.0069e-03 + SC3 = 0.1311 + W_density = 0.7033 + H_density = 0.6000 iter = 6 - elapsed = 0.044260 - tTotal = 0.3107 - subIterW = 6 - subIterH = 4 - relError = 0.2491 + elapsed = 0.023493 + tTotal = 0.1947 + subIterW = 4 + subIterH = 5 + relError = 0.2684 toshow = scalar structure containing the fields: - initGrNormW = 370.94 - initGrNormH = 483.38 - initNorm = 14.465 - SC1 = 2.5214e-03 - SC2 = 1.7600e-03 - SC3 = 0.1220 - W_density = 0.7700 - H_density = 0.5700 + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 2.6560e-03 + SC2 = 1.7795e-03 + SC3 = 0.1136 + W_density = 0.6967 + H_density = 0.5950 iter = 7 - elapsed = 0.037552 - tTotal = 0.3483 - subIterW = 5 - subIterH = 5 - relError = 0.2467 + elapsed = 0.028667 + tTotal = 0.2234 + subIterW = 7 + subIterH = 2 + relError = 0.2663 toshow = scalar structure containing the fields: - initGrNormW = 370.94 - initGrNormH = 483.38 - initNorm = 14.465 - SC1 = 2.1632e-03 - SC2 = 1.4926e-03 - SC3 = 0.1084 - W_density = 0.7633 - H_density = 0.5800 + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 2.5449e-03 + SC2 = 1.7102e-03 + SC3 = 0.1136 + W_density = 0.7033 + H_density = 0.6050 iter = 8 - elapsed = 0.042783 - tTotal = 0.3911 + elapsed = 0.023540 + tTotal = 0.2469 subIterW = 2 subIterH = 5 - relError = 0.2448 + relError = 0.2646 toshow = scalar structure containing the fields: - initGrNormW = 370.94 - initGrNormH = 483.38 - initNorm = 14.465 - SC1 = 1.9158e-03 - SC2 = 1.3142e-03 - SC3 = 0.095452 - W_density = 0.7633 - H_density = 0.5650 + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 2.4953e-03 + SC2 = 1.6668e-03 + SC3 = 0.1165 + W_density = 0.7067 + H_density = 0.6000 iter = 9 - elapsed = 0.036713 - tTotal = 0.4278 - subIterW = 4 - subIterH = 3 - relError = 0.2435 + elapsed = 0.024407 + tTotal = 0.2713 + subIterW = 7 + subIterH = 5 + relError = 0.2632 toshow = scalar structure containing the fields: - initGrNormW = 370.94 - initGrNormH = 483.38 - initNorm = 14.465 - SC1 = 1.6972e-03 - SC2 = 1.1677e-03 - SC3 = 0.081861 - W_density = 0.7667 - H_density = 0.5700 + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 2.2560e-03 + SC2 = 1.5025e-03 + SC3 = 0.1146 + W_density = 0.7133 + H_density = 0.5750 iter = 10 - elapsed = 0.034476 - tTotal = 0.4623 - subIterW = 3 - subIterH = 0 - relError = 0.2425 + elapsed = 0.029578 + tTotal = 0.3009 + subIterW = 7 + subIterH = 7 + relError = 0.2619 toshow = scalar structure containing the fields: - initGrNormW = 370.94 - initGrNormH = 483.38 - initNorm = 14.465 - SC1 = 1.5352e-03 - SC2 = 1.0378e-03 - SC3 = 0.074573 - W_density = 0.7600 - H_density = 0.5500 + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 2.2634e-03 + SC2 = 1.5029e-03 + SC3 = 0.1114 + W_density = 0.7100 + H_density = 0.5750 iter = 11 - elapsed = 0.037410 - tTotal = 0.4997 - subIterW = 6 - subIterH = 4 - relError = 0.2418 + elapsed = 0.024476 + tTotal = 0.3254 + subIterW = 4 + subIterH = 2 + relError = 0.2607 toshow = scalar structure containing the fields: - initGrNormW = 370.94 - initGrNormH = 483.38 - initNorm = 14.465 - SC1 = 1.3956e-03 - SC2 = 9.4060e-04 - SC3 = 0.064266 - W_density = 0.7467 - H_density = 0.5550 + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 2.2886e-03 + SC2 = 1.5059e-03 + SC3 = 0.1116 + W_density = 0.7100 + H_density = 0.5800 iter = 12 - elapsed = 0.034967 - tTotal = 0.5346 - subIterW = 6 - subIterH = 0 - relError = 0.2412 + elapsed = 0.022866 + tTotal = 0.3483 + subIterW = 3 + subIterH = 3 + relError = 0.2595 toshow = scalar structure containing the fields: - initGrNormW = 370.94 - initGrNormH = 483.38 - initNorm = 14.465 - SC1 = 1.3101e-03 - SC2 = 8.8039e-04 - SC3 = 0.062124 - W_density = 0.7467 - H_density = 0.5500 + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 2.1699e-03 + SC2 = 1.4322e-03 + SC3 = 0.057181 + W_density = 0.7100 + H_density = 0.5750 iter = 13 - elapsed = 0.036167 - tTotal = 0.5708 - subIterW = 4 - subIterH = 3 - relError = 0.2407 + elapsed = 0.017404 + tTotal = 0.3657 + subIterW = 0 + subIterH = 2 + relError = 0.2584 toshow = scalar structure containing the fields: - initGrNormW = 370.94 - initGrNormH = 483.38 - initNorm = 14.465 - SC1 = 1.2850e-03 - SC2 = 8.7126e-04 - SC3 = 0.059912 - W_density = 0.7533 - H_density = 0.5600 + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 2.1133e-03 + SC2 = 1.4159e-03 + SC3 = 0.1025 + W_density = 0.7167 + H_density = 0.6000 iter = 14 - elapsed = 0.035048 - tTotal = 0.6058 - subIterW = 2 - subIterH = 2 - relError = 0.2403 + elapsed = 0.023937 + tTotal = 0.3896 + subIterW = 6 + subIterH = 4 + relError = 0.2573 toshow = scalar structure containing the fields: - initGrNormW = 370.94 - initGrNormH = 483.38 - initNorm = 14.465 - SC1 = 1.2159e-03 - SC2 = 8.2437e-04 - SC3 = 0.058603 - W_density = 0.7600 - H_density = 0.5500 + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 2.1573e-03 + SC2 = 1.4411e-03 + SC3 = 0.1040 + W_density = 0.7167 + H_density = 0.5950 iter = 15 - elapsed = 0.035891 - tTotal = 0.6417 - subIterW = 4 - subIterH = 3 - relError = 0.2399 + elapsed = 0.021818 + tTotal = 0.4114 + subIterW = 6 + subIterH = 0 + relError = 0.2563 toshow = scalar structure containing the fields: - initGrNormW = 370.94 - initGrNormH = 483.38 - initNorm = 14.465 - SC1 = 1.2209e-03 - SC2 = 8.2044e-04 - SC3 = 0.029520 - W_density = 0.7600 - H_density = 0.5400 + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 2.1723e-03 + SC2 = 1.4294e-03 + SC3 = 0.1023 + W_density = 0.7033 + H_density = 0.5900 iter = 16 - elapsed = 0.028280 - tTotal = 0.6700 - subIterW = 0 - subIterH = 4 - relError = 0.2395 + elapsed = 0.022879 + tTotal = 0.4343 + subIterW = 6 + subIterH = 3 + relError = 0.2552 toshow = scalar structure containing the fields: - initGrNormW = 370.94 - initGrNormH = 483.38 - initNorm = 14.465 - SC1 = 1.2712e-03 - SC2 = 8.5425e-04 - SC3 = 0.064073 - W_density = 0.7567 - H_density = 0.5450 + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 1.9911e-03 + SC2 = 1.3061e-03 + SC3 = 0.098163 + W_density = 0.7067 + H_density = 0.5750 iter = 17 - elapsed = 0.032509 - tTotal = 0.7025 - subIterW = 0 - subIterH = 0 - relError = 0.2391 + elapsed = 0.023159 + tTotal = 0.4575 + subIterW = 7 + subIterH = 3 + relError = 0.2542 toshow = scalar structure containing the fields: - initGrNormW = 370.94 - initGrNormH = 483.38 - initNorm = 14.465 - SC1 = 1.3195e-03 - SC2 = 8.7617e-04 - SC3 = 0.065782 - W_density = 0.7467 - H_density = 0.5400 + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 1.9436e-03 + SC2 = 1.2711e-03 + SC3 = 0.093017 + W_density = 0.7033 + H_density = 0.5800 iter = 18 - elapsed = 0.035508 - tTotal = 0.7380 - subIterW = 5 + elapsed = 0.022342 + tTotal = 0.4798 + subIterW = 7 subIterH = 0 - relError = 0.2387 + relError = 0.2532 toshow = scalar structure containing the fields: - initGrNormW = 370.94 - initGrNormH = 483.38 - initNorm = 14.465 - SC1 = 1.2951e-03 - SC2 = 8.5736e-04 - SC3 = 0.065696 - W_density = 0.7467 - H_density = 0.5300 + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 1.9007e-03 + SC2 = 1.2544e-03 + SC3 = 0.085313 + W_density = 0.7067 + H_density = 0.5900 iter = 19 - elapsed = 0.035068 - tTotal = 0.7731 - subIterW = 2 + elapsed = 0.023059 + tTotal = 0.5029 + subIterW = 5 subIterH = 2 - relError = 0.2382 + relError = 0.2522 toshow = scalar structure containing the fields: - initGrNormW = 370.94 - initGrNormH = 483.38 - initNorm = 14.465 - SC1 = 1.1608e-03 - SC2 = 7.6615e-04 - SC3 = 0.062257 - W_density = 0.7433 - H_density = 0.5300 + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 1.8720e-03 + SC2 = 1.2393e-03 + SC3 = 0.043546 + W_density = 0.7067 + H_density = 0.5900 iter = 20 - elapsed = 0.033781 - tTotal = 0.8069 + elapsed = 0.017761 + tTotal = 0.5206 subIterW = 0 subIterH = 2 - relError = 0.2378 + relError = 0.2513 toshow = scalar structure containing the fields: - initGrNormW = 370.94 - initGrNormH = 483.38 - initNorm = 14.465 - SC1 = 1.1228e-03 - SC2 = 7.3880e-04 - SC3 = 0.060149 - W_density = 0.7400 - H_density = 0.5350 + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 1.8002e-03 + SC2 = 1.1953e-03 + SC3 = 0.078816 + W_density = 0.7167 + H_density = 0.5850 iter = 21 - elapsed = 0.033165 - tTotal = 0.8400 - subIterW = 0 - subIterH = 0 - relError = 0.2374 + elapsed = 0.020981 + tTotal = 0.5416 + subIterW = 3 + subIterH = 3 + relError = 0.2504 toshow = scalar structure containing the fields: - initGrNormW = 370.94 - initGrNormH = 483.38 - initNorm = 14.465 - SC1 = 1.0318e-03 - SC2 = 6.7479e-04 - SC3 = 0.058887 - W_density = 0.7400 - H_density = 0.5250 + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 1.7282e-03 + SC2 = 1.1406e-03 + SC3 = 0.077945 + W_density = 0.7167 + H_density = 0.5700 iter = 22 - elapsed = 0.037291 - tTotal = 0.8773 + elapsed = 0.022477 + tTotal = 0.5641 subIterW = 4 - subIterH = 2 - relError = 0.2371 + subIterH = 5 + relError = 0.2496 toshow = scalar structure containing the fields: - initGrNormW = 370.94 - initGrNormH = 483.38 - initNorm = 14.465 - SC1 = 9.8074e-04 - SC2 = 6.3944e-04 - SC3 = 0.055774 - W_density = 0.7367 - H_density = 0.5250 + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 1.6191e-03 + SC2 = 1.0589e-03 + SC3 = 0.074924 + W_density = 0.7133 + H_density = 0.5650 iter = 23 - elapsed = 0.027505 - tTotal = 0.9048 - subIterW = 0 - subIterH = 0 - relError = 0.2368 + elapsed = 0.025905 + tTotal = 0.5900 + subIterW = 3 + subIterH = 3 + relError = 0.2489 toshow = scalar structure containing the fields: - initGrNormW = 370.94 - initGrNormH = 483.38 - initNorm = 14.465 - SC1 = 9.0326e-04 - SC2 = 5.9073e-04 - SC3 = 0.024360 - W_density = 0.7367 - H_density = 0.5250 + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 1.5288e-03 + SC2 = 1.0090e-03 + SC3 = 0.070108 + W_density = 0.7200 + H_density = 0.5650 iter = 24 - elapsed = 0.021187 - tTotal = 0.9260 + elapsed = 0.022345 + tTotal = 0.6123 + subIterW = 4 + subIterH = 2 + relError = 0.2482 + +toshow = + + scalar structure containing the fields: + + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 1.3834e-03 + SC2 = 9.1580e-04 + SC3 = 0.066612 + W_density = 0.7233 + H_density = 0.5700 + iter = 25 + elapsed = 0.021807 + tTotal = 0.6341 subIterW = 0 + subIterH = 4 + relError = 0.2476 + +toshow = + + scalar structure containing the fields: + + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 1.3107e-03 + SC2 = 8.6505e-04 + SC3 = 0.063716 + W_density = 0.7233 + H_density = 0.5600 + iter = 26 + elapsed = 0.020814 + tTotal = 0.6550 + subIterW = 2 + subIterH = 2 + relError = 0.2470 + +toshow = + + scalar structure containing the fields: + + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 1.2746e-03 + SC2 = 8.3616e-04 + SC3 = 0.060279 + W_density = 0.7133 + H_density = 0.5650 + iter = 27 + elapsed = 0.021437 + tTotal = 0.6764 + subIterW = 3 subIterH = 0 - relError = 0.2365 + relError = 0.2465 toshow = scalar structure containing the fields: - initGrNormW = 370.94 - initGrNormH = 483.38 - initNorm = 14.465 - SC1 = 8.4199e-04 - SC2 = 5.4898e-04 - SC3 = 0.023238 - W_density = 0.7367 - H_density = 0.5250 - iter = 25 - elapsed = 0.028835 - tTotal = 0.9549 + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 1.1780e-03 + SC2 = 7.7984e-04 + SC3 = 0.056971 + W_density = 0.7233 + H_density = 0.5700 + iter = 28 + elapsed = 0.022198 + tTotal = 0.6986 + subIterW = 5 + subIterH = 3 + relError = 0.2461 + +toshow = + + scalar structure containing the fields: + + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 1.1040e-03 + SC2 = 7.3306e-04 + SC3 = 0.053831 + W_density = 0.7267 + H_density = 0.5650 + iter = 29 + elapsed = 0.024847 + tTotal = 0.7234 + subIterW = 3 + subIterH = 0 + relError = 0.2457 + +toshow = + + scalar structure containing the fields: + + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 1.0428e-03 + SC2 = 6.8825e-04 + SC3 = 0.051936 + W_density = 0.7233 + H_density = 0.5650 + iter = 30 + elapsed = 0.021758 + tTotal = 0.7452 subIterW = 0 + subIterH = 4 + relError = 0.2453 + +toshow = + + scalar structure containing the fields: + + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 1.0159e-03 + SC2 = 6.7252e-04 + SC3 = 0.050619 + W_density = 0.7200 + H_density = 0.5750 + iter = 31 + elapsed = 0.019749 + tTotal = 0.7649 + subIterW = 0 + subIterH = 2 + relError = 0.2450 + +toshow = + + scalar structure containing the fields: + + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 9.9242e-04 + SC2 = 6.5499e-04 + SC3 = 0.049033 + W_density = 0.7133 + H_density = 0.5800 + iter = 32 + elapsed = 0.021115 + tTotal = 0.7861 + subIterW = 2 + subIterH = 0 + relError = 0.2447 + +toshow = + + scalar structure containing the fields: + + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 9.7165e-04 + SC2 = 6.4323e-04 + SC3 = 0.048606 + W_density = 0.7133 + H_density = 0.5800 + iter = 33 + elapsed = 0.022351 + tTotal = 0.8084 + subIterW = 2 + subIterH = 2 + relError = 0.2444 + +toshow = + + scalar structure containing the fields: + + initGrNormW = 364.77 + initGrNormH = 430.32 + initNorm = 14.111 + SC1 = 9.9590e-04 + SC2 = 6.5730e-04 + SC3 = 0.050761 + W_density = 0.7167 + H_density = 0.5700 + iter = 34 + elapsed = 0.020574 + tTotal = 0.8290 + subIterW = 1 subIterH = 2 - relError = 0.2363 + relError = 0.2441 Stop: tolerance reached. final = scalar structure containing the fields: - iterations = 25 - elapsed_total = 0.9549 - relative_error = 0.2363 - W_density = 0.7367 - H_density = 0.5250 + iterations = 34 + elapsed_total = 0.8290 + relative_error = 0.2441 + W_density = 0.7167 + H_density = 0.5700 ***** demo m = 300; @@ -5323,14 +5518,16 @@ 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/funm.m] ->>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/funm.m -***** function b = fsin (a) - b = sin (a); -***** endfunction +[inst/smwsolve.m] +>>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/smwsolve.m ***** test - mtx = randn (100); - assert (funm (mtx, "fsin"), thfm (mtx, "sin"), 1e-9) + 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/vec_projection.m] >>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/vec_projection.m @@ -5342,26 +5539,36 @@ 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/cod.m] ->>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/cod.m -***** test - a = rand (5, 10); - [q, r, z] = cod (a); - assert (norm (q*r*z' - a) / norm (a) < 1e-10); -***** test - a = rand (5, 10) + i * rand (5, 10); - [q, r, z] = cod (a); - assert (norm (q*r*z' - a) / norm (a) < 1e-10); +[inst/@blksparse/ctranspose.m] +>>>>> /build/reproducible-path/octave-linear-algebra-2.2.3/inst/@blksparse/ctranspose.m ***** test - a = rand (5, 10); - [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_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); + 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_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/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_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 Checking C++ files ... Done running the unit tests. @@ -5393,12 +5600,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/12084/tmp/hooks/B01_cleanup starting +I: user script /srv/workspace/pbuilder/12084/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/22840 and its subdirectories -I: Current time: Wed Dec 11 05:53:42 -12 2024 -I: pbuilder-time-stamp: 1733939622 +I: removing directory /srv/workspace/pbuilder/12084 and its subdirectories +I: Current time: Thu Dec 12 07:59:41 +14 2024 +I: pbuilder-time-stamp: 1733939981