// reference:
// https://github.com/xianyi/OpenBLAS/blob/develop/reference/dtrsmf.f
/// <summary>
/// [dtrsm] Solves one of the matrix equations:
/// <list type="bullet">
/// <item>
/// <term><c>(side, transa) = (Left, None)</c></term>
/// <description><c>A * X = alpha * B</c></description>
/// </item>
/// <item>
/// <term><c>(side, transa) = (Right, None)</c></term>
/// <description><c>X * A = alpha * B</c></description>
/// </item>
/// </list>
/// </summary>
/// <typeparam name="T"></typeparam>
/// <param name="side"></param>
/// <param name="uplo"></param>
/// <param name="diag"> <c>true</c> if <paramref name="a"/> is a unit triangular; otherwise, <c>false</c>. </param>
/// <param name="alpha"></param>
/// <param name="a"></param>
/// <param name="b"></param>
public static void SolveTriangleMatrix <T>(OperandSide side, TriangleKind uplo, bool diag, T alpha, INdArray <T> a, MutableNdArray <T> b)
{
if (ValueTrait.Equals(alpha, Zero <T>()))
{
VectorOperation.Identity(NdArray.Zeros <T>(b.Shape), b);
return;
}
switch ((side, uplo))
{
case (OperandSide.Left, TriangleKind.Upper): SolveTriangleMatrixLU(diag, alpha, a, b); break;
case (OperandSide.Left, TriangleKind.Lower): SolveTriangleMatrixLL(diag, alpha, a, b); break;
case (OperandSide.Right, TriangleKind.Upper): SolveTriangleMatrixRU(diag, alpha, a, b); break;
case (OperandSide.Right, TriangleKind.Lower): SolveTriangleMatrixRL(diag, alpha, a, b); break;
default:
Guard.ThrowArgumentError("Invalid configs.");
break;
}
}
/// <summary>
/// Creates a new UnaryOperator
/// </summary>
/// <param name="strOp">the string representation of the operator</param>
/// <param name="doOp">the method that processes the operand</param>
/// <param name="side">the side that the operand is on</param>
public UnaryOperator(string strOp, UnaryOpDelegate doOp, OperandSide side = OperandSide.Right)
{
OperatorString = strOp;
ExecuteOperator = doOp;
Side = side;
}
