refactoring

conrhs.tpl
View Options

			!>
			!> @file conrhs.tpl
			!>
			!> @brief
			!>
			!> @copyright
			!> Copyright (©) 2021 EPFL (Ecole Polytechnique Fédérale de Lausanne)
			!> SPC (Swiss Plasma Center)
			!>
			!> spclibs 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.
			!>
			!> spclibs 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 General Public License for more details.
			!>
			!> You should have received a copy of the GNU Lesser General Public License
			!> along with this program. If not, see <https://www.gnu.org/licenses/>.
			!>
			!> @authors
			!> (in alphabetical order)
			!> @author Trach-Minh Tran <trach-minh.tran@epfl.ch>
			!>
			DOUBLE PRECISION, POINTER :: xg(:,:), wg(:,:)
			DOUBLE PRECISION, ALLOCATABLE :: fun(:,:)
			INTEGER :: ndim, n, nidbas, ng
			INTEGER :: i, ig, it, irow
			LOGICAL :: nlper
			!
			CALL get_dim(spl, ndim, n, nidbas)
			nlper = spl%period
			xg => spl%gausx ! xg(ng,n)
			wg => spl%gausw ! wg(ng,n)
			ng = SIZE(xg,1)
			ALLOCATE(fun(0:nidbas,1))
			!
			farr = 0.0d0
			DO i=1,n
			DO ig=1,ng
			CALL basfun(xg(ig,i), spl, fun, i)
			contrib = wg(ig,i)*frhs(xg(ig,i))
			DO it=0,nidbas
			irow = i+it
			IF(nlper) irow = MODULO(irow-1,n) +1
			farr(irow) = farr(irow)+contrib*fun(it,1)
			END DO
			END DO
			END DO
			!
			DEALLOCATE(fun)