gpamela/gpfield/flin3.F

*
* $Id$
*
* $Log$
*
*CMZ :  2.01/00 05/04/2000  14.35.18  by  Marialuigia Ambriola
*CMZ :  2.00/00 03/03/2000  15.39.05  by  Francesco Cafagna
*CMZ :  1.02/00 21/10/99  16.14.45  by  Francesco Cafagna
*CMZ :  1.00/02 25/03/96  15.18.51  by  Francesco Cafagna
*-- Author :    Paolo Papini   16/02/96
      SUBROUTINE FLIN3(V0,V1,V2,V3,F0,F1,F2,F3,V,FF)
************************************************************************
*                                                                      *
*                                 *
* A bit of housekeeping: delete print statement etc. etc., by F. Cafagna
* Called by: GUFLD                                                     *
* Author: Paolo Papini   16/02/96                                      *
*                                                                      *
************************************************************************
#include "gcunit.inc"
*
      REAL V0(3),V1(3),V2(3),V3(3),V(3),FF
      REAL F0,F1,F2,F3,DX,DY,DZ,X1,X2,X3,Y1,Y2,Y3,Z1,Z2,Z3,DET,
     +     S1,S3,S4
*
      F1=F1-F0
      F2=F2-F0
      F3=F3-F0
      DX=V(1)-V0(1)
      DY=V(2)-V0(2)
      DZ=V(3)-V0(3)
      X1=V1(1)-V0(1)
      X2=V2(1)-V0(1)
      X3=V3(1)-V0(1)
      Y1=V1(2)-V0(2)
      Y2=V2(2)-V0(2)
      Y3=V3(2)-V0(2)
      Z1=V1(3)-V0(3)
      Z2=V2(3)-V0(3)
      Z3=V3(3)-V0(3)
      DET=X1*Y2*Z3+Y1*Z2*X3+Z1*X2*Y3-Z1*Y2*X3-Z2*Y3*X1-Z3*Y1*X2
      S1=F1*(Y2*Z3-Z2*Y3)+F2*(Z1*Y3-Z3*Y1)+F3*(Y1*Z2-Z1*Y2)
      S2=F1*(Z2*X3-Z3*X2)+F2*(X1*Z3-Z1*X3)+F3*(Z1*X2-Z2*X1)
      S3=F1*(X2*Y3-Y2*X3)+F2*(Y1*X3-Y3*X1)+F3*(X1*Y2-Y1*X2)
      IF(DET.EQ.0.) THEN
*
* GMAIL routine used
*
         CHMAIL=' flin3 ERROR: MATRIX DETERMINANT EQUAL ZERO'
         CALL GMAIL(1,0)
         FF=F0
         RETURN
      END IF
      FF=F0+(S1*DX+S2*DY+S3*DZ)/DET
      RETURN
      END
1	*
2	* $Id$
3	*
4	* $Log$
5	*
6	*CMZ : 2.01/00 05/04/2000 14.35.18 by Marialuigia Ambriola
7	*CMZ : 2.00/00 03/03/2000 15.39.05 by Francesco Cafagna
8	*CMZ : 1.02/00 21/10/99 16.14.45 by Francesco Cafagna
9	*CMZ : 1.00/02 25/03/96 15.18.51 by Francesco Cafagna
10	*-- Author : Paolo Papini 16/02/96
11	SUBROUTINE FLIN3(V0,V1,V2,V3,F0,F1,F2,F3,V,FF)
12	************************************************************************
13	* *
14	* *
15	* A bit of housekeeping: delete print statement etc. etc., by F. Cafagna
16	* Called by: GUFLD *
17	* Author: Paolo Papini 16/02/96 *
18	* *
19	************************************************************************
20	#include "gcunit.inc"
21	*
22	REAL V0(3),V1(3),V2(3),V3(3),V(3),FF
23	REAL F0,F1,F2,F3,DX,DY,DZ,X1,X2,X3,Y1,Y2,Y3,Z1,Z2,Z3,DET,
24	+ S1,S3,S4
25	*
26	F1=F1-F0
27	F2=F2-F0
28	F3=F3-F0
29	DX=V(1)-V0(1)
30	DY=V(2)-V0(2)
31	DZ=V(3)-V0(3)
32	X1=V1(1)-V0(1)
33	X2=V2(1)-V0(1)
34	X3=V3(1)-V0(1)
35	Y1=V1(2)-V0(2)
36	Y2=V2(2)-V0(2)
37	Y3=V3(2)-V0(2)
38	Z1=V1(3)-V0(3)
39	Z2=V2(3)-V0(3)
40	Z3=V3(3)-V0(3)
41	DET=X1Y2Z3+Y1Z2X3+Z1X2Y3-Z1Y2X3-Z2Y3X1-Z3Y1X2
42	S1=F1(Y2Z3-Z2Y3)+F2(Z1Y3-Z3Y1)+F3(Y1Z2-Z1*Y2)
43	S2=F1(Z2X3-Z3X2)+F2(X1Z3-Z1X3)+F3(Z1X2-Z2*X1)
44	S3=F1(X2Y3-Y2X3)+F2(Y1X3-Y3X1)+F3(X1Y2-Y1*X2)
45	IF(DET.EQ.0.) THEN
46	*
47	* GMAIL routine used
48	*
49	CHMAIL=' flin3 ERROR: MATRIX DETERMINANT EQUAL ZERO'
50	CALL GMAIL(1,0)
51	FF=F0
52	RETURN
53	END IF
54	FF=F0+(S1DX+S2DY+S3*DZ)/DET
55	RETURN
56	END