using System;
using System.Collections.Generic;
using System.Runtime.Serialization;
using System.Text;
using System.IO;
using DataPerformer;
using CategoryTheory;
using DiagramUI;
namespace Grav36
{
public class Gravity : IObjectTransformer, ICategoryObject, IPropertiesEditor,
IObjectFactory
{
#region Fields
private System.Drawing.Icon icon;
private double[] R = new double[3];
private double[] C = new double[700];
private double[] S = new double[700];
private double[] HP = new double[37];
private double[] CO = new double[37];
private double[] SI = new double[37];
private double[] AR = new double[37];
private double[] CF = new double[37];
private double[] PNK = new double[37];
private double[] ANAI = new double[478];
private int n0;
private int nk;
const Double ret = 0;
char point = '.';
static private readonly string[] inps = new string[] { "x", "y", "z" };
static private readonly string[] outs = new string[] { "Gx", "Gy", "Gz" };
private object obj;
StreamReader reader;
static private readonly char[] sep = " ".ToCharArray();
string[] inp;
string current;
int pos;
static private readonly double[] SK = new double[]
{
1.732050807568877E0, 1.936491673103709E0,2.091650066335189E0,
2.218529918662356E0, 2.326813808623286E0,2.421824596249695E0,
2.506826616960176E0, 2.583977731709147E0,2.654784752117980E0,
2.720344864917320E0, 2.781483843970261E0,2.838840060634283E0,
2.892918063839265E0, 2.944124128779573E0,2.992790634483277E0,
3.039193256447120E0, 3.083563388216997E0,3.126097306274296E0,
3.166963057815222E0, 3.206305722292480E0,3.244251489527417E0,
3.280910862053330E0, 3.316381199514726E0,3.350748761981671E0,
3.384090366884451E0, 3.416474744628894E0,3.447963656780269E0,
3.478612825366963E0, 3.508472710600489E0,3.537589165949811E0,
3.566003993231045E0, 3.593755415611321E0,3.620878482777667E0,
3.647405419702514E0, 3.673365928240249E0,3.698787449063569E0
};
#endregion
#region Constructors
public Gravity()
{
icon = Grav36.Properties.Resources.EARTH;
System.Globalization.CultureInfo culture = System.Globalization.CultureInfo.CurrentCulture;
if (culture.TwoLetterISOLanguageName.ToLower().Equals("ru"))
{
this.point = ',';
}
}
#endregion
#region ICategoryObject Members
ICategory ICategoryObject.Category
{
get { throw new Exception("The method or operation is not implemented."); }
}
ICategoryArrow ICategoryObject.Id
{
get { throw new Exception("The method or operation is not implemented."); }
}
#endregion
#region IAssociatedObject Members
object IAssociatedObject.Object
{
get
{
return obj;
}
set
{
obj = value;
}
}
#endregion
#region IObjectTransformer Members
string[] IObjectTransformer.Input
{
get { return inps; }
}
string[] IObjectTransformer.Output
{
get { return outs; }
}
object IObjectTransformer.GetInputType(int i)
{
return ret;
}
object IObjectTransformer.GetOutputType(int i)
{
return ret;
}
void IObjectTransformer.Calculate(object[] input, object[] output)
{
double x = (double)input[0];
double y = (double)input[1];
double z = (double)input[2];
double fx, fy, fz;
forces(n0, nk, x, y, z, out fx, out fy, out fz);
output[0] = fx;
output[1] = fy;
output[2] = fz;
}
#endregion
#region IPropertiesEditor Members
object IPropertiesEditor.Editor
{
get
{
// Returns a property editor and icon
return new object[] { new FormEditor(this), icon };
}
}
object IPropertiesEditor.Properties
{
get
{
return saver;
}
set
{
saver = value as List<object>;
}
}
#endregion
#region IObjectFactory Members
string[] IObjectFactory.Names
{
get { return new string[] { "Gravity 36 * 36" }; }
}
ICategoryObject IObjectFactory.this[string name]
{
get { return this; }
}
#endregion
#region Specific Members
public int N0
{
get
{
return n0;
}
set
{
n0 = value;
}
}
public int NK
{
get
{
return nk;
}
set
{
nk = value;
}
}
public double[] MUR
{
get
{
return R;
}
}
public double[] Cnm
{
get
{
return C;
}
}
public double[] Snm
{
get
{
return S;
}
}
public void Load(string filename)
{
double[] ANAi = new double[37];
double[] Hp = new double[700];
int N, N1;
int N0 = 36;
int N01 = N0 - 1;
int NP = N01 * (N0 + 4) / 2;
for (int i = 0; i < 700; i++)
{
C[i] = 0; S[i] = 0;
}
reader = new StreamReader(filename);
reader.ReadLine();
read(out R[1]);
read(out R[2]);
read(out R[0]);
m200:
int I, J; double CC, SS;
if (!(read(out I))) goto m300;
if (!(read(out J))) goto m300;
if (!(read(out CC))) goto m300;
if (!(read(out SS))) goto m300;
if (I <= N0)
{
if (J < 2)
{
N = N01 * J + I - 1;
}
else
{
N = N0 - J + 1;
N = NP - N * (N + 1) / 2 + I - J + 1;
}
C[N - 1] = CC;
S[N - 1] = SS;
}
goto m200;
m300:
reader.Close();
N01 = N0 + 1;
for (N = 1; N <= N0; N++) Hp[N - 1] = 1 + (double)N;
for (N1 = 1; N1 <= N01; N1++)
{
N = N1 - 1;
if (N > 3)
ANAi[N1 - 1] = N * (N0 - 1) + 1 - (N - 3) * (N - 2) / 2;
else
ANAi[N1 - 1] = 1 + N * (N0 - 1);
}
for (int i = 36; i < 700; i++) Hp[i] = 0;
for (int i = 0; i < 37; i++) Hp[222 + i] = ANAi[i];
for (int i = 0; i < 37; i++)
{
HP[i] = Hp[i];
CO[i] = Hp[i + 37];
SI[i] = Hp[i + 2 * 37];
AR[i] = Hp[i + 3 * 37];
CF[i] = Hp[i + 4 * 37];
PNK[i] = Hp[i + 5 * 37];
}
for (int i = 0; i < 478; i++) ANAI[i] = Hp[i + 6 * 37];
}
private bool read(out int i)
{
if (!next())
{
i = 0;
return false;
}
try
{
i = Int32.Parse(current);
return true;
}
catch (Exception)
{
}
i = 0;
return false;
}
private bool read(out double a)
{
if (!next())
{
a = 0;
return false;
}
try
{
a = Double.Parse(current.Replace('.', point));
return true;
}
catch (Exception)
{
a = 0;
return false;
}
}
private void prepare()
{
if (inp == null)
{
string s = reader.ReadLine();
inp = s.Split(sep);
pos = 0;
}
current = inp[pos].Trim();
if (current.Length == 0)
{
}
++pos;
if (pos >= inp.Length)
{
inp = null;
}
}
void nextLine()
{
string s = reader.ReadLine();
inp = s.Split(sep);
}
bool next()
{
if (inp == null)
{
string s = reader.ReadLine();
if (s == null)
{
return false;
}
inp = s.Split(sep);
}
while (true)
{
if (pos == inp.Length)
{
pos = 0;
string s = reader.ReadLine();
if (s == null)
{
return false;
}
inp = s.Split(sep);
}
current = inp[pos];
++pos;
if (current.Length > 0)
{
break;
}
}
return true;
}
private double sqrt(double x)
{
return Math.Sqrt(x);
}
private double sqrt(int i)
{
return Math.Sqrt((double)i);
}
private void forces(int N0, int NK, double X, double Y, double Z, out double FX,
out double FY, out double FZ)
{
bool LOG;
double FR, FF, FL, R2, R3, R1, SF, GR, A, P20 = 0, P30 = 0, PN0, CK1, CK2, TG, AN;
int N, N5, N4, N3, N2, N1, J, M;
FR = 0.0;
FF = 0;
FL = 0;
R2 = X * X + Y * Y;
R3 = 1 / (R2 + Z * Z);
R1 = sqrt(R3);
R2 = sqrt(R2);
SF = Z * R1;
CF[0] = R2 * R1;
R2 = 1 / R2;
CO[0] = X * R2;
SI[0] = Y * R2;
GR = R[0] * R3;
if (N0 != 0 || NK != 0)
{
CF[1] = CF[0] * CF[0];
AR[0] = R[1] * R1;
for (N = 1; N < N0; N++)
{ //1
N3 = N - 1;
AR[N] = AR[0] * AR[N3];
A = C[N3] * AR[N];
if (N == 1)
{ //1.1
P20 = /*SQ[4]*/sqrt(5) * (1 - 1.5 * CF[1]);
PNK[0] =/*SQ[14]*/sqrt(15) * CF[0] * SF;
FR = 3 * A * P20;
FF = A * PNK[0] * sqrt(3);//SQ[2];
} //1.1
else if (N == 2)
{ //1.2
P30 = /*SQ[6]*/sqrt(7) * (1 - 2.5 * CF[1]) * SF;
PNK[1] =/*SQ[62]/SQ[23]*/sqrt(63) / sqrt(24) * CF[0] * (4 - 5 * CF[1]);
FR += 4 * A * P30;
FF += A * PNK[1] * sqrt(6);//SQ[5];
} //1.2
else
{ //1.3
N1 = N + N3 + 1;
N2 = N1 + 2;
N4 = N + 1;
AN = (double)(N + 1);
PN0 = sqrt(N2 + 1)/*SQ[N2]*// AN * (/*SQ[N1]*/sqrt(N1 + 1) * SF * P30 - (double)(N3 + 1) / /*SQ[N1-2]*/sqrt(N1 - 1) * P20);
PNK[N3] = sqrt(N2 + 1) / (sqrt(N3 + 1) * sqrt(N4 + 1)) * (sqrt(N1 + 1) * SF * PNK[N - 2] - sqrt(N + 1) *
sqrt(N - 1) / sqrt(N1 - 1) * PNK[N - 3]);
FR += (double)(N4 + 1) * A * PN0;
FF += A */*SQ[N]*SQ[N4]/SQ[1]*/sqrt(N + 1) * sqrt(N4 + 1) / sqrt(2) * PNK[N3];
P20 = P30;
P30 = PN0;
} //1.3
} //1 CONTINUE
if (NK == 0) goto m5;
LOG = (NK >= 3);
A = CO[0] + CO[0];
CO[1] = A * CO[0] - 1;
SI[1] = A * SI[0];
TG = Z * R2;
if (LOG)
for (N = 2; N < NK; N++)
{ //2
N1 = N - 1;
N2 = N - 2;
CF[N] = CF[0] * CF[N1];
CO[N] = A * CO[N1] - CO[N2];
SI[N] = A * SI[N1] - SI[N2];
} //2
CK1 = (C[35] * CO[0] + S[35] * SI[0]) * AR[1];
CK2 = (C[35] * SI[0] - S[35] * CO[0]) * AR[1];
A = PNK[0];
PNK[0] = SK[1] * CF[1];
FR += CK1 * 3 * A;
FF += CK1 * (PNK[0] + PNK[0] - TG * A);
FL += CK2 * A;
J = 35;
if (LOG)
for (N = 2; N < NK; N++) //commain
{ //3
J++;
N1 = N - 1;
N2 = N + N + 1;
N3 = N + 1;
A = PNK[N1];
CK1 = (C[J] * CO[0] + S[J] * SI[0]) * AR[N];
CK2 = (C[J] * SI[0] - S[J] * CO[0]) * AR[N];
if (N == 2)
PNK[1] = /*SQ[6]*/sqrt(7) * SF * PNK[0];
else
PNK[N1] = sqrt(N2 + 2) / (sqrt(N - 1) * sqrt(N + 3)) * (sqrt(N2) *
SF * PNK[N - 2] - sqrt(N3 + 1) * sqrt(N - 2) / sqrt(N2 - 2) * PNK[N - 3]);
FR += (N3 + 1) * CK1 * A;
FF += CK1 * (PNK[N1] */*SQ[N1]*/sqrt(N1 + 1) */*SQ[N+2]*/sqrt(N + 3) - TG * A);
FL += CK2 * A;
} //ENDIF //3
for (M = 1; M < NK; M++)
{ // 4
J = (int)(ANAI[1 + M]) - 1;
for (N = M; N < NK; N++)
{ //4.1
N1 = N - M;
N2 = N + M + 1;
N3 = N + N + 2;
N4 = N1 - 2;
N5 = N1 - 3;
A = PNK[N1];
AN = (double)(M + 1) * A;
CK1 = AR[N] * (C[J] * CO[M] + S[J] * SI[M]);
CK2 = AR[N] * (C[J] * SI[M] - S[J] * CO[M]);
if (N1 > 2)
PNK[N1 - 1] = sqrt(N3 + 1) / (sqrt(N4 + 1) * sqrt(N2 + 2)) * (sqrt(N3 - 1) *
SF * PNK[N4] - sqrt(N2 + 1) * sqrt(N5 + 1) / sqrt(N3 - 3) * PNK[N5]);
else if (N1 == 0)
{ //4.1.1
FR += HP[N] * CK1 * A;
FF -= CK1 * AN * TG;
FL += CK2 * AN;
goto m12;
} //4.1.1
else if (N1 == 1)
PNK[0] = SK[N] * CF[N];
else if (N1 == 2)
PNK[1] = /*SQ[N3]*/sqrt(N3 + 1) * SF * PNK[0];
FR += HP[N] * CK1 * A;
FF += CK1 * (PNK[N1 - 1] */*SQ[N1-1]*SQ[N2+1]*/sqrt(N1) * sqrt(N2 + 2) - TG * AN);
FL += CK2 * AN;
m12: J++;
} //4.1
} //4
// 5
} //0
m5:
FR = -GR * FR - R[2] * R3;
FF = GR * FF;
FL = -GR / CF[0] * FL;
A = FF * SF;
FX = FR * (X * R1) - A * CO[0] - FL * SI[0];
FY = FR * (Y * R1) - A * SI[0] + FL * CO[0];
FZ = FR * SF + FF * CF[0];
}
private List<object> saver
{
get
{
List<object> l = new List<object>();
l.Add(R);
l.Add(C);
l.Add(S);
l.Add(HP);
l.Add(CO);
l.Add(SI);
l.Add(AR);
l.Add(CF);
l.Add(PNK);
l.Add(ANAI);
l.Add(n0);
l.Add(nk);
return l;
}
set
{
int k = 0;
R = value[k] as double[]; ++k;
C = value[k] as double[]; ++k;
S = value[k] as double[]; ++k;
HP = value[k] as double[]; ++k;
CO = value[k] as double[]; ++k;
SI = value[k] as double[]; ++k;
AR = value[k] as double[]; ++k;
CF = value[k] as double[]; ++k;
PNK = value[k] as double[]; ++k;
ANAI = value[k] as double[]; ++k;
n0 = (int) value[k]; ++k;
nk = (int) value[k]; ++k;
GC.Collect();
}
}
#endregion
}
}
Submit an article or tip