source: ThirdParty/FilterSQP/trunk/filtersqp.patch

Last change on this file was 1554, checked in by andreasw, 10 years ago

merging patch correction from stable/1.0 into trunk

File size: 1.4 KB
  • src/bqpd.f

    diff -Naur source-orig/src/bqpd.f source/src/bqpd.f
    old new  
    434434        if(lcyc(i).ne.lcyc(i+ncyc))goto16
    435435      enddo
    436436      gtol=1.D1*gtol
     437      aaaa = 1.D0
     438      bbbb = 0.D0
     439      if (gtol .eq. aaaa/bbbb) then
     440        ifail = 8
     441        return
     442      endif
    437443   16 continue
    438444c     if(info(1).ge.117)iprint=3
    439445c     if(info(1).ge.90)stop
  • src/filter.f

    diff -Naur source-orig/src/filter.f source/src/filter.f
    old new  
    447447c     ... ensure that the initial point is feasible wrt linear c/s
    448448      call linear_feasible (n, m, x, c, blo, bup, cstype, lin_feas)
    449449      if (.not. lin_feas) then
     450         do i=1,n
     451            d(i) = 0.d0
     452         enddo
    450453         call objfun1(x,n,f,s,user,iuser,flag)
    451454         hc = h_one (m, c, blo(n1), bup(n1))
    452455         phi = f + hc
  • src/sparseL.f

    diff -Naur source-orig/src/sparseL.f source/src/sparseL.f
    old new  
    17011701      common/epsc/eps,tol,emin
    17021702      common/noutc/nout
    17031703      parameter (thresh=1.D-1,growth=1.D1)
     1704      qr_ = 0
    17041705c     write(nout,*)'update_L:  p,q =',pp,qq
    17051706      nup=nup+1
    17061707      if(qq.gt.n)then
Note: See TracBrowser for help on using the repository browser.