R/simulatorDT.R

In saezlab/CNORdt: Add-on to CellNOptR: Discretized time treatments

#
#  This file is part of the CNO software
#
#  Copyright (c) 2011-2012 - EBI
#
#  File author(s): CNO developers (cno-dev@ebi.ac.uk)
#
#  Distributed under the GPLv2 License.
#  See accompanying file LICENSE.txt or copy at
#      http://www.gnu.org/licenses/gpl-2.0.html
#
#  CNO website: http://www.ebi.ac.uk/saezrodriguez/software.html
#
##############################################################################
# $Id: $

simulatorDT <- function(CNOlist, model, simList, indices, boolUpdates, prevSim=NULL) {
	
	if ((class(CNOlist) == "CNOlist") == FALSE) {
    	CNOlist = CellNOptR::CNOlist(CNOlist)
	}
	
	# CONSTANTS
	# cnolist
	nStimuli <- as.integer(length(indices$stimulated))
	nInhibitors <- as.integer(length(indices$inhibited))
	nSignals <- as.integer(length(indices$signals))
	nCond <- as.integer(dim(CNOlist@signals[[1]])[1])
	nTimes <- as.integer(length(CNOlist@timepoints))

	# model
	nReacs <- as.integer(length(model$reacID))
	nSpecies <- as.integer(length(model$namesSpecies))

	# simList
	nMaxInputs <- as.integer(dim(simList$finalCube)[2])
	finalCube = as.integer(as.vector(t(simList$finalCube))-1)
	ixNeg = as.integer(as.vector(t(simList$ixNeg)))
	ignoreCube = as.integer(as.vector(t(simList$ignoreCube)))
	maxIx = as.integer(simList$maxIx-1)
	
	# index
	indexSignals <- as.integer(as.vector(indices$signals)-1)
	indexStimuli <- as.integer(as.vector(indices$stimulated)-1)
	indexInhibitors <- as.integer(as.vector(indices$inhibited)-1)

	# STRUCTURES
	# cnolist
	valueInhibitors <- as.integer(t(CNOlist@inhibitors))
	valueStimuli <- as.integer(t(CNOlist@stimuli))
	
	# additional to DT
	boolUpdates <- as.integer(boolUpdates)
	# set any NA to 2 (if present)
	if(is.null(prevSim)) {
		prevSim = matrix(2, nrow=dim(CNOlist@signals[[1]])[1],
		ncol=length(model$reacID))
	} else {
		prevSim[is.na(prevSim)] = 2
	}
	prevSim <- as.integer(t(prevSim))

	
	res = .Call("simulatorDT",
		# variables	
		nStimuli,
		nInhibitors,
		nSignals,
		nCond,
		nTimes,
		nReacs,
		nSpecies,
		nMaxInputs,
		boolUpdates,
		# vectors
		maxIx,		
		indexSignals, 
		indexStimuli, 
		indexInhibitors,
		# matrices
		finalCube,
		ixNeg,
		ignoreCube,
		valueInhibitors,
		valueStimuli,
		prevSim
	)
	
	return(res)
}