| 64 |
|
Shape* shape; |
| 65 |
|
}; |
| 66 |
|
void registerHydrodynamicsModels(); |
| 67 |
+ |
void writeHydroProps(std::ostream& os); |
| 68 |
|
|
| 69 |
|
int main(int argc, char* argv[]){ |
| 70 |
|
//register force fields |
| 113 |
|
identMat(2,2) = 1.0; |
| 114 |
|
|
| 115 |
|
Globals* simParams = info->getSimParams(); |
| 116 |
< |
double temperature; |
| 117 |
< |
double viscosity; |
| 116 |
> |
RealType temperature; |
| 117 |
> |
RealType viscosity; |
| 118 |
|
|
| 119 |
|
if (simParams->haveViscosity()) { |
| 120 |
|
viscosity = simParams->getViscosity(); |
| 127 |
|
if (simParams->haveTargetTemp()) { |
| 128 |
|
temperature = simParams->getTargetTemp(); |
| 129 |
|
} else { |
| 130 |
< |
sprintf(painCave.errMsg, "target temperaturemust be set\n"); |
| 130 |
> |
sprintf(painCave.errMsg, "target temperature must be set\n"); |
| 131 |
|
painCave.isFatal = 1; |
| 132 |
|
simError(); |
| 133 |
|
} |
| 182 |
|
|
| 183 |
|
//if beads option is turned on, skip the calculation |
| 184 |
|
if (!args_info.beads_flag) { |
| 185 |
< |
shape->calcHydroProps(model, viscosity, temperature); |
| 186 |
< |
model->writeHydroProps(outputDiff); |
| 185 |
> |
model->calcHydroProps(shape, viscosity, temperature); |
| 186 |
> |
model->writeHydroProps(outputDiff); |
| 187 |
|
} |
| 188 |
< |
|
| 188 |
> |
|
| 189 |
|
delete model; |
| 190 |
|
} |
| 191 |
|
|
| 201 |
|
HydrodynamicsModelFactory::getInstance()->registerHydrodynamicsModel(new HydrodynamicsModelBuilder<AnalyticalModel>("AnalyticalModel")); |
| 202 |
|
|
| 203 |
|
} |
| 203 |
– |
|