mpskyl.f90 File Reference

Skyline storage. More...

Go to the source code of this file.

Modules
module	mpskyl
	QL data.

Functions/Subroutines
subroutine	mpslan (n, lrows, npgrp, nsparr)
	Analyse skyline stucture. More...
subroutine	mpslsp (npgrp, nsparr, nsparc)
	Analyse skyline stucture. More...
subroutine	mpsldc (npgrp, nsparr, nsparc, gmat, nsize, nrank, mon)
	Cholesky decomposition. More...
subroutine	mpslsv (gij, gvec, nfit)
	Backward, forward substitution. More...
subroutine	mpslsv2 (npgrp, nsparr, nsparc, gmat, gvec, nsize)
	Backward, forward substitution. More...

Variables
integer(mpi)	mpskyl::npar
	number of parameters More...
integer(mpi), dimension(:), allocatable	mpskyl::lastrowincol
	last (non-zero) row in column More...
integer(mpi), dimension(:), allocatable	mpskyl::pargroupstoregions
	mapping of parameter groups to regions (in row) More...
integer(mpi), dimension(:), allocatable	mpskyl::regionoffsets
	offsets for parameter group regions (in row) More...
integer(mpi)	mpskyl::nthrd
	number of threads More...

Detailed Description

Skyline storage.

Author

Claus Kleinwort, DESY, 2023 (Claus.nosp@m..Kle.nosp@m.inwor.nosp@m.t@de.nosp@m.sy.de)

Copyright

Copyright (c) 2023 Deutsches Elektronen-Synchroton, Member of the Helmholtz Association, (DESY), HAMBURG, GERMANY

This library is free software; you can redistribute it and/or modify it under the terms of the GNU Library General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version.

This library 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 Library General Public License for more details.

You should have received a copy of the GNU Library General Public License along with this program (see the file COPYING.LIB for more details); if not, write to the Free Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.

Sparse storage for (column based) skyline matrix with Cholesky decomposition.

Definition in file mpskyl.f90.

Function/Subroutine Documentation

mpslan()

subroutine mpslan	(	integer(mpi), intent(in)	n,
		integer(mpi), dimension(:), intent(in)	lrows,
		integer(mpi), dimension(:), intent(in)	npgrp,
		integer(mpl), dimension(:,:), intent(out)	nsparr
	)

Analyse skyline stucture.

Parameters

[in]	n	number of rows (parameters)
[in]	lrows	last row in columns
[in]	npgrp	parameter groups
[out]	nsparr	row offsets

Definition at line 45 of file mpskyl.f90.

References mpskyl::lastrowincol, mpskyl::npar, mpskyl::nthrd, mpskyl::pargroupstoregions, and mpskyl::regionoffsets.

mpsldc()

subroutine mpsldc	(	integer(mpi), dimension(:), intent(in)	npgrp,
		integer(mpl), dimension(:,:), intent(in)	nsparr,
		integer(mpi), dimension(:), intent(in)	nsparc,
		real(mpd), dimension(:), intent(inout)	gmat,
		integer(mpi), intent(in)	nsize,
		integer(mpi), intent(out)	nrank,
		integer(mpi), intent(in)	mon
	)

Cholesky decomposition.

Parameters

[in]	npgrp	parameter groups
[in]	nsparr	row offsets
[in]	nsparc	column offsets
[in,out]	gmat	(skyline) matrix, replaced by decomposition
[in]	nsize	size of matrix
[out]	nrank	rank of matrix
[in]	mon	monitoring flag

Definition at line 187 of file mpskyl.f90.

References monpgs(), mpskyl::npar, mpskyl::pargroupstoregions, and mpskyl::regionoffsets.

mpslsp()

subroutine mpslsp	(	integer(mpi), dimension(:), intent(in)	npgrp,
		integer(mpl), dimension(:,:), intent(in)	nsparr,
		integer(mpi), dimension(:), intent(out)	nsparc
	)

Analyse skyline stucture.

Parameters

[in]	npgrp	parameter groups
[in]	nsparr	row offsets
[out]	nsparv	column offsets

Definition at line 120 of file mpskyl.f90.

References mpskyl::lastrowincol, mpskyl::npar, and mpskyl::nthrd.

mpslsv()

subroutine mpslsv	(	real(mpd)	gij,
		real(mpd), dimension(:), intent(inout)	gvec,
		integer(mpi), intent(in)	nfit
	)

Backward, forward substitution.

Parameters

[in]	gij	external procedure access global matrix element (i,j)
[in,out]	gvec	r.h.s, replaced by solution
[in]	nfit	number of parameters

Definition at line 390 of file mpskyl.f90.

mpslsv2()

subroutine mpslsv2	(	integer(mpi), dimension(:), intent(in)	npgrp,
		integer(mpl), dimension(:,:), intent(in)	nsparr,
		integer(mpi), dimension(:), intent(in)	nsparc,
		real(mpd), dimension(:), intent(in)	gmat,
		real(mpd), dimension(:), intent(inout)	gvec,
		integer(mpi), intent(in)	nsize
	)

Backward, forward substitution.

Parameters

[in]	npgrp	parameter groups
[in]	nsparr	row offsets
[in]	nsparc	column offsets
[in]	gmat	(skyline) matrix, replaced by decomposition
[in,out]	gvec	r.h.s, replaced by solution
[in]	nsize	size of matrix

Definition at line 437 of file mpskyl.f90.

References mpskyl::npar, mpskyl::pargroupstoregions, and mpskyl::regionoffsets.