File: parallel_part.c

package info (click to toggle)
garlic 1.6-3
  • links: PTS, VCS
  • area: main
  • in suites: bullseye, buster, sid
  • size: 4,516 kB
  • sloc: ansic: 52,465; makefile: 2,254
file content (77 lines) | stat: -rw-r--r-- 2,085 bytes parent folder | download | duplicates (5)
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
/* Copyright (C) 2001 Damir Zucic */

/*=============================================================================

				parallel_part.c

Purpose:
	Find the part of the vector which is parallel with the reference
	vector.

Input:
	(1) Pointer to VectorS (where the result will be stored).
	(2) Pointer to VectorS structure (the reference vector).
	(3) Pointer to VectorS structure (the vector which should be
	    projected).

Output:
	(1) Parallel part of the specified vector.
	(2) Return value.

Return value:
	(1) Positive on success.
	(2) Negative on failure.

========includes:============================================================*/

#include <stdio.h>

#include <math.h>

#include <X11/Xlib.h>
#include <X11/Xutil.h>
#include <X11/Xos.h>
#include <X11/Xatom.h>

#include "defines.h"
#include "typedefs.h"

/*======function prototypes:=================================================*/

double		ScalarProduct_ (VectorS *, VectorS *);

/*======find parallel part of input vector:==================================*/

int ParallelPart_ (VectorS *parallel_vectorSP,
		   VectorS *reference_vectorSP, VectorS *input_vectorSP)
{
double		abs_value_squared;
double		reciprocal_denominator;
double		scalar_product;
double		scale_factor;

/* Check the squared absolute value of the reference vector: */
abs_value_squared = ScalarProduct_ (reference_vectorSP, reference_vectorSP);
if (abs_value_squared == 0.0) return -1;

/* Prepare the reciprocal denominator: */
reciprocal_denominator = 1.0 / abs_value_squared;

/* Prepare the scalar product of two vectors: */
scalar_product = ScalarProduct_ (input_vectorSP, reference_vectorSP);

/* Prepare the scale factor: */
scale_factor = scalar_product * reciprocal_denominator;

/* Prepare the result: */
parallel_vectorSP->x = scale_factor * reference_vectorSP->x;
parallel_vectorSP->y = scale_factor * reference_vectorSP->y;
parallel_vectorSP->z = scale_factor * reference_vectorSP->z;

/* Return positive value (success indicator): */
return 1;
}

/*===========================================================================*/