From 4f9511516904395d34e478ef3e7aec16ec7c4542 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Lisa=20Sch=C3=A4tzle?= <lisa.schaetzle@kit.edu> Date: Thu, 11 Jan 2024 08:41:22 +0000 Subject: [PATCH] Update evaluateFarfieldSecondOrder.m --- src/evaluateFarfieldSecondOrder.m | 33 +++++++++++++++++++------------ 1 file changed, 20 insertions(+), 13 deletions(-) diff --git a/src/evaluateFarfieldSecondOrder.m b/src/evaluateFarfieldSecondOrder.m index 62b0d5a..848fe3b 100644 --- a/src/evaluateFarfieldSecondOrder.m +++ b/src/evaluateFarfieldSecondOrder.m @@ -1,16 +1,23 @@ -function [F] = evaluateFarfieldSecondOrder(k, sampling, objects, par, R,q) -% Calculates a discretized version of the Born farfield operator by using -% the composite trapezoidal rule for the outer integral and calculations as suggested -% by Vainikko in order to approximate the individual Born farfields. -% INPUT: kappa - wave number -% sampling - cell array that contains all information about the -% discretization of the two spheres -% (detectors, illumination directions) -% object - String, that determines the shape of the refractive -% index n=1-q -% par - Vector, that contains all parameters for calculating q -% R,z - Ass.: supp(q) is contained in B_R(z) -% OUTPUT: F - Born farfield operator, nxhat*nd array +function [F] = evaluateFarfieldSecondOrder(k, sampling, objects, par, R, q) + +% Calculates a discretized version of the Born farfield operator of second order +% by using the composite trapezoidal rule for the outer integral and calculations +% as suggested by Vainikko in order to approximate the individual Born farfields. +% +% INPUT: kappa Wave number, >0. +% sampling Structure containing information about the discretization. +% objects Structure including two strings that set the shape-types of +% the two scatterers. +% par Vector containing all further information about the scatterers +% shape, 2*6-array. +% R Sizes of the two scatterers, i.e. radii of balls containing them, +% vector of 2. +% q Values of the contrast function within the individual scatterers, +% >-1. +% +% OUTPUT: F Born farfield operator of second order, nxhat*nd-array. +% +% ************************************************************************************ kR = k * R; % kz = k * z; -- GitLab