current_source_ac.h

Go to the documentation of this file.

/*
 * Copyright (c) 2017 The University of Manchester
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
 * You may obtain a copy of the License at
 *
 *     https://www.apache.org/licenses/LICENSE-2.0
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 */
 
#ifndef _CURRENT_SOURCE_AC_H_
#define _CURRENT_SOURCE_AC_H_
 
#include <sincos.h>
 
// Structure for AC current source
typedef struct ac_source_t {
    uint32_t start;
    uint32_t stop;
    REAL amplitude;
    REAL offset;
    REAL frequency;
    REAL phase;
} ac_source_t;
 
static ac_source_t **ac_source;
 
static bool current_source_ac_init(uint32_t n_ac_sources, uint32_t *next) {
    ac_source = spin1_malloc(n_ac_sources * sizeof(uint32_t*));
    for (uint32_t n_ac=0; n_ac < n_ac_sources; n_ac++) {
        ac_source[n_ac] = spin1_malloc(sizeof(ac_source_t));
        if (ac_source[n_ac] == NULL) {
            log_error("Unable to allocate DC source parameters - out of DTCM");
            return false;
        }
        *next += sizeof(ac_source_t) / 4;
    }
    return true;
}
 
static bool current_source_ac_load_parameters(
        address_t cs_address, uint32_t n_ac_sources, uint32_t *next) {
    for (uint32_t n_ac=0; n_ac < n_ac_sources; n_ac++) {
        spin1_memcpy(ac_source[n_ac], &cs_address[*next], sizeof(ac_source_t));
        *next += sizeof(ac_source_t) / 4;
    }
    return true;
}
 
static REAL current_source_ac_get_offset(uint32_t cs_index, uint32_t time) {
    if ((time >= ac_source[cs_index]->start) && (time < ac_source[cs_index]->stop)) {
        REAL time_value = kbits((time - ac_source[cs_index]->start) << 15);
        REAL sin_value = sink((time_value * ac_source[cs_index]->frequency) +
                ac_source[cs_index]->phase);
        REAL ac_current_offset = ac_source[cs_index]->offset + (
                ac_source[cs_index]->amplitude * sin_value);
        return ac_current_offset;
    }
    return ZERO;
}
 
#endif // _CURRENT_SOURCE_AC_H_