nurbs_curve – ForCAD

type, public :: nurbs_curve

Inherits

Help

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 control points (1D array: [nc])
real(kind=rk),	private,	allocatable	::	Xc(:,:)	Control points (2D array: [nc, dim])
real(kind=rk),	private,	allocatable	::	Xg(:,:)	Geometry points (2D array: [ng, dim])
real(kind=rk),	private,	allocatable	::	Xt(:)	Evaluation points (1D array: [ng])
integer,	private		::	degree	Degree (order) of the curve
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	::	knot(:)	Knot vector (1D array)
integer,	private		::	nc	Number of control points
integer,	private		::	ng	Number of geometry points

Type-Bound Procedures

procedure, public :: ansatz

Compute the shape functions, derivative of shape functions and dL

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

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_curve),	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)			::	dL
integer,	intent(in),	optional		::	ngauss

generic, public :: basis => basis_vector, basis_scalar

Compute the basis functions of the NURBS curve

private pure subroutine basis_vector(this, res, Xt, Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

private pure subroutine basis_scalar(this, Xt, Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

procedure, public :: cmp_Xg

Compute geometry points

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_curve),	intent(in)			::	this
real(kind=rk),	intent(in)			::	Xt

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

procedure, public :: cmp_degree

Compute degree of the NURBS curve

private pure subroutine cmp_degree(this)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

procedure, public :: cmp_elem

Generate IGA element connectivity

private pure function cmp_elem(this) result(elemConn)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

Return Value integer, allocatable, (:,:)

procedure, public :: cmp_elem_Xc_vis

Generate connectivity for control points

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_curve),	intent(in)			::	this
integer,	intent(in),	optional		::	p

Return Value integer, allocatable, (:,:)

procedure, public :: cmp_elem_Xg_vis

Generate connectivity for geometry points

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_curve),	intent(in)			::	this
integer,	intent(in),	optional		::	p

Return Value integer, allocatable, (:,:)

procedure, public :: cmp_elem_Xth

Generate connectivity for parameter points

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_curve),	intent(in)			::	this
integer,	intent(in),	optional		::	p

Return Value integer, allocatable, (:,:)

procedure, public :: cmp_length

Compute the length of the NURBS curve

private pure subroutine cmp_length(this, length, ngauss)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

procedure, public :: cmp_nc

Compute number of required control points

private pure subroutine cmp_nc(this)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

procedure, public :: create

Generate geometry points

private pure subroutine create(this, res, Xt)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_curve),	intent(inout)			::	this
integer,	intent(in),	optional		::	res
real(kind=rk),	intent(in),	optional,	contiguous	::	Xt(:)

generic, public :: derivative => derivative_vector, derivative_scalar

Compute the derivative of the NURBS curve

private pure subroutine derivative_vector(this, res, Xt, dTgc, Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_curve),	intent(inout)			::	this
integer,	intent(in),	optional		::	res
real(kind=rk),	intent(in),	optional,	contiguous	::	Xt(:)
real(kind=rk),	intent(out),		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_curve),	intent(inout)			::	this
real(kind=rk),	intent(in)			::	Xt
real(kind=rk),	intent(out),		allocatable	::	dTgc(:)
real(kind=rk),	intent(out),	optional,	allocatable	::	Tgc(:)
integer,	intent(in),	optional,	contiguous	::	elem(:)

generic, public :: derivative2 => derivative2_vector, derivative2_scalar

Compute the second derivative of the NURBS curve

private pure subroutine derivative2_vector(this, res, Xt, d2Tgc, dTgc, Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_curve),	intent(inout)			::	this
integer,	intent(in),	optional		::	res
real(kind=rk),	intent(in),	optional,	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_scalar(this, Xt, d2Tgc, dTgc, Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

procedure, public :: elevate_degree

Elevate the degree of the curve

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

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

procedure, public :: export_Xc

Export control points to VTK file

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_curve),	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

procedure, public :: export_Xg

Export geometry points to VTK file

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_curve),	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

procedure, public :: export_Xth

Export parameter space to VTK file

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_curve),	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

procedure, public :: export_iges

Export the NURBS curve to an IGES file

private impure subroutine export_iges(this, filename)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_curve),	intent(inout)			::	this
character(len=*),	intent(in)			::	filename

procedure, public :: finalize

Finalize the NURBS curve object

private pure elemental subroutine finalize(this)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

generic, public :: get_Wc => get_Wc_all, get_Wci

Get weights

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_curve),	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_curve),	intent(in)			::	this
integer,	intent(in)			::	n

