59 for (
int j = 0; j < nBlock(); ++j) {
60 monomerId = block(j).monomerId();
61 block(j).setupSolver(wFields[monomerId]);
80 for (
int i = 0; i < nParams_; ++i) {
84 for (
int i = 0; i < nBlock(); ++i) {
85 prefactor = exp(mu_)/length();
86 block(i).computeStress(prefactor);
88 for (
int j=0; j < nParams_; ++j){
89 stress_ [j] += block(i).stress(j);
Field of real double precision values on an FFT mesh.
void computeStress()
Compute stress from a polymer chain, needs a pointer to basis.
void setClassName(const char *className)
Set class name string.
void setMu(double mu)
Set the mu (chemical potential) for this species.
void compute(DArray< RField< D > > const &wFields, double phiTot=1.0)
Compute solution to MDE and block concentrations.
void setPhi(double phi)
Set the phi (volume fraction) for this species.
void setNParams(int nParams)
Store the number of lattice parameters in the unit cell.
Dynamically allocatable contiguous array template.
#define UTIL_CHECK(condition)
Assertion macro suitable for serial or parallel production code.
PSCF package top-level namespace.