forcad_nurbs_volume

This module defines the 'nurbs_volume' type for representing a Non-Uniform Rational B-Spline (NURBS) volume.

Uses

- forcad_kinds
- forcad_utils
- fordebug
Help

Graph Key

Nodes of different colours represent the following:

Solid arrows point from a submodule to the (sub)module which it is descended from. Dashed arrows point from a module or program unit to modules which it uses.
Where possible, edges connecting nodes are given different colours to make them easier to distinguish in large graphs.

Interfaces

public interface compute_Tgc

private pure function compute_Tgc_nurbs_3d_vector(Xt, knot1, knot2, knot3, degree, nc, ng, Wc) result(Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=rk),	intent(in),		contiguous	::	Xt(:,:)
real(kind=rk),	intent(in),		contiguous	::	knot1(:)
real(kind=rk),	intent(in),		contiguous	::	knot2(:)
real(kind=rk),	intent(in),		contiguous	::	knot3(:)
integer,	intent(in)			::	degree(3)
integer,	intent(in)			::	nc(3)
integer,	intent(in),	optional		::	ng(3)
real(kind=rk),	intent(in),		contiguous	::	Wc(:)

Return Value real(kind=rk), allocatable, (:,:)

private pure function compute_Tgc_bspline_3d_vector(Xt, knot1, knot2, knot3, degree, nc, ng) result(Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=rk),	intent(in),		contiguous	::	Xt(:,:)
real(kind=rk),	intent(in),		contiguous	::	knot1(:)
real(kind=rk),	intent(in),		contiguous	::	knot2(:)
real(kind=rk),	intent(in),		contiguous	::	knot3(:)
integer,	intent(in)			::	degree(3)
integer,	intent(in)			::	nc(3)
integer,	intent(in),	optional		::	ng(3)

Return Value real(kind=rk), allocatable, (:,:)

private pure function compute_Tgc_nurbs_3d_scalar(Xt, knot1, knot2, knot3, degree, nc, Wc) result(Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Attributes		Name
real(kind=rk),	intent(in),	contiguous	::	Xt(:)
real(kind=rk),	intent(in),	contiguous	::	knot1(:)
real(kind=rk),	intent(in),	contiguous	::	knot2(:)
real(kind=rk),	intent(in),	contiguous	::	knot3(:)
integer,	intent(in)		::	degree(3)
integer,	intent(in)		::	nc(3)
real(kind=rk),	intent(in),	contiguous	::	Wc(:)

Return Value real(kind=rk), allocatable, (:)

private pure function compute_Tgc_bspline_3d_scalar(Xt, knot1, knot2, knot3, degree, nc) result(Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Attributes		Name
real(kind=rk),	intent(in),	contiguous	::	Xt(:)
real(kind=rk),	intent(in),	contiguous	::	knot1(:)
real(kind=rk),	intent(in),	contiguous	::	knot2(:)
real(kind=rk),	intent(in),	contiguous	::	knot3(:)
integer,	intent(in)		::	degree(3)
integer,	intent(in)		::	nc(3)

Return Value real(kind=rk), allocatable, (:)

public interface compute_dTgc

private pure subroutine compute_dTgc_nurbs_3d_vector(Xt, knot1, knot2, knot3, degree, nc, ng, Wc, dTgc, Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=rk),	intent(in),		contiguous	::	Xt(:,:)
real(kind=rk),	intent(in),		contiguous	::	knot1(:)
real(kind=rk),	intent(in),		contiguous	::	knot2(:)
real(kind=rk),	intent(in),		contiguous	::	knot3(:)
integer,	intent(in)			::	degree(3)
integer,	intent(in)			::	nc(3)
integer,	intent(in),	optional		::	ng(3)
real(kind=rk),	intent(in),		contiguous	::	Wc(:)
real(kind=rk),	intent(out),		allocatable	::	dTgc(:,:,:)
real(kind=rk),	intent(out),		allocatable	::	Tgc(:,:)

private pure subroutine compute_dTgc_bspline_3d_vector(Xt, knot1, knot2, knot3, degree, nc, ng, dTgc, Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=rk),	intent(in),		contiguous	::	Xt(:,:)
real(kind=rk),	intent(in),		contiguous	::	knot1(:)
real(kind=rk),	intent(in),		contiguous	::	knot2(:)
real(kind=rk),	intent(in),		contiguous	::	knot3(:)
integer,	intent(in)			::	degree(3)
integer,	intent(in)			::	nc(3)
integer,	intent(in),	optional		::	ng(3)
real(kind=rk),	intent(out),		allocatable	::	dTgc(:,:,:)
real(kind=rk),	intent(out),		allocatable	::	Tgc(:,:)

private pure subroutine compute_dTgc_nurbs_3d_scalar(Xt, knot1, knot2, knot3, degree, nc, Wc, dTgc, Tgc, elem)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

If elem is not present: Wc refers to the full weight vector. If elem is present: Wc refers to the element-local weight vector (Wce).

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=rk),	intent(in),		contiguous	::	Xt(:)
real(kind=rk),	intent(in),		contiguous	::	knot1(:)
real(kind=rk),	intent(in),		contiguous	::	knot2(:)
real(kind=rk),	intent(in),		contiguous	::	knot3(:)
integer,	intent(in)			::	degree(3)
integer,	intent(in)			::	nc(3)
real(kind=rk),	intent(in),		contiguous	::	Wc(:)
real(kind=rk),	intent(out),		allocatable	::	dTgc(:,:)
real(kind=rk),	intent(out),		allocatable	::	Tgc(:)
integer,	intent(in),	optional,	contiguous	::	elem(:)

private pure subroutine compute_dTgc_bspline_3d_scalar(Xt, knot1, knot2, knot3, degree, nc, dTgc, Tgc, elem)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=rk),	intent(in),		contiguous	::	Xt(:)
real(kind=rk),	intent(in),		contiguous	::	knot1(:)
real(kind=rk),	intent(in),		contiguous	::	knot2(:)
real(kind=rk),	intent(in),		contiguous	::	knot3(:)
integer,	intent(in)			::	degree(3)
integer,	intent(in)			::	nc(3)
real(kind=rk),	intent(out),		allocatable	::	dTgc(:,:)
real(kind=rk),	intent(out),		allocatable	::	Tgc(:)
integer,	intent(in),	optional,	contiguous	::	elem(:)

private interface compute_Xg

private pure function compute_Xg_nurbs_3d(Xt, knot1, knot2, knot3, degree, nc, ng, Xc, Wc) result(Xg)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=rk),	intent(in),		contiguous	::	Xt(:,:)
real(kind=rk),	intent(in),		contiguous	::	knot1(:)
real(kind=rk),	intent(in),		contiguous	::	knot2(:)
real(kind=rk),	intent(in),		contiguous	::	knot3(:)
integer,	intent(in)			::	degree(3)
integer,	intent(in)			::	nc(3)
integer,	intent(in),	optional		::	ng(3)
real(kind=rk),	intent(in),		contiguous	::	Xc(:,:)
real(kind=rk),	intent(in),		contiguous	::	Wc(:)

Return Value real(kind=rk), allocatable, (:,:)

