FLA_Apply_pivots_rn.h File Reference

(r)

Go to the source code of this file.

Functions
FLA_Error	FLA_Apply_pivots_rn_opt_var1 (FLA_Obj p, FLA_Obj A)
FLA_Error	FLA_Apply_pivots_rn_ops_var1 (int n, float *a, int a_rs, int a_cs, int k1, int k2, int *p, int incp)
FLA_Error	FLA_Apply_pivots_rn_opd_var1 (int n, double *a, int a_rs, int a_cs, int k1, int k2, int *p, int incp)
FLA_Error	FLA_Apply_pivots_rn_opc_var1 (int n, scomplex *a, int a_rs, int a_cs, int k1, int k2, int *p, int incp)
FLA_Error	FLA_Apply_pivots_rn_opz_var1 (int n, dcomplex *a, int a_rs, int a_cs, int k1, int k2, int *p, int incp)

---

Function Documentation

FLA_Error FLA_Apply_pivots_rn_opc_var1	(	int	n,
		scomplex *	a,
		int	a_rs,
		int	a_cs,
		int	k1,
		int	k2,
		int *	p,
		int	incp
	)

FLA_Error FLA_Apply_pivots_rn_opd_var1	(	int	n,
		double *	a,
		int	a_rs,
		int	a_cs,
		int	k1,
		int	k2,
		int *	p,
		int	incp
	)

FLA_Error FLA_Apply_pivots_rn_ops_var1	(	int	n,
		float *	a,
		int	a_rs,
		int	a_cs,
		int	k1,
		int	k2,
		int *	p,
		int	incp
	)

FLA_Error FLA_Apply_pivots_rn_opt_var1	(	FLA_Obj	p,
		FLA_Obj	A
	)

References FLA_Apply_pivots_ln_opc_var1(), FLA_Apply_pivots_ln_opd_var1(), FLA_Apply_pivots_ln_opi_var1(), FLA_Apply_pivots_ln_ops_var1(), FLA_Apply_pivots_ln_opz_var1(), FLA_Obj_col_stride(), FLA_Obj_datatype(), FLA_Obj_length(), FLA_Obj_row_stride(), FLA_Obj_vector_dim(), and FLA_Obj_vector_inc().

Referenced by FLA_Apply_pivots_rn().

{
  FLA_Datatype datatype;
  int          m_A;
  int          rs_A, cs_A;
  int          inc_p;
  int          k1_0, k2_0;

  datatype = FLA_Obj_datatype( A );

  m_A      = FLA_Obj_length( A );

  // Swap the stride; FLA_Apply_pivots_ln_ops_var1 already consider the memory access pattern.
  cs_A     = FLA_Obj_row_stride( A );
  rs_A     = FLA_Obj_col_stride( A );

  inc_p    = -FLA_Obj_vector_inc( p );

  // Use zero-based indices.
  k1_0     = 0;
  k2_0     = ( int ) FLA_Obj_vector_dim( p ) - 1;

  switch ( datatype )
  {
    case FLA_INT:
    {
      int*   buff_A = FLA_INT_PTR( A );
      int*   buff_p = FLA_INT_PTR( p );

      FLA_Apply_pivots_ln_opi_var1( m_A,
                                    buff_A, rs_A, cs_A,
                                    k1_0,
                                    k2_0,
                                    buff_p, inc_p );

      break;
    }

    case FLA_FLOAT:
    {
      float* buff_A = FLA_FLOAT_PTR( A );
      int*   buff_p = FLA_INT_PTR( p );

      FLA_Apply_pivots_ln_ops_var1( m_A,
                                    buff_A, rs_A, cs_A,
                                    k1_0,
                                    k2_0,
                                    buff_p, inc_p );

      break;
    }

    case FLA_DOUBLE:
    {
      double* buff_A = FLA_DOUBLE_PTR( A );
      int*    buff_p = FLA_INT_PTR( p );

      FLA_Apply_pivots_ln_opd_var1( m_A,
                                    buff_A, rs_A, cs_A,
                                    k1_0,
                                    k2_0,
                                    buff_p, inc_p );

      break;
    }

    case FLA_COMPLEX:
    {
      scomplex* buff_A = FLA_COMPLEX_PTR( A );
      int*      buff_p = FLA_INT_PTR( p );

      FLA_Apply_pivots_ln_opc_var1( m_A,
                                    buff_A, rs_A, cs_A,
                                    k1_0,
                                    k2_0,
                                    buff_p, inc_p );

      break;
    }

    case FLA_DOUBLE_COMPLEX:
    {
      dcomplex* buff_A = FLA_DOUBLE_COMPLEX_PTR( A );
      int*      buff_p = FLA_INT_PTR( p );

      FLA_Apply_pivots_ln_opz_var1( m_A,
                                    buff_A, rs_A, cs_A,
                                    k1_0,
                                    k2_0,
                                    buff_p, inc_p );

      break;
    }
  }

  return FLA_SUCCESS;
}

FLA_Error FLA_Apply_pivots_rn_opz_var1	(	int	n,
		dcomplex *	a,
		int	a_rs,
		int	a_cs,
		int	k1,
		int	k2,
		int *	p,
		int	incp
	)