Updating docs for year 2025

README.md

@@ -58,7 +58,7 @@ First we need to configure our project by modyfing `docusaurus.config.js` file.
 
 - `organizationName`, The GitHub user or organization name.
 - `projectName`, The name of the deployment repository.
-- `deploymentBranch`, The name of the deployment branch. Bu default it is `gh-pages` for non-organization GitHub pages repositories (that is, `projectName` not ending in `.github.io`). Otherwise, it needs to be explicit as a config field or environment variable.
+- `deploymentBranch`, The name of the deployment branch. But default it is `gh-pages` for non-organization GitHub pages repositories (that is, `projectName` not ending in `.github.io`). Otherwise, it needs to be explicit as a config field or environment variable.
 
 In this repository we have set the following parameters in `docusaurus.config.js` file:
 

docs/blog/DirichletBC/code/mesh/mesh.F90

@@ -0,0 +1,88 @@
+!> author: Vikas Sharma, Ph. D.
+! date: 2025-10-02
+! summary: A Simple mesh generator
+
+PROGRAM main
+USE Gmsh_Class
+USE GlobalData
+USE MSHFile_Class
+USE HDF5File_Class
+USE Display_Method
+
+TYPE(Gmsh_) :: gmsh
+
+CHARACTER(LEN=*), PARAMETER :: prefix = "square"
+REAL(DFP), PARAMETER :: lx = 1.0_DFP
+REAL(DFP), PARAMETER :: ly = 1.0_DFP
+INTEGER(I4B), PARAMETER :: order = 1
+
+REAL(DFP), PARAMETER :: meshSize = 1.0
+REAL(DFP), PARAMETER :: corner(3) = 0.0_DFP
+INTEGER(I4B), PARAMETER :: nnod(2) = [3, 3]
+INTEGER(I4B) :: ierr
+REAL(DFP) :: x, y, z, lc
+CHARACTER(:), ALLOCATABLE :: suffix, title
+
+TYPE(MSHFile_) :: mshFile
+TYPE(HDF5File_) :: hdf5file
+
+suffix = ToString(nnod(1))//"x"//ToString(nnod(2))
+
+title = prefix//"_"//suffix
+
+ierr = gmsh%Initialize()
+ierr = gmsh%model%add(title)
+
+x = corner(1); y = corner(2); z = corner(3); lc = meshSize
+ierr = gmsh%model%geo%addPoint(x=x, y=y, z=z, meshSize=lc, tag=1)
+
+x = x + lx; y = y; z = z; lc = lc
+ierr = gmsh%model%geo%addPoint(x=x, y=y, z=z, meshSize=lc, tag=2)
+
+x = x; y = y + ly; z = z; lc = lc
+ierr = gmsh%model%geo%addPoint(x=x, y=y, z=z, meshSize=lc, tag=3)
+
+x = corner(1); y = y; z = z; lc = lc
+ierr = gmsh%model%geo%addPoint(x=x, y=y, z=z, meshSize=lc, tag=4)
+
+ierr = gmsh%model%geo%addLine(1, 2, 1)
+ierr = gmsh%model%geo%addLine(2, 3, 2)
+ierr = gmsh%model%geo%addLine(3, 4, 3)
+ierr = gmsh%model%geo%addLine(4, 1, 4)
+
+ierr = gmsh%model%geo%addCurveLoop([1, 2, 3, 4], tag=1)
+
+ierr = gmsh%model%geo%addPlaneSurface([1], 1)
+
+ierr = gmsh%model%geo%mesh%setTransfiniteCurve(1, nnod(1), &
+                                               "Progression", 1.0_DFP)
+ierr = gmsh%model%geo%mesh%setTransfiniteCurve(3, nnod(1), &
+                                               "Progression", 1.0_DFP)
+ierr = gmsh%model%geo%mesh%setTransfiniteCurve(2, nnod(2), &
+                                               "Progression", 1.0_DFP)
+ierr = gmsh%model%geo%mesh%setTransfiniteCurve(4, nnod(2), &
+                                               "Progression", 1.0_DFP)
+
+ierr = gmsh%model%geo%mesh%setTransfiniteSurface(1)
+ierr = gmsh%model%geo%mesh%setRecombine(2, 1)
+
+ierr = gmsh%model%geo%synchronize()
+
+ierr = gmsh%model%mesh%generate(2)
+
+ierr = gmsh%model%mesh%setOrder(order)
+
+ierr = gmsh%WRITE(title//'.msh')
+
+ierr = gmsh%Finalize()
+
+CALL mshFile%Initiate(filename=title//'.msh', STATUS="OLD", ACTION="READ")
+CALL mshFile%OPEN()
+CALL mshFile%READ()
+CALL hdf5file%Initiate(title//'.h5', MODE="NEW")
+CALL hdf5file%OPEN()
+CALL mshFile%Export(hdf5=hdf5file, group="")
+CALL mshFile%DEALLOCATE()
+CALL hdf5file%DEALLOCATE()
+
+END PROGRAM main

docs/blog/DirichletBC/code/mesh/runner.toml

@@ -0,0 +1,3 @@
+BuildType = "Debug"
+BuildDir = "/tmp/easifem-tests/Blog/understanding-fedof/mesh/"
+TargetLibs = ["easifemClasses"]

docs/blog/DirichletBC/code/runner.toml

@@ -0,0 +1,3 @@
+BuildType = "Debug"
+BuildDir = "/tmp/easifem-tests/Blog/DirichletBC/"
+TargetLibs = ["easifemClasses"]

docs/blog/DirichletBC/code/test1.F90

@@ -0,0 +1,104 @@
+!> author: Vikas Sharma, Ph. D.
+! date: 2025-10-23
+! summary: Study of DirichletBC class
+
+PROGRAM main
+USE GlobalData, ONLY: DFP, I4B, LGT
+USE ExceptionHandler_Class, ONLY: e, EXCEPTION_INFORMATION
+USE DirichletBC_Class
+USE FEDomain_Class
+USE AbstractMesh_Class
+USE FEDOF_Class
+USE ScalarField_Class
+USE Display_Method
+USE ReallocateUtility
+
+IMPLICIT NONE
+
+!----------------------------------------------------------------------------
+! Parameters
+!----------------------------------------------------------------------------
+CHARACTER(*), PARAMETER :: tomlFileName = "./test1.toml"
+
+!----------------------------------------------------------------------------
+! Types and variables
+!----------------------------------------------------------------------------
+TYPE(DirichletBC_) :: obj
+TYPE(FEDomain_) :: dom
+TYPE(ScalarField_) :: u
+TYPE(FEDOF_) :: fedof, geofedof
+INTEGER(I4B) :: tsize, iNodeOnNode, iNodeOnEdge, iNodeOnFace, nrow, ncol
+
+!----------------------------------------------------------------------------
+! Allocatables and pointers
+!----------------------------------------------------------------------------
+CLASS(AbstractMesh_), POINTER :: cellMesh, boundaryMesh
+INTEGER(I4B), ALLOCATABLE :: nodeNum(:)
+REAL(DFP), ALLOCATABLE :: nodalValue(:, :)
+
+!----------------------------------------------------------------------------
+! Setting the verbosity
+! In Debug mode, there will many messages printed to the screen
+! Following code suppress information-exception.
+!----------------------------------------------------------------------------
+
+CALL e%setQuietMode(EXCEPTION_INFORMATION, .TRUE.)
+
+!----------------------------------------------------------------------------
+! Domain: Initiate FEDomain and print info
+!----------------------------------------------------------------------------
+
+CALL dom%ImportFromToml(tomlName="domain", filename=tomlFileName)
+! CALL dom%DisplayDomainInfo(msg="DomainInfo: ")
+
+!----------------------------------------------------------------------------
+! ScalarField and FEDOF: Initiate scalarField and fedof
+!----------------------------------------------------------------------------
+
+CALL u%ImportFromToml(tomlName="u", fedof=fedof, geofedof=geofedof, dom=dom, &
+                      filename=tomlFileName)
+! CALL fedof%Display(msg="FEDOF info: ")
+
+!----------------------------------------------------------------------------
+! DirichletBC: Import from toml
+!----------------------------------------------------------------------------
+
+CALL obj%ImportFromToml(filename=tomlFileName, dom=dom, tomlName="dbc")
+! CALL obj%Display("DirichletBC Info: ")
+
+!----------------------------------------------------------------------------
+! DirichletBC: Get the total node numbers
+!----------------------------------------------------------------------------
+tsize = obj%GetTotalNodeNum(fedof=fedof)
+CALL Display(tsize, "Total Node Num: ")
+
+!----------------------------------------------------------------------------
+! DirichletBC: Get the node numbers
+!----------------------------------------------------------------------------
+
+CALL Reallocate(nodeNum, tsize, isExpand=.TRUE., expandFactor=2)
+CALL obj%Get(fedof=fedof, nodeNum=nodeNum, tsize=tsize, &
+             iNodeOnNode=iNodeOnNode, iNodeOnEdge=iNodeOnEdge, &
+             iNodeOnFace=iNodeOnFace)
+CALL Display(tsize, "tsize = ")
+CALL Display(nodeNum(1:tsize), "nodeNum", full=.TRUE., orient="ROW")
+CALL Display(iNodeOnNode, "iNodeOnNode = ")
+CALL Display(iNodeOnFace, "iNodeOnFace = ")
+CALL Display(iNodeOnEdge, "iNodeOnEdge = ")
+
+!----------------------------------------------------------------------------
+! DirichletBC: Get the nodal values
+!----------------------------------------------------------------------------
+
+CALL Reallocate(nodalValue, tsize, 2, isExpand=.TRUE., expandFactor=2)
+CALL obj%Get(fedof=fedof, nodeNum=nodeNum, nodalValue=nodalValue, &
+             nrow=nrow, ncol=ncol)
+CALL Display(nodalValue(1:nrow, 1:ncol), "nodalValue: ", &
+             full=.TRUE.)
+CALL Display([nrow, ncol], "[nrow, ncol]: ")
+
+!----------------------------------------------------------------------------
+! Cleanup
+!----------------------------------------------------------------------------
+CALL dom%DEALLOCATE()
+END PROGRAM main

docs/blog/DirichletBC/code/test1.toml

@@ -0,0 +1,43 @@
+[domain]
+filename = "./mesh/square_3x3.h5"
+group = "/"
+totalMedium = 1
+
+[u]
+name = "u"
+engine = "NATIVE_SERIAL"
+fieldType = "Normal"
+spaceCompo = 1
+timeCompo = 1
+fedofName = "space"
+geofedofName = "geofedof"
+
+[u.geofedof]
+baseContinuity = "H1"
+baseInterpolation = "Lagrange"
+ipType = "Equidistance"
+baseType = "Monomial"
+
+[u.space]
+baseContinuity = "H1"
+baseInterpolation = "Hierarchical"
+# ipType = "Equidistance"
+# baseType = "Monomial"
+order = 4
+scaleForQuadOrder = 2
+quadOptName = "quadOpt"
+
+[u.space.quadOpt]
+isHomogeneous = true
+quadratureType = "GaussLegendre"
+# order = 2
+
+[dbc]
+name = "left_bottom_fix"
+idof = 1
+nodalValueType = "Constant"
+value = 10.0
+[dbc.boundary]
+isSelectionByMeshID = true
+[dbc.boundary.meshID]
+line = [1, 4]

docs/blog/DirichletBC/code/test2.F90

@@ -0,0 +1,111 @@
+!> author: Vikas Sharma, Ph. D.
+! date: 2025-10-23
+! summary: Study of DirichletBC class
+
+PROGRAM main
+USE GlobalData, ONLY: DFP, I4B, LGT
+USE ExceptionHandler_Class, ONLY: e, EXCEPTION_INFORMATION
+USE DirichletBC_Class
+USE FEDomain_Class
+USE AbstractMesh_Class
+USE FEDOF_Class
+USE ScalarField_Class
+USE Display_Method
+USE ReallocateUtility
+
+IMPLICIT NONE
+
+!----------------------------------------------------------------------------
+! Parameters
+!----------------------------------------------------------------------------
+CHARACTER(*), PARAMETER :: tomlFileName = "./test2.toml"
+
+!----------------------------------------------------------------------------
+! Types and variables
+!----------------------------------------------------------------------------
+TYPE(DirichletBC_) :: obj
+TYPE(FEDomain_) :: dom
+TYPE(ScalarField_) :: u
+TYPE(FEDOF_) :: fedof, geofedof
+INTEGER(I4B) :: tsize, iNodeOnNode, iNodeOnEdge, iNodeOnFace, nrow, ncol
+
+!----------------------------------------------------------------------------
+! Allocatables and pointers
+!----------------------------------------------------------------------------
+CLASS(AbstractMesh_), POINTER :: cellMesh, boundaryMesh
+INTEGER(I4B), ALLOCATABLE :: nodeNum(:)
+REAL(DFP), ALLOCATABLE :: nodalValue(:, :)
+
+!----------------------------------------------------------------------------
+! Setting the verbosity
+! In Debug mode, there will many messages printed to the screen
+! Following code suppress information-exception.
+!----------------------------------------------------------------------------
+
+! CALL e%setQuietMode(EXCEPTION_INFORMATION, .TRUE.)
+
+!----------------------------------------------------------------------------
+! Domain: Initiate FEDomain and print info
+!----------------------------------------------------------------------------
+
+CALL dom%ImportFromToml(tomlName="domain", filename=tomlFileName)
+! CALL dom%DisplayDomainInfo(msg="DomainInfo: ")
+
+!----------------------------------------------------------------------------
+! ScalarField and FEDOF: Initiate scalarField and fedof
+!----------------------------------------------------------------------------
+
+CALL u%ImportFromToml(tomlName="u", fedof=fedof, geofedof=geofedof, dom=dom, &
+                      filename=tomlFileName)
+! CALL fedof%Display(msg="FEDOF info: ")
+
+!----------------------------------------------------------------------------
+! DirichletBC: Import from toml
+!----------------------------------------------------------------------------
+
+CALL obj%ImportFromToml(filename=tomlFileName, dom=dom, tomlName="dbc")
+! CALL obj%Display("DirichletBC Info: ")
+
+!----------------------------------------------------------------------------
+! DirichletBC: Get the total node numbers
+!----------------------------------------------------------------------------
+tsize = obj%GetTotalNodeNum(fedof=fedof)
+CALL Display(tsize, "Total Node Num: ")
+
+!----------------------------------------------------------------------------
+! DirichletBC: Get the node numbers
+!----------------------------------------------------------------------------
+
+CALL Reallocate(nodeNum, tsize, isExpand=.TRUE., expandFactor=2)
+CALL obj%Get(fedof=fedof, nodeNum=nodeNum, tsize=tsize, &
+             iNodeOnNode=iNodeOnNode, iNodeOnEdge=iNodeOnEdge, &
+             iNodeOnFace=iNodeOnFace)
+CALL Display(tsize, "tsize = ")
+CALL Display(nodeNum(1:tsize), "nodeNum", full=.TRUE., orient="ROW")
+CALL Display(iNodeOnNode, "iNodeOnNode = ")
+CALL Display(iNodeOnFace, "iNodeOnFace = ")
+CALL Display(iNodeOnEdge, "iNodeOnEdge = ")
+
+!----------------------------------------------------------------------------
+! DirichletBC: Get the nodal values
+!----------------------------------------------------------------------------
+
+CALL Reallocate(nodalValue, tsize, 2, isExpand=.TRUE., expandFactor=2)
+CALL Display(SHAPE(nodalValue), "shape(nodalValue) before Get: ")
+nrow = 0; ncol = 0
+CALL obj%Get(fedof=fedof, geofedof=geofedof, nodeNum=nodeNum, &
+             nodalValue=nodalValue, nrow=nrow, ncol=ncol)
+
+CALL Display(nodeNum(1:nrow), "nodeNum(1:nrow): ", &
+             full=.TRUE., advance="NO")
+CALL Display(nodalValue(1:nrow, 1:ncol), "nodalValue: ", &
+             full=.TRUE.)
+CALL Display([nrow, ncol], "[nrow, ncol]: ")
+
+!----------------------------------------------------------------------------
+! Cleanup
+!----------------------------------------------------------------------------
+
+CALL dom%DEALLOCATE()
+
+END PROGRAM main

docs/blog/DirichletBC/code/test2.lua

@@ -0,0 +1,3 @@
+function BottomFunc(x, y, z, t)
+	return x ^ 2 + y
+end