additional_input_ca2_adaptive_impl.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 _ADDITIONAL_INPUT_CA2_ADAPTIVE_H_
#define _ADDITIONAL_INPUT_CA2_ADAPTIVE_H_
 
#include "additional_input.h"
 
struct additional_input_params_t {
    REAL    tau_ca2;
    REAL    i_ca2;
    REAL    i_alpha;
    REAL    time_step;
};
 
struct additional_input_t {
    REAL    exp_tau_ca2;
    REAL    i_ca2;
    REAL    i_alpha;
};
 
static inline void additional_input_initialise(
        additional_input_t *state, additional_input_params_t *params,
        uint32_t n_steps_per_timestep) {
    REAL ts = kdivui(params->time_step, n_steps_per_timestep);
    state->exp_tau_ca2 = expk(-kdivk(ts, params->tau_ca2));
    state->i_ca2 = params->i_ca2;
    state->i_alpha = params->i_alpha;
}
 
static inline void additional_input_save_state(additional_input_t *state,
        additional_input_params_t *params) {
    params->i_ca2 = state->i_ca2;
}
 
static inline input_t additional_input_get_input_value_as_current(
        additional_input_t *additional_input,
        UNUSED state_t membrane_voltage) {
    // Decay Ca2 trace
    additional_input->i_ca2 *= additional_input->exp_tau_ca2;
 
    // Return the Ca2
    return -additional_input->i_ca2;
}
 
static inline void additional_input_has_spiked(
        additional_input_t *additional_input) {
    // Apply influx of calcium to trace
    additional_input->i_ca2 += additional_input->i_alpha;
}
 
#endif // _ADDITIONAL_INPUT_CA2_ADAPTIVE_H_