1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
|
subroutine jacl2(neq,t,tq,ml,mu,pd,nrowpd)
c!but
c jacl2 cree la matrice jacobienne necessaire a Lsoda,
c qui correspond en fait a la hessienne du probleme
c d'approximation L2.
c!liste d'appel
c entree :
c - neq. tableau entier de taille 3+(nq+1)*(nq+2)
c neq(1)=nq est le degre effectif du polynome q
c neq(2)=ng est le nombre de coefficient de fourier
c neq(3)=dgmax degre maximum pour q (l'adresse des coeff de
c fourier dans tq est neq(3)+2
c neq(4:(nq+1)*(nq+2)) tableau de travail entier
c - t est une variable parametrique necessaire a Lsoda.
c - tq. tableau reel de taille au moins
c 7+dgmax+5*nq+6*ng+nq*ng+nq**2*(ng+1)
c tq(1:nq+1) est le tableau des coefficients du polynome q.
c tq(dgmax+2:dgmax+ng+2) est le tableau des coefficients
c de fourier
c tq(dgmax+ng+3:) est un tableau de travail de taille au moins
c 5+5*nq+5*ng+nq*ng+nq**2*(ng+1)
c - ml et mu sont les parametres du stockage par bande
c de la matrice qui n a pas lieu ici ,ils donc ignores.
c
c sortie :
c - pd. est le tableau ou l on range la matrice pleine
c dont les elements sont etablis par la sub. Hessien
c - nrowpd. est le nombre de ligne du tableau pd
c!
c Copyright INRIA
implicit double precision (a-h,o-y)
dimension tq(*),pd(nrowpd,*)
dimension neq(*)
c
call hessl2(neq,tq,pd,nrowpd)
nq=neq(1)
c write(6,'(''jac='')')
c do 10 i=1,nq
c write(6,'(5(e10.3,2x))') (pd(i,j),j=1,nq)
c 10 continue
c
return
end
subroutine jacl2n(neq,t,tq,ml,mu,pd,nrowpd)
c!but
c jacl2 cree la matrice jacobienne necessaire a Lsoda,
c qui correspond en fait a la hessienne du probleme
c d'approximation L2.
c!liste d'appel
c entree :
c - neq. tableau entier de taille 3+(nq+1)*(nq+2)
c neq(1)=nq est le degre effectif du polynome q
c neq(2)=ng est le nombre de coefficient de fourier
c neq(3)=dgmax degre maximum pour q (l'adresse des coeff de
c fourier dans tq est neq(3)+2
c neq(4:(nq+1)*(nq+2)) tableau de travail entier
c - t est une variable parametrique necessaire a Lsoda.
c - tq. tableau reel de taille au moins
c 7+dgmax+5*nq+6*ng+nq*ng+nq**2*(ng+1)
c tq(1:nq+1) est le tableau des coefficients du polynome q.
c tq(dgmax+2:dgmax+ng+2) est le tableau des coefficients
c de fourier
c tq(dgmax+ng+3:) est un tableau de travail de taille au moins
c 5+5*nq+5*ng+nq*ng+nq**2*(ng+1)
c - ml et mu sont les parametres du stockage par bande
c de la matrice qui n a pas lieu ici ,ils donc ignores.
c
c sortie :
c - pd. est le tableau ou l on range la matrice pleine
c dont les elements sont etablis par la sub. Hessien
c - nrowpd. est le nombre de ligne du tableau pd
c!
implicit double precision (a-h,o-y)
dimension tq(*),pd(nrowpd,*)
dimension neq(*)
c
call hessl2(neq,tq,pd,nrowpd)
nq=neq(1)
do 20 i=1,nq
do 10 j=1,nq
pd(i,j)=-pd(i,j)
10 continue
20 continue
c write(6,'(''jac='')')
c do 10 i=1,nq
c write(6,'(5(e10.3,2x))') (pd(i,j),j=1,nq)
c 10 continue
c
return
end
|
