Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F75855943
deltaPolarMechanicalForces.hpp
No One
Temporary
Actions
Download File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Subscribers
None
File Metadata
Details
File Info
Storage
Attached
Created
Sun, Aug 4, 17:23
Size
1 KB
Mime Type
text/x-c++
Expires
Tue, Aug 6, 17:23 (2 d)
Engine
blob
Format
Raw Data
Handle
19619372
Attached To
R9411 tisue modeling
deltaPolarMechanicalForces.hpp
View Options
#ifndef __DELTA_POLAR_MECHANICAL_FORCES_H__
#define __DELTA_POLAR_MECHANICAL_FORCES_H__
#include <iostream>
#include <boost/array.hpp>
#include "deltaMechanicalForces.hpp"
class
DeltaPolarMechanicalForces
:
public
DeltaMechanicalForces
{
protected
:
Lattice
*
lattice
;
Interactor
interaction
;
double
gamma_d
;
double
h_dc
;
double
dprod
;
double
hfdp
;
double
del_sig
;
double
tau_sig
;
double
tau_mech
;
protected
:
double
deltacisfunc
(
double
delta
)
{
double
a
=
pow
(
h_dc
,
hill_co
);
double
b
=
pow
(
delta
,
hill_co
);
return
a
/
(
a
+
b
);}
virtual
bool
isHC
(
Cell
*
c
){
return
m_signal
[
c
->
id
()]
<
0.05
;}
virtual
bool
isSC
(
Cell
*
c
){
return
m_signal
[
c
->
id
()]
>=
2
;}
// double weightfunc(double s)
// {double a=pow(h_we,hill_co);double b=pow(s,hill_co);return we_min+we_r*a/(a+b);}
public
:
DeltaMechanicalForces
(){
initParameters
();}
void
initParameters
();
void
operator
()(
const
state_type
&
y
,
state_type
&
dydt
,
const
double
t
);
virtual
float
getWeight
(
Boundary
*
b
);
// void setState(state_type& s){MechanicalForces::setState(s);CellInteraction::setState(s);}
virtual
void
initState
(
state_type
&
s
);
virtual
void
initState
(
state_type
&
s
,
double
*
delta
);
virtual
double
getDeltaProduction
(
double
si
){
return
dprod
/
(
1
+
pow
(
del_sig
*
si
,
hfdp
));}
virtual
uint
cellStateDim
(){
return
1
;}
virtual
uint
BoundaryStateDim
(){
return
0
;}
virtual
uint
setStateVectorIndex
();
virtual
void
setLattice
(
Lattice
*
lat
);
void
fillDeltaVector
(
double
*
v
,
state_type
&
s
);
// void fillSignalVector(double *v); // in Mechanical forces
void
updateDelta
(
const
state_type
&
y
,
state_type
&
dydt
);
virtual
void
updateSignal
(
const
state_type
&
y
);
// void findT1Transitions();
// virtual void checkT1Transition(Boundary *b);
// bool fromState(const state_type &y);
};
// void myWrapper(Boundary *b, void *args){ checkT1TransitionWrapper<DeltaMechanicalForces>(b, args);}
//checkT1TransitionWrapper<DeltaMechanicalForces>
#endif
Event Timeline
Log In to Comment