summaryrefslogtreecommitdiff
path: root/support/dktools/f2lpara.c
blob: aa372f0df85707fc69e4b63f6c9125c61d6ccb85 (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
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
/*
Copyright (C) 2014-2020, Dirk Krause
SPDX-License-Identifier: BSD-3-Clause
*/

/*
	WARNING: This file was generated by the dkct program (see
	http://dktools.sourceforge.net/ for details).
	Changes you make here will be lost if dkct is run again!
	You should modify the original source and run dkct on it.
	Original source: f2lpara.ctr
*/

/**	@file f2lpara.c The f2lpara module.
*/


#line 54 "f2lpara.ctr"


#include "dk3ma.h"
#include "f2lpara.h"
#include "dk3unused.h"





#line 63 "f2lpara.ctr"


double
f2lpara_derived_left(double yneg, double y0, double ypos, int *ec)
{
  double		back	= 0.0;
  int			myec	= 0;
  /* back = 2 * y0 - 0.5 * ypos - 1.5 * yneg */
  back = dk3ma_d_sub_ok(
    dk3ma_d_mul_ok(2.0, y0, &myec),
    dk3ma_d_add_ok(
      (0.5 * ypos),
      dk3ma_d_mul_ok(1.5, yneg, &myec),
      &myec
    ),
    &myec
  );
  if (myec) { if (ec) { *ec = myec; } }
  return back;
}


double
f2lpara_derived_center(
  double yneg, double DK3_SILENCE_ARG_UNUSED(y0), double ypos, int *ec
)
{
  double		back	= 0.0;
  int			myec	= 0;
  /* back = 0.5 * (ypos - yneg)  */
  DK3_UNUSED_ARG(y0)
  back = 0.5 * dk3ma_d_sub_ok(ypos, yneg, &myec);
  if (myec) { if (ec) { *ec = myec; } }
  return back;
}


double
f2lpara_derived_right(double yneg, double y0, double ypos, int *ec)
{
  double		back	= 0.0;
  int			myec	= 0;
  /* back = 1.5 * ypos + 0.5 * yneg - 2 * y0 */
  back = dk3ma_d_sub_ok(
    dk3ma_d_add_ok(
      dk3ma_d_mul_ok(1.5, ypos, &myec),
      (0.5 * yneg),
      &myec
    ),
    dk3ma_d_mul_ok(2.0, y0, &myec),
    &myec
  );
  if (myec) { if (ec) { *ec = myec; } }
  return back;
}