← Back to branch summary

~ubuntu-branches/ubuntu/saucy/python-scipy/saucy

« back to all changes in this revision

Viewing changes to Lib/integrate/odepack/adjlr.f

  • Committer: Bazaar Package Importer
  • Author(s): Ondrej Certik
  • Date: 2008-06-16 22:58:01 UTC
  • mfrom: (2.1.24 intrepid)
  • Revision ID: james.westby@ubuntu.com-20080616225801-irdhrpcwiocfbcmt
Tags: 0.6.0-12
http://bugs.debian.org/489149
* The description updated to match the current SciPy (Closes: #489149).
* Standards-Version bumped to 3.8.0 (no action needed)
* Build-Depends: netcdf-dev changed to libnetcdf-dev

Show diffs side-by-side

added added

removed removed

Lines of Context:
Lib/integrate/odepack/adjlr.f
1
 
      subroutine adjlr (n, isp, ldif)
2
 
      integer n, isp, ldif
3
 
      dimension isp(1)
4
 
c-----------------------------------------------------------------------
5
 
c this routine computes an adjustment, ldif, to the required
6
 
c integer storage space in iwk (sparse matrix work space).
7
 
c it is called only if the word length ratio is lrat = 1.
8
 
c this is to account for the possibility that the symbolic lu phase
9
 
c may require more storage than the numerical lu and solution phases.
10
 
c-----------------------------------------------------------------------
11
 
      integer ip, jlmax, jumax, lnfc, lsfc, nzlu
12
 
c
13
 
      ip = 2*n + 1
14
 
c get jlmax = ijl(n) and jumax = iju(n) (sizes of jl and ju). ----------
15
 
      jlmax = isp(ip)
16
 
      jumax = isp(ip+ip)
17
 
c nzlu = (size of l) + (size of u) = (il(n+1)-il(1)) + (iu(n+1)-iu(1)).
18
 
      nzlu = isp(n+1) - isp(1) + isp(ip+n+1) - isp(ip+1)
19
 
      lsfc = 12*n + 3 + 2*max0(jlmax,jumax)
20
 
      lnfc = 9*n + 2 + jlmax + jumax + nzlu
21
 
      ldif = max0(0, lsfc - lnfc)
22
 
      return
23
 
c----------------------- end of subroutine adjlr -----------------------
24
 
      end