private pure function compute_Xg_bspline_3d(Xt, knot1, knot2, knot3, degree, nc, ng, Xc) result(Xg)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=rk),	intent(in),		contiguous	::	Xt(:,:)
real(kind=rk),	intent(in),		contiguous	::	knot1(:)
real(kind=rk),	intent(in),		contiguous	::	knot2(:)
real(kind=rk),	intent(in),		contiguous	::	knot3(:)
integer,	intent(in)			::	degree(3)
integer,	intent(in)			::	nc(3)
integer,	intent(in),	optional		::	ng(3)
real(kind=rk),	intent(in),		contiguous	::	Xc(:,:)

Return Value real(kind=rk), allocatable, (:,:)

private pure function compute_Xg_nurbs_3d_1point(Xt, knot1, knot2, knot3, degree, nc, Xc, Wc) result(Xg)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Attributes		Name
real(kind=rk),	intent(in),	contiguous	::	Xt(:)
real(kind=rk),	intent(in),	contiguous	::	knot1(:)
real(kind=rk),	intent(in),	contiguous	::	knot2(:)
real(kind=rk),	intent(in),	contiguous	::	knot3(:)
integer,	intent(in)		::	degree(3)
integer,	intent(in)		::	nc(3)
real(kind=rk),	intent(in),	contiguous	::	Xc(:,:)
real(kind=rk),	intent(in),	contiguous	::	Wc(:)

Return Value real(kind=rk), (size(Xc,2))

private pure function compute_Xg_bspline_3d_1point(Xt, knot1, knot2, knot3, degree, nc, Xc) result(Xg)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Attributes		Name
real(kind=rk),	intent(in),	contiguous	::	Xt(:)
real(kind=rk),	intent(in),	contiguous	::	knot1(:)
real(kind=rk),	intent(in),	contiguous	::	knot2(:)
real(kind=rk),	intent(in),	contiguous	::	knot3(:)
integer,	intent(in)		::	degree(3)
integer,	intent(in)		::	nc(3)
real(kind=rk),	intent(in),	contiguous	::	Xc(:,:)

Return Value real(kind=rk), (size(Xc,2))

private interface compute_d2Tgc

private pure subroutine compute_d2Tgc_nurbs_3d_vector(Xt, knot1, knot2, knot3, degree, nc, ng, Wc, d2Tgc, dTgc, Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=rk),	intent(in),		contiguous	::	Xt(:,:)
real(kind=rk),	intent(in),		contiguous	::	knot1(:)
real(kind=rk),	intent(in),		contiguous	::	knot2(:)
real(kind=rk),	intent(in),		contiguous	::	knot3(:)
integer,	intent(in)			::	degree(3)
integer,	intent(in)			::	nc(3)
integer,	intent(in),	optional		::	ng(3)
real(kind=rk),	intent(in),		contiguous	::	Wc(:)
real(kind=rk),	intent(out),		allocatable	::	d2Tgc(:,:,:)
real(kind=rk),	intent(out),		allocatable	::	dTgc(:,:,:)
real(kind=rk),	intent(out),		allocatable	::	Tgc(:,:)

private pure subroutine compute_d2Tgc_bspline_3d_vector(Xt, knot1, knot2, knot3, degree, nc, ng, d2Tgc, dTgc, Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=rk),	intent(in),		contiguous	::	Xt(:,:)
real(kind=rk),	intent(in),		contiguous	::	knot1(:)
real(kind=rk),	intent(in),		contiguous	::	knot2(:)
real(kind=rk),	intent(in),		contiguous	::	knot3(:)
integer,	intent(in)			::	degree(3)
integer,	intent(in)			::	nc(3)
integer,	intent(in),	optional		::	ng(3)
real(kind=rk),	intent(out),		allocatable	::	d2Tgc(:,:,:)
real(kind=rk),	intent(out),		allocatable	::	dTgc(:,:,:)
real(kind=rk),	intent(out),		allocatable	::	Tgc(:,:)

