diff --git a/math/Makefile b/math/Makefile
index 3845593a0c53..710586b5faca 100644
--- a/math/Makefile
+++ b/math/Makefile
@@ -476,6 +476,7 @@
     SUBDIR += mbasecalc
     SUBDIR += mcsim
     SUBDIR += mdal
+    SUBDIR += mesquite
     SUBDIR += metis
     SUBDIR += mfem
     SUBDIR += mingw32-libgmp
diff --git a/math/mesquite/Makefile b/math/mesquite/Makefile
new file mode 100644
index 000000000000..90a13dcce027
--- /dev/null
+++ b/math/mesquite/Makefile
@@ -0,0 +1,35 @@
+PORTNAME=	mesquite
+DISTVERSION=	2.3.0
+CATEGORIES=	math
+MASTER_SITES=	https://software.sandia.gov/mesquite/
+PKGNAMESUFFIX=	-mesh-toolkit
+
+MAINTAINER=	yuri@FreeBSD.org
+COMMENT=	Mesh Quality Improvement Toolkit is a suite of mesh algorithms
+WWW=		https://software.sandia.gov/mesquite
+
+LICENSE=	LGPL21
+LICENSE_FILE=	${WRKSRC}/lgpl.txt
+
+USES=		gmake libtool
+USE_LDCONFIG=	yes
+
+GNU_CONFIGURE=	yes
+CONFIGURE_ARGS=	--enable-shared --disable-static
+
+INSTALL_TARGET=	install-strip
+TEST_TARGET=	check
+
+OPTIONS_DEFINE=	EXAMPLES
+
+PORTEXAMPLES=	*
+
+post-install-EXAMPLES-on:
+	@${MKDIR} ${STAGEDIR}${EXAMPLESDIR}
+	@${MV} ${STAGEDIR}${DOCSDIR}/examples/* ${STAGEDIR}${EXAMPLESDIR}
+	@${RM} -r ${STAGEDIR}${DOCSDIR}
+
+post-install-EXAMPLES-off:
+	@${RM} -r ${STAGEDIR}${DOCSDIR}
+
+.include <bsd.port.mk>
diff --git a/math/mesquite/distinfo b/math/mesquite/distinfo
new file mode 100644
index 000000000000..eb607f23e1ff
--- /dev/null
+++ b/math/mesquite/distinfo
@@ -0,0 +1,3 @@
+TIMESTAMP = 1666841766
+SHA256 (mesquite-2.3.0.tar.gz) = 4ab4ceadfa596e16c00dbb0e8b830a9112fa1b73291ca07633ec379a39b8bb28
+SIZE (mesquite-2.3.0.tar.gz) = 12450450
diff --git a/math/mesquite/pkg-descr b/math/mesquite/pkg-descr
new file mode 100644
index 000000000000..9149e52ba841
--- /dev/null
+++ b/math/mesquite/pkg-descr
@@ -0,0 +1,7 @@
+Mesquite (Mesh Quality Improvement Toolkit) is designed to provide a
+stand-alone, portable, comprehensive suite of mesh quality improvement
+algorithms and components that can be used to construct custom quality
+improvement algorithms. Mesquite provides a robust and effective mesh
+improvement toolkit that allows both meshing researchers application
+scientists to benefit from the latest developments in mesh quality control
+and improvement.
diff --git a/math/mesquite/pkg-plist b/math/mesquite/pkg-plist
new file mode 100644
index 000000000000..f01ed9867c72
--- /dev/null
+++ b/math/mesquite/pkg-plist
@@ -0,0 +1,235 @@
+bin/msqquality
+bin/msqshape
+bin/vtkrandom
+bin/vtkxform
+include/AWMetric.hpp
+include/AWMetricBarrier.hpp
+include/AWMetricNonBarrier.hpp
+include/AWQualityMetric.hpp
+include/AWShape2DB1.hpp
+include/AWShape2DNB1.hpp
+include/AWShape2DNB2.hpp
+include/AWShapeOrientNB1.hpp
+include/AWShapeSizeB1.hpp
+include/AWShapeSizeOrientNB1.hpp
+include/AWSizeB1.hpp
+include/AWSizeNB1.hpp
+include/AWUntangleBeta.hpp
+include/AddQualityMetric.hpp
+include/AffineMapMetric.hpp
+include/ArrayMesh.hpp
+include/AspectRatioGammaQualityMetric.hpp
+include/AveragingQM.hpp
+include/Bits.hpp
+include/BoundedCylinderDomain.hpp
+include/CachingTargetCalculator.hpp
+include/CompareQM.hpp
+include/CompositeOFAdd.hpp
+include/CompositeOFMultiply.hpp
+include/CompositeOFScalarAdd.hpp
+include/CompositeOFScalarMultiply.hpp
+include/ConditionNumberFunctions.hpp
+include/ConditionNumberQualityMetric.hpp
+include/ConicDomain.hpp
+include/ConjugateGradient.hpp
+include/CurveDomain.hpp
+include/CylinderDomain.hpp
+include/DeformingDomainWrapper.hpp
+include/DomainClassifier.hpp
+include/EdgeIterator.hpp
+include/EdgeLengthMetric.hpp
+include/EdgeLengthQualityMetric.hpp
+include/EdgeLengthRangeQualityMetric.hpp
+include/EdgeQM.hpp
+include/ElemSampleQM.hpp
+include/ElementAvgQM.hpp
+include/ElementMaxQM.hpp
+include/ElementPMeanP.hpp
+include/ElementPatches.hpp
+include/ElementQM.hpp
+include/Exponent.hpp
+include/ExtraData.hpp
+include/ExtraDataUser.hpp
+include/FeasibleNewton.hpp
+include/FileTokenizer.hpp
+include/GlobalPatch.hpp
+include/HexLagrangeShape.hpp
+include/IQInterface.hpp
+include/IdealElements.hpp
+include/IdealShapeTarget.hpp
+include/IdealWeightInverseMeanRatio.hpp
+include/IdealWeightMeanRatio.hpp
+include/Instruction.hpp
+include/InstructionQueue.hpp
+include/InvTransBarrier.hpp
+include/InverseMetricWeight.hpp
+include/JacobianCalculator.hpp
+include/LInfTemplate.hpp
+include/LPtoPTemplate.hpp
+include/LVQDTargetCalculator.hpp
+include/LambdaConstant.hpp
+include/LambdaTarget.hpp
+include/LaplaceWrapper.hpp
+include/LaplacianSmoother.hpp
+include/LinearHexahedron.hpp
+include/LinearPrism.hpp
+include/LinearPyramid.hpp
+include/LinearQuadrilateral.hpp
+include/LinearTetrahedron.hpp
+include/LinearTriangle.hpp
+include/LocalSizeQualityMetric.hpp
+include/MappingFunction.hpp
+include/Matrix3D.hpp
+include/MaxTemplate.hpp
+include/MeanRatioFunctions.hpp
+include/MeshDecorator.hpp
+include/MeshDomain1D.hpp
+include/MeshImpl.hpp
+include/MeshImplData.hpp
+include/MeshImplTags.hpp
+include/MeshInterface.hpp
+include/MeshTransform.hpp
+include/MeshUtil.hpp
+include/MeshWriter.hpp
+include/Mesquite.hpp
+include/Mesquite_all_headers.hpp
+include/MetricWeight.hpp
+include/MsqDebug.hpp
+include/MsqError.hpp
+include/MsqFPE.hpp
+include/MsqFreeVertexIndexIterator.hpp
+include/MsqGeomPrim.hpp
+include/MsqHessian.hpp
+include/MsqIGeom.hpp
+include/MsqIMesh.hpp
+include/MsqIMeshP.hpp
+include/MsqIRel.hpp
+include/MsqInterrupt.hpp
+include/MsqMatrix.hpp
+include/MsqMeshEntity.hpp
+include/MsqTimer.hpp
+include/MsqVertex.hpp
+include/MultiplyQualityMetric.hpp
+include/NodeSet.hpp
+include/NonGradient.hpp
+include/NonSmoothDescent.hpp
+include/NullImprover.hpp
+include/NumericalQM.hpp
+include/OFEvaluator.hpp
+include/ObjectiveFunction.hpp
+include/ObjectiveFunctionTemplate.hpp
+include/PMeanPMetric.hpp
+include/PMeanPTemplate.hpp
+include/ParallelHelperInterface.hpp
+include/ParallelHelperNoMPIStub.hpp
+include/ParallelMeshImpl.hpp
+include/ParallelMeshInterface.hpp
+include/PatchData.hpp
+include/PatchIterator.hpp
+include/PatchPowerMeanP.hpp
+include/PatchSet.hpp
+include/PatchSetUser.hpp
+include/PaverMinEdgeLengthWrapper.hpp
+include/PlanarDomain.hpp
+include/PowerQualityMetric.hpp
+include/QuadLagrangeShape.hpp
+include/QualityAssessor.hpp
+include/QualityImprover.hpp
+include/QualityMetric.hpp
+include/QuasiNewton.hpp
+include/Randomize.hpp
+include/RefMeshTargetCalculator.hpp
+include/RefSizeTargetCalculator.hpp
+include/ReferenceMesh.hpp
+include/RelaxationSmoother.hpp
+include/RemainingWeight.hpp
+include/Sample.hpp
+include/ScalarAddQualityMetric.hpp
+include/ScalarMultiplyQualityMetric.hpp
+include/Settings.hpp
+include/ShapeImprovementWrapper.hpp
+include/ShapeImprover.hpp
+include/SimpleStats.hpp
+include/SizeAdaptShapeWrapper.hpp
+include/SizeMetric.hpp
+include/SlaveBoundaryVertices.hpp
+include/SmartLaplacianSmoother.hpp
+include/SphericalDomain.hpp
+include/StdDevTemplate.hpp
+include/SteepestDescent.hpp
+include/SymMatrix3D.hpp
+include/TInverseMeanRatio.hpp
+include/TMPQualityMetric.hpp
+include/TMetric.hpp
+include/TMetricBarrier.hpp
+include/TMetricNonBarrier.hpp
+include/TMixed.hpp
+include/TOffset.hpp
+include/TPower2.hpp
+include/TQualityMetric.hpp
+include/TScale.hpp
+include/TShape2DNB2.hpp
+include/TShape3DB2.hpp
+include/TShapeB1.hpp
+include/TShapeNB1.hpp
+include/TShapeOrientB1.hpp
+include/TShapeOrientB2.hpp
+include/TShapeOrientNB1.hpp
+include/TShapeOrientNB2.hpp
+include/TShapeSize2DB2.hpp
+include/TShapeSize2DNB1.hpp
+include/TShapeSize2DNB2.hpp
+include/TShapeSize3DB2.hpp
+include/TShapeSize3DB4.hpp
+include/TShapeSize3DNB1.hpp
+include/TShapeSizeB1.hpp
+include/TShapeSizeB3.hpp
+include/TShapeSizeNB3.hpp
+include/TShapeSizeOrientB1.hpp
+include/TShapeSizeOrientB2.hpp
+include/TShapeSizeOrientNB1.hpp
+include/TSizeB1.hpp
+include/TSizeNB1.hpp
+include/TSquared.hpp
+include/TSum.hpp
+include/TTau.hpp
+include/TUntangle1.hpp
+include/TUntangleBeta.hpp
+include/TUntangleMu.hpp
+include/TagVertexMesh.hpp
+include/TargetCalculator.hpp
+include/TargetReader.hpp
+include/TargetWriter.hpp
+include/TerminationCriterion.hpp
+include/TetDihedralWeight.hpp
+include/TetLagrangeShape.hpp
+include/TopologyInfo.hpp
+include/TriLagrangeShape.hpp
+include/TrustRegion.hpp
+include/UntangleBetaQualityMetric.hpp
+include/UntangleWrapper.hpp
+include/VarianceTemplate.hpp
+include/Vector3D.hpp
+include/VertexConditionNumberQualityMetric.hpp
+include/VertexMaxQM.hpp
+include/VertexMover.hpp
+include/VertexPMeanP.hpp
+include/VertexPatches.hpp
+include/VertexQM.hpp
+include/VertexSlaver.hpp
+include/ViscousCFDTetShapeWrapper.hpp
+include/VtkTypeInfo.hpp
+include/WeightCalculator.hpp
+include/WeightReader.hpp
+include/Wrapper.hpp
+include/XYPlanarDomain.hpp
+include/XYRectangle.hpp
+include/mesquite_config.h
+include/mesquite_version.h
+lib/libmesquite.so
+lib/libmesquite.so.0
+lib/libmesquite.so.0.0.0
+man/man1/msqquality.1.gz
+man/man1/msqshape.1.gz
+man/man1/vtkrandom.1.gz
+man/man1/vtkxform.1.gz