diff --git a/science/getdp/Makefile b/science/getdp/Makefile
index 17ce626ea0e4..ea2ab883a61d 100644
--- a/science/getdp/Makefile
+++ b/science/getdp/Makefile
@@ -1,32 +1,31 @@
 PORTNAME=	getdp
-PORTVERSION=	3.5.0
-PORTREVISION=	2
+PORTVERSION=	3.6.0
 CATEGORIES=	science
 MASTER_SITES=	http://www.geuz.org/getdp/src/
 DISTNAME=	${PORTNAME}-${PORTVERSION}-source
 
 MAINTAINER=	stephen@FreeBSD.org
 COMMENT=	Rather general finite element solver using mixed finite elements
 WWW=		https://www.geuz.org/getdp/
 
 LICENSE=	GPLv2+
 LICENSE_FILE=	${WRKSRC}/LICENSE.txt
 
 LIB_DEPENDS=	libgsl.so:math/gsl	\
 		liblapack.so:math/lapack
 
 USES=		fortran tar:tgz cmake
 
 OPTIONS_DEFINE=	DOCS EXAMPLES
 
 .include <bsd.port.pre.mk>
 
 .if ${GCC_DEFAULT} >= 10
 # workaround for Rank mismatch in argument 'lines' at (1) (rank-1 and scalar)
 FFLAGS+=	-fallow-argument-mismatch
 .endif
 
 do-test:
 	cd ${WRKSRC}/demos && yes "" | ${STAGEDIR}${PREFIX}/bin/${PORTNAME} magnet.pro -solve MagSta_phi
 
 .include <bsd.port.post.mk>
diff --git a/science/getdp/distinfo b/science/getdp/distinfo
index 4f4d49029d3f..21133c0a3622 100644
--- a/science/getdp/distinfo
+++ b/science/getdp/distinfo
@@ -1,3 +1,3 @@
-TIMESTAMP = 1652456807
-SHA256 (getdp-3.5.0-source.tgz) = 0bf76c49efa9b88401a457c12f7aabd97597ad49ef632d3fb9308aa8ea437bdc
-SIZE (getdp-3.5.0-source.tgz) = 1901571
+TIMESTAMP = 1709878592
+SHA256 (getdp-3.6.0-source.tgz) = 37571bef65bbec3afe4f3787a4b82fe45b856fe489ab219c7e7fa439f12f6600
+SIZE (getdp-3.6.0-source.tgz) = 1900780
diff --git a/science/getdp/files/patch-src_kernel_LinAlg__SPARSKIT.cpp b/science/getdp/files/patch-src_kernel_LinAlg__SPARSKIT.cpp
new file mode 100644
index 000000000000..cf21fca26281
--- /dev/null
+++ b/science/getdp/files/patch-src_kernel_LinAlg__SPARSKIT.cpp
@@ -0,0 +1,11 @@
+--- src/kernel/LinAlg_SPARSKIT.cpp.orig	2024-03-08 06:48:58 UTC
++++ src/kernel/LinAlg_SPARSKIT.cpp
+@@ -117,7 +117,7 @@ void LinAlg_CreateVector(gVector *V, gSolver *Solver, 
+   V->N = n;
+ }
+ 
+-void LinAlg_CreateMatrix(gMatrix *M, gSolver *Solver, int n, int m)
++void LinAlg_CreateMatrix(gMatrix *M, gSolver *Solver, int n, int m, bool silent)
+ {
+   init_matrix(n, &M->M, &Solver->Params);
+ }