modified installation guide, added ubuntu installation hints and arpack... (d22be54b) · Commits · xerus / xerus

doc/jekyll/_posts/2000-12-20-building_xerus.md

+36 −0

Original line number	Diff line number	Diff line
		@@ -27,6 +27,42 @@ E.g. to install all dependencies on a fedora system execute
		dnf install gcc-c++ openblas-devel suitesparse-devel lapack-devel boost-devel binutils-devel
		~~~

		Or to install all dependencies for Ubuntu (tested for 18.04) execute:
		~~~ bash
		apt-get update
		apt install make g++ libboost-all-dev binutils-dev libsuitesparse-dev libz-dev libiberty-dev libopenblas-dev liblapack-dev liblapacke-dev gfortran
		~~~

		To build the optional ARPACK extension, you need to install ARPACK and the c++ headers by hand. First clone the git project to an appropriate place, enter the project directory, configure the sources and make the library:
		~~~ bash
		git clone https://github.com/opencollab/arpack-ng.git
		cd arpack-ng
		apt-get install autoconf
		sh bootstrap
		./configure F77=gfortran --enable-icb
		make
		make check
		~~~
		You can add your own installation path with --prefix=<installation-path> for the configure command.
		Unfortunately, the C wrapper is broken. You need to modify the file_exists arpack.hpp:
		~~~ bash
		cd ICB/
		vi arpack.hpp
		~~~
		In arpack.hpp comment out all 4 instances of the function neupd.

		After that, to install the packages type:
		~~~ bash
		make install
		~~~

		Once ARPACK is installed you need to modify the xerus config file. Just remove the hashtags in front of ARPACK_LIBRARIES and the addtional flag below:
		~~~ bash
		# (optional) Uncomment if needed, iterative eigenvalue solver ARPACK-ng, see https://github.com/opencollab/arpack-ng
		ARPACK_LIBRARIES = -larpack
		OTHER += -DARPACK_LIBRARIES
		~~~

		To build the python bindings you will furthermore need the python development headers, `numpy` as well as `boost-python` or
		`boost-python3` depending on the python version you wish to use. E.g. for a fedora system and if you want to use python 2 simply execute
		~~~ bash

src/xerus/htNetwork.cpp

+10 −4

Original line number	Diff line number	Diff line
		@@ -168,6 +168,7 @@ namespace xerus {
		}
		dimensions = _tensor.dimensions;
		Tensor remains;
		auto epsPerSite = numComp < 2 ? _eps : misc::hard_equal(_eps, 0.0) ? EPSILON : _eps / std::sqrt(double(2*numComp-3));



		@@ -213,7 +214,8 @@ namespace xerus {

		xerus::reshuffle(remains, remains, ithmode);

		calculate_svd(newNode, singularValues, remains, remains, N, _maxRanks[pos - 1], _eps);

		calculate_svd(newNode, singularValues, remains, remains, N, _maxRanks[pos - 1], epsPerSite);
		if (isOperator){
		xerus::reshuffle(newNode, newNode, {1,2,0});
		} else {
		@@ -247,7 +249,7 @@ namespace xerus {
		}
		xerus::reshuffle(remains, remains, ithmode);

		calculate_svd(newNode, singularValues, remains, remains, 2, _maxRanks[numCompOnLvl + pos - 2], _eps); // TODO fix maxRanks
		calculate_svd(newNode, singularValues, remains, remains, 2, _maxRanks[numCompOnLvl + pos - 2], epsPerSite); // TODO fix maxRanks
		xerus::reshuffle(newNode, newNode, {1,2,0}); // first parent then children
		set_component(numCompOnLvl + pos - 1, std::move(newNode));
		newNode.reset();
		@@ -727,8 +729,12 @@ namespace xerus {

		canonicalize_root();

		if(numComponents > 1) {
		auto epsPerSite = misc::hard_equal(_eps, 0.0) ? EPSILON : _eps / std::sqrt(double(2*numComponents-3)); // Taken from HIERARCHICAL SINGULAR VALUE DECOMPOSITION OF TENSORS
		// from Lars Grasedyck
		for (size_t n = numComponents - 1; n > 0; --n) {
		round_edge(n, (n + 1) / 2 - 1, _maxRanks[n - 1], _eps, 0.0);
		round_edge(n, (n + 1) / 2 - 1, _maxRanks[n - 1], epsPerSite, 0.0);
		}
		}

		assume_core_position(0);