shape_C_2d.f90 Source File

program shape_C_2d

    use forcad, only: rk, nurbs_surface

    implicit none
    type(nurbs_surface) :: shape

    !-----------------------------------------------------------------------------
    ! Setting up NURBS C-shape
    !-----------------------------------------------------------------------------

    !> Set a C-shape with radius 2.0 and center at [0.0, 0.0, 0.0]
    call shape%set_C(center = [0.0_rk, 0.0_rk, 0.0_rk], radius1 = 1.0_rk, radius2 = 2.0_rk)

    !> Export control points to a VTK file
    call shape%export_Xc('vtk/shape_C_2d_Xc.vtk')

    !-----------------------------------------------------------------------------
    ! Creating C-shape
    !-----------------------------------------------------------------------------

    !> Generate the NURBS C-shape with a resolution of 100
    call shape%create(100, 60)

    !> Export the generated cirlce to a VTK file
    call shape%export_Xg('vtk/shape_C_2d_Xg.vtk')

    !-----------------------------------------------------------------------------
    ! Visualization using PyVista
    ! Note: PyVista is required for visualization. Install it using `pip install pyvista`
    !-----------------------------------------------------------------------------

    !> Show the control geometry and geometry using PyVista
    call shape%show('vtk/shape_C_2d_Xc.vtk','vtk/shape_C_2d_Xg.vtk')

    !-----------------------------------------------------------------------------
    ! Finalizing
    !-----------------------------------------------------------------------------

    !> Finalize the NURBS curve object
    call shape%finalize()

end program