Return Value real(kind=rk)

generic, public :: get_Xc => get_Xc_all, get_Xci, get_Xcid

Get control points

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_curve),	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_curve),	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_curve),	intent(in)	::	this
integer,	intent(in)	::	n
integer,	intent(in)	::	dir

Return Value real(kind=rk)

generic, public :: get_Xg => get_Xg_all, get_Xgi, get_Xgid

Get geometry points

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_curve),	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_curve),	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_curve),	intent(in)	::	this
integer,	intent(in)	::	n
integer,	intent(in)	::	dir

Return Value real(kind=rk)

procedure, public :: get_Xt

Get parameter values

private pure function get_Xt(this) result(Xt)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

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

procedure, public :: get_continuity

Compute and return the continuity of the curve

private pure function get_continuity(this) result(c)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

Return Value integer, allocatable, (:)

procedure, public :: get_degree

Get degree of the NURBS curve

private pure function get_degree(this) result(degree)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

Return Value integer

procedure, public :: get_elem

Get IGA element connectivity

private pure function get_elem(this) result(elemConn)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

Return Value integer, allocatable, (:,:)

procedure, public :: get_elem_Xc_vis

Get connectivity for control points

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_curve),	intent(in)			::	this

Return Value integer, allocatable, (:,:)

procedure, public :: get_elem_Xg_vis

Get connectivity for geometry points

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_curve),	intent(in)			::	this

Return Value integer, allocatable, (:,:)

generic, public :: get_knot => get_knoti, get_knot_all

Get knot vector

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

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_curve),	intent(in)			::	this
integer,	intent(in)			::	i

Return Value real(kind=rk)

private pure function get_knot_all(this) result(knot)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

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

procedure, public :: get_multiplicity

Compute and return the multiplicity of the knots

private pure function get_multiplicity(this) result(m)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

Return Value integer, allocatable, (:)

procedure, public :: get_nc

Get number of control points

private pure function get_nc(this) result(nc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

Return Value integer

procedure, public :: get_ng

Get number of geometry points

private pure function get_ng(this) result(ng)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

Return Value integer

procedure, public :: insert_knots

Insert knots into the knot vector

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

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_curve),	intent(inout)			::	this
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(:,:)

procedure, public :: is_rational

Check if the NURBS curve is rational

private pure function is_rational(this) result(r)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

Return Value logical

procedure, public :: lsq_fit_bspline

Fit B-spline curve to structured data points using least squares

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_curve),	intent(inout)		::	this
real(kind=rk),	intent(in),	contiguous	::	Xt(:)
real(kind=rk),	intent(in),	contiguous	::	Xdata(:,:)
integer,	intent(in)		::	ndata

procedure, public :: lsq_fit_nurbs

Fit NURBS curve to structured data points using least squares

private pure subroutine lsq_fit_nurbs(this, Xt, Xdata, ndata, maxit, tol, lambda_xc, mu0, reg_logw)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_curve),	intent(inout)			::	this
real(kind=rk),	intent(in),		contiguous	::	Xt(:)
real(kind=rk),	intent(in),		contiguous	::	Xdata(:,:)
integer,	intent(in)			::	ndata
integer,	intent(in),	optional		::	maxit
real(kind=rk),	intent(in),	optional		::	tol
real(kind=rk),	intent(in),	optional		::	lambda_xc
real(kind=rk),	intent(in),	optional		::	mu0
real(kind=rk),	intent(in),	optional		::	reg_logw

procedure, public :: modify_Wc

Modify weights

private pure subroutine modify_Wc(this, W, num)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

procedure, public :: modify_Xc

Modify control points

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

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

procedure, public :: nearest_point

Find the nearest point on the NURBS curve (Approximation)

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_curve),	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
integer,	intent(out),	optional		::	id

procedure, public :: nearest_point2

Find the nearest point on the NURBS curve (Minimization - Newtons method)

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_curve),	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
real(kind=rk),	intent(out),	optional		::	nearest_Xg(size(this%Xc,2))

procedure, public :: remove_knots

Remove knots from the knot vector

private pure subroutine remove_knots(this, Xth, r)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

procedure, public :: rotate_Xc

Rotate control points

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

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

procedure, public :: rotate_Xg

Rotate geometry points

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

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

generic, public :: set => set1, set1a, set2, set3, set4

Set NURBS curve

