test_nurbs_volume Program

Uses

  • program~~test_nurbs_volume~~UsesGraph program~test_nurbs_volume test_nurbs_volume forunittest forunittest program~test_nurbs_volume->forunittest module~forcad forcad program~test_nurbs_volume->module~forcad module~forcad_utils forcad_utils program~test_nurbs_volume->module~forcad_utils module~forcad_kinds forcad_kinds module~forcad->module~forcad_kinds module~forcad_nurbs_curve forcad_nurbs_curve module~forcad->module~forcad_nurbs_curve module~forcad_nurbs_surface forcad_nurbs_surface module~forcad->module~forcad_nurbs_surface module~forcad_nurbs_volume forcad_nurbs_volume module~forcad->module~forcad_nurbs_volume module~forcad_utils->module~forcad_kinds module~forcad_nurbs_curve->module~forcad_utils module~forcad_nurbs_curve->module~forcad_kinds module~forcad_nurbs_surface->module~forcad_utils module~forcad_nurbs_surface->module~forcad_kinds module~forcad_nurbs_volume->module~forcad_utils module~forcad_nurbs_volume->module~forcad_kinds

Calls

program~~test_nurbs_volume~~CallsGraph program~test_nurbs_volume test_nurbs_volume check check program~test_nurbs_volume->check interface~ndgrid ndgrid program~test_nurbs_volume->interface~ndgrid n n program~test_nurbs_volume->n none~basis~3 nurbs_volume%basis program~test_nurbs_volume->none~basis~3 none~derivative2~3 nurbs_volume%derivative2 program~test_nurbs_volume->none~derivative2~3 none~derivative~3 nurbs_volume%derivative program~test_nurbs_volume->none~derivative~3 none~get_degree~3 nurbs_volume%get_degree program~test_nurbs_volume->none~get_degree~3 none~get_knot~3 nurbs_volume%get_knot program~test_nurbs_volume->none~get_knot~3 none~get_nc~3 nurbs_volume%get_nc program~test_nurbs_volume->none~get_nc~3 none~get_wc~3 nurbs_volume%get_Wc program~test_nurbs_volume->none~get_wc~3 none~get_xc~3 nurbs_volume%get_Xc program~test_nurbs_volume->none~get_xc~3 none~get_xg~3 nurbs_volume%get_Xg program~test_nurbs_volume->none~get_xg~3 none~set~3 nurbs_volume%set program~test_nurbs_volume->none~set~3 proc~cmp_elem_xc_vis~3 nurbs_volume%cmp_elem_Xc_vis program~test_nurbs_volume->proc~cmp_elem_xc_vis~3 proc~cmp_elem_xg_vis~3 nurbs_volume%cmp_elem_Xg_vis program~test_nurbs_volume->proc~cmp_elem_xg_vis~3 proc~cmp_elem~3 nurbs_volume%cmp_elem program~test_nurbs_volume->proc~cmp_elem~3 proc~cmp_nc~3 nurbs_volume%cmp_nc program~test_nurbs_volume->proc~cmp_nc~3 proc~cmp_volume nurbs_volume%cmp_volume program~test_nurbs_volume->proc~cmp_volume proc~create~3 nurbs_volume%create program~test_nurbs_volume->proc~create~3 proc~elevate_degree~3 nurbs_volume%elevate_degree program~test_nurbs_volume->proc~elevate_degree~3 proc~export_xc~3 nurbs_volume%export_Xc program~test_nurbs_volume->proc~export_xc~3 proc~export_xg~3 nurbs_volume%export_Xg program~test_nurbs_volume->proc~export_xg~3 proc~export_xth~3 nurbs_volume%export_Xth program~test_nurbs_volume->proc~export_xth~3 proc~finalize~3 nurbs_volume%finalize program~test_nurbs_volume->proc~finalize~3 proc~get_continuity~3 nurbs_volume%get_continuity program~test_nurbs_volume->proc~get_continuity~3 proc~get_elem_xc_vis~3 nurbs_volume%get_elem_Xc_vis program~test_nurbs_volume->proc~get_elem_xc_vis~3 proc~get_elem_xg_vis~3 nurbs_volume%get_elem_Xg_vis program~test_nurbs_volume->proc~get_elem_xg_vis~3 proc~get_elem~3 nurbs_volume%get_elem program~test_nurbs_volume->proc~get_elem~3 proc~get_multiplicity~3 nurbs_volume%get_multiplicity program~test_nurbs_volume->proc~get_multiplicity~3 proc~get_xt~3 nurbs_volume%get_Xt program~test_nurbs_volume->proc~get_xt~3 proc~insert_knots~3 nurbs_volume%insert_knots program~test_nurbs_volume->proc~insert_knots~3 proc~lsq_fit_bspline~3 nurbs_volume%lsq_fit_bspline program~test_nurbs_volume->proc~lsq_fit_bspline~3 proc~modify_wc~3 nurbs_volume%modify_Wc program~test_nurbs_volume->proc~modify_wc~3 proc~modify_xc~3 nurbs_volume%modify_Xc program~test_nurbs_volume->proc~modify_xc~3 proc~nearest_point2~3 nurbs_volume%nearest_point2 program~test_nurbs_volume->proc~nearest_point2~3 proc~nearest_point~3 nurbs_volume%nearest_point program~test_nurbs_volume->proc~nearest_point~3 proc~remove_knots~3 nurbs_volume%remove_knots program~test_nurbs_volume->proc~remove_knots~3 proc~rotate_xc~3 nurbs_volume%rotate_Xc program~test_nurbs_volume->proc~rotate_xc~3 proc~rotate_xg~3 nurbs_volume%rotate_Xg program~test_nurbs_volume->proc~rotate_xg~3 proc~set_c~3 nurbs_volume%set_C program~test_nurbs_volume->proc~set_c~3 proc~set_elem_xc_vis~3 nurbs_volume%set_elem_Xc_vis program~test_nurbs_volume->proc~set_elem_xc_vis~3 proc~set_elem_xg_vis~3 nurbs_volume%set_elem_Xg_vis program~test_nurbs_volume->proc~set_elem_xg_vis~3 proc~set_elem~3 nurbs_volume%set_elem program~test_nurbs_volume->proc~set_elem~3 proc~set_half_ring~2 nurbs_volume%set_half_ring program~test_nurbs_volume->proc~set_half_ring~2 proc~set_hexahedron nurbs_volume%set_hexahedron program~test_nurbs_volume->proc~set_hexahedron proc~set_ring~2 nurbs_volume%set_ring program~test_nurbs_volume->proc~set_ring~2 proc~translate_xc~3 nurbs_volume%translate_Xc program~test_nurbs_volume->proc~translate_xc~3 proc~translate_xg~3 nurbs_volume%translate_Xg program~test_nurbs_volume->proc~translate_xg~3 xdata xdata program~test_nurbs_volume->xdata xg_eval xg_eval program~test_nurbs_volume->xg_eval xt xt program~test_nurbs_volume->xt xt1 xt1 program~test_nurbs_volume->xt1 xt2 xt2 program~test_nurbs_volume->xt2 xt3 xt3 program~test_nurbs_volume->xt3 proc~ndgrid2 ndgrid2 interface~ndgrid->proc~ndgrid2 proc~ndgrid3 ndgrid3 interface~ndgrid->proc~ndgrid3 proc~basis_scalar~3 nurbs_volume%basis_scalar none~basis~3->proc~basis_scalar~3 proc~basis_vector~3 nurbs_volume%basis_vector none~basis~3->proc~basis_vector~3 proc~derivative2_scalar~3 nurbs_volume%derivative2_scalar none~derivative2~3->proc~derivative2_scalar~3 proc~derivative2_vector~3 nurbs_volume%derivative2_vector none~derivative2~3->proc~derivative2_vector~3 proc~derivative_scalar~3 nurbs_volume%derivative_scalar none~derivative~3->proc~derivative_scalar~3 proc~derivative_vector~3 nurbs_volume%derivative_vector none~derivative~3->proc~derivative_vector~3 proc~get_degree_all~2 nurbs_volume%get_degree_all none~get_degree~3->proc~get_degree_all~2 proc~get_degree_dir~2 nurbs_volume%get_degree_dir none~get_degree~3->proc~get_degree_dir~2 proc~get_knot_all~3 nurbs_volume%get_knot_all none~get_knot~3->proc~get_knot_all~3 proc~get_knoti~3 nurbs_volume%get_knoti none~get_knot~3->proc~get_knoti~3 proc~get_nc_all~2 nurbs_volume%get_nc_all none~get_nc~3->proc~get_nc_all~2 proc~get_nc_dir~2 nurbs_volume%get_nc_dir none~get_nc~3->proc~get_nc_dir~2 proc~get_wc_all~3 nurbs_volume%get_Wc_all none~get_wc~3->proc~get_wc_all~3 proc~get_wci~3 nurbs_volume%get_Wci none~get_wc~3->proc~get_wci~3 proc~get_xc_all~3 nurbs_volume%get_Xc_all none~get_xc~3->proc~get_xc_all~3 proc~get_xcid~3 nurbs_volume%get_Xcid none~get_xc~3->proc~get_xcid~3 proc~get_xci~3 nurbs_volume%get_Xci none~get_xc~3->proc~get_xci~3 proc~get_xg_all~3 nurbs_volume%get_Xg_all none~get_xg~3->proc~get_xg_all~3 proc~get_xgid~3 nurbs_volume%get_Xgid none~get_xg~3->proc~get_xgid~3 proc~get_xgi~3 nurbs_volume%get_Xgi none~get_xg~3->proc~get_xgi~3 proc~set1~3 nurbs_volume%set1 none~set~3->proc~set1~3 proc~set2~3 nurbs_volume%set2 none~set~3->proc~set2~3 proc~set3~3 nurbs_volume%set3 none~set~3->proc~set3~3 proc~set4~3 nurbs_volume%set4 none~set~3->proc~set4~3 interface~elemconn_c0 elemConn_C0 proc~cmp_elem_xc_vis~3->interface~elemconn_c0 proc~cmp_elem_xg_vis~3->interface~elemconn_c0 proc~cmp_elem~3->proc~get_multiplicity~3 interface~elemconn_cn elemConn_Cn proc~cmp_elem~3->interface~elemconn_cn interface~unique unique proc~cmp_elem~3->interface~unique interface~compute_multiplicity compute_multiplicity proc~cmp_nc~3->interface~compute_multiplicity proc~cmp_volume->proc~cmp_elem~3 proc~ansatz~3 nurbs_volume%ansatz proc~cmp_volume->proc~ansatz~3 reduce reduce proc~cmp_volume->reduce proc~create~3->interface~ndgrid interface~compute_xg~3 compute_Xg proc~create~3->interface~compute_xg~3 proc~is_rational~3 nurbs_volume%is_rational proc~create~3->proc~is_rational~3 proc~elevate_degree~3->none~get_knot~3 proc~elevate_degree~3->none~set~3 proc~elevate_degree_a_5_9 elevate_degree_A_5_9 proc~elevate_degree~3->proc~elevate_degree_a_5_9 proc~export_xc~3->proc~cmp_elem_xc_vis~3 proc~export_vtk_legacy export_vtk_legacy proc~export_xc~3->proc~export_vtk_legacy proc~export_xg~3->proc~cmp_elem_xg_vis~3 proc~export_xg~3->proc~export_vtk_legacy proc~export_xth~3->interface~ndgrid proc~export_xth~3->none~set~3 proc~export_xth~3->proc~cmp_elem~3 proc~export_xth~3->interface~unique proc~export_xth~3->proc~export_vtk_legacy proc~get_continuity~3->interface~compute_multiplicity proc~get_multiplicity~3->interface~compute_multiplicity proc~insert_knots~3->none~get_knot~3 proc~insert_knots~3->none~set~3 proc~insert_knots~3->interface~compute_multiplicity proc~findspan findspan proc~insert_knots~3->proc~findspan proc~insert_knot_a_5_1 insert_knot_A_5_1 proc~insert_knots~3->proc~insert_knot_a_5_1 proc~basis_bspline basis_bspline proc~lsq_fit_bspline~3->proc~basis_bspline proc~kron kron proc~lsq_fit_bspline~3->proc~kron proc~solve solve proc~lsq_fit_bspline~3->proc~solve proc~modify_wc~3->none~get_knot~3 proc~modify_wc~3->none~get_wc~3 proc~modify_wc~3->none~get_xc~3 proc~modify_wc~3->none~set~3 proc~modify_xc~3->none~get_knot~3 proc~modify_xc~3->none~get_wc~3 proc~modify_xc~3->none~get_xc~3 proc~modify_xc~3->none~set~3 proc~nearest_point2~3->none~derivative2~3 proc~nearest_point2~3->proc~create~3 proc~nearest_point2~3->proc~finalize~3 proc~nearest_point2~3->proc~nearest_point~3 proc~cmp_xg~3 nurbs_volume%cmp_Xg proc~nearest_point2~3->proc~cmp_xg~3 proc~inv inv proc~nearest_point2~3->proc~inv proc~nearest_point_help_3d nearest_point_help_3d proc~nearest_point~3->proc~nearest_point_help_3d proc~remove_knots~3->none~get_knot~3 proc~remove_knots~3->none~set~3 proc~remove_knots~3->interface~compute_multiplicity proc~remove_knots~3->proc~findspan proc~remove_knots_a_5_8 remove_knots_A_5_8 proc~remove_knots~3->proc~remove_knots_a_5_8 proc~rotation rotation proc~rotate_xc~3->proc~rotation proc~rotate_xg~3->proc~rotation proc~set_c~3->none~set~3 proc~set_half_ring~2->none~set~3 proc~set_hexahedron->none~set~3 proc~hexahedron_xc hexahedron_Xc proc~set_hexahedron->proc~hexahedron_xc proc~set_ring~2->none~set~3 proc~compute_multiplicity1 compute_multiplicity1 interface~compute_multiplicity->proc~compute_multiplicity1 proc~compute_multiplicity2 compute_multiplicity2 interface~compute_multiplicity->proc~compute_multiplicity2 proc~compute_xg_bspline_3d compute_Xg_bspline_3d interface~compute_xg~3->proc~compute_xg_bspline_3d proc~compute_xg_bspline_3d_1point compute_Xg_bspline_3d_1point interface~compute_xg~3->proc~compute_xg_bspline_3d_1point proc~compute_xg_nurbs_3d compute_Xg_nurbs_3d interface~compute_xg~3->proc~compute_xg_nurbs_3d proc~compute_xg_nurbs_3d_1point compute_Xg_nurbs_3d_1point interface~compute_xg~3->proc~compute_xg_nurbs_3d_1point proc~cmp_elemconn_c0_l cmp_elemConn_C0_L interface~elemconn_c0->proc~cmp_elemconn_c0_l proc~cmp_elemconn_c0_s cmp_elemConn_C0_S interface~elemconn_c0->proc~cmp_elemconn_c0_s proc~cmp_elemconn_c0_v cmp_elemConn_C0_V interface~elemconn_c0->proc~cmp_elemconn_c0_v proc~cmp_elemconn_cn_l cmp_elemConn_Cn_L interface~elemconn_cn->proc~cmp_elemconn_cn_l proc~cmp_elemconn_cn_s cmp_elemConn_Cn_S interface~elemconn_cn->proc~cmp_elemconn_cn_s proc~cmp_elemconn_cn_v cmp_elemConn_Cn_V interface~elemconn_cn->proc~cmp_elemconn_cn_v proc~unique_integer unique_integer interface~unique->proc~unique_integer proc~unique_real unique_real interface~unique->proc~unique_real proc~ansatz~3->interface~ndgrid proc~ansatz~3->none~derivative~3 proc~ansatz~3->none~set~3 proc~ansatz~3->proc~cmp_elem~3 proc~ansatz~3->interface~unique proc~ansatz~3->proc~inv interface~gauss_leg gauss_leg proc~ansatz~3->interface~gauss_leg proc~det det proc~ansatz~3->proc~det proc~basis_scalar~3->proc~is_rational~3 interface~compute_tgc~3 compute_Tgc proc~basis_scalar~3->interface~compute_tgc~3 proc~basis_vector~3->interface~ndgrid proc~basis_vector~3->proc~is_rational~3 proc~basis_vector~3->interface~compute_tgc~3 proc~cmp_xg~3->interface~compute_xg~3 proc~cmp_xg~3->proc~is_rational~3 proc~derivative2_scalar~3->proc~is_rational~3 interface~compute_d2tgc~3 compute_d2Tgc proc~derivative2_scalar~3->interface~compute_d2tgc~3 proc~derivative2_vector~3->interface~ndgrid proc~derivative2_vector~3->proc~is_rational~3 proc~derivative2_vector~3->interface~compute_d2tgc~3 proc~derivative_scalar~3->proc~is_rational~3 interface~compute_dtgc~3 compute_dTgc proc~derivative_scalar~3->interface~compute_dtgc~3 proc~derivative_vector~3->interface~ndgrid proc~derivative_vector~3->proc~is_rational~3 proc~derivative_vector~3->interface~compute_dtgc~3 proc~elevate_degree_a_5_9->interface~compute_multiplicity proc~bincoeff bincoeff proc~elevate_degree_a_5_9->proc~bincoeff proc~get_nc_dir~2->interface~compute_multiplicity proc~inv->proc~inv proc~inv->proc~det cosd cosd proc~rotation->cosd sind sind proc~rotation->sind proc~set1~3->proc~cmp_nc~3 proc~cmp_degree~3 nurbs_volume%cmp_degree proc~set1~3->proc~cmp_degree~3 proc~set2~3->proc~cmp_nc~3 proc~compute_knot_vector compute_knot_vector proc~set2~3->proc~compute_knot_vector proc~set3~3->proc~cmp_degree~3 proc~compute_d2tgc_bspline_3d_scalar compute_d2Tgc_bspline_3d_scalar interface~compute_d2tgc~3->proc~compute_d2tgc_bspline_3d_scalar proc~compute_d2tgc_bspline_3d_vector compute_d2Tgc_bspline_3d_vector interface~compute_d2tgc~3->proc~compute_d2tgc_bspline_3d_vector proc~compute_d2tgc_nurbs_3d_scalar compute_d2Tgc_nurbs_3d_scalar interface~compute_d2tgc~3->proc~compute_d2tgc_nurbs_3d_scalar proc~compute_d2tgc_nurbs_3d_vector compute_d2Tgc_nurbs_3d_vector interface~compute_d2tgc~3->proc~compute_d2tgc_nurbs_3d_vector proc~compute_dtgc_bspline_3d_scalar compute_dTgc_bspline_3d_scalar interface~compute_dtgc~3->proc~compute_dtgc_bspline_3d_scalar proc~compute_dtgc_bspline_3d_vector compute_dTgc_bspline_3d_vector interface~compute_dtgc~3->proc~compute_dtgc_bspline_3d_vector proc~compute_dtgc_nurbs_3d_scalar compute_dTgc_nurbs_3d_scalar interface~compute_dtgc~3->proc~compute_dtgc_nurbs_3d_scalar proc~compute_dtgc_nurbs_3d_vector compute_dTgc_nurbs_3d_vector interface~compute_dtgc~3->proc~compute_dtgc_nurbs_3d_vector proc~compute_tgc_bspline_3d_scalar compute_Tgc_bspline_3d_scalar interface~compute_tgc~3->proc~compute_tgc_bspline_3d_scalar proc~compute_tgc_bspline_3d_vector compute_Tgc_bspline_3d_vector interface~compute_tgc~3->proc~compute_tgc_bspline_3d_vector proc~compute_tgc_nurbs_3d_scalar compute_Tgc_nurbs_3d_scalar interface~compute_tgc~3->proc~compute_tgc_nurbs_3d_scalar proc~compute_tgc_nurbs_3d_vector compute_Tgc_nurbs_3d_vector interface~compute_tgc~3->proc~compute_tgc_nurbs_3d_vector proc~gauss_legendre_1d gauss_legendre_1D interface~gauss_leg->proc~gauss_legendre_1d proc~gauss_legendre_2d gauss_legendre_2D interface~gauss_leg->proc~gauss_legendre_2d proc~gauss_legendre_3d gauss_legendre_3D interface~gauss_leg->proc~gauss_legendre_3d proc~factln factln proc~bincoeff->proc~factln proc~cmp_degree~3->proc~get_multiplicity~3 proc~repelem repelem proc~compute_knot_vector->proc~repelem proc~compute_xg_bspline_3d->proc~basis_bspline proc~compute_xg_bspline_3d->proc~kron proc~compute_xg_bspline_3d_1point->proc~basis_bspline proc~compute_xg_bspline_3d_1point->proc~kron proc~cmp_tgc_3d cmp_Tgc_3d proc~compute_xg_nurbs_3d->proc~cmp_tgc_3d proc~compute_xg_nurbs_3d_1point->proc~basis_bspline proc~compute_xg_nurbs_3d_1point->proc~kron proc~cmp_tgc_3d->proc~basis_bspline proc~cmp_tgc_3d->proc~kron proc~compute_d2tgc_bspline_3d_scalar->proc~kron interface~basis_bspline_2der basis_bspline_2der proc~compute_d2tgc_bspline_3d_scalar->interface~basis_bspline_2der proc~compute_d2tgc_bspline_3d_vector->proc~kron proc~compute_d2tgc_bspline_3d_vector->interface~basis_bspline_2der proc~compute_d2tgc_nurbs_3d_scalar->proc~kron proc~compute_d2tgc_nurbs_3d_scalar->interface~basis_bspline_2der proc~compute_d2tgc_nurbs_3d_vector->proc~kron proc~compute_d2tgc_nurbs_3d_vector->interface~basis_bspline_2der proc~compute_dtgc_bspline_3d_scalar->proc~kron interface~basis_bspline_der basis_bspline_der proc~compute_dtgc_bspline_3d_scalar->interface~basis_bspline_der proc~compute_dtgc_bspline_3d_vector->proc~kron proc~compute_dtgc_bspline_3d_vector->interface~basis_bspline_der proc~compute_dtgc_nurbs_3d_scalar->proc~kron proc~compute_dtgc_nurbs_3d_scalar->interface~basis_bspline_der proc~compute_dtgc_nurbs_3d_vector->proc~kron proc~compute_dtgc_nurbs_3d_vector->interface~basis_bspline_der proc~compute_tgc_bspline_3d_scalar->proc~basis_bspline proc~compute_tgc_bspline_3d_scalar->proc~kron proc~compute_tgc_bspline_3d_vector->proc~basis_bspline proc~compute_tgc_bspline_3d_vector->proc~kron proc~compute_tgc_nurbs_3d_scalar->proc~basis_bspline proc~compute_tgc_nurbs_3d_scalar->proc~kron proc~compute_tgc_nurbs_3d_vector->proc~basis_bspline proc~compute_tgc_nurbs_3d_vector->proc~kron proc~gauss_legendre gauss_legendre proc~gauss_legendre_1d->proc~gauss_legendre proc~gauss_legendre_2d->interface~ndgrid proc~gauss_legendre_2d->proc~kron proc~gauss_legendre_2d->proc~gauss_legendre proc~gauss_legendre_3d->interface~ndgrid proc~gauss_legendre_3d->proc~kron proc~gauss_legendre_3d->proc~gauss_legendre proc~basis_bspline_2der_a basis_bspline_2der_A interface~basis_bspline_2der->proc~basis_bspline_2der_a proc~basis_bspline_2der_b basis_bspline_2der_B interface~basis_bspline_2der->proc~basis_bspline_2der_b proc~basis_bspline_2der_c basis_bspline_2der_C interface~basis_bspline_2der->proc~basis_bspline_2der_c proc~basis_bspline_der_a basis_bspline_der_A interface~basis_bspline_der->proc~basis_bspline_der_a proc~basis_bspline_der_b basis_bspline_der_B interface~basis_bspline_der->proc~basis_bspline_der_b

