omnet_helpers.R

#
# Copyright (C) 2016-2022 Stefan Joerer <joerer@ccs-labs.org>
# Copyright (C) 2016-2022 Bastian Bloessl <bloessl@ccs-labs.org>
# Copyright (C) 2016-2022 Michele Segata <segata@ccs-labs.org>
#
# SPDX-License-Identifier: GPL-2.0-or-later
#
# This program is free software; you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation; either version 2 of the License, or
# (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program; if not, write to the Free Software
# Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
#

library(omnetpp)
if (!suppressWarnings(require(reshape, warn.conflicts = F, quietly=T))) {
    library(reshape2, warn.conflicts = F)
}

rVersion4OrGreater <- function() {
    return(as.numeric(version$major) >= 4)
}

getVectorNames <- function(dataset) {
    if (rVersion4OrGreater()) {
        return(unique(dataset$vectors$name))
    } else {
        return(levels(dataset$vectors$name))
    }
}

getScalars <- function(scaFiles, module, scalarName, attributeName) {

    selectString <- paste('module(', module, ') AND name("', scalarName,'")', sep='')
    dataset <- loadDataset(scaFiles, add('scalar', select=selectString))

    # attrvalues is a factor with levels
    attrvalues <- unique(subset(dataset$runattrs$attrvalue, dataset$runattrs$attrname == attributeName))
    if (rVersion4OrGreater()) {
        attrvalues <- sort(as.integer(attrvalues[as.integer(attrvalues)]))
    } else {
        attrvalues <- sort(as.integer(levels(attrvalues)[as.integer(attrvalues)]))
    }

    res <- matrix(0, nrow=0, ncol=2)

    # iterate through distances
    for (i in 1:length(attrvalues)) {
        a <- attrvalues[i]
        runids <- subset(dataset$runattrs$runid, dataset$runattrs$attrname == attributeName & dataset$runattrs$attrvalue == a)
        p <- subset(dataset$scalars$value, is.element(dataset$scalars$runid, runids))
        res <- rbind(res, cbind(a, p))
    }

    return(res)
}

loadVectorsShaped <- function(vecFiles, ...){
    d <- loadDataset(vecFiles, ...)

    vectornames <- getVectorNames(d)
    if(length(vectornames) == 0) {
        warning("no vectors loaded!")
        return
    }
    if(length(d$scalars$name)){
        warning("scalars loaded, but will not be considered in this function")
    }

    #preprocess first vector
    reskeys <- subset(d$vectors$resultkey, d$vectors$name == vectornames[1])
    v <- loadVectors(d,reskeys)
    result <- v$vectordata

    #process remaining vectors
    if(length(vectornames) > 1) {
        for(name in vectornames[2:length(vectornames)]) {
            reskeys <- subset(d$vectors$resultkey, d$vectors$name == name)
            v <- loadVectors(d,reskeys)
            result <- merge(result, v$vectordata, by=c("eventno", "resultkey", "x"), all.x = TRUE, all.y = TRUE, suffixes = c("",name))
        }
    }
    names(result) <- c(names(result)[1:2], "time", vectornames)
    return (result)
}

listVectors <- function(vecFiles, ...){
    d <- loadDataset(vecFiles, ...)
    return (getVectorNames(d))
}

print.runs <- function(vecFiles) {
    d <- loadDataset(vecFiles)
    print(cast(d$runattrs, runid~attrname, value='attrvalue'))
}