private pure subroutine set1(this, knot, Xc, Wc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Set knot vector, control points and weights for the NURBS curve object.

Arguments

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

private pure subroutine set1a(this, knot, Xc, Wc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Set knot vector, control points and weights for the NURBS curve object.

Arguments

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

private pure subroutine set2(this, Xth_dir, degree, continuity, Xc, Wc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Set NURBS curve using nodes of parameter space (Xth), degree, continuity, control points and weights.

Arguments

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

private pure subroutine set3(this, Xc, Wc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

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

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_curve),	intent(inout)			::	this
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_curve),	intent(inout)			::	this
integer,	intent(in)			::	degree
integer,	intent(in)			::	nc
real(kind=rk),	intent(in),		contiguous	::	Xc(:,:)
real(kind=rk),	intent(in),	optional,	contiguous	::	Wc(:)

procedure, public :: set_C

Set a C-shape

private pure subroutine set_C(this, center, radius)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Attributes		Name
class(nurbs_curve),	intent(inout)		::	this
real(kind=rk),	intent(in),	contiguous	::	center(:)
real(kind=rk),	intent(in)		::	radius

procedure, public :: set_circle

Set a circle

private pure subroutine set_circle(this, center, radius)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Attributes		Name
class(nurbs_curve),	intent(inout)		::	this
real(kind=rk),	intent(in),	contiguous	::	center(:)
real(kind=rk),	intent(in)		::	radius

procedure, public :: set_elem

Set IGA element connectivity

private pure subroutine set_elem(this, elemConn)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

procedure, public :: set_elem_Xc_vis

Set connectivity for control points

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_curve),	intent(inout)			::	this
integer,	intent(in),		contiguous	::	elemConn(:,:)

procedure, public :: set_elem_Xg_vis

Set connectivity for geometry points

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_curve),	intent(inout)			::	this
integer,	intent(in),		contiguous	::	elemConn(:,:)

procedure, public :: set_half_circle

Set a half circle

private pure subroutine set_half_circle(this, center, radius)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Attributes		Name
class(nurbs_curve),	intent(inout)		::	this
real(kind=rk),	intent(in),	contiguous	::	center(:)
real(kind=rk),	intent(in)		::	radius

procedure, public :: show

Show the NURBS object using PyVista

private impure subroutine show(this, vtkfile_Xc, vtkfile_Xg)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent		Name
class(nurbs_curve),	intent(inout)	::	this
character(len=*),	intent(in)	::	vtkfile_Xc
character(len=*),	intent(in)	::	vtkfile_Xg

procedure, public :: translate_Xc

Translate control points

private pure subroutine translate_Xc(this, vec)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

procedure, public :: translate_Xg

Translate geometry points

private pure subroutine translate_Xg(this, vec)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

procedure, private :: basis_scalar

Compute the basis functions of the NURBS curve

private pure subroutine basis_scalar(this, Xt, Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

procedure, private :: basis_vector

Compute the basis functions of the NURBS curve

private pure subroutine basis_vector(this, res, Xt, Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

procedure, private :: derivative2_scalar

Compute the second derivative of the NURBS curve

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_curve),	intent(inout)			::	this
real(kind=rk),	intent(in)			::	Xt
real(kind=rk),	intent(out),		allocatable	::	d2Tgc(:)
real(kind=rk),	intent(out),	optional,	allocatable	::	dTgc(:)
real(kind=rk),	intent(out),	optional,	allocatable	::	Tgc(:)

procedure, private :: derivative2_vector

Compute the second derivative of the NURBS curve

private pure subroutine derivative2_vector(this, res, Xt, d2Tgc, dTgc, Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_curve),	intent(inout)			::	this
integer,	intent(in),	optional		::	res
real(kind=rk),	intent(in),	optional,	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(:,:)

procedure, private :: derivative_scalar

Compute the derivative of the NURBS curve

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_curve),	intent(inout)			::	this
real(kind=rk),	intent(in)			::	Xt
real(kind=rk),	intent(out),		allocatable	::	dTgc(:)
real(kind=rk),	intent(out),	optional,	allocatable	::	Tgc(:)
integer,	intent(in),	optional,	contiguous	::	elem(:)

procedure, private :: derivative_vector

Compute the derivative of the NURBS curve

private pure subroutine derivative_vector(this, res, Xt, dTgc, Tgc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

procedure, private :: get_Wc_all

Get all weights

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_curve),	intent(in)			::	this

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

procedure, private :: get_Wci

Get i-th weight

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_curve),	intent(in)			::	this
integer,	intent(in)			::	n

Return Value real(kind=rk)

procedure, private :: get_Xc_all

