			#include<stdio.h>
			#include<math.h>
			#include<fonction.h>
			#include<constant.h>
			#include<dimension.h>
			#include<structure.h>

			/****************************************************************/
			/* nom: e_amp */
			/* auteur: Jean-Paul Kneib */
			/* date: 10/02/92 */
			/* place: Toulouse *
			****************************************************************
			* Return the determinant of the amplification matrix ie the amplification
			* at the position position.
			*
			* - position is the position where to compute the amplification
			* - dl0s is the distance between lens[0] and zs
			* - dos is the normalised distance distcosmo1(zs)
			* - zs is the redshift of the source
			*
			* Global variables read :
			* - in e_grad2() : G, lens, lens_table
			*/
			double e_amp(const struct point *position, double dl0s, double dos, double zs)
			{
			struct matrix M;

			M = e_grad2(position, dl0s, zs);
			M.a /= dos;
			M.b /= dos;
			M.c /= dos;

			return(((1. - M.a)(1. - M.c) - M.bM.b));
			}

			double e_amp_gal(struct galaxie image, double np_b0)
			{
			struct matrix M;

			M = e_grad2_gal(image, np_b0);
			M.a /= image->dos;
			M.b /= image->dos;
			M.c /= image->dos;

			return(((1. - M.a)(1. - M.c) - M.bM.b));
			}

R1448/ChiBenchmark/lenstoolchi/src/e_amp.c
c99f2e455a51RC1

e_amp.c
View Options