Variables

Type Attributes Name Initial
real(kind=rk), allocatable :: Tgc(:,:)
real(kind=rk), allocatable :: Tgc1(:)
real(kind=rk), allocatable :: Tgc1b(:)
real(kind=rk), allocatable :: Tgcb(:,:)
real(kind=rk), allocatable :: Wc(:)
real(kind=rk), allocatable :: Xc(:,:)
real(kind=rk), allocatable :: Xg(:,:)
real(kind=rk), allocatable :: Xgb(:,:)
type(nurbs_volume) :: bsp
real(kind=rk), allocatable :: d2Tgc(:,:,:)
real(kind=rk), allocatable :: d2Tgc1(:,:)
real(kind=rk), allocatable :: d2Tgc1b(:,:)
real(kind=rk), allocatable :: d2Tgcb(:,:,:)
real(kind=rk), allocatable :: dTgc(:,:,:)
real(kind=rk), allocatable :: dTgc1(:,:)
real(kind=rk), allocatable :: dTgc1b(:,:)
real(kind=rk), allocatable :: dTgcb(:,:,:)
integer, allocatable :: elemConn(:,:)
integer :: i
integer :: id
real(kind=rk) :: knot1(4)
real(kind=rk) :: knot2(4)
real(kind=rk) :: knot3(4)
real(kind=rk), allocatable :: nearest_Xg(:)
real(kind=rk), allocatable :: nearest_Xt(:)
type(nurbs_volume) :: nurbs
type(unit_test) :: ut
real(kind=rk) :: volume
real(kind=rk) :: volumeb