shape_half_ring_3d.f90 Source File

program shape_half_ring_3d

    use forcad, only: rk, nurbs_volume

    implicit none
    type(nurbs_volume) :: shape


    !> Set up a half ring centered at 0,0,0 with inner radius 1, outer radius 2, and length 1.
    call shape%set_half_ring([0.0_rk, 0.0_rk, 0.0_rk], 1.0_rk, 2.0_rk, 1.0_rk)

    !> Export the control points to a VTK file for visualization.
    call shape%export_Xc('vtk/shape_half_ring_3d_Xc.vtk')

    !> Create the shape using the specified number of elements in each direction.
    call shape%create(60, 15, 10)

    !> Export the geometry to a VTK file for visualization.
    call shape%export_Xg('vtk/shape_half_ring_3d_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_half_ring_3d_Xc.vtk','vtk/shape_half_ring_3d_Xg.vtk')

    !> Finalize and clean up the shape object.
    call shape%finalize()

end program