|
|
@ -71,7 +71,8 @@ INDsetup(SMPmatrix *matrix, GENmodel *inModel, CKTcircuit *ckt, int *states) |
|
|
} |
|
|
} |
|
|
|
|
|
|
|
|
/* Lundin's geometry correction factor */ |
|
|
/* Lundin's geometry correction factor */ |
|
|
model->INDspecInd *= Lundin(model->INDlength, model->INDcsect); |
|
|
|
|
|
|
|
|
if(model->INDlengthGiven && (model->INDdiaGiven || model->INDcsectGiven)) |
|
|
|
|
|
model->INDspecInd *= Lundin(model->INDlength, model->INDcsect); |
|
|
|
|
|
|
|
|
/* |
|
|
/* |
|
|
double kl = Lundin(model->INDlength, model->INDcsect); |
|
|
double kl = Lundin(model->INDlength, model->INDcsect); |
|
|
@ -145,6 +146,8 @@ static double Lundin(double l, double csec) |
|
|
/* x = solenoid diam. / length */ |
|
|
/* x = solenoid diam. / length */ |
|
|
double num, den, kk, x, xx, xxxx; |
|
|
double num, den, kk, x, xx, xxxx; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
if (csec < 1e-12 || l < 1e-6) { |
|
|
if (csec < 1e-12 || l < 1e-6) { |
|
|
fprintf(stderr, "Warning: coil geometries too small (< 1um length dimensions),\n"); |
|
|
fprintf(stderr, "Warning: coil geometries too small (< 1um length dimensions),\n"); |
|
|
fprintf(stderr, " Lundin's correction factor will not be calculated\n"); |
|
|
fprintf(stderr, " Lundin's correction factor will not be calculated\n"); |
|
|
|