blob: 94f0db9c73b6270ebd6b053083ee1dbde3cf5959 (
plain)
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
|
*DECK DETPAR2
SUBROUTINE DETPAR2(V1,V2,V3,U1,U2,U3,AS,BS,CS)
*
*----------------------------------------------------------------------
*Purpose: routine de HQSIMEX
*
*Author(s):
* M. Beaudet
*
*Parameters:
* V1
* V2
* V3
* U1
* U2
* U3
* AS
* BS
* CS
*
*----------------------------------------------------------------------
*
CHARACTER*6 CLNAME
CLNAME = 'PAR'
ANUM = U1*(V2-V3)+U3*(V1-V2)+U2*(V3-V1)
ADEN = (V1-V2)*(V1-V3)*(V2-V3)
AS = ANUM/ADEN
BS = (U2-U3-AS*(V2*V2-V3*V3))/(V2-V3)
CS = U1-BS*V1-AS*V1*V1
RETURN
END
|
