Merge pull request #17 from gha3mi/dev

Added New Features

CHANGELOG.md

@@ -4,6 +4,21 @@ All notable changes to [ForCAD](https://github.com/gha3mi/forcad) will be docume
 
 The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.1.0/),
 
+## [0.6.0]
+
+### Added
+
+- Added generic `get_nc()` method to `nurbs_surface` and `nurbs_volume` derived types.
+- Added `det`, `inv`, `dyad` and `gauss_leg` procedures in the `forcad_utils` module.
+- Added procedure `set1a` to the `nurbs_curve` derived type.
+- Added procedure `set4` to `nurbs_curve`,  `nurbs_area` and `nurbs_volume` derived types.
+- Added optional input variable `elem` to `derivative_scalar` procedures.
+- Added `ansatz` procedures to compute shape functions, derivatives of shape functions and (dV, dA, dL).
+- Added `cmp_length()` to compute the length of a NURBS curve.
+- Added `cmp_area()` to compute the area of a NURBS surface.
+- Added `cmp_volume()` to compute the volume of a NURBS volume.
+- Added examples for `cmp_length()`, `cmp_area()`, and `cmp_volume()`.
+
 ## [0.5.1]
 
 ### Changed

ROADMAP.md

@@ -27,6 +27,7 @@ This roadmap outlines upcoming features and enhancements for ForCAD. Contributio
     - [x] Improvements and bug fixes.
 
 - Additional Ideas:
+    - [ ] Improve `ROADMAP.md`.
     - [ ] Add `reduce_degree()` method for curves, surfaces and volumes.
     - [ ] Add `morph()` method for morphing box.
     - [ ] Add support binary `vtk` files.
@@ -48,5 +49,5 @@ This roadmap outlines upcoming features and enhancements for ForCAD. Contributio
     - [ ] Improve documentation.
     - [ ] Add offset method.
     - [ ] Add curve, surface and volume fitting.
-    - [ ] Calculate length, area and volume.
+    - [x] Calculate length, area and volume.
     - [ ] Think about a GUI.

VERSION

@@ -1 +1 @@
-0.4.1
+0.6.0

example/cmp_area.f90

@@ -0,0 +1,23 @@
+program compute_area
+
+    use forcad
+
+    implicit none
+
+    type(nurbs_surface) :: shape
+    real(rk) :: area
+    real(rk) :: Xc(4,3)
+
+    Xc(1,:) = [0.0_rk, 0.0_rk, 0.0_rk]
+    Xc(2,:) = [2.0_rk, 0.0_rk, 0.0_rk]
+    Xc(3,:) = [0.0_rk, 2.0_rk, 0.0_rk]
+    Xc(4,:) = [2.0_rk, 2.0_rk, 0.0_rk]
+
+    call shape%set(&
+        knot1=[0.0_rk, 0.0_rk, 1.0_rk, 1.0_rk],&
+        knot2=[0.0_rk, 0.0_rk, 1.0_rk, 1.0_rk],&
+        Xc=Xc)
+
+    call shape%cmp_area(area)
+    print*, area
+end program

example/cmp_length.f90

@@ -0,0 +1,20 @@
+program compute_length
+
+    use forcad
+
+    implicit none
+
+    type(nurbs_curve) :: shape
+    real(rk) :: length
+    real(rk) :: Xc(2,3)
+
+    Xc(1,:) = [0.0_rk, 0.0_rk, 0.0_rk]
+    Xc(2,:) = [2.0_rk, 0.0_rk, 0.0_rk]
+
+    call shape%set(&
+        knot=[0.0_rk, 0.0_rk, 1.0_rk, 1.0_rk],&
+        Xc=Xc)
+
+    call shape%cmp_length(length)
+    print*, length
+end program

example/cmp_volume.f90

@@ -0,0 +1,28 @@
+program compute_volume
+
+    use forcad
+
+    implicit none
+
+    type(nurbs_volume) :: shape
+    real(rk) :: volume
+    real(rk) :: Xc(8,3)
+
+    Xc(1,:) = [0.0_rk, 0.0_rk, 0.0_rk]
+    Xc(2,:) = [2.0_rk, 0.0_rk, 0.0_rk]
+    Xc(3,:) = [0.0_rk, 2.0_rk, 0.0_rk]
+    Xc(4,:) = [2.0_rk, 2.0_rk, 0.0_rk]
+    Xc(5,:) = [0.0_rk, 0.0_rk, 2.0_rk]
+    Xc(6,:) = [2.0_rk, 0.0_rk, 2.0_rk]
+    Xc(7,:) = [0.0_rk, 2.0_rk, 2.0_rk]
+    Xc(8,:) = [2.0_rk, 2.0_rk, 2.0_rk]
+
+    call shape%set(&
+        knot1=[0.0_rk, 0.0_rk, 1.0_rk, 1.0_rk],&
+        knot2=[0.0_rk, 0.0_rk, 1.0_rk, 1.0_rk],&
+        knot3=[0.0_rk, 0.0_rk, 1.0_rk, 1.0_rk],&
+        Xc=Xc)
+
+    call shape%cmp_volume(volume)
+    print*, volume
+end program

ford.yml

@@ -1,5 +1,5 @@
 project:             ForCAD
-version:             0.5.1
+version:             0.6.0
 year:                2024
 summary:             <p style="text-align: center"><img alt="ForCAD" src="https://github.com/gha3mi/forcad/raw/main/logo/logo.png" style="max-width: 100%; height: auto;"></p>
 project_github:      https://github.com/gha3mi/forcad

fpm.toml

@@ -1,5 +1,5 @@
 name       = "forcad"
-version    = "0.5.1"
+version    = "0.6.0"
 license    = "BSD 3-Clause License"
 author     = "Seyed Ali Ghasemi"
 maintainer = "info@gha3mi.com"
@@ -26,6 +26,7 @@ forimage    = {git="https://github.com/gha3mi/forimage.git"}
 forcolormap = {git="https://github.com/vmagnin/forcolormap.git"}
 fortime     = {git="https://github.com/gha3mi/fortime.git"}
 openmp      = "*"
+stdlib      = "*"
 
 [dev-dependencies]
 forunittest = {git="https://github.com/gha3mi/forunittest"}
@@ -158,4 +159,19 @@ main       = "nearest_point_3d.f90"
 [[example]]
 name       = "nearest_point_2d_bench"
 source-dir = "example"
-main       = "nearest_point_2d_bench.f90"
+main       = "nearest_point_2d_bench.f90"
+
+[[example]]
+name       = "cmp_length"
+source-dir = "example"
+main       = "cmp_length.f90"
+
+[[example]]
+name       = "cmp_area"
+source-dir = "example"
+main       = "cmp_area.f90"
+
+[[example]]
+name       = "cmp_volume"
+source-dir = "example"
+main       = "cmp_volume.f90"