Skip to content

Commit 42cdeb8

Browse files
committed
Fix: correct interface for sb10ad.
The fix is a bit of a guess, but after this the signature matches that expected by synthesis.py's use of _wrapper.sb10ad, and a simple test case works.
1 parent 76950ba commit 42cdeb8

File tree

1 file changed

+4
-4
lines changed

1 file changed

+4
-4
lines changed

slycot/src/synthesis.pyf

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -410,13 +410,13 @@ subroutine sb10ad(job,n,m,np,ncon,nmeas,gamma,a,lda,b,ldb,c,ldc,d,ldd,ak,ldak,bk
410410
double precision intent(out), dimension(ncon,nmeas), depend(ncon,nmeas) :: dk
411411
integer intent(hide), depend(dk) :: lddk=shape(dk,0)
412412
double precision intent(out), dimension(2*n,2*n), depend(n) :: ac
413-
integer intent(hide), depend(ac) :: ldak=shape(ac,0)
413+
integer intent(hide), depend(ac) :: ldac=shape(ac,0)
414414
double precision intent(out), dimension(2*n,m-ncon), depend(n,m,ncon) :: bc
415-
integer intent(hide), depend(bc) :: ldbk=shape(bc,0)
415+
integer intent(hide), depend(bc) :: ldbc=shape(bc,0)
416416
double precision intent(out), dimension(np-nmeas,2*n), depend(np,nmeas,n) :: cc
417-
integer intent(hide), depend(cc) :: ldck=shape(cc,0)
417+
integer intent(hide), depend(cc) :: ldcc=shape(cc,0)
418418
double precision intent(out), dimension(np-nmeas,m-ncon), depend(np,nmeas,m,ncon) :: dc
419-
integer intent(hide), depend(dc) :: lddk=shape(dc,0)
419+
integer intent(hide), depend(dc) :: lddc=shape(dc,0)
420420
double precision intent(out), dimension(4) :: rcond
421421
double precision optional :: gtol=0.0
422422
double precision optional :: actol=0.0

0 commit comments

Comments
 (0)