コード例 #1
0
ファイル: matrix.cpp プロジェクト: sisl/Smile.jl
double dmatrix_GetAtCoord( void * void_dmat, int * theCoordinates, unsigned len )
{
	DSL_Dmatrix * dmat = reinterpret_cast<DSL_Dmatrix*>(void_dmat);
	
	DSL_intArray coords = DSL_intArray(len);
	for (unsigned i = 0; i < len; i++)
		coords[i] = theCoordinates[i];

	double retval = dmat->Subscript(coords);
	return retval;
}
コード例 #2
0
ファイル: DSLPNLConverter.cpp プロジェクト: JacobCWard/PyPNL
void DSLPNLConverter::CreateFactors(DSL_network& dslNet, CBNet* pnlBNet)
{
    int i,j,k;
    
    // Read number of nodes in the net
    int numberOfNodes = dslNet.GetNumberOfNodes();
    
    // This is a way PNL likes it	
    pnlBNet->AllocFactors();
    
    DSL_Dmatrix* dslMatrix;
    CCPD* pnlCPD;
    
    for (i=0;i<numberOfNodes;i++)
    {
        
        // Get parents of the ith node 
        // IMPORTANT -- we should preserve order from DSL_network, since
        // probabilities will be according DSL ordering
        
        DSL_intArray dslParents;
        dslParents = dslNet.GetParents(dslNet.FindNode(theIds[i]));
        
        // establish sizes and allocate memory
        int numberOfNodesInDomain = dslParents.NumItems() + 1;
        int* domain = new int[numberOfNodesInDomain];
        CNodeType** nodeTypes = new CNodeType*[numberOfNodesInDomain];
        
        // establish members of the domain
        for (j=0;j<numberOfNodesInDomain-1;j++)
            domain[j] = dslParents[j];
        
        domain[numberOfNodesInDomain-1] = i;
        
        // Fill up node types
        for (j=0;j<numberOfNodesInDomain;j++)
            nodeTypes[j] = const_cast <CNodeType*> (pnlBNet->GetNodeType(domain[j]));
        
        // Read CPT from SMILE
        dslNet.GetNode(dslNet.FindNode(theIds[i]))->Definition()->GetDefinition(&dslMatrix);
        
        // Alloc space for CPT
        int sizeOfCPT = dslMatrix->GetSize();
        float* flatCPT = new float[sizeOfCPT];
        
        // Here we convert 'copy' numbers from SMILE to PNL
        // The painful part is convert double to float.
        // Additionally we check if after conversion they sum-up to 1
        int numberOfMyStates = nodeTypes[numberOfNodesInDomain-1]->GetNodeSize();
        
        int iterations  = sizeOfCPT/numberOfMyStates;
        for (j=0;j<iterations;j++)
        {
            float sum = 0.0f;
            for (k=0;k<numberOfMyStates;k++)
            {
                flatCPT[j*numberOfMyStates+k] = static_cast <float> (dslMatrix->Subscript(j*numberOfMyStates+k));
                sum += flatCPT[j*numberOfMyStates+k];
            }
            if (sum!=1.0f)
            {
                for (k=0;k<numberOfMyStates;k++)
                    flatCPT[j*numberOfMyStates+k] /= sum;      
            }
        }
        
#ifdef DSLPNL_DEBUG
        std::cerr << "Node "<< i << " domain : ";
        for (j=0;j<numberOfNodesInDomain;j++)
        {
            std::cerr << domain[j] << " ";
        }
        std::cerr <<  std::endl;
        for (j=0;j<sizeOfCPT;j++)
            std::cerr << flatCPT[j] << " ";
        std::cerr <<  std::endl;
#endif
        
        CModelDomain* pMD = pnlBNet->GetModelDomain();
        pnlCPD = CTabularCPD::Create(domain, numberOfNodesInDomain, pMD, flatCPT);
        if (pnlCPD==NULL)
        {
            std::cout << "We got a problem with creating CPD" << std::endl;
            return;
        }
        pnlBNet->AttachFactor(pnlCPD);
        
        delete[] nodeTypes;
        delete[] domain;
        delete[] flatCPT; 
    }
}
コード例 #3
0
ファイル: matrix.cpp プロジェクト: sisl/Smile.jl
int dmatrix_CoordinatesToIndex( void * void_dmat, void * void_intarr )
{
	DSL_Dmatrix * dmat = reinterpret_cast<DSL_Dmatrix*>(void_dmat);
	DSL_intArray * arr = reinterpret_cast<DSL_intArray*>(void_intarr);
	return dmat->CoordinatesToIndex(*arr);
}
コード例 #4
0
ファイル: matrix.cpp プロジェクト: sisl/Smile.jl
int dmatrix_GetSizeOfDimension( void * void_dmat, int aDimension )
{
	DSL_Dmatrix * dmat = reinterpret_cast<DSL_Dmatrix*>(void_dmat);
	return dmat->GetSizeOfDimension( aDimension );
}
コード例 #5
0
ファイル: matrix.cpp プロジェクト: sisl/Smile.jl
int dmatrix_GetSize( void * void_dmat )
{
	DSL_Dmatrix * dmat = reinterpret_cast<DSL_Dmatrix*>(void_dmat);
	return dmat->GetSize();
}
コード例 #6
0
ファイル: matrix.cpp プロジェクト: sisl/Smile.jl
int dmatrix_GetNumberOfDimensions( void * void_dmat )
{
	DSL_Dmatrix * dmat = reinterpret_cast<DSL_Dmatrix*>(void_dmat);
	return dmat->GetNumberOfDimensions();
}
コード例 #7
0
ファイル: matrix.cpp プロジェクト: sisl/Smile.jl
void dmatrix_SetAtInd( void * void_dmat, int index, double value )
{
	DSL_Dmatrix * dmat = reinterpret_cast<DSL_Dmatrix*>(void_dmat);
	dmat->Subscript(index) = value;
}
コード例 #8
0
ファイル: matrix.cpp プロジェクト: sisl/Smile.jl
double dmatrix_GetAtInd( void * void_dmat, int index )
{
	DSL_Dmatrix * dmat = reinterpret_cast<DSL_Dmatrix*>(void_dmat);
	double retval = dmat->Subscript(index);
	return retval;	
}