Get all control points

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_curve),	intent(in)			::	this

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

procedure, private :: get_Xci

Get i-th control point

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_curve),	intent(in)			::	this
integer,	intent(in)			::	n

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

procedure, private :: get_Xcid

Get i-th control point in a specific direction

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

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

Return Value real(kind=rk)

procedure, private :: get_Xg_all

Get all geometry points

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_curve),	intent(in)			::	this

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

procedure, private :: get_Xgi

Get i-th geometry point

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_curve),	intent(in)			::	this
integer,	intent(in)			::	n

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

procedure, private :: get_Xgid

Get i-th geometry point in a specific direction

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

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

Return Value real(kind=rk)

procedure, private :: get_knot_all

Get all knot vectors

private pure function get_knot_all(this) result(knot)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

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

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

procedure, private :: get_knoti

Get i-th knot value

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

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Arguments

Type	Intent	Optional	Attributes		Name
class(nurbs_curve),	intent(in)			::	this
integer,	intent(in)			::	i

Return Value real(kind=rk)

procedure, private :: set1

Set knot vector, control points and weights for the NURBS curve object

private pure subroutine set1(this, knot, Xc, Wc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Set knot vector, control points and weights for the NURBS curve object.

Arguments

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

procedure, private :: set1a

private pure subroutine set1a(this, knot, Xc, Wc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Set knot vector, control points and weights for the NURBS curve object.

Arguments

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

procedure, private :: set2

Set NURBS curve using nodes of parameter space, degree, continuity, control points and weights

private pure subroutine set2(this, Xth_dir, degree, continuity, Xc, Wc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

Set NURBS curve using nodes of parameter space (Xth), degree, continuity, control points and weights.

Arguments

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

procedure, private :: set3

Set Bezier or Rational Bezier curve using control points and weights

private pure subroutine set3(this, Xc, Wc)

Author: Seyed Ali Ghasemi
License: BSD 3-Clause

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

Arguments

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

procedure, private :: set4

Set NURBS curve using degree, number of control points, control points and weights

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_curve),	intent(inout)			::	this
integer,	intent(in)			::	degree
integer,	intent(in)			::	nc
real(kind=rk),	intent(in),		contiguous	::	Xc(:,:)
real(kind=rk),	intent(in),	optional,	contiguous	::	Wc(:)

nurbs_curve Derived Type

Contents

Variables

Type-Bound Procedures

type, public :: nurbs_curve

Inherits

Components

Type-Bound Procedures

procedure, public :: ansatz

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

Arguments

generic, public :: basis => basis_vector, basis_scalar

private pure subroutine basis_vector(this, res, Xt, Tgc)

Arguments

private pure subroutine basis_scalar(this, Xt, Tgc)

Arguments

procedure, public :: cmp_Xg

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

Arguments

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

procedure, public :: cmp_degree

private pure subroutine cmp_degree(this)

Arguments

procedure, public :: cmp_elem

private pure function cmp_elem(this) result(elemConn)

Arguments

Return Value integer, allocatable, (:,:)

procedure, public :: cmp_elem_Xc_vis

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

Arguments

Return Value integer, allocatable, (:,:)

procedure, public :: cmp_elem_Xg_vis

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

Arguments

Return Value integer, allocatable, (:,:)

procedure, public :: cmp_elem_Xth

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

Arguments

Return Value integer, allocatable, (:,:)

procedure, public :: cmp_length

private pure subroutine cmp_length(this, length, ngauss)

Arguments

procedure, public :: cmp_nc

private pure subroutine cmp_nc(this)

Arguments

procedure, public :: create

private pure subroutine create(this, res, Xt)

Arguments

generic, public :: derivative => derivative_vector, derivative_scalar

private pure subroutine derivative_vector(this, res, Xt, dTgc, Tgc)

Arguments

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

Arguments

generic, public :: derivative2 => derivative2_vector, derivative2_scalar

private pure subroutine derivative2_vector(this, res, Xt, d2Tgc, dTgc, Tgc)

Arguments

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

Arguments

procedure, public :: elevate_degree

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

Arguments

procedure, public :: export_Xc

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

Arguments

procedure, public :: export_Xg

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

Arguments

procedure, public :: export_Xth

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

Arguments

procedure, public :: export_iges

private impure subroutine export_iges(this, filename)

Arguments

procedure, public :: finalize

private pure elemental subroutine finalize(this)

Arguments

generic, public :: get_Wc => get_Wc_all, get_Wci

private pure function get_Wc_all(this) result(Wc)

Arguments

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