laplace.py
a2f123cbd58emaster

laplace.py
View Options

			# Copyright (C) 2018 by the RROMPy authors
			#
			# This file is part of RROMPy.
			#
			# RROMPy is free software: you can redistribute it and/or modify
			# it under the terms of the GNU Lesser General Public License as published by
			# the Free Software Foundation, either version 3 of the License, or
			# (at your option) any later version.
			#
			# RROMPy is distributed in the hope that it will be useful,
			# but WITHOUT ANY WARRANTY; without even the implied warranty of
			# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
			# GNU Lesser General Public License for more details.
			#
			# You should have received a copy of the GNU Lesser General Public License
			# along with RROMPy. If not, see <http://www.gnu.org/licenses/>.
			#

			import numpy as np
			from rrompy.hfengines.linear_problem import (LaplaceDiskGaussian,
			LaplaceDiskGaussian2)

			def test_laplace_disk():
			solver = LaplaceDiskGaussian(n = 20, verbosity = 0)
			mu = 1.5
			solver.setSolver("BICG", {"tol" : 1e-15})
			uh = solver.solve(mu)
			assert np.isclose(solver.norm(uh), 1.0534030774205372, rtol = 1e-5)
			assert np.isclose(solver.norm(solver.residual(uh, mu)), 5.27345e-13,
			rtol = 1e-1)

			def test_laplace_disk_2():
			solver = LaplaceDiskGaussian2(n = 20, verbosity = 0)
			mu = tuple([0., 1.5])
			uh = solver.solve(mu)
			assert np.isclose(solver.norm(uh), 1.0534030774205372, rtol = 1e-5)
			assert np.isclose(solver.norm(solver.residual(uh, mu)), 5.27345e-13,
			rtol = 1e-1)