private pure subroutine compute_d2Tgc_nurbs_3d_scalar(Xt, knot1, knot2, knot3, degree, nc, Wc, d2Tgc, dTgc, Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Attributes		Name
real(kind=rk),	intent(in),	contiguous	::	Xt(:)
real(kind=rk),	intent(in),	contiguous	::	knot1(:)
real(kind=rk),	intent(in),	contiguous	::	knot2(:)
real(kind=rk),	intent(in),	contiguous	::	knot3(:)
integer,	intent(in)		::	degree(3)
integer,	intent(in)		::	nc(3)
real(kind=rk),	intent(in),	contiguous	::	Wc(:)
real(kind=rk),	intent(out),	allocatable	::	d2Tgc(:,:)
real(kind=rk),	intent(out),	allocatable	::	dTgc(:,:)
real(kind=rk),	intent(out),	allocatable	::	Tgc(:)

private pure subroutine compute_d2Tgc_bspline_3d_scalar(Xt, knot1, knot2, knot3, degree, nc, d2Tgc, dTgc, Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Attributes		Name
real(kind=rk),	intent(in),	contiguous	::	Xt(:)
real(kind=rk),	intent(in),	contiguous	::	knot1(:)
real(kind=rk),	intent(in),	contiguous	::	knot2(:)
real(kind=rk),	intent(in),	contiguous	::	knot3(:)
integer,	intent(in)		::	degree(3)
integer,	intent(in)		::	nc(3)
real(kind=rk),	intent(out),	allocatable	::	d2Tgc(:,:)
real(kind=rk),	intent(out),	allocatable	::	dTgc(:,:)
real(kind=rk),	intent(out),	allocatable	::	Tgc(:)

Derived Types

type, public :: nurbs_volume

Components

Type	Visibility	Attributes		Name
type(debug),	public		::	err	101: size mismatch (weights vs control points), 102: missing control points, 103: missing knot vector, 104: missing geometry points, 105: missing weights, 106: lsq fit underdetermined
real(kind=rk),	private,	allocatable	::	Wc(:)	Weights for the control points (1D array: [nc(1)nc(2)nc(3)])
real(kind=rk),	private,	allocatable	::	Xc(:,:)	Control points (2D array: [nc(1)nc(2)nc(3), dim])
real(kind=rk),	private,	allocatable	::	Xg(:,:)	Geometry points (2D array: [ng(1)ng(2)ng(3), dim])
real(kind=rk),	private,	allocatable	::	Xt(:,:)	Evaluation parameter values (2D array: [ng(1)ng(2)ng(3), dim]
real(kind=rk),	private,	allocatable	::	Xt1(:)	Evaluation parameter values in the first direction (1D array: [ng(1)])
real(kind=rk),	private,	allocatable	::	Xt2(:)	Evaluation parameter values in the second direction (1D array: [ng(2)])
real(kind=rk),	private,	allocatable	::	Xt3(:)	Evaluation parameter values in the third direction (1D array: [ng(3)])
integer,	private		::	degree(3)	Degree (order) of the volume
integer,	private,	allocatable	::	elemConn(:,:)	IGA element connectivity
integer,	private,	allocatable	::	elemConn_Xc_vis(:,:)	Connectivity for visualization of control points
integer,	private,	allocatable	::	elemConn_Xg_vis(:,:)	Connectivity for visualization of geometry points
real(kind=rk),	private,	allocatable	::	knot1(:)	Knot vector in the first direction (1D array)
real(kind=rk),	private,	allocatable	::	knot2(:)	Knot vector in the second direction (1D array)
real(kind=rk),	private,	allocatable	::	knot3(:)	Knot vector in the third direction (1D array)
integer,	private		::	nc(3)	Number of control points in each direction
integer,	private		::	ng(3)	Number of geometry points in each direction

Type-Bound Procedures

procedure, public :: ansatz	Compute the shape functions, derivative of shape functions and dV Read more…
generic, public :: basis => basis_vector, basis_scalar	Compute the basis functions of the NURBS volume Read more…
procedure, public :: cmp_Xg	Compute geometry points Read more…
procedure, public :: cmp_degree	Compute degree of the NURBS volume Read more…
procedure, public :: cmp_degreeFace	Compute degrees of the faces Read more…
procedure, public :: cmp_elem	Generate IGA element connectivity Read more…
procedure, public :: cmp_elemFace	Compute faces of the IGA elements Read more…
procedure, public :: cmp_elemFace_Xc_vis	Compute faces of the control points Read more…
procedure, public :: cmp_elemFace_Xg_vis	Compute faces of the geometry points Read more…
procedure, public :: cmp_elem_Xc_vis	Generate connectivity for control points Read more…
procedure, public :: cmp_elem_Xg_vis	Generate connectivity for geometry points Read more…
procedure, public :: cmp_elem_Xth	Generate connectivity for parameter points Read more…
procedure, public :: cmp_nc	Compute number of required control points Read more…
procedure, public :: cmp_volume	Compute the volume of the NURBS volume Read more…
procedure, public :: create	Generate geometry points Read more…
generic, public :: derivative => derivative_vector, derivative_scalar	Compute the derivative of the NURBS volume Read more…
generic, public :: derivative2 => derivative2_vector, derivative2_scalar	Compute the second derivative of the NURBS volume Read more…
procedure, public :: elevate_degree	Elevate the degree of the NURBS volume Read more…
procedure, public :: export_Xc	Export control points to VTK file Read more…
procedure, public :: export_Xg	Export geometry points to VTK file Read more…
procedure, public :: export_Xth	Export parameter space to VTK file Read more…
procedure, public :: export_Xth_in_Xg	Export parameter space in geometry points to VTK file Read more…
procedure, public :: finalize	Finalize the NURBS volume object Read more…
generic, public :: get_Wc => get_Wc_all, get_Wci	Get weights Read more…
generic, public :: get_Xc => get_Xc_all, get_Xci, get_Xcid	Get control points Read more…
generic, public :: get_Xg => get_Xg_all, get_Xgi, get_Xgid	Get geometry points Read more…
procedure, public :: get_Xt	Get parameter values Read more…
procedure, public :: get_continuity	Compute and return the continuity of the NURBS volume Read more…
generic, public :: get_degree => get_degree_all, get_degree_dir	Get degree of the NURBS volume Read more…
procedure, public :: get_elem	Get IGA element connectivity Read more…
procedure, public :: get_elem_Xc_vis	Get connectivity for control points Read more…
procedure, public :: get_elem_Xg_vis	Get connectivity for geometry points Read more…
generic, public :: get_knot => get_knoti, get_knot_all	Get knot vector Read more…
procedure, public :: get_multiplicity	Compute and return the multiplicity of the knots Read more…
generic, public :: get_nc => get_nc_all, get_nc_dir	Get number of control points Read more…
procedure, public :: get_ng	Get number of geometry points Read more…
procedure, public :: insert_knots	Insert knots into the knot vector Read more…
procedure, public :: is_rational	Check if the NURBS volume is rational Read more…
procedure, public :: lsq_fit_bspline	Fit B-spline volume to structured data points using least squares Read more…
procedure, public :: modify_Wc	Modify weights Read more…
procedure, public :: modify_Xc	Modify control points Read more…
procedure, public :: nearest_point	Find the nearest point on the NURBS volume (Approximation) Read more…
procedure, public :: nearest_point2	Find the nearest point on the NURBS volume (Minimization - Newtons method) Read more…
procedure, public :: put_to_nurbs	Put a shape to a NURBS volume Read more…
procedure, public :: remove_knots	Remove knots from the knot vector Read more…
procedure, public :: rotate_Xc	Rotate control points Read more…
procedure, public :: rotate_Xg	Rotate geometry points Read more…
generic, public :: set => set1, set2, set3, set4	Set NURBS volume Read more…
procedure, public :: set_C	Set a C-shape Read more…
procedure, public :: set_elem	Set IGA element connectivity Read more…
procedure, public :: set_elem_Xc_vis	Set connectivity for control points Read more…
procedure, public :: set_elem_Xg_vis	Set connectivity for geometry points Read more…
procedure, public :: set_half_ring	Set a half ring Read more…
procedure, public :: set_hexahedron	Set a hexahedron Read more…
procedure, public :: set_ring	Set a ring Read more…
procedure, public :: show	Show the NURBS object using PyVista Read more…
procedure, public :: translate_Xc	Translate control points Read more…
procedure, public :: translate_Xg	Translate geometry points Read more…
procedure, private :: basis_scalar	Compute the basis functions of the NURBS volume Read more…
procedure, private :: basis_vector	Compute the basis functions of the NURBS volume Read more…
procedure, private :: derivative2_scalar	Compute the second derivative of the NURBS volume Read more…
procedure, private :: derivative2_vector	Compute the second derivative of the NURBS volume Read more…
procedure, private :: derivative_scalar	Compute the derivative of the NURBS volume Read more…
procedure, private :: derivative_vector	Compute the derivative of the NURBS volume Read more…
procedure, private :: get_Wc_all	Get all weights Read more…
procedure, private :: get_Wci	Get i-th weight Read more…
procedure, private :: get_Xc_all	Get all control points Read more…
procedure, private :: get_Xci	Get i-th control point Read more…
procedure, private :: get_Xcid	Get i-th control point in a specific direction Read more…
procedure, private :: get_Xg_all	Get all geometry points Read more…
procedure, private :: get_Xgi	Get i-th geometry point Read more…
procedure, private :: get_Xgid	Get i-th geometry point in a specific direction Read more…
procedure, private :: get_degree_all	Get degree of the NURBS volume in all directions Read more…
procedure, private :: get_degree_dir	Get degree of the NURBS volume in a specific direction Read more…
procedure, private :: get_knot_all	Get all knot vectors Read more…
procedure, private :: get_knoti	Get i-th knot value Read more…
procedure, private :: get_nc_all	Get number of control points in all directions Read more…
procedure, private :: get_nc_dir	Get number of control points in a specific direction Read more…
procedure, private :: set1	Set knot vectors, control points and weights for the NURBS volume object Read more…
procedure, private :: set2	Set NURBS volume using nodes of parameter space, degree, continuity, control points and weights Read more…
procedure, private :: set3	Set Bezier or Rational Bezier volume using control points and weights Read more…
procedure, private :: set4	Set NURBS volume using degree, number of control points, control points and weights Read more…

Functions

private pure function cmp_Tgc_3d(Xti, knot1, knot2, knot3, nc, degree, Wc) result(Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Attributes		Name
real(kind=rk),	intent(in),	contiguous	::	Xti(:)
real(kind=rk),	intent(in),	contiguous	::	knot1(:)
real(kind=rk),	intent(in),	contiguous	::	knot2(:)
real(kind=rk),	intent(in),	contiguous	::	knot3(:)
integer,	intent(in)		::	nc(3)
integer,	intent(in)		::	degree(3)
real(kind=rk),	intent(in),	contiguous	::	Wc(:)

Return Value real(kind=rk), (nc(1)nc(2)nc(3))

private pure function cmp_Xg(this, Xt) result(Xg)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this
real(kind=rk),	intent(in),		contiguous	::	Xt(:)

Return Value real(kind=rk), allocatable, (:)

private pure function cmp_degreeFace(this, face) result(degree)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this
integer,	intent(in)			::	face

Return Value integer, (3)

private pure function cmp_elem(this) result(elemConn)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this

Return Value integer, allocatable, (:,:)

private pure function cmp_elemFace(this, elem, face) result(elemConn)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

number of nodes in each direction

Arguments

Type	Intent		Name
class(nurbs_volume),	intent(in)	::	this
integer,	intent(in)	::	elem
integer,	intent(in)	::	face

Return Value integer, allocatable, (:)

private pure function cmp_elemFace_Xc_vis(this, elem, face) result(elemConn)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

number of nodes in each direction

Arguments

Type	Intent		Name
class(nurbs_volume),	intent(in)	::	this
integer,	intent(in)	::	elem
integer,	intent(in)	::	face

Return Value integer, allocatable, (:)

private pure function cmp_elemFace_Xg_vis(this, elem, face) result(elemConn)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

number of nodes in each direction

Arguments

Type	Intent		Name
class(nurbs_volume),	intent(in)	::	this
integer,	intent(in)	::	elem
integer,	intent(in)	::	face

Return Value integer, allocatable, (:)

private pure function cmp_elem_Xc_vis(this, p) result(elemConn)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this
integer,	intent(in),	optional,	contiguous	::	p(:)

Return Value integer, allocatable, (:,:)

private pure function cmp_elem_Xg_vis(this, p) result(elemConn)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this
integer,	intent(in),	optional,	contiguous	::	p(:)

Return Value integer, allocatable, (:,:)

private pure function cmp_elem_Xth(this, p) result(elemConn)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this
integer,	intent(in),	optional,	contiguous	::	p(:)

Return Value integer, allocatable, (:,:)

private pure function compute_Tgc_bspline_3d_scalar(Xt, knot1, knot2, knot3, degree, nc) result(Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Attributes		Name
real(kind=rk),	intent(in),	contiguous	::	Xt(:)
real(kind=rk),	intent(in),	contiguous	::	knot1(:)
real(kind=rk),	intent(in),	contiguous	::	knot2(:)
real(kind=rk),	intent(in),	contiguous	::	knot3(:)
integer,	intent(in)		::	degree(3)
integer,	intent(in)		::	nc(3)

Return Value real(kind=rk), allocatable, (:)

private pure function compute_Tgc_bspline_3d_vector(Xt, knot1, knot2, knot3, degree, nc, ng) result(Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=rk),	intent(in),		contiguous	::	Xt(:,:)
real(kind=rk),	intent(in),		contiguous	::	knot1(:)
real(kind=rk),	intent(in),		contiguous	::	knot2(:)
real(kind=rk),	intent(in),		contiguous	::	knot3(:)
integer,	intent(in)			::	degree(3)
integer,	intent(in)			::	nc(3)
integer,	intent(in),	optional		::	ng(3)

Return Value real(kind=rk), allocatable, (:,:)

private pure function compute_Tgc_nurbs_3d_scalar(Xt, knot1, knot2, knot3, degree, nc, Wc) result(Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Attributes		Name
real(kind=rk),	intent(in),	contiguous	::	Xt(:)
real(kind=rk),	intent(in),	contiguous	::	knot1(:)
real(kind=rk),	intent(in),	contiguous	::	knot2(:)
real(kind=rk),	intent(in),	contiguous	::	knot3(:)
integer,	intent(in)		::	degree(3)
integer,	intent(in)		::	nc(3)
real(kind=rk),	intent(in),	contiguous	::	Wc(:)

Return Value real(kind=rk), allocatable, (:)

private pure function compute_Tgc_nurbs_3d_vector(Xt, knot1, knot2, knot3, degree, nc, ng, Wc) result(Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=rk),	intent(in),		contiguous	::	Xt(:,:)
real(kind=rk),	intent(in),		contiguous	::	knot1(:)
real(kind=rk),	intent(in),		contiguous	::	knot2(:)
real(kind=rk),	intent(in),		contiguous	::	knot3(:)
integer,	intent(in)			::	degree(3)
integer,	intent(in)			::	nc(3)
integer,	intent(in),	optional		::	ng(3)
real(kind=rk),	intent(in),		contiguous	::	Wc(:)

Return Value real(kind=rk), allocatable, (:,:)

private pure function compute_Xg_bspline_3d(Xt, knot1, knot2, knot3, degree, nc, ng, Xc) result(Xg)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=rk),	intent(in),		contiguous	::	Xt(:,:)
real(kind=rk),	intent(in),		contiguous	::	knot1(:)
real(kind=rk),	intent(in),		contiguous	::	knot2(:)
real(kind=rk),	intent(in),		contiguous	::	knot3(:)
integer,	intent(in)			::	degree(3)
integer,	intent(in)			::	nc(3)
integer,	intent(in),	optional		::	ng(3)
real(kind=rk),	intent(in),		contiguous	::	Xc(:,:)

Return Value real(kind=rk), allocatable, (:,:)

private pure function compute_Xg_bspline_3d_1point(Xt, knot1, knot2, knot3, degree, nc, Xc) result(Xg)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Attributes		Name
real(kind=rk),	intent(in),	contiguous	::	Xt(:)
real(kind=rk),	intent(in),	contiguous	::	knot1(:)
real(kind=rk),	intent(in),	contiguous	::	knot2(:)
real(kind=rk),	intent(in),	contiguous	::	knot3(:)
integer,	intent(in)		::	degree(3)
integer,	intent(in)		::	nc(3)
real(kind=rk),	intent(in),	contiguous	::	Xc(:,:)

Return Value real(kind=rk), (size(Xc,2))

private pure function compute_Xg_nurbs_3d(Xt, knot1, knot2, knot3, degree, nc, ng, Xc, Wc) result(Xg)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=rk),	intent(in),		contiguous	::	Xt(:,:)
real(kind=rk),	intent(in),		contiguous	::	knot1(:)
real(kind=rk),	intent(in),		contiguous	::	knot2(:)
real(kind=rk),	intent(in),		contiguous	::	knot3(:)
integer,	intent(in)			::	degree(3)
integer,	intent(in)			::	nc(3)
integer,	intent(in),	optional		::	ng(3)
real(kind=rk),	intent(in),		contiguous	::	Xc(:,:)
real(kind=rk),	intent(in),		contiguous	::	Wc(:)

Return Value real(kind=rk), allocatable, (:,:)

private pure function compute_Xg_nurbs_3d_1point(Xt, knot1, knot2, knot3, degree, nc, Xc, Wc) result(Xg)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Attributes		Name
real(kind=rk),	intent(in),	contiguous	::	Xt(:)
real(kind=rk),	intent(in),	contiguous	::	knot1(:)
real(kind=rk),	intent(in),	contiguous	::	knot2(:)
real(kind=rk),	intent(in),	contiguous	::	knot3(:)
integer,	intent(in)		::	degree(3)
integer,	intent(in)		::	nc(3)
real(kind=rk),	intent(in),	contiguous	::	Xc(:,:)
real(kind=rk),	intent(in),	contiguous	::	Wc(:)

Return Value real(kind=rk), (size(Xc,2))

private pure function get_Wc_all(this) result(Wc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this

Return Value real(kind=rk), allocatable, (:)

private pure function get_Wci(this, n) result(Wc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this
integer,	intent(in)			::	n

Return Value real(kind=rk)

private pure function get_Xc_all(this) result(Xc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this

Return Value real(kind=rk), allocatable, (:,:)

private pure function get_Xci(this, n) result(Xc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this
integer,	intent(in)			::	n

Return Value real(kind=rk), allocatable, (:)

private pure function get_Xcid(this, n, dir) result(Xc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent		Name
class(nurbs_volume),	intent(in)	::	this
integer,	intent(in)	::	n
integer,	intent(in)	::	dir

Return Value real(kind=rk)

private pure function get_Xg_all(this) result(Xg)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this

Return Value real(kind=rk), allocatable, (:,:)

private pure function get_Xgi(this, n) result(Xg)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this
integer,	intent(in)			::	n

Return Value real(kind=rk), allocatable, (:)

private pure function get_Xgid(this, n, dir) result(Xg)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent		Name
class(nurbs_volume),	intent(in)	::	this
integer,	intent(in)	::	n
integer,	intent(in)	::	dir

Return Value real(kind=rk)

private pure function get_Xt(this, dir) result(Xt)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this
integer,	intent(in)			::	dir

Return Value real(kind=rk), allocatable, (:)

private pure function get_continuity(this, dir) result(c)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this
integer,	intent(in)			::	dir

Return Value integer, allocatable, (:)

private pure function get_degree_all(this) result(degree)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this

Return Value integer, (3)

private pure function get_degree_dir(this, dir) result(degree)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this
integer,	intent(in)			::	dir

Return Value integer

private pure function get_elem(this) result(elemConn)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this

Return Value integer, allocatable, (:,:)

private pure function get_elem_Xc_vis(this) result(elemConn)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this

Return Value integer, allocatable, (:,:)

private pure function get_elem_Xg_vis(this) result(elemConn)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this

Return Value integer, allocatable, (:,:)

private pure function get_knot_all(this, dir) result(knot)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this
integer,	intent(in)			::	dir

Return Value real(kind=rk), allocatable, (:)

private pure function get_knoti(this, dir, i) result(knot)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent		Name
class(nurbs_volume),	intent(in)	::	this
integer,	intent(in)	::	dir
integer,	intent(in)	::	i

Return Value real(kind=rk)

private pure function get_multiplicity(this, dir) result(m)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this
integer,	intent(in)			::	dir

Return Value integer, allocatable, (:)

private pure function get_nc_all(this) result(nc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this

Return Value integer, (3)

private pure function get_nc_dir(this, dir) result(nc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this
integer,	intent(in)			::	dir

Return Value integer

private pure function get_ng(this) result(ng)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this

Return Value integer, (3)

private pure function is_rational(this) result(r)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this

Return Value logical

Subroutines

private pure subroutine ansatz(this, ie, ig, Tgc, dTgc_dXg, dV, ngauss)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(inout)			::	this
integer,	intent(in)			::	ie
integer,	intent(in)			::	ig
real(kind=rk),	intent(out),		allocatable	::	Tgc(:)
real(kind=rk),	intent(out),		allocatable	::	dTgc_dXg(:,:)
real(kind=rk),	intent(out)			::	dV
integer,	intent(in),	optional		::	ngauss(3)

private pure subroutine basis_scalar(this, Xt, Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Attributes		Name
class(nurbs_volume),	intent(inout)		::	this
real(kind=rk),	intent(in),	contiguous	::	Xt(:)
real(kind=rk),	intent(out),	allocatable	::	Tgc(:)

private pure subroutine basis_vector(this, res1, res2, res3, Xt1, Xt2, Xt3, Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(inout)			::	this
integer,	intent(in),	optional		::	res1
integer,	intent(in),	optional		::	res2
integer,	intent(in),	optional		::	res3
real(kind=rk),	intent(in),	optional,	contiguous	::	Xt1(:)
real(kind=rk),	intent(in),	optional,	contiguous	::	Xt2(:)
real(kind=rk),	intent(in),	optional,	contiguous	::	Xt3(:)
real(kind=rk),	intent(out),		allocatable	::	Tgc(:,:)

private pure subroutine cmp_degree(this, dir)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(inout)			::	this
integer,	intent(in),	optional		::	dir

private pure subroutine cmp_nc(this, dir)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(inout)			::	this
integer,	intent(in),	optional		::	dir

private pure subroutine cmp_volume(this, volume, ngauss)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional		Name
class(nurbs_volume),	intent(inout)		::	this
real(kind=rk),	intent(out)		::	volume
integer,	intent(in),	optional	::	ngauss(3)

private pure subroutine compute_d2Tgc_bspline_3d_scalar(Xt, knot1, knot2, knot3, degree, nc, d2Tgc, dTgc, Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Attributes		Name
real(kind=rk),	intent(in),	contiguous	::	Xt(:)
real(kind=rk),	intent(in),	contiguous	::	knot1(:)
real(kind=rk),	intent(in),	contiguous	::	knot2(:)
real(kind=rk),	intent(in),	contiguous	::	knot3(:)
integer,	intent(in)		::	degree(3)
integer,	intent(in)		::	nc(3)
real(kind=rk),	intent(out),	allocatable	::	d2Tgc(:,:)
real(kind=rk),	intent(out),	allocatable	::	dTgc(:,:)
real(kind=rk),	intent(out),	allocatable	::	Tgc(:)

private pure subroutine compute_d2Tgc_bspline_3d_vector(Xt, knot1, knot2, knot3, degree, nc, ng, d2Tgc, dTgc, Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=rk),	intent(in),		contiguous	::	Xt(:,:)
real(kind=rk),	intent(in),		contiguous	::	knot1(:)
real(kind=rk),	intent(in),		contiguous	::	knot2(:)
real(kind=rk),	intent(in),		contiguous	::	knot3(:)
integer,	intent(in)			::	degree(3)
integer,	intent(in)			::	nc(3)
integer,	intent(in),	optional		::	ng(3)
real(kind=rk),	intent(out),		allocatable	::	d2Tgc(:,:,:)
real(kind=rk),	intent(out),		allocatable	::	dTgc(:,:,:)
real(kind=rk),	intent(out),		allocatable	::	Tgc(:,:)

private pure subroutine compute_d2Tgc_nurbs_3d_scalar(Xt, knot1, knot2, knot3, degree, nc, Wc, d2Tgc, dTgc, Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Attributes		Name
real(kind=rk),	intent(in),	contiguous	::	Xt(:)
real(kind=rk),	intent(in),	contiguous	::	knot1(:)
real(kind=rk),	intent(in),	contiguous	::	knot2(:)
real(kind=rk),	intent(in),	contiguous	::	knot3(:)
integer,	intent(in)		::	degree(3)
integer,	intent(in)		::	nc(3)
real(kind=rk),	intent(in),	contiguous	::	Wc(:)
real(kind=rk),	intent(out),	allocatable	::	d2Tgc(:,:)
real(kind=rk),	intent(out),	allocatable	::	dTgc(:,:)
real(kind=rk),	intent(out),	allocatable	::	Tgc(:)

private pure subroutine compute_d2Tgc_nurbs_3d_vector(Xt, knot1, knot2, knot3, degree, nc, ng, Wc, d2Tgc, dTgc, Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=rk),	intent(in),		contiguous	::	Xt(:,:)
real(kind=rk),	intent(in),		contiguous	::	knot1(:)
real(kind=rk),	intent(in),		contiguous	::	knot2(:)
real(kind=rk),	intent(in),		contiguous	::	knot3(:)
integer,	intent(in)			::	degree(3)
integer,	intent(in)			::	nc(3)
integer,	intent(in),	optional		::	ng(3)
real(kind=rk),	intent(in),		contiguous	::	Wc(:)
real(kind=rk),	intent(out),		allocatable	::	d2Tgc(:,:,:)
real(kind=rk),	intent(out),		allocatable	::	dTgc(:,:,:)
real(kind=rk),	intent(out),		allocatable	::	Tgc(:,:)

private pure subroutine compute_dTgc_bspline_3d_scalar(Xt, knot1, knot2, knot3, degree, nc, dTgc, Tgc, elem)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=rk),	intent(in),		contiguous	::	Xt(:)
real(kind=rk),	intent(in),		contiguous	::	knot1(:)
real(kind=rk),	intent(in),		contiguous	::	knot2(:)
real(kind=rk),	intent(in),		contiguous	::	knot3(:)
integer,	intent(in)			::	degree(3)
integer,	intent(in)			::	nc(3)
real(kind=rk),	intent(out),		allocatable	::	dTgc(:,:)
real(kind=rk),	intent(out),		allocatable	::	Tgc(:)
integer,	intent(in),	optional,	contiguous	::	elem(:)

private pure subroutine compute_dTgc_bspline_3d_vector(Xt, knot1, knot2, knot3, degree, nc, ng, dTgc, Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=rk),	intent(in),		contiguous	::	Xt(:,:)
real(kind=rk),	intent(in),		contiguous	::	knot1(:)
real(kind=rk),	intent(in),		contiguous	::	knot2(:)
real(kind=rk),	intent(in),		contiguous	::	knot3(:)
integer,	intent(in)			::	degree(3)
integer,	intent(in)			::	nc(3)
integer,	intent(in),	optional		::	ng(3)
real(kind=rk),	intent(out),		allocatable	::	dTgc(:,:,:)
real(kind=rk),	intent(out),		allocatable	::	Tgc(:,:)

private pure subroutine compute_dTgc_nurbs_3d_scalar(Xt, knot1, knot2, knot3, degree, nc, Wc, dTgc, Tgc, elem)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

If elem is not present: Wc refers to the full weight vector. If elem is present: Wc refers to the element-local weight vector (Wce).

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=rk),	intent(in),		contiguous	::	Xt(:)
real(kind=rk),	intent(in),		contiguous	::	knot1(:)
real(kind=rk),	intent(in),		contiguous	::	knot2(:)
real(kind=rk),	intent(in),		contiguous	::	knot3(:)
integer,	intent(in)			::	degree(3)
integer,	intent(in)			::	nc(3)
real(kind=rk),	intent(in),		contiguous	::	Wc(:)
real(kind=rk),	intent(out),		allocatable	::	dTgc(:,:)
real(kind=rk),	intent(out),		allocatable	::	Tgc(:)
integer,	intent(in),	optional,	contiguous	::	elem(:)

private pure subroutine compute_dTgc_nurbs_3d_vector(Xt, knot1, knot2, knot3, degree, nc, ng, Wc, dTgc, Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=rk),	intent(in),		contiguous	::	Xt(:,:)
real(kind=rk),	intent(in),		contiguous	::	knot1(:)
real(kind=rk),	intent(in),		contiguous	::	knot2(:)
real(kind=rk),	intent(in),		contiguous	::	knot3(:)
integer,	intent(in)			::	degree(3)
integer,	intent(in)			::	nc(3)
integer,	intent(in),	optional		::	ng(3)
real(kind=rk),	intent(in),		contiguous	::	Wc(:)
real(kind=rk),	intent(out),		allocatable	::	dTgc(:,:,:)
real(kind=rk),	intent(out),		allocatable	::	Tgc(:,:)

private pure subroutine create(this, res1, res2, res3, Xt1, Xt2, Xt3, Xt)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(inout)			::	this
integer,	intent(in),	optional		::	res1
integer,	intent(in),	optional		::	res2
integer,	intent(in),	optional		::	res3
real(kind=rk),	intent(in),	optional,	contiguous	::	Xt1(:)
real(kind=rk),	intent(in),	optional,	contiguous	::	Xt2(:)
real(kind=rk),	intent(in),	optional,	contiguous	::	Xt3(:)
real(kind=rk),	intent(in),	optional,	contiguous	::	Xt(:,:)

private pure subroutine derivative2_scalar(this, Xt, d2Tgc, dTgc, Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(inout)			::	this
real(kind=rk),	intent(in),		contiguous	::	Xt(:)
real(kind=rk),	intent(out),		allocatable	::	d2Tgc(:,:)
real(kind=rk),	intent(out),	optional,	allocatable	::	dTgc(:,:)
real(kind=rk),	intent(out),	optional,	allocatable	::	Tgc(:)

private pure subroutine derivative2_vector(this, res1, res2, res3, Xt1, Xt2, Xt3, d2Tgc, dTgc, Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(inout)			::	this
integer,	intent(in),	optional		::	res1
integer,	intent(in),	optional		::	res2
integer,	intent(in),	optional		::	res3
real(kind=rk),	intent(in),	optional,	contiguous	::	Xt1(:)
real(kind=rk),	intent(in),	optional,	contiguous	::	Xt2(:)
real(kind=rk),	intent(in),	optional,	contiguous	::	Xt3(:)
real(kind=rk),	intent(out),		allocatable	::	d2Tgc(:,:,:)
real(kind=rk),	intent(out),	optional,	allocatable	::	dTgc(:,:,:)
real(kind=rk),	intent(out),	optional,	allocatable	::	Tgc(:,:)

private pure subroutine derivative_scalar(this, Xt, dTgc, Tgc, elem)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(inout)			::	this
real(kind=rk),	intent(in),		contiguous	::	Xt(:)
real(kind=rk),	intent(out),		allocatable	::	dTgc(:,:)
real(kind=rk),	intent(out),	optional,	allocatable	::	Tgc(:)
integer,	intent(in),	optional,	contiguous	::	elem(:)

private pure subroutine derivative_vector(this, res1, res2, res3, Xt1, Xt2, Xt3, dTgc, Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(inout)			::	this
integer,	intent(in),	optional		::	res1
integer,	intent(in),	optional		::	res2
integer,	intent(in),	optional		::	res3
real(kind=rk),	intent(in),	optional,	contiguous	::	Xt1(:)
real(kind=rk),	intent(in),	optional,	contiguous	::	Xt2(:)
real(kind=rk),	intent(in),	optional,	contiguous	::	Xt3(:)
real(kind=rk),	intent(out),		allocatable	::	dTgc(:,:,:)
real(kind=rk),	intent(out),	optional,	allocatable	::	Tgc(:,:)

private pure subroutine elevate_degree(this, dir, t, B, Bs)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(inout)			::	this
integer,	intent(in)			::	dir
integer,	intent(in)			::	t
real(kind=rk),	intent(out),	optional,	allocatable	::	B(:,:)
real(kind=rk),	intent(out),	optional,	allocatable	::	Bs(:,:)

private impure subroutine export_Xc(this, filename, point_data, field_names, encoding)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(inout)			::	this
character(len=*),	intent(in)			::	filename
real(kind=rk),	intent(in),	optional,	contiguous	::	point_data(:,:)
character(len=*),	intent(in),	optional,	contiguous	::	field_names(:)
character(len=*),	intent(in),	optional		::	encoding

private impure subroutine export_Xg(this, filename, point_data, field_names, encoding)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(inout)			::	this
character(len=*),	intent(in)			::	filename
real(kind=rk),	intent(in),	optional,	contiguous	::	point_data(:,:)
character(len=*),	intent(in),	optional,	contiguous	::	field_names(:)
character(len=*),	intent(in),	optional		::	encoding

private impure subroutine export_Xth(this, filename, point_data, field_names, encoding)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this
character(len=*),	intent(in)			::	filename
real(kind=rk),	intent(in),	optional,	contiguous	::	point_data(:,:)
character(len=*),	intent(in),	optional,	contiguous	::	field_names(:)
character(len=*),	intent(in),	optional		::	encoding

private impure subroutine export_Xth_in_Xg(this, filename, res, encoding)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional		Name
class(nurbs_volume),	intent(inout)		::	this
character(len=*),	intent(in)		::	filename
integer,	intent(in),	optional	::	res
character(len=*),	intent(in),	optional	::	encoding

private pure subroutine finalize(this)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(inout)			::	this

private pure subroutine insert_knots(this, dir, Xth, r, B, Bs)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(inout)			::	this
integer,	intent(in)			::	dir
real(kind=rk),	intent(in),		contiguous	::	Xth(:)
integer,	intent(in),		contiguous	::	r(:)
real(kind=rk),	intent(out),	optional,	allocatable	::	B(:,:)
real(kind=rk),	intent(out),	optional,	allocatable	::	Bs(:,:)

private pure subroutine lsq_fit_bspline(this, Xt, Xdata, ndata)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Attributes		Name
class(nurbs_volume),	intent(inout)		::	this
real(kind=rk),	intent(in),	contiguous	::	Xt(:,:)
real(kind=rk),	intent(in),	contiguous	::	Xdata(:,:)
integer,	intent(in)		::	ndata(3)

private pure subroutine modify_Wc(this, W, num)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent		Name
class(nurbs_volume),	intent(inout)	::	this
real(kind=rk),	intent(in)	::	W
integer,	intent(in)	::	num

private pure subroutine modify_Xc(this, X, num, dir)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent		Name
class(nurbs_volume),	intent(inout)	::	this
real(kind=rk),	intent(in)	::	X
integer,	intent(in)	::	num
integer,	intent(in)	::	dir

private pure subroutine nearest_point(this, point_Xg, nearest_Xg, nearest_Xt, id)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(in)			::	this
real(kind=rk),	intent(in),		contiguous	::	point_Xg(:)
real(kind=rk),	intent(out),	optional		::	nearest_Xg(size(point_Xg))
real(kind=rk),	intent(out),	optional		::	nearest_Xt(3)
integer,	intent(out),	optional		::	id

private impure subroutine nearest_point2(this, point_Xg, tol, maxit, nearest_Xt, nearest_Xg)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(inout)			::	this
real(kind=rk),	intent(in),		contiguous	::	point_Xg(:)
real(kind=rk),	intent(in)			::	tol
integer,	intent(in)			::	maxit
real(kind=rk),	intent(out)			::	nearest_Xt(3)
real(kind=rk),	intent(out),	optional		::	nearest_Xg(size(this%Xc,2))

private pure subroutine put_to_nurbs(this, X, elemConn)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(inout)			::	this
real(kind=rk),	intent(in),		contiguous	::	X(:,:)
integer,	intent(in),	optional,	contiguous	::	elemConn(:,:)

private pure subroutine remove_knots(this, dir, Xth, r)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Attributes		Name
class(nurbs_volume),	intent(inout)		::	this
integer,	intent(in)		::	dir
real(kind=rk),	intent(in),	contiguous	::	Xth(:)
integer,	intent(in),	contiguous	::	r(:)

private pure subroutine rotate_Xc(this, alpha, beta, theta)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent		Name
class(nurbs_volume),	intent(inout)	::	this
real(kind=rk),	intent(in)	::	alpha
real(kind=rk),	intent(in)	::	beta
real(kind=rk),	intent(in)	::	theta

private pure subroutine rotate_Xg(this, alpha, beta, theta)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent		Name
class(nurbs_volume),	intent(inout)	::	this
real(kind=rk),	intent(in)	::	alpha
real(kind=rk),	intent(in)	::	beta
real(kind=rk),	intent(in)	::	theta

private pure subroutine set1(this, knot1, knot2, knot3, Xc, Wc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Set control points and weights for the NURBS volume object.

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(inout)			::	this
real(kind=rk),	intent(in),		contiguous	::	knot1(:)
real(kind=rk),	intent(in),		contiguous	::	knot2(:)
real(kind=rk),	intent(in),		contiguous	::	knot3(:)
real(kind=rk),	intent(in),		contiguous	::	Xc(:,:)
real(kind=rk),	intent(in),	optional,	contiguous	::	Wc(:)

private pure subroutine set2(this, Xth_dir1, Xth_dir2, Xth_dir3, degree, continuity1, continuity2, continuity3, Xc, Wc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Set control points and weights for the NURBS volume object.

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(inout)			::	this
real(kind=rk),	intent(in),		contiguous	::	Xth_dir1(:)
real(kind=rk),	intent(in),		contiguous	::	Xth_dir2(:)
real(kind=rk),	intent(in),		contiguous	::	Xth_dir3(:)
integer,	intent(in),		contiguous	::	degree(:)
integer,	intent(in),		contiguous	::	continuity1(:)
integer,	intent(in),		contiguous	::	continuity2(:)
integer,	intent(in),		contiguous	::	continuity3(:)
real(kind=rk),	intent(in),	optional,	contiguous	::	Xc(:,:)
real(kind=rk),	intent(in),	optional,	contiguous	::	Wc(:)

private pure subroutine set3(this, nc, Xc, Wc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Set Bezier or Rational Bezier volume using control points and weights.

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(inout)			::	this
integer,	intent(in),		contiguous	::	nc(:)
real(kind=rk),	intent(in),		contiguous	::	Xc(:,:)
real(kind=rk),	intent(in),	optional,	contiguous	::	Wc(:)

private pure subroutine set4(this, degree, nc, Xc, Wc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(inout)			::	this
integer,	intent(in),		contiguous	::	degree(:)
integer,	intent(in),		contiguous	::	nc(:)
real(kind=rk),	intent(in),		contiguous	::	Xc(:,:)
real(kind=rk),	intent(in),	optional,	contiguous	::	Wc(:)

private pure subroutine set_C(this, center, radius1, radius2, length)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Attributes		Name
class(nurbs_volume),	intent(inout)		::	this
real(kind=rk),	intent(in),	contiguous	::	center(:)
real(kind=rk),	intent(in)		::	radius1
real(kind=rk),	intent(in)		::	radius2
real(kind=rk),	intent(in)		::	length

private pure subroutine set_elem(this, elemConn)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(inout)			::	this
integer,	intent(in),		contiguous	::	elemConn(:,:)

private pure subroutine set_elem_Xc_vis(this, elemConn)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(inout)			::	this
integer,	intent(in),		contiguous	::	elemConn(:,:)

private pure subroutine set_elem_Xg_vis(this, elemConn)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(inout)			::	this
integer,	intent(in),		contiguous	::	elemConn(:,:)

private pure subroutine set_half_ring(this, center, radius1, radius2, length)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Attributes		Name
class(nurbs_volume),	intent(inout)		::	this
real(kind=rk),	intent(in),	contiguous	::	center(:)
real(kind=rk),	intent(in)		::	radius1
real(kind=rk),	intent(in)		::	radius2
real(kind=rk),	intent(in)		::	length

private pure subroutine set_hexahedron(this, L, nc, Wc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(inout)			::	this
real(kind=rk),	intent(in),		contiguous	::	L(:)
integer,	intent(in),		contiguous	::	nc(:)
real(kind=rk),	intent(in),	optional,	contiguous	::	Wc(:)

private pure subroutine set_ring(this, center, radius1, radius2, length)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Attributes		Name
class(nurbs_volume),	intent(inout)		::	this
real(kind=rk),	intent(in),	contiguous	::	center(:)
real(kind=rk),	intent(in)		::	radius1
real(kind=rk),	intent(in)		::	radius2
real(kind=rk),	intent(in)		::	length

private impure subroutine show(this, vtkfile_Xc, vtkfile_Xg, vtkfile_Xth_in_Xg)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional		Name
class(nurbs_volume),	intent(inout)		::	this
character(len=*),	intent(in)		::	vtkfile_Xc
character(len=*),	intent(in)		::	vtkfile_Xg
character(len=*),	intent(in),	optional	::	vtkfile_Xth_in_Xg

private pure subroutine translate_Xc(this, vec)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(inout)			::	this
real(kind=rk),	intent(in),		contiguous	::	vec(:)

private pure subroutine translate_Xg(this, vec)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_volume),	intent(inout)			::	this
real(kind=rk),	intent(in),		contiguous	::	vec(:)

forcad_nurbs_volume Module

Contents

Interfaces

Derived Types

Functions

Subroutines

Uses

Used by

Interfaces

public interface compute_Tgc

private pure function compute_Tgc_nurbs_3d_vector(Xt, knot1, knot2, knot3, degree, nc, ng, Wc) result(Tgc)

Arguments

Return Value real(kind=rk), allocatable, (:,:)

private pure function compute_Tgc_bspline_3d_vector(Xt, knot1, knot2, knot3, degree, nc, ng) result(Tgc)

Arguments

Return Value real(kind=rk), allocatable, (:,:)

private pure function compute_Tgc_nurbs_3d_scalar(Xt, knot1, knot2, knot3, degree, nc, Wc) result(Tgc)

Arguments

Return Value real(kind=rk), allocatable, (:)

private pure function compute_Tgc_bspline_3d_scalar(Xt, knot1, knot2, knot3, degree, nc) result(Tgc)

Arguments

Return Value real(kind=rk), allocatable, (:)

public interface compute_dTgc

private pure subroutine compute_dTgc_nurbs_3d_vector(Xt, knot1, knot2, knot3, degree, nc, ng, Wc, dTgc, Tgc)

Arguments

private pure subroutine compute_dTgc_bspline_3d_vector(Xt, knot1, knot2, knot3, degree, nc, ng, dTgc, Tgc)

Arguments

private pure subroutine compute_dTgc_nurbs_3d_scalar(Xt, knot1, knot2, knot3, degree, nc, Wc, dTgc, Tgc, elem)

Arguments

private pure subroutine compute_dTgc_bspline_3d_scalar(Xt, knot1, knot2, knot3, degree, nc, dTgc, Tgc, elem)

Arguments

private interface compute_Xg

private pure function compute_Xg_nurbs_3d(Xt, knot1, knot2, knot3, degree, nc, ng, Xc, Wc) result(Xg)

Arguments

Return Value real(kind=rk), allocatable, (:,:)

private pure function compute_Xg_bspline_3d(Xt, knot1, knot2, knot3, degree, nc, ng, Xc) result(Xg)

Arguments

Return Value real(kind=rk), allocatable, (:,:)

private pure function compute_Xg_nurbs_3d_1point(Xt, knot1, knot2, knot3, degree, nc, Xc, Wc) result(Xg)

Arguments

Return Value real(kind=rk), (size(Xc,2))

private pure function compute_Xg_bspline_3d_1point(Xt, knot1, knot2, knot3, degree, nc, Xc) result(Xg)

Arguments

Return Value real(kind=rk), (size(Xc,2))

private interface compute_d2Tgc

private pure subroutine compute_d2Tgc_nurbs_3d_vector(Xt, knot1, knot2, knot3, degree, nc, ng, Wc, d2Tgc, dTgc, Tgc)

Arguments

private pure subroutine compute_d2Tgc_bspline_3d_vector(Xt, knot1, knot2, knot3, degree, nc, ng, d2Tgc, dTgc, Tgc)

Arguments

private pure subroutine compute_d2Tgc_nurbs_3d_scalar(Xt, knot1, knot2, knot3, degree, nc, Wc, d2Tgc, dTgc, Tgc)

Arguments

private pure subroutine compute_d2Tgc_bspline_3d_scalar(Xt, knot1, knot2, knot3, degree, nc, d2Tgc, dTgc, Tgc)

Arguments

Derived Types

type, public :: nurbs_volume

Components

Type-Bound Procedures

Functions

private pure function cmp_Tgc_3d(Xti, knot1, knot2, knot3, nc, degree, Wc) result(Tgc)

Arguments

Return Value real(kind=rk), (nc(1)*nc(2)*nc(3))

private pure function cmp_Xg(this, Xt) result(Xg)

Arguments

Return Value real(kind=rk), allocatable, (:)

private pure function cmp_degreeFace(this, face) result(degree)

Arguments

Return Value integer, (3)

private pure function cmp_elem(this) result(elemConn)

Arguments

Return Value integer, allocatable, (:,:)

private pure function cmp_elemFace(this, elem, face) result(elemConn)

Arguments

Return Value integer, allocatable, (:)

private pure function cmp_elemFace_Xc_vis(this, elem, face) result(elemConn)

Arguments

Return Value integer, allocatable, (:)

private pure function cmp_elemFace_Xg_vis(this, elem, face) result(elemConn)

Arguments

Return Value integer, allocatable, (:)

private pure function cmp_elem_Xc_vis(this, p) result(elemConn)

Return Value real(kind=rk), (nc(1)nc(2)nc(3))