129 double const m1 = 4.0 *
lch / acos(-1.) / ls;
134 double const m2 = std::pow(2., 5. / 3.) - 3.;
135 double const n = 2.5;
136 double const a1 = std::pow(d_ip, n) +
137 m1 * (1. - d_ip) * (1. + m2 * (1. - d_ip));
138 double const a2 = n * std::pow(d_ip, n - 1.) -
139 2. * m1 * m2 * (1. - d_ip) - m1;
140 return (n * std::pow(d_ip, n - 1.) * a1 -
141 std::pow(d_ip, n) * a2) /
146 double const m2 = -0.5;
148 double const a1 = std::pow(d_ip, n) +
149 m1 * (1. - d_ip) * (1. + m2 * (1. - d_ip));
150 double const a2 = n * std::pow(d_ip, n - 1.) -
151 2. * m1 * m2 * (1. - d_ip) - m1;
152 return (n * std::pow(d_ip, n - 1.) * a1 -
153 std::pow(d_ip, n) * a2) /
160 double const m1 = 4.0 *
lch / acos(-1.) / ls;
165 double const m2 = std::pow(2., 5. / 3.) - 3.;
166 double const n = 2.5;
167 double a1 = std::pow(d_ip, n) +
168 m1 * (1. - d_ip) * (1. + m2 * (1. - d_ip));
169 double a2 = n * std::pow(d_ip, n - 1.) -
170 2. * m1 * m2 * (1. - d_ip) - m1;
171 return (2. * a2 * a2 * std::pow(d_ip, n) -
172 a1 * std::pow(d_ip, n) *
174 n * std::pow(d_ip, n - 2.) * (n - 1.)) -
175 2. * a1 * a2 * n * std::pow(d_ip, n - 1.) +
176 a1 * a1 * n * std::pow(d_ip, n - 2.) * (n - 1.)) /
181 double const m2 = -0.5;
183 double const a1 = std::pow(d_ip, n) +
184 m1 * (1. - d_ip) * (1. + m2 * (1. - d_ip));
185 double const a2 = n * std::pow(d_ip, n - 1.) -
186 2. * m1 * m2 * (1. - d_ip) - m1;
187 return (2. * a2 * a2 * std::pow(d_ip, n) -
188 a1 * std::pow(d_ip, n) *
190 n * std::pow(d_ip, n - 2.) * (n - 1.)) -
191 2. * a1 * a2 * n * std::pow(d_ip, n - 1.) +
192 a1 * a1 * n * std::pow(d_ip, n - 2.) * (n - 1.)) /