← Back to branch summary

~ubuntu-branches/ubuntu/karmic/scilab/karmic

« back to all changes in this revision

Viewing changes to routines/optim/n1fc1.f

  • Committer: Bazaar Package Importer
  • Author(s): Torsten Werner
  • Date: 2002-03-21 16:57:43 UTC
  • Revision ID: james.westby@ubuntu.com-20020321165743-e9mv12c1tb1plztg
Tags: upstream-2.6
Import upstream version 2.6

Show diffs side-by-side

added added

removed removed

Lines of Context:
routines/optim/n1fc1.f
 
1
      subroutine n1fc1(simul,prosca,n,xn,fn,g,dxmin,df1,epsf,zero,imp,
 
2
     &                 io,mode,iter,nsim,memax,iz,rz,dz,izs,rzs,dzs)
 
3
C          dimension iz=2*(memax+1)
 
4
C          dimension rz=5*n+(n+4)*memax
 
5
C          dimension dz=(memax+9)*memax+8
 
6
c     Copyright INRIA
 
7
 
 
8
      implicit double precision (a-h,o-z)
 
9
      external simul, prosca
 
10
      dimension iz(*), rz(*), dz(*), xn(n), g(n), izs(*), dzs(*)
 
11
      real rzs(*)
 
12
      dimension i5(1), d3(1), d4(1)
 
13
C
 
14
      if (n.gt.0 .and. df1.gt.0.d0 .and. epsf.ge.0.d0 .and. zero.ge.0.d0
 
15
     &    .and. iter.ge.0 .and. nsim.ge.0 .and. memax.ge.1 .and.
 
16
     &    dxmin.gt.0.d0) goto 10
 
17
      mode = 2
 
18
C     appel incoherent
 
19
      call n1fc1o(io,1,i1,i2,i3,i4,i5,d1,d2,d3,d4)
 
20
      goto 999
 
21
 10   ns = 1
 
22
      ngd = ns + n
 
23
      nx = ngd + n
 
24
      nsa = nx + n
 
25
      ngg = nsa + n
 
26
      nal = ngg + n
 
27
      naps = nal + memax
 
28
      nanc = naps + memax
 
29
      npoids = nanc + memax
 
30
      nq = npoids + memax
 
31
      njc = 1
 
32
      nic = njc + memax + 1
 
33
      nr = 1
 
34
      na = nr + (memax+1)*(memax+1)
 
35
      ne = na + memax + 1
 
36
      nrr = ne + memax + 1
 
37
      nxga = nrr + memax + 1
 
38
      ny = nxga + memax + 1
 
39
      nw1 = ny + memax + 1
 
40
      nw2 = nw1 + memax + 1
 
41
C
 
42
      niz = 2 * (memax+1)
 
43
      nrz = nq + n*memax - 1
 
44
      ndz = nw2 + memax
 
45
      if (imp .gt. 0) call n1fc1o(io,2,n,memax,niz,nrz,ndz,d1,d2,d3,d4)
 
46
      do 110 i = 1,niz
 
47
 110  iz(i) = 0
 
48
      do 120 i = 1,nrz
 
49
 120  rz(i) = 0.d0
 
50
      do 130 i = 1,ndz
 
51
 130  dz(i) = 0.d0
 
52
      call n1fc1a(simul,prosca,n,mode,xn,fn,g,df1,epsf,dxmin,imp,zero,
 
53
     &            io,ntot,iter,nsim,memax,rz(ns),rz(ngd),rz(nx),rz(nsa),
 
54
     &            rz(ngg),rz(nal),rz(naps),rz(nanc),rz(npoids),rz(nq),
 
55
     &            iz(njc),iz(nic),dz(nr),dz(na),dz(ne),dz(nrr),dz(nxga),
 
56
     &            dz(ny),dz(nw1),dz(nw2),izs,rzs,dzs)
 
57
      iz(1) = ntot
 
58
 999  return
 
59
      end