mrdivide_helper.c

Go to the documentation of this file.

/*
 * Academic License - for use in teaching, academic research, and meeting
 * course requirements at degree granting institutions only.  Not for
 * government, commercial, or other organizational use.
 * File: mrdivide_helper.c
 *
 * MATLAB Coder version            : 24.2
 * C/C++ source code generated on  : 05-May-2025 19:53:20
 */
 
/* Include Files */
#include "mrdivide_helper.h"
#include <math.h>
 
/* Function Definitions */
/*
 * Arguments    : const double A[4]
 *                const double B[4]
 *                double Y[4]
 * Return Type  : void
 */
void mrdiv(const double A[4], const double B[4], double Y[4])
{
  double a21;
  double a22;
  double a22_tmp;
  int Y_tmp;
  int r1;
  int r2;
  if (fabs(B[1]) > fabs(B[0])) {
    r1 = 1;
    r2 = 0;
  } else {
    r1 = 0;
    r2 = 1;
  }
  a21 = B[r2] / B[r1];
  a22_tmp = B[r1 + 2];
  a22 = B[r2 + 2] - a21 * a22_tmp;
  Y_tmp = r1 << 1;
  Y[Y_tmp] = A[0] / B[r1];
  r2 <<= 1;
  Y[r2] = (A[2] - Y[Y_tmp] * a22_tmp) / a22;
  Y[Y_tmp] -= Y[r2] * a21;
  Y[Y_tmp + 1] = A[1] / B[r1];
  Y[r2 + 1] = (A[3] - Y[Y_tmp + 1] * a22_tmp) / a22;
  Y[Y_tmp + 1] -= Y[r2 + 1] * a21;
}
 
 
/* Function Definitions */
/*
 * Arguments    : const double B_data[]
 *                int B_size
 *                double Y_data[]
 *                int Y_size[2]
 * Return Type  : void
 */
void mrdiv_2(const double B_data[], int B_size, double Y_data[], int Y_size[2])
{
  if (B_size == 0) {
    Y_size[0] = 1;
    Y_size[1] = 0;
  } else {
    Y_size[0] = 1;
    Y_size[1] = 1;
    Y_data[0] = 1.0 / B_data[0];
  }
}
 
/*
 * File trailer for mrdivide_helper.c
 *
 * [EOF]
 */