solve.hpp

Go to the documentation of this file.

#pragma once
 
#include "NumCpp/Core/Internal/Error.hpp"
#include "NumCpp/Core/Internal/StaticAsserts.hpp"
#include "NumCpp/Functions/dot.hpp"
#include "NumCpp/Linalg/inv.hpp"
#include "NumCpp/NdArray.hpp"
 
namespace nc::linalg
{
    //============================================================================
    // Method Description:
    template<typename dtype>
    NdArray<double> solve(const NdArray<dtype>& inA, const NdArray<dtype>& inB)
    {
        STATIC_ASSERT_ARITHMETIC(dtype);
 
        if (!inA.issquare())
        {
            THROW_INVALID_ARGUMENT_ERROR("input array a must be square.");
        }
 
        if (!inB.isflat())
        {
            THROW_INVALID_ARGUMENT_ERROR("input array b must be flat.");
        }
 
        if (inA.numCols() != inB.size())
        {
            THROW_INVALID_ARGUMENT_ERROR("input array b size must be the same as the square size of a.");
        }
 
        return dot(inv(inA), inB.template astype<double>().reshape(inB.size(), 1)).reshape(inB.shape());
    }
} // namespace nc::linalg