gaussquadrature Module Reference

Module for computation of 1d and tensorized 2d Gauss points. More...

Data Types
type	gaussq

Functions/Subroutines
subroutine	init_gaussq (this, lun_err, ngauss)
	Static constructor for gaussquadrature::gaussq
subroutine	kill_gaussq (this, lun_err)
	Static destructor for gaussquadrature::gaussq
subroutine	info_gaussq (this, lun_out)
	Info procedure for gaussquadrature::gaussq
subroutine	on_interval (this, a, b)
	Compute Gauss points on the interval \([a,b]\).
subroutine	on_cell (this, ax, bx, ay, by)
	Compute Gauss points of \([a_x,b_x] \times [a_y,b_y]\).

Detailed Description

Module for computation of 1d and tensorized 2d Gauss points.

Function/Subroutine Documentation

init_gaussq()

subroutine gaussquadrature::init_gaussq	(	class(gaussq), intent(inout)	this,
		integer, intent(in)	lun_err,
		integer, intent(in)	ngauss )

Static constructor for gaussquadrature::gaussq

---

Allocate all the variables. Define variable gaussq::coeff1d and gaussq::weight1d.

Parameters

[in]	lun_err	unit number for error messages
[in]	ngauss	Number of Gauss points

kill_gaussq()

subroutine gaussquadrature::kill_gaussq ( class(gaussq), intent(inout) this, integer, intent(in) lun_err )

private

Static destructor for gaussquadrature::gaussq

---

Deallocate all the variables.

Parameters

[in]

lun_err

unit number for error messagges

info_gaussq()

subroutine gaussquadrature::info_gaussq ( class(gaussq), intent(in) this, integer, intent(in) lun_out )

private

Info procedure for gaussquadrature::gaussq

---

Parameters

[in]

lun_out

unit number for output messagges

on_interval()

subroutine gaussquadrature::on_interval ( class(gaussq), intent(inout) this, real(double), intent(in) a, real(double), intent(in) b )

private

Compute Gauss points on the interval \([a,b]\).

---

Define variables gaussq::coord_ab and gaussq::weight_ab.

Parameters

[in]	a	left endpoint of interval
[in]	b	right endpoint of interval

on_cell()

subroutine gaussquadrature::on_cell ( class(gaussq), intent(inout) this, real(double), intent(in) ax, real(double), intent(in) bx, real(double), intent(in) ay, real(double), intent(in) by )

private

Compute Gauss points of \([a_x,b_x] \times [a_y,b_y]\).

---

Gauss points in \(\mathbb{R}^{2}\) are computed by tensorization of 1d Gauss points. Define variables gaussq::coord_cell and gaussq::weight_cell.

Parameters

[in]	ax	left endpoint interval variable \(x\).
[in]	bx	right endpoint interval variable \(x\).
[in]	ay	left endpoint interval variable \(y\).
[in]	by	right endpoint interval variable y.