C user function for scaling parameters definition for electric model.
if (ielarc >= 1) {
else
}
double econs = 1.5e5;
double coepot = 0.;
double coepoa = 1.;
double amex = 1.e30;
double aiex = -1.e30;
double emax = 0.;
double *w1;
for (
int iel = 0; iel <
ncel; iel++) {
w1[iel] = pow(xelec * xelec + yelec * yelec + zelec * zelec, 0.5);
}
bft_printf(
"min and max for E : %14.5E %15.4E\n", amex, aiex);
if (aiex > econs) {
double diff = 0.;
double xyzmax[3] = {-1.e10, -1.e10, -1.e10};
for (
int iel = 0; iel <
ncel; iel++) {
diff = aiex - w1[iel];
if (diff < 1.e-6) {
emax = w1[iel];
xyzmax[0] = xyzcen[3 * iel ];
xyzmax[1] = xyzcen[3 * iel + 1];
xyzmax[2] = xyzcen[3 * iel + 2];
}
}
bft_printf(
"restrike point : %14.5E %14.5E %14.5E\n",
}
if (z1 < 0.)
z1 = 0.;
if (z2 > 2.e-2)
z2 = 2.e-2;
for (
int iel = 0; iel <
ncel; iel++) {
if (xyzcen[3 * iel + 2] > z1 && xyzcen[3 * iel + 2] < z2) {
rayo /= denom;
rayo = pow(rayo, 0.5);
if (rayo < 5.e-4 && posi <= 0.)
}
}
}
else {
}
double somje = 0.;
for (
int iel = 0; iel <
ncel; iel++) {
}
if (fabs(somje) > 1.-20)
bft_printf(
"imposed current %14.5E, Dpot %14.5E, Somje %14.5E\n",
somje);
for (
int ifac = 0; ifac <
nfac; ifac++) {
if (fabs(surfac[ifac][0]) < 1.e-8 && fabs(surfac[ifac][1]) < 1.e-8 &&
int iel =
mesh->i_face_cells[ifac][0];
}
}
if (fabs(elcou) > 1.e-6)
elcou = fabs(elcou);
else
elcou = 0.;
if (fabs(elcou) > 1.e20)
coepot = coepoa;
double dtj = 1.e15;
double dtjm = dtj;
double delhsh = 0.;
double cdtj = 20.;
for (
int iel = 0; iel <
ncel; iel++) {
if (fabs(
CS_F_(rho)->val[iel]) > 1.e-20)
if (fabs(delhsh) > 1.e-20)
dtjm =
CS_F_(
h)->val[iel] / delhsh;
else
dtjm = dtj;
dtjm = fabs(dtjm);
}
double cpmx = pow(cdtj * dtj, 0.5);
coepot = cpmx;
if (coepoa > 1.05)
coepot = cpmx;
else
coepot = coepoa;
}
for (
int iel = 0; iel <
ncel; iel++)
if (ielarc > 0)
for (int i = 0; i < 3 ; i++)
for (int iel = 0; iel < 3 ; iel++)
for (int iel = 0; iel < 3 ; iel++)
}