1: //
   2: // Heccer : a compartmental solver that implements efficient Crank-Nicolson
   3: // integration for neuronal models.
   4: //
   5: 
   6: //////////////////////////////////////////////////////////////////////////////
   7: //'
   8: //' Heccer : testbed C implementation
   9: //'
  10: //' Copyright (C) 2006-2008 Hugo Cornelis
  11: //'
  12: //' functional ideas .. Hugo Cornelis, hugo.cornelis@gmail.com
  13: //'
  14: //' coding ............ Hugo Cornelis, hugo.cornelis@gmail.com
  15: //'
  16: //////////////////////////////////////////////////////////////////////////////
  17: 
  18: 
  19: #include "../../heccer/compartment.h"
  20: #include "../../heccer/heccer.h"
  21: 
  22: 
  23: struct Compartment comp[] =
  24: {
  25:     {
  26:         //m administrative overhead
  27: 
  28:         {
  29:             //m type of structure
  30: 
  31:             MATH_TYPE_Compartment,
  32:         },
  33: 
  34:         //m index of parent compartment, -1 for none
  35: 
  36:         -1,
  37: 
  38: /*      //m first mechanism */
  39: 
  40: /*      NULL, */
  41: 
  42: /*      //m number of mechanisms */
  43: 
  44: /*      0, */
  45: 
  46:         //m descriptive values, alphabetical order
  47: 
  48: /*     double dCm; */
  49: 
  50:         4.57537e-11, // unscaled 0.0164,
  51: 
  52: /*     double dEm; */
  53: 
  54:         -0.08,
  55: 
  56: /*     double dInitVm; */
  57: 
  58:         -0.068,
  59: 
  60: /*     double dInject;           */
  61: 
  62:         0,
  63: 
  64: /*     double dRa; */
  65: 
  66:         360502, // unscaled 2.5,
  67: 
  68: /*     double dRm; */
  69: 
  70:         3.58441e+08, // unscaled 1
  71:     },
  72: 
  73:     {
  74:         //m administrative overhead
  75: 
  76:         {
  77:             //m type of structure
  78: 
  79:             MATH_TYPE_Compartment,
  80:         },
  81: 
  82:         //m index of parent compartment, -1 for none
  83: 
  84:         0,
  85: 
  86: /*      //m first mechanism */
  87: 
  88: /*      NULL, */
  89: 
  90: /*      //m number of mechanisms */
  91: 
  92: /*      0, */
  93: 
  94:         //m descriptive values, alphabetical order
  95: 
  96: /*     double dCm; */
  97: 
  98:         4.57537e-11, // unscaled 0.0164,
  99: 
 100: /*     double dEm; */
 101: 
 102:         -0.08,
 103: 
 104: /*     double dInitVm; */
 105: 
 106:         -0.068,
 107: 
 108: /*     double dInject;           */
 109: 
 110:         0,
 111: 
 112: /*     double dRa; */
 113: 
 114:         360502, // unscaled 2.5,
 115: 
 116: /*     double dRm; */
 117: 
 118:         3.58441e+08, // unscaled 1
 119:     },
 120: 
 121:     {
 122:         //m administrative overhead
 123: 
 124:         {
 125:             //m type of structure
 126: 
 127:             MATH_TYPE_Compartment,
 128:         },
 129: 
 130:         //m index of parent compartment, -1 for none
 131: 
 132:         1,
 133: 
 134: /*      //m first mechanism */
 135: 
 136: /*      NULL, */
 137: 
 138: /*      //m number of mechanisms */
 139: 
 140: /*      0, */
 141: 
 142:         //m descriptive values, alphabetical order
 143: 
 144: /*     double dCm; */
 145: 
 146:         4.57537e-11, // unscaled 0.0164,
 147: 
 148: /*     double dEm; */
 149: 
 150:         -0.08,
 151: 
 152: /*     double dInitVm; */
 153: 
 154:         -0.068,
 155: 
 156: /*     double dInject;           */
 157: 
 158:         0,
 159: 
 160: /*     double dRa; */
 161: 
 162:         360502, // unscaled 2.5,
 163: 
 164: /*     double dRm; */
 165: 
 166:         3.58441e+08, // unscaled 1
 167:     },
 168: 
 169:     {
 170:         //m administrative overhead
 171: 
 172:         {
 173:             //m type of structure
 174: 
 175:             MATH_TYPE_Compartment,
 176:         },
 177: 
 178:         //m index of parent compartment, -1 for none
 179: 
 180:         0,
 181: 
 182: /*      //m first mechanism */
 183: 
 184: /*      NULL, */
 185: 
 186: /*      //m number of mechanisms */
 187: 
 188: /*      0, */
 189: 
 190:         //m descriptive values, alphabetical order
 191: 
 192: /*     double dCm; */
 193: 
 194:         4.57537e-11, // unscaled 0.0164,
 195: 
 196: /*     double dEm; */
 197: 
 198:         -0.08,
 199: 
 200: /*     double dInitVm; */
 201: 
 202:         -0.068,
 203: 
 204: /*     double dInject;           */
 205: 
 206:         0,
 207: 
 208: /*     double dRa; */
 209: 
 210:         360502, // unscaled 2.5,
 211: 
 212: /*     double dRm; */
 213: 
 214:         3.58441e+08, // unscaled 1
 215:     },
 216: };
 217: 
 218: 
 219: int piC2m[] =
 220: {
 221:     0,
 222:     0,
 223:     0,
 224:     0,
 225:     -1,
 226: };
 227: 
 228: 
 229: struct Intermediary inter =
 230: {
 231:     //m compartment array
 232: 
 233:     4,
 234: 
 235:     comp,
 236: 
 237:     //m all other mathematical components
 238: 
 239:     NULL,
 240: 
 241:     //m compartment 2 first mechanism number
 242: 
 243:     piC2m,
 244: };
 245: 
 246: 
 247: #include "main.c"
 248: 
 249: 
 250: