1 /*
2 * Copyright (c) 2015, The University of Oxford
3 * All rights reserved.
4 *
5 * Redistribution and use in source and binary forms, with or without
6 * modification, are permitted provided that the following conditions are met:
7 * 1. Redistributions of source code must retain the above copyright notice,
8 * this list of conditions and the following disclaimer.
9 * 2. Redistributions in binary form must reproduce the above copyright notice,
10 * this list of conditions and the following disclaimer in the documentation
11 * and/or other materials provided with the distribution.
12 * 3. Neither the name of the University of Oxford nor the names of its
13 * contributors may be used to endorse or promote products derived from this
14 * software without specific prior written permission.
15 *
16 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
17 * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
18 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
19 * ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE
20 * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
21 * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
22 * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
23 * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
24 * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
25 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
26 * POSSIBILITY OF SUCH DAMAGE.
27 */
28
29 #ifndef OSKAR_PRIVATE_COND2_2X2_H_
30 #define OSKAR_PRIVATE_COND2_2X2_H_
31
32 #include <math.h>
33
34 #include "oskar_global.h"
35 #include "math/define_multiply.h"
36 #include "utility/oskar_vector_types.h"
37
38 #ifdef __cplusplus
39 extern "C" {
40 #endif
41
42 /* Private inline functions. */
43
44 OSKAR_INLINE
oskar_cond2_2x2_inline_f(const float4c * RESTRICT in)45 float oskar_cond2_2x2_inline_f(const float4c* RESTRICT in)
46 {
47 float sum, diff, t1, t2, a, b;
48 float4c p, q;
49 p = *in;
50 q = p;
51 OSKAR_MUL_COMPLEX_MATRIX_CONJUGATE_TRANSPOSE_IN_PLACE(float2, p, q)
52 sum = p.a.x + p.d.x;
53 t1 = p.a.x - p.d.x;
54 t1 = t1 * t1;
55 t2 = 4.0f * ((p.b.x * p.b.x) + (p.b.y * p.b.y));
56 diff = sqrtf(t1 + t2);
57 a = sqrtf(0.5f * (sum + diff));
58 b = sqrtf(0.5f * (sum - diff));
59 return (a > b) ? a / b : b / a;
60 }
61
62 OSKAR_INLINE
oskar_cond2_2x2_inline_d(const double4c * RESTRICT in)63 double oskar_cond2_2x2_inline_d(const double4c* RESTRICT in)
64 {
65 double sum, diff, t1, t2, a, b;
66 double4c p, q;
67 p = *in;
68 q = p;
69 OSKAR_MUL_COMPLEX_MATRIX_CONJUGATE_TRANSPOSE_IN_PLACE(double2, p, q)
70 sum = p.a.x + p.d.x;
71 t1 = p.a.x - p.d.x;
72 t1 = t1 * t1;
73 t2 = 4.0 * ((p.b.x * p.b.x) + (p.b.y * p.b.y));
74 diff = sqrt(t1 + t2);
75 a = sqrt(0.5 * (sum + diff));
76 b = sqrt(0.5 * (sum - diff));
77 return (a > b) ? a / b : b / a;
78 }
79
80 #ifdef __cplusplus
81 }
82 #endif
83
84 #endif /* OSKAR_PRIVATE_COND2_2X2_H_ */
85