/// <summary>
/// Performs training iteration.
/// </summary>
public void Iteration()
{
//Next epoch
++_epoch;
//Noise intensity
double intensity = _alphas[Math.Min(_maxEpoch, _epoch) - 1];
//Adjusted predictors
Matrix predictors = PreparePredictors((double)intensity);
//Decomposition
QRD decomposition = new QRD(predictors);
//New waights
double[] newWaights = new double[_net.NumOfWeights];
//Regression for each output neuron
for (int outputIdx = 0; outputIdx < _net.NumOfOutputValues; outputIdx++)
{
//Regression
Matrix solution = decomposition.Solve(_outputSingleColMatrixCollection[outputIdx]);
//Store weights
//Input weights
for (int i = 0; i < solution.NumOfRows - 1; i++)
{
newWaights[outputIdx * _net.NumOfInputValues + i] = solution.Data[i][0];
}
//Bias weight
newWaights[_net.NumOfOutputValues * _net.NumOfInputValues + outputIdx] = solution.Data[solution.NumOfRows - 1][0];
}
//Set new weights and compute error
_net.SetWeights(newWaights);
_mse = _net.ComputeBatchErrorStat(_inputVectorCollection, _outputVectorCollection).MeanSquare;
return;
}
/// <summary>
/// Get WHATDATACODEVALUEQUAL Records from QRD
/// </summary>
public static IEnumerable GetWHATDATACODEVALUEQUALRecords(this QRD message)
{
object[] result = message.GetRecords("WHATDATACODEVALUEQUALRepetitionsUsed", "GetWHATDATACODEVALUEQUAL");
if ((result != null) && (result.Count() > 0))
{
for (int i = 0; i < result.Count(); i++)
{
yield return(result[i]);
}
}
}
/// <summary>
/// Get WHOSUBJECTFILTER Records from QRD
/// </summary>
public static IEnumerable GetWHOSUBJECTFILTERRecords(this QRD message)
{
object[] result = message.GetRecords("WHOSUBJECTFILTERRepetitionsUsed", "GetWHOSUBJECTFILTER");
if ((result != null) && (result.Count() > 0))
{
for (int i = 0; i < result.Count(); i++)
{
yield return(result[i]);
}
}
}
/// <summary>
/// Performs training iteration.
/// </summary>
public void Iteration()
{
//Next epoch
++_epoch;
//Noise intensity
double intensity = _alphas[Math.Min(_maxEpoch, _epoch) - 1];
//Adjusted predictors
Matrix predictors = PreparePredictors((double)intensity);
//Decomposition
QRD decomposition = null;
bool useableQRD = true;
try
{
//Try to create QRD. Any exception signals numerical unstability
decomposition = new QRD(predictors);
}
catch
{
//Creation of QRD object throws exception. QRD object is not ready for use.
useableQRD = false;
if (_epoch == 1)
{
//No previous successful epoch so stop training
throw;
}
}
if (useableQRD)
{
//QRD is ready for use (low probability of numerical unstability)
//New weights
double[] newWeights = new double[_net.NumOfWeights];
//Regression for each output neuron
for (int outputIdx = 0; outputIdx < _net.NumOfOutputValues; outputIdx++)
{
//Regression
Matrix solution = decomposition.Solve(_outputSingleColMatrixCollection[outputIdx]);
//Store weights
//Input weights
for (int i = 0; i < solution.NumOfRows - 1; i++)
{
newWeights[outputIdx * _net.NumOfInputValues + i] = solution.Data[i][0];
}
//Bias weight
newWeights[_net.NumOfOutputValues * _net.NumOfInputValues + outputIdx] = solution.Data[solution.NumOfRows - 1][0];
}
//Set new weights and compute error
_net.SetWeights(newWeights);
_mse = _net.ComputeBatchErrorStat(_inputVectorCollection, _outputVectorCollection).MeanSquare;
}
return;
}
/// <summary>
/// Constructs an HL7 segment for the specified Segments enum object
/// </summary>
/// <param name="seg">The Segments enum object to construct for</param>
public Segment(Segments seg)
{
switch (seg)
{
case Segments.ABS:
ABS abs = new ABS();
Name = abs.Name;
Description = abs.Description;
Fields = abs.Fields;
break;
case Segments.ACC:
ACC acc = new ACC();
Name = acc.Name;
Description = acc.Description;
Fields = acc.Fields;
break;
case Segments.ADD:
ADD add = new ADD();
Name = add.Name;
Description = add.Description;
Fields = add.Fields;
break;
case Segments.AFF:
AFF aff = new AFF();
Name = aff.Name;
Description = aff.Description;
Fields = aff.Fields;
break;
case Segments.AIG:
AIG aig = new AIG();
Name = aig.Name;
Description = aig.Description;
Fields = aig.Fields;
break;
case Segments.AIL:
AIL ail = new AIL();
Name = ail.Name;
Description = ail.Description;
Fields = ail.Fields;
break;
case Segments.AIP:
AIP aip = new AIP();
Name = aip.Name;
Description = aip.Description;
Fields = aip.Fields;
break;
case Segments.AIS:
AIS ais = new AIS();
Name = ais.Name;
Description = ais.Description;
Fields = ais.Fields;
break;
case Segments.AL1:
AL1 al1 = new AL1();
Name = al1.Name;
Description = al1.Description;
Fields = al1.Fields;
break;
case Segments.APR:
APR apr = new APR();
Name = apr.Name;
Description = apr.Description;
Fields = apr.Fields;
break;
case Segments.ARQ:
ARQ arq = new ARQ();
Name = arq.Name;
Description = arq.Description;
Fields = arq.Fields;
break;
case Segments.AUT:
AUT aut = new AUT();
Name = aut.Name;
Description = aut.Description;
Fields = aut.Fields;
break;
case Segments.BHS:
BHS bhs = new BHS();
Name = bhs.Name;
Description = bhs.Description;
Fields = bhs.Fields;
break;
case Segments.BLC:
BLC blc = new BLC();
Name = blc.Name;
Description = blc.Description;
Fields = blc.Fields;
break;
case Segments.BLG:
BLG blg = new BLG();
Name = blg.Name;
Description = blg.Description;
Fields = blg.Fields;
break;
case Segments.BPO:
BPO bpo = new BPO();
Name = bpo.Name;
Description = bpo.Description;
Fields = bpo.Fields;
break;
case Segments.BPX:
BPX bpx = new BPX();
Name = bpx.Name;
Description = bpx.Description;
Fields = bpx.Fields;
break;
case Segments.BTS:
BTS bts = new BTS();
Name = bts.Name;
Description = bts.Description;
Fields = bts.Fields;
break;
case Segments.BTX:
BTX btx = new BTX();
Name = btx.Name;
Description = btx.Description;
Fields = btx.Fields;
break;
case Segments.CDM:
CDM cdm = new CDM();
Name = cdm.Name;
Description = cdm.Description;
Fields = cdm.Fields;
break;
case Segments.CER:
CER cer = new CER();
Name = cer.Name;
Description = cer.Description;
Fields = cer.Fields;
break;
case Segments.CM0:
CM0 cm0 = new CM0();
Name = cm0.Name;
Description = cm0.Description;
Fields = cm0.Fields;
break;
case Segments.CM1:
CM1 cm1 = new CM1();
Name = cm1.Name;
Description = cm1.Description;
Fields = cm1.Fields;
break;
case Segments.CM2:
CM2 cm2 = new CM2();
Name = cm2.Name;
Description = cm2.Description;
Fields = cm2.Fields;
break;
case Segments.CNS:
CNS cns = new CNS();
Name = cns.Name;
Description = cns.Description;
Fields = cns.Fields;
break;
case Segments.CON:
CON con = new CON();
Name = con.Name;
Description = con.Description;
Fields = con.Fields;
break;
case Segments.CSP:
CSP csp = new CSP();
Name = csp.Name;
Description = csp.Description;
Fields = csp.Fields;
break;
case Segments.CSR:
CSR csr = new CSR();
Name = csr.Name;
Description = csr.Description;
Fields = csr.Fields;
break;
case Segments.CSS:
CSS css = new CSS();
Name = css.Name;
Description = css.Description;
Fields = css.Fields;
break;
case Segments.CTD:
CTD ctd = new CTD();
Name = ctd.Name;
Description = ctd.Description;
Fields = ctd.Fields;
break;
case Segments.CTI:
CTI cti = new CTI();
Name = cti.Name;
Description = cti.Description;
Fields = cti.Fields;
break;
case Segments.DB1:
DB1 db1 = new DB1();
Name = db1.Name;
Description = db1.Description;
Fields = db1.Fields;
break;
case Segments.DG1:
DG1 dg1 = new DG1();
Name = dg1.Name;
Description = dg1.Description;
Fields = dg1.Fields;
break;
case Segments.DRG:
DRG drg = new DRG();
Name = drg.Name;
Description = drg.Description;
Fields = drg.Fields;
break;
case Segments.DSC:
DSC dsc = new DSC();
Name = dsc.Name;
Description = dsc.Description;
Fields = dsc.Fields;
break;
case Segments.DSP:
DSP dsp = new DSP();
Name = dsp.Name;
Description = dsp.Description;
Fields = dsp.Fields;
break;
case Segments.ECD:
ECD ecd = new ECD();
Name = ecd.Name;
Description = ecd.Description;
Fields = ecd.Fields;
break;
case Segments.ECR:
ECR ecr = new ECR();
Name = ecr.Name;
Description = ecr.Description;
Fields = ecr.Fields;
break;
case Segments.EDU:
EDU edu = new EDU();
Name = edu.Name;
Description = edu.Description;
Fields = edu.Fields;
break;
case Segments.EQL:
EQL eql = new EQL();
Name = eql.Name;
Description = eql.Description;
Fields = eql.Fields;
break;
case Segments.EQP:
EQP eqp = new EQP();
Name = eqp.Name;
Description = eqp.Description;
Fields = eqp.Fields;
break;
case Segments.EQU:
EQU equ = new EQU();
Name = equ.Name;
Description = equ.Description;
Fields = equ.Fields;
break;
case Segments.ERQ:
ERQ erq = new ERQ();
Name = erq.Name;
Description = erq.Description;
Fields = erq.Fields;
break;
case Segments.ERR:
ERR err = new ERR();
Name = err.Name;
Description = err.Description;
Fields = err.Fields;
break;
case Segments.EVN:
EVN evn = new EVN();
Name = evn.Name;
Description = evn.Description;
Fields = evn.Fields;
break;
case Segments.FAC:
FAC fac = new FAC();
Name = fac.Name;
Description = fac.Description;
Fields = fac.Fields;
break;
case Segments.FHS:
FHS fhs = new FHS();
Name = fhs.Name;
Description = fhs.Description;
Fields = fhs.Fields;
break;
case Segments.FT1:
FT1 ft1 = new FT1();
Name = ft1.Name;
Description = ft1.Description;
Fields = ft1.Fields;
break;
case Segments.FTS:
FTS fts = new FTS();
Name = fts.Name;
Description = fts.Description;
Fields = fts.Fields;
break;
case Segments.GOL:
GOL gol = new GOL();
Name = gol.Name;
Description = gol.Description;
Fields = gol.Fields;
break;
case Segments.GP1:
GP1 gp1 = new GP1();
Name = gp1.Name;
Description = gp1.Description;
Fields = gp1.Fields;
break;
case Segments.GP2:
GP2 gp2 = new GP2();
Name = gp2.Name;
Description = gp2.Description;
Fields = gp2.Fields;
break;
case Segments.GT1:
GT1 gt1 = new GT1();
Name = gt1.Name;
Description = gt1.Description;
Fields = gt1.Fields;
break;
case Segments.IAM:
IAM iam = new IAM();
Name = iam.Name;
Description = iam.Description;
Fields = iam.Fields;
break;
case Segments.IIM:
IIM iim = new IIM();
Name = iim.Name;
Description = iim.Description;
Fields = iim.Fields;
break;
case Segments.IN1:
IN1 in1 = new IN1();
Name = in1.Name;
Description = in1.Description;
Fields = in1.Fields;
break;
case Segments.IN2:
IN2 in2 = new IN2();
Name = in2.Name;
Description = in2.Description;
Fields = in2.Fields;
break;
case Segments.IN3:
IN3 in3 = new IN3();
Name = in3.Name;
Description = in3.Description;
Fields = in3.Fields;
break;
case Segments.INV:
INV inv = new INV();
Name = inv.Name;
Description = inv.Description;
Fields = inv.Fields;
break;
case Segments.IPC:
IPC ipc = new IPC();
Name = ipc.Name;
Description = ipc.Description;
Fields = ipc.Fields;
break;
case Segments.ISD:
ISD isd = new ISD();
Name = isd.Name;
Description = isd.Description;
Fields = isd.Fields;
break;
case Segments.LAN:
LAN lan = new LAN();
Name = lan.Name;
Description = lan.Description;
Fields = lan.Fields;
break;
case Segments.LCC:
LCC lcc = new LCC();
Name = lcc.Name;
Description = lcc.Description;
Fields = lcc.Fields;
break;
case Segments.LCH:
LCH lch = new LCH();
Name = lch.Name;
Description = lch.Description;
Fields = lch.Fields;
break;
case Segments.LDP:
LDP ldp = new LDP();
Name = ldp.Name;
Description = ldp.Description;
Fields = ldp.Fields;
break;
case Segments.LOC:
LOC loc = new LOC();
Name = loc.Name;
Description = loc.Description;
Fields = loc.Fields;
break;
case Segments.LRL:
LRL lrl = new LRL();
Name = lrl.Name;
Description = lrl.Description;
Fields = lrl.Fields;
break;
case Segments.MFA:
MFA mfa = new MFA();
Name = mfa.Name;
Description = mfa.Description;
Fields = mfa.Fields;
break;
case Segments.MFE:
MFE mfe = new MFE();
Name = mfe.Name;
Description = mfe.Description;
Fields = mfe.Fields;
break;
case Segments.MFI:
MFI mfi = new MFI();
Name = mfi.Name;
Description = mfi.Description;
Fields = mfi.Fields;
break;
case Segments.MRG:
MRG mrg = new MRG();
Name = mrg.Name;
Description = mrg.Description;
Fields = mrg.Fields;
break;
case Segments.MSA:
MSA msa = new MSA();
Name = msa.Name;
Description = msa.Description;
Fields = msa.Fields;
break;
case Segments.MSH:
MSH msh = new MSH();
Name = msh.Name;
Description = msh.Description;
Fields = msh.Fields;
break;
case Segments.NCK:
NCK nck = new NCK();
Name = nck.Name;
Description = nck.Description;
Fields = nck.Fields;
break;
case Segments.NDS:
NDS nds = new NDS();
Name = nds.Name;
Description = nds.Description;
Fields = nds.Fields;
break;
case Segments.NK1:
NK1 nk1 = new NK1();
Name = nk1.Name;
Description = nk1.Description;
Fields = nk1.Fields;
break;
case Segments.NPU:
NPU npu = new NPU();
Name = npu.Name;
Description = npu.Description;
Fields = npu.Fields;
break;
case Segments.NSC:
NSC nsc = new NSC();
Name = nsc.Name;
Description = nsc.Description;
Fields = nsc.Fields;
break;
case Segments.NST:
NST nst = new NST();
Name = nst.Name;
Description = nst.Description;
Fields = nst.Fields;
break;
case Segments.NTE:
NTE nte = new NTE();
Name = nte.Name;
Description = nte.Description;
Fields = nte.Fields;
break;
case Segments.OBR:
OBR obr = new OBR();
Name = obr.Name;
Description = obr.Description;
Fields = obr.Fields;
break;
case Segments.OBX:
OBX obx = new OBX();
Name = obx.Name;
Description = obx.Description;
Fields = obx.Fields;
break;
case Segments.ODS:
ODS ods = new ODS();
Name = ods.Name;
Description = ods.Description;
Fields = ods.Fields;
break;
case Segments.ODT:
ODT odt = new ODT();
Name = odt.Name;
Description = odt.Description;
Fields = odt.Fields;
break;
case Segments.OM1:
OM1 om1 = new OM1();
Name = om1.Name;
Description = om1.Description;
Fields = om1.Fields;
break;
case Segments.OM2:
OM2 om2 = new OM2();
Name = om2.Name;
Description = om2.Description;
Fields = om2.Fields;
break;
case Segments.OM3:
OM3 om3 = new OM3();
Name = om3.Name;
Description = om3.Description;
Fields = om3.Fields;
break;
case Segments.OM4:
OM4 om4 = new OM4();
Name = om4.Name;
Description = om4.Description;
Fields = om4.Fields;
break;
case Segments.OM5:
OM5 om5 = new OM5();
Name = om5.Name;
Description = om5.Description;
Fields = om5.Fields;
break;
case Segments.OM6:
OM6 om6 = new OM6();
Name = om6.Name;
Description = om6.Description;
Fields = om6.Fields;
break;
case Segments.OM7:
OM7 om7 = new OM7();
Name = om7.Name;
Description = om7.Description;
Fields = om7.Fields;
break;
case Segments.ORC:
ORC orc = new ORC();
Name = orc.Name;
Description = orc.Description;
Fields = orc.Fields;
break;
case Segments.ORG:
ORG org = new ORG();
Name = org.Name;
Description = org.Description;
Fields = org.Fields;
break;
case Segments.OVR:
OVR ovr = new OVR();
Name = ovr.Name;
Description = ovr.Description;
Fields = ovr.Fields;
break;
case Segments.PCR:
PCR pcr = new PCR();
Name = pcr.Name;
Description = pcr.Description;
Fields = pcr.Fields;
break;
case Segments.PD1:
PD1 pd1 = new PD1();
Name = pd1.Name;
Description = pd1.Description;
Fields = pd1.Fields;
break;
case Segments.PDA:
PDA pda = new PDA();
Name = pda.Name;
Description = pda.Description;
Fields = pda.Fields;
break;
case Segments.PDC:
PDC pdc = new PDC();
Name = pdc.Name;
Description = pdc.Description;
Fields = pdc.Fields;
break;
case Segments.PEO:
PEO peo = new PEO();
Name = peo.Name;
Description = peo.Description;
Fields = peo.Fields;
break;
case Segments.PES:
PES pes = new PES();
Name = pes.Name;
Description = pes.Description;
Fields = pes.Fields;
break;
case Segments.PID:
PID pid = new PID();
Name = pid.Name;
Description = pid.Description;
Fields = pid.Fields;
break;
case Segments.PR1:
PR1 pr1 = new PR1();
Name = pr1.Name;
Description = pr1.Description;
Fields = pr1.Fields;
break;
case Segments.PRA:
PRA pra = new PRA();
Name = pra.Name;
Description = pra.Description;
Fields = pra.Fields;
break;
case Segments.PRB:
PRB prb = new PRB();
Name = prb.Name;
Description = prb.Description;
Fields = prb.Fields;
break;
case Segments.PRC:
PRC prc = new PRC();
Name = prc.Name;
Description = prc.Description;
Fields = prc.Fields;
break;
case Segments.PRD:
PRD prd = new PRD();
Name = prd.Name;
Description = prd.Description;
Fields = prd.Fields;
break;
case Segments.PSH:
PSH psh = new PSH();
Name = psh.Name;
Description = psh.Description;
Fields = psh.Fields;
break;
case Segments.PTH:
PTH pth = new PTH();
Name = pth.Name;
Description = pth.Description;
Fields = pth.Fields;
break;
case Segments.PV1:
PV1 pv1 = new PV1();
Name = pv1.Name;
Description = pv1.Description;
Fields = pv1.Fields;
break;
case Segments.PV2:
PV2 pv2 = new PV2();
Name = pv2.Name;
Description = pv2.Description;
Fields = pv2.Fields;
break;
case Segments.QAK:
QAK qak = new QAK();
Name = qak.Name;
Description = qak.Description;
Fields = qak.Fields;
break;
case Segments.QID:
QID qid = new QID();
Name = qid.Name;
Description = qid.Description;
Fields = qid.Fields;
break;
case Segments.QPD:
QPD qpd = new QPD();
Name = qpd.Name;
Description = qpd.Description;
Fields = qpd.Fields;
break;
case Segments.QRD:
QRD qrd = new QRD();
Name = qrd.Name;
Description = qrd.Description;
Fields = qrd.Fields;
break;
case Segments.QRF:
QRF qrf = new QRF();
Name = qrf.Name;
Description = qrf.Description;
Fields = qrf.Fields;
break;
case Segments.QRI:
QRI qri = new QRI();
Name = qri.Name;
Description = qri.Description;
Fields = qri.Fields;
break;
case Segments.RCP:
RCP rcp = new RCP();
Name = rcp.Name;
Description = rcp.Description;
Fields = rcp.Fields;
break;
case Segments.RDF:
RDF rdf = new RDF();
Name = rdf.Name;
Description = rdf.Description;
Fields = rdf.Fields;
break;
case Segments.RF1:
RF1 rf1 = new RF1();
Name = rf1.Name;
Description = rf1.Description;
Fields = rf1.Fields;
break;
case Segments.RGS:
RGS rgs = new RGS();
Name = rgs.Name;
Description = rgs.Description;
Fields = rgs.Fields;
break;
case Segments.RMI:
RMI rmi = new RMI();
Name = rmi.Name;
Description = rmi.Description;
Fields = rmi.Fields;
break;
case Segments.ROL:
ROL rol = new ROL();
Name = rol.Name;
Description = rol.Description;
Fields = rol.Fields;
break;
case Segments.RQ1:
RQ1 rq1 = new RQ1();
Name = rq1.Name;
Description = rq1.Description;
Fields = rq1.Fields;
break;
case Segments.RQD:
RQD rqd = new RQD();
Name = rqd.Name;
Description = rqd.Description;
Fields = rqd.Fields;
break;
case Segments.RXA:
RXA rxa = new RXA();
Name = rxa.Name;
Description = rxa.Description;
Fields = rxa.Fields;
break;
case Segments.RXC:
RXC rxc = new RXC();
Name = rxc.Name;
Description = rxc.Description;
Fields = rxc.Fields;
break;
case Segments.RXD:
RXD rxd = new RXD();
Name = rxd.Name;
Description = rxd.Description;
Fields = rxd.Fields;
break;
case Segments.RXE:
RXE rxe = new RXE();
Name = rxe.Name;
Description = rxe.Description;
Fields = rxe.Fields;
break;
case Segments.RXG:
RXG rxg = new RXG();
Name = rxg.Name;
Description = rxg.Description;
Fields = rxg.Fields;
break;
case Segments.RXO:
RXO rxo = new RXO();
Name = rxo.Name;
Description = rxo.Description;
Fields = rxo.Fields;
break;
case Segments.RXR:
RXR rxr = new RXR();
Name = rxr.Name;
Description = rxr.Description;
Fields = rxr.Fields;
break;
case Segments.SAC:
SAC sac = new SAC();
Name = sac.Name;
Description = sac.Description;
Fields = sac.Fields;
break;
case Segments.SCH:
SCH sch = new SCH();
Name = sch.Name;
Description = sch.Description;
Fields = sch.Fields;
break;
case Segments.SFT:
SFT sft = new SFT();
Name = sft.Name;
Description = sft.Description;
Fields = sft.Fields;
break;
case Segments.SID:
SID sid = new SID();
Name = sid.Name;
Description = sid.Description;
Fields = sid.Fields;
break;
case Segments.SPM:
SPM spm = new SPM();
Name = spm.Name;
Description = spm.Description;
Fields = spm.Fields;
break;
case Segments.SPR:
SPR spr = new SPR();
Name = spr.Name;
Description = spr.Description;
Fields = spr.Fields;
break;
case Segments.STF:
STF stf = new STF();
Name = stf.Name;
Description = stf.Description;
Fields = stf.Fields;
break;
case Segments.TCC:
TCC tcc = new TCC();
Name = tcc.Name;
Description = tcc.Description;
Fields = tcc.Fields;
break;
case Segments.TCD:
TCD tcd = new TCD();
Name = tcd.Name;
Description = tcd.Description;
Fields = tcd.Fields;
break;
case Segments.TQ1:
TQ1 tq1 = new TQ1();
Name = tq1.Name;
Description = tq1.Description;
Fields = tq1.Fields;
break;
case Segments.TQ2:
TQ2 tq2 = new TQ2();
Name = tq2.Name;
Description = tq2.Description;
Fields = tq2.Fields;
break;
case Segments.TXA:
TXA txa = new TXA();
Name = txa.Name;
Description = txa.Description;
Fields = txa.Fields;
break;
case Segments.UB1:
UB1 ub1 = new UB1();
Name = ub1.Name;
Description = ub1.Description;
Fields = ub1.Fields;
break;
case Segments.UB2:
UB2 ub2 = new UB2();
Name = ub2.Name;
Description = ub2.Description;
Fields = ub2.Fields;
break;
case Segments.URD:
URD urd = new URD();
Name = urd.Name;
Description = urd.Description;
Fields = urd.Fields;
break;
case Segments.URS:
URS urs = new URS();
Name = urs.Name;
Description = urs.Description;
Fields = urs.Fields;
break;
case Segments.VAR:
VAR var = new VAR();
Name = var.Name;
Description = var.Description;
Fields = var.Fields;
break;
case Segments.VTQ:
VTQ vtq = new VTQ();
Name = vtq.Name;
Description = vtq.Description;
Fields = vtq.Fields;
break;
}
}
public static double[,] Solve(this double[,] matrix, double[,] b)
{
return(matrix.IsSquare() ? (LUD.ComputeLUD(matrix)).Solve(b) : (QRD.ComputeQRD(matrix)).Solve(b));
}
/// <summary>
/// Add a new WHATDATACODEVALUEQUAL to QRD
/// </summary>
public static ST AddWHATDATACODEVALUEQUAL(this QRD message)
{
return(message.GetWHATDATACODEVALUEQUAL(message.WHATDATACODEVALUEQUALRepetitionsUsed));
}
/// <summary>
/// Get all WHATDATACODEVALUEQUAL Records from QRD
/// </summary>
public static List <ST> GetAllWHATDATACODEVALUEQUALRecords(this QRD message)
{
return(message.GetAllRecords <ST>("WHATDATACODEVALUEQUALRepetitionsUsed", "GetWHATDATACODEVALUEQUAL"));
}
/// <summary>
/// Add a new WHATDEPARTMENTDATACODE to QRD
/// </summary>
public static ST AddWHATDEPARTMENTDATACODE(this QRD message)
{
return(message.GetWHATDEPARTMENTDATACODE(message.WHATDEPARTMENTDATACODERepetitionsUsed));
}
/// <summary>
/// Get all WHATDEPARTMENTDATACODE Records from QRD
/// </summary>
public static List <ST> GetAllWHATDEPARTMENTDATACODERecords(this QRD message)
{
return(message.GetAllRecords <ST>("WHATDEPARTMENTDATACODERepetitionsUsed", "GetWHATDEPARTMENTDATACODE"));
}
/// <summary>
/// Add a new WHATSUBJECTFILTER to QRD
/// </summary>
public static ID AddWHATSUBJECTFILTER(this QRD message)
{
return(message.GetWHATSUBJECTFILTER(message.WHATSUBJECTFILTERRepetitionsUsed));
}
/// <summary>
/// Get all WHATSUBJECTFILTER Records from QRD
/// </summary>
public static List <ID> GetAllWHATSUBJECTFILTERRecords(this QRD message)
{
return(message.GetAllRecords <ID>("WHATSUBJECTFILTERRepetitionsUsed", "GetWHATSUBJECTFILTER"));
}
/// <summary>
/// Performs training iteration.
/// </summary>
public bool Iteration()
{
//Fetch next noise intensity
double newMaxNoise = _maxNoiseSeeker.Next;
//Check continue conditions
if (AttemptEpoch == MaxAttemptEpoch || Math.Abs(_currMaxNoise - newMaxNoise) < StopNoiseDifference)
{
//Try new attempt
if (!NextAttempt())
{
//Next attempt is not available
return(false);
}
else
{
newMaxNoise = _maxNoiseSeeker.Next;
}
}
//Next epoch
++AttemptEpoch;
_currMaxNoise = newMaxNoise;
InfoMessage = $"maxNoise={_currMaxNoise.ToString(CultureInfo.InvariantCulture)}";
//Adjusted predictors
Matrix predictors = PreparePredictors(_currMaxNoise);
//Decomposition
QRD decomposition = null;
bool useableQRD = true;
try
{
//Try to create QRD. Any exception signals numerical unstability
decomposition = new QRD(predictors);
}
catch
{
//Creation of QRD object throws exception. QRD object is not ready for use.
useableQRD = false;
if (AttemptEpoch == 1)
{
//No previous successful epoch so stop training
throw;
}
}
if (useableQRD)
{
//QRD is ready for use (low probability of numerical unstability)
//New weights
double[] newWeights = new double[_net.NumOfWeights];
//Regression for each output neuron
for (int outputIdx = 0; outputIdx < _net.NumOfOutputValues; outputIdx++)
{
//Regression
Matrix solution = decomposition.Solve(_outputSingleColMatrixCollection[outputIdx]);
//Store weights
//Input weights
for (int i = 0; i < solution.NumOfRows - 1; i++)
{
newWeights[outputIdx * _net.NumOfInputValues + i] = solution.Data[i][0];
}
//Bias weight
newWeights[_net.NumOfOutputValues * _net.NumOfInputValues + outputIdx] = solution.Data[solution.NumOfRows - 1][0];
}
//Set new weights and compute error
_net.SetWeights(newWeights);
MSE = _net.ComputeBatchErrorStat(_inputVectorCollection, _outputVectorCollection).MeanSquare;
}
_maxNoiseSeeker.ProcessError(MSE);
return(true);
}
//Methods
public void Run()
{
//Linear algebra test
double[] flatData =
{
0.2, 5, 17.3, 1.01, 54, 7,
2.2, 5.5, 12.13, 11.57, 5.71, -85,
-70.1, 15, -18.3, 0.3, 42, -6.25,
0.042, 1, 7.75, -81.01, -21.29, 5.44,
0.1, 4, -4.3, 18.01, 7.12, -3.14,
-80.1, 24.4, 4.3, 12.03, 2.789, -13
};
Matrix testM = new Matrix(6, 6, flatData);
/*
* //Inversion test
* Matrix resultM = new Matrix(testM);
* resultM.SingleThreadInverse();
*/
/*
* //Transpose test
* Matrix resultM = testM.Transpose();
*/
/*
* //Multiply test
* Matrix resultM = Matrix.Multiply(testM, testM);
* for (int i = 0; i < resultM.NumOfRows; i++)
* {
* Console.WriteLine($"{resultM.Data[i][0]}; {resultM.Data[i][1]}; {resultM.Data[i][2]}; {resultM.Data[i][3]}; {resultM.Data[i][4]}; {resultM.Data[i][5]}");
* }
*/
;
int numOfweights = 4;
int xIdx, dIdx = 0;
double[][] data = new double[5][];
data[dIdx] = new double[numOfweights];
xIdx = -1;
data[dIdx][++xIdx] = 1; //Bias
data[dIdx][++xIdx] = 2;
data[dIdx][++xIdx] = 1;
data[dIdx][++xIdx] = 3;
++dIdx;
data[dIdx] = new double[numOfweights];
xIdx = -1;
data[dIdx][++xIdx] = 1; //Bias
data[dIdx][++xIdx] = 1;
data[dIdx][++xIdx] = 3;
data[dIdx][++xIdx] = -3;
++dIdx;
data[dIdx] = new double[numOfweights];
xIdx = -1;
data[dIdx][++xIdx] = 1; //Bias
data[dIdx][++xIdx] = -2;
data[dIdx][++xIdx] = 4;
data[dIdx][++xIdx] = 4;
++dIdx;
data[dIdx] = new double[numOfweights];
xIdx = -1;
data[dIdx][++xIdx] = 1; //Bias
data[dIdx][++xIdx] = -5;
data[dIdx][++xIdx] = 7;
data[dIdx][++xIdx] = 6;
++dIdx;
data[dIdx] = new double[numOfweights];
xIdx = -1;
data[dIdx][++xIdx] = 1; //Bias
data[dIdx][++xIdx] = 1;
data[dIdx][++xIdx] = 12;
data[dIdx][++xIdx] = 5;
Matrix M = new Matrix(data, true);
Vector desired = new Vector(5);
dIdx = -1;
desired.Data[++dIdx] = 8;
desired.Data[++dIdx] = 13;
desired.Data[++dIdx] = 5;
desired.Data[++dIdx] = 7;
desired.Data[++dIdx] = 10;
Vector weights = Matrix.RidgeRegression(M, desired, 0);
//Display results
for (int i = 0; i < data.Length; i++)
{
double result = 0;
for (int j = 0; j < weights.Length; j++)
{
result += data[i][j] * weights[j];
}
Console.WriteLine($"Computed {result}, Desired {desired.Data[i]}");
}
for (int i = 0; i < weights.Length; i++)
{
Console.WriteLine($"Weight[{i}] = {weights[i]}");
}
Console.WriteLine();
//QRD
QRD decomposition = new QRD(M);
Matrix B = new Matrix(desired.Length, 1, desired.Data);
Matrix W = decomposition.Solve(B);
//Display results
for (int i = 0; i < data.Length; i++)
{
double result = 0;
for (int j = 0; j < W.Data.Length; j++)
{
result += data[i][j] * W.Data[j][0];
}
Console.WriteLine($"Computed {result}, Desired {desired.Data[i]}");
}
for (int i = 0; i < W.Data.Length; i++)
{
Console.WriteLine($"Weight[{i}] = {W.Data[i][0]}");
}
;
//TimeSeriesGenerator.SaveTimeSeriesToCsvFile("MackeyGlass_big.csv", "Value", TimeSeriesGenerator.GenMackeyGlassTimeSeries(16000), CultureInfo.InvariantCulture);
MackeyGlassGeneratorSettings modSettings = new MackeyGlassGeneratorSettings(18, 0.1, 0.2);
IGenerator generator = new MackeyGlassGenerator(modSettings);
int steps = 100;
for (int i = 0; i < steps; i++)
{
Console.WriteLine(generator.Next());
}
Console.ReadLine();
generator.Reset();
for (int i = 0; i < steps; i++)
{
Console.WriteLine(generator.Next());
}
Console.ReadLine();
///*
SimpleIFSettings settings = new SimpleIFSettings(1,
new RandomValueSettings(15, 15),
new RandomValueSettings(0.05, 0.05),
new RandomValueSettings(5, 5),
new RandomValueSettings(20, 20),
0
);
IActivationFunction af = ActivationFactory.Create(settings, new Random(0));
//*/
TestActivation(af, 800, 0.15, 10, 600);
return;
}