/* nag_complex_create (a02bac) Example Program.
 *
 * Copyright 2023 Numerical Algorithms Group.
 *
 * Mark 29.1, 2023.
 */

#include <nag.h>
#include <stdio.h>

int main(void) {
  Integer exit_status = 0;
  Complex v, w, z;
  double r, theta, x, y;
  Nag_Boolean equal, not_equal;

  printf("nag_complex_create (a02bac) Example Program Results\n");

  x = 2.0;
  y = -3.0;
  /* nag_complex_create (a02bac).
   * Complex number from real and imaginary parts
   */
  z = nag_complex_create(x, y);

  printf("  %-21s  %s    %8s =  %7.4f, %7.4f\n", "", "", "x, y", x, y);
  printf("  %-21s: %s = %8s = (%7.4f, %7.4f)\n", "nag_complex_create", "z",
         "(x,y)", z.re, z.im);
  /* nag_complex_real (a02bbc).
   * Real part of a complex number
   */
  printf("  %-21s: %s    %8s =  %7.4f\n", "nag_complex_real", "", "real(z)",
         nag_complex_real(z));
  /* nag_complex_imag (a02bcc).
   * Imaginary part of a complex number
   */
  printf("  %-21s: %s    %8s =  %7.4f\n", "nag_complex_imag", "", "imag(z)",
         nag_complex_imag(z));
  /* nag_complex_create (a02bac), see above. */
  v = nag_complex_create(3.0, 1.25);
  /* nag_complex_create (a02bac), see above. */
  w = nag_complex_create(2.5, -1.75);
  printf("  %-21s: %s    %8s = (%7.4f, %7.4f)\n", "nag_complex_create", "", "v",
         v.re, v.im);
  printf("  %-21s: %s    %8s = (%7.4f, %7.4f)\n", "nag_complex_create", "", "w",
         w.re, w.im);
  /* nag_complex_add (a02cac).
   * Addition of two complex numbers
   */
  z = nag_complex_add(v, w);
  printf("  %-21s: %s = %8s = (%7.4f, %7.4f)\n", "nag_complex_add", "z", "v+w",
         z.re, z.im);
  /* nag_complex_subtract (a02cbc).
   * Subtraction of two complex numbers
   */
  z = nag_complex_subtract(v, w);
  printf("  %-21s: %s = %8s = (%7.4f, %7.4f)\n", "nag_complex_subtract", "z",
         "v-w", z.re, z.im);
  /* nag_complex_multiply (a02ccc).
   * Multiplication of two complex numbers
   */
  z = nag_complex_multiply(v, w);
  printf("  %-21s: %s = %8s = (%7.4f, %7.4f)\n", "nag_complex_multiply", "z",
         "v*w", z.re, z.im);
  /* nag_complex_divide (a02cdc).
   * Quotient of two complex numbers
   */
  z = nag_complex_divide(v, w);
  printf("  %-21s: %s = %8s = (%7.4f, %7.4f)\n", "nag_complex_divide", "z",
         "v/w", z.re, z.im);
  /* nag_complex_negate (a02cec).
   * Negation of a complex number
   */
  z = nag_complex_negate(w);
  printf("  %-21s: %s = %8s = (%7.4f, %7.4f)\n", "nag_complex_negate", "z",
         "-w", z.re, z.im);
  /* nag_complex_conjg (a02cfc).
   * Conjugate of a complex number
   */
  z = nag_complex_conjg(w);
  printf("  %-21s: %s = %8s = (%7.4f, %7.4f)\n", "nag_complex_conjg", "z",
         "conjg(w)", z.re, z.im);
  /* nag_complex_equal (a02cgc).
   * Equality of two complex numbers
   */
  equal = nag_complex_equal(v, w);
  if (equal)
    printf("  %-21s: %s == %s\n", "nag_complex_equal", "v", "w");
  else
    printf("  %-21s: %s != %s\n", "nag_complex_equal", "v", "w");
  /* nag_complex_not_equal (a02chc).
   * Inequality of two complex numbers
   */
  not_equal = nag_complex_not_equal(w, z);
  if (not_equal)
    printf("  %-21s: %s != %s\n\n", "nag_complex_not_equal", "w", "z");
  else
    printf("  %-21s: %s == %s\n\n", "nag_complex_not_equal", "w", "z");

  /* nag_complex_arg (a02dac).
   * Argument of a complex number
   */
  theta = nag_complex_arg(z);
  printf("  %-21s: %s    %8s =  %7.4f\n", "nag_complex_arg", "", "arg(z)",
         theta);
  /* nag_complex_abs (a02dbc).
   * Modulus of a complex number
   */
  r = nag_complex_abs(z);
  printf("  %-21s: %s = %8s =  %7.4f\n", "nag_complex_abs", "r", "abs(z)", r);
  /* nag_complex_sqrt (a02dcc).
   * Square root of a complex number
   */
  v = nag_complex_sqrt(z);
  printf("  %-21s: %s = %8s = (%7.4f, %7.4f)\n", "nag_complex_sqrt", "v",
         "sqrt(z)", v.re, v.im);
  /* nag_complex_i_power (a02ddc).
   * Complex number raised to integer power
   */
  v = nag_complex_i_power(z, (Integer)3);
  printf("  %-21s: %s = %8s = (%7.4f, %7.4f)\n", "nag_complex_i_power", "v",
         "z**3", v.re, v.im);
  /* nag_complex_r_power (a02dec).
   * Complex number raised to real power
   */
  v = nag_complex_r_power(z, 2.5);
  printf("  %-21s: %s = %8s = (%7.4f, %7.4f)\n", "nag_complex_r_power", "v",
         "z**2.5", v.re, v.im);
  /* nag_complex_c_power (a02dfc).
   * Complex number raised to complex power
   */
  v = nag_complex_c_power(z, w);
  printf("  %-21s: %s = %8s = (%7.4f,%8.4f)\n", "nag_complex_c_power", "v",
         "z**w", v.re, v.im);
  /* nag_complex_log (a02dgc).
   * Complex logarithm
   */
  v = nag_complex_log(z);
  printf("  %-21s: %s = %8s = (%7.4f, %7.4f)\n", "nag_complex_log", "v",
         "log(z)", v.re, v.im);
  /* nag_complex_exp (a02dhc).
   * Complex exponential
   */
  z = nag_complex_exp(v);
  printf("  %-21s: %s = %8s = (%7.4f, %7.4f)\n", "nag_complex_exp", "z",
         "exp(v)", z.re, z.im);
  /* nag_complex_sin (a02djc).
   * Complex sine
   */
  v = nag_complex_sin(z);
  printf("  %-21s: %s = %8s = (%7.4f, %7.4f)\n", "nag_complex_sin", "v",
         "sin(z)", v.re, v.im);
  /* nag_complex_cos (a02dkc).
   * Complex cosine
   */
  v = nag_complex_cos(z);
  printf("  %-21s: %s = %8s = (%7.4f, %7.4f)\n", "nag_complex_cos", "v",
         "cos(z)", v.re, v.im);
  /* nag_complex_tan (a02dlc).
   * Complex tangent
   */
  v = nag_complex_tan(z);
  printf("  %-21s: %s = %8s = (%7.4f, %7.4f)\n", "nag_complex_tan", "v",
         "tan(z)", v.re, v.im);
  /* nag_complex_divide (a02cdc), see above. */
  v = nag_complex_divide(nag_complex_sin(z), nag_complex_cos(z));
  printf("  %-21s:%13s = (%7.4f, %7.4f)\n", "nag_complex_divide",
         "sin(z)/cos(z)", v.re, v.im);

  return exit_status;
}