/*
* MATLAB Compiler: 3.0
* Date: Wed Jul 06 02:27:15 2005
* Arguments: "-B" "macro_default" "-O" "all" "-O" "fold_scalar_mxarrays:on"
* "-O" "fold_non_scalar_mxarrays:on" "-O" "optimize_integer_for_loops:on" "-O"
* "array_indexing:on" "-O" "optimize_conditionals:on" "-m" "-W" "main" "-L"
* "C" "-t" "-T" "link:exe" "-h" "libmmfile.mlib" "-B" "sgl" "-m" "-W" "main"
* "-L" "C" "-t" "-T" "link:exe" "-h" "libmmfile.mlib" "-W" "mainhg"
* "libmwsglm.mlib" "micampo.m" "matrices1lazo.m" "leeActivo.m" "leePasivo.m"
* "matActivo.m"
*/
#include "matactivo.h"
#include "libmatlbm.h"
static mxArray * _mxarray0_;
static mxArray * _mxarray1_;
static mxArray * _mxarray2_;
static mxArray * _mxarray3_;
static mxArray * _mxarray4_;
static mxArray * _mxarray5_;
static mxArray * _mxarray6_;
static mxArray * _mxarray7_;
static mxArray * _mxarray8_;
static mxArray * _mxarray9_;
static mxArray * _mxarray10_;
void InitializeModule_matactivo(void) {
_mxarray0_ = mclInitializeDouble(0.0);
_mxarray1_ = mclInitializeDouble(1.0);
_mxarray2_ = mclInitializeDoubleVector(0, 0, (double *)NULL);
_mxarray3_ = mclInitializeDouble(2.0);
_mxarray4_ = mclInitializeDouble(3.141592653589793);
_mxarray5_ = mclInitializeDouble(180.0);
_mxarray6_ = mclInitializeComplex(0.0, 1.0);
_mxarray7_ = mclInitializeDouble(6.0);
_mxarray8_ = mclInitializeDouble(3.0);
_mxarray9_ = mclInitializeDouble(8.0);
_mxarray10_ = mclInitializeDouble(4.0);
}
void TerminateModule_matactivo(void) {
mxDestroyArray(_mxarray10_);
mxDestroyArray(_mxarray9_);
mxDestroyArray(_mxarray8_);
mxDestroyArray(_mxarray7_);
mxDestroyArray(_mxarray6_);
mxDestroyArray(_mxarray5_);
mxDestroyArray(_mxarray4_);
mxDestroyArray(_mxarray3_);
mxDestroyArray(_mxarray2_);
mxDestroyArray(_mxarray1_);
mxDestroyArray(_mxarray0_);
}
static mxArray * Mmatactivo(mxArray * * mislazos,
int nargout_,
mxArray * coordlazos,
mxArray * mislazoss);
_mexLocalFunctionTable _local_function_table_matactivo
= { 0, (mexFunctionTableEntry *)NULL };
/*
* The function "mlfMatactivo" contains the normal interface for the
* "matactivo" M-function from file "c:\matlab6p5\work\matactivo.m" (lines
* 1-109). This function processes any input arguments and passes them to the
* implementation version of the function, appearing above.
*/
mxArray * mlfMatactivo(mxArray * * mislazos,
mxArray * coordlazos,
mxArray * mislazoss) {
int nargout = 1;
mxArray * I = NULL;
mxArray * mislazos__ = NULL;
mlfEnterNewContext(1, 2, mislazos, coordlazos, mislazoss);
if (mislazos != NULL) {
++nargout;
}
I = Mmatactivo(&mislazos__, nargout, coordlazos, mislazoss);
mlfRestorePreviousContext(1, 2, mislazos, coordlazos, mislazoss);
if (mislazos != NULL) {
mclCopyOutputArg(mislazos, mislazos__);
} else {
mxDestroyArray(mislazos__);
}
return mlfReturnValue(I);
}
/*
* The function "mlxMatactivo" contains the feval interface for the "matactivo"
* M-function from file "c:\matlab6p5\work\matactivo.m" (lines 1-109). The
* feval function calls the implementation version of matactivo through this
* function. This function processes any input arguments and passes them to the
* implementation version of the function, appearing above.
*/
void mlxMatactivo(int nlhs, mxArray * plhs[], int nrhs, mxArray * prhs[]) {
mxArray * mprhs[2];
mxArray * mplhs[2];
int i;
if (nlhs > 2) {
mlfError(
mxCreateString(
"Run-time Error: File: matactivo Line: 19 Column"
": 1 The function \"matactivo\" was called with "
"more than the declared number of outputs (2)."),
NULL);
}
if (nrhs > 2) {
mlfError(
mxCreateString(
"Run-time Error: File: matactivo Line: 19 Column"
": 1 The function \"matactivo\" was called with "
"more than the declared number of inputs (2)."),
NULL);
}
for (i = 0; i < 2; ++i) {
mplhs[i] = NULL;
}
for (i = 0; i < 2 && i < nrhs; ++i) {
mprhs[i] = prhs[i];
}
for (; i < 2; ++i) {
mprhs[i] = NULL;
}
mlfEnterNewContext(0, 2, mprhs[0], mprhs[1]);
mplhs[0] = Mmatactivo(&mplhs[1], nlhs, mprhs[0], mprhs[1]);
mlfRestorePreviousContext(0, 2, mprhs[0], mprhs[1]);
plhs[0] = mplhs[0];
for (i = 1; i < 2 && i < nlhs; ++i) {
plhs[i] = mplhs[i];
}
for (; i < 2; ++i) {
mxDestroyArray(mplhs[i]);
}
}
/*
* The function "Mmatactivo" is the implementation version of the "matactivo"
* M-function from file "c:\matlab6p5\work\matactivo.m" (lines 1-109). It
* contains the actual compiled code for that M-function. It is a static
* function and must only be called from one of the interface functions,
* appearing below.
*/
/*
* % Calcula la descomposicion en lazos simples de una serie
* % de lazos.
* %
* % Sintaxis: [I,mislazos]=matActivo(lazos)
* %
* % lazos es una matriz de dimensiones (m,4) donde cada fila representa
* % a un lazo y las columnas son, por este orden:
* %
* % columna 1 -> Coordenada x del conductor 1 del lazo.
* % columna 2 -> Coordenada y del conductor 1 del lazo.
* % columna 3 -> Coordenada x del conductor 2 del lazo.
* % columna 4 -> Coordenada y del conductor 2 del lazo.
* % columna 5 -> Intensidad por el lazo 1.
* % columna 6 -> Intensidad por el lazo 2.
* %
* % mislazos es una matriz que contiene los lazos originales
* % descompuestos en lazos simples de dos conductores.
*
* function [I,mislazos]=matActivo(coordlazos,mislazoss)
*/
static mxArray * Mmatactivo(mxArray * * mislazos,
int nargout_,
mxArray * coordlazos,
mxArray * mislazoss) {
mexLocalFunctionTable save_local_function_table_
= mclSetCurrentLocalFunctionTable(&_local_function_table_matactivo);
mxArray * I = NULL;
mxArray * t = NULL;
mxArray * u = NULL;
mxArray * iiii = NULL;
mxArray * numero = NULL;
mclCopyArray(&coordlazos);
mclCopyArray(&mislazoss);
/*
*
* % Convertimos la matriz de lazos que nos pasan a otra de lazos simples.
*
* % Calculamos fila inicial. El valor sera 1 si no habia lazos activos.
* if mislazoss==0
*/
if (mclEqBool(mclVa(mislazoss, "mislazoss"), _mxarray0_)) {
/*
*
* numero=1
*/
mlfAssign(&numero, _mxarray1_);
mclPrintArray(mclVv(numero, "numero"), "numero");
/*
* else
*/
} else {
/*
* [numero,iiii]=size(mislazoss)
*/
mlfSize(
mlfVarargout(&numero, &iiii, NULL),
mclVa(mislazoss, "mislazoss"),
NULL);
mclPrintArray(mclVv(numero, "numero"), "numero");
mclPrintArray(mclVv(iiii, "iiii"), "iiii");
/*
* numero=numero+1;
*/
mlfAssign(&numero, mclPlus(mclVv(numero, "numero"), _mxarray1_));
/*
* end
*/
}
/*
*
* if numero>1
*/
if (mclGtBool(mclVv(numero, "numero"), _mxarray1_)) {
/*
* mislazos=mislazoss;
*/
mlfAssign(mislazos, mclVa(mislazoss, "mislazoss"));
/*
* end
*/
}
/*
*
* for u=1:size(coordlazos,1)
*/
{
int v_ = mclForIntStart(1);
int e_
= mclForIntEnd(
mlfSize(
mclValueVarargout(),
mclVa(coordlazos, "coordlazos"),
_mxarray1_));
if (v_ > e_) {
mlfAssign(&u, _mxarray2_);
} else {
/*
*
* if(coordlazos(u,1)==2) % Si es simple se copian los datos tal cual.
*
* mislazos(numero,1)=2;
* for(t=2:5)
* mislazos(numero,t)=coordlazos(u,t);
* end
*
* % Intensidad
* mislazos(numero,6)=coordlazos(u,10)*(cos(coordlazos(u,11)*pi/180)+i*sin(coordlazos(u,11)*pi/180));
* I(numero,1)=mislazos(numero,6);
*
* numero=numero+1;
*
* end
*
* if(coordlazos(u,1)==3) % Si es de 3 conductores salen 2 lazos simples.
*
* % PRIMER LAZO SIMPLE
* mislazos(numero,1)=2;
* for(t=2:5)
* mislazos(numero,t)=coordlazos(u,t);
* end
*
* % Intensidad
* mislazos(numero,8)=coordlazos(u,10)*(cos(coordlazos(u,11)*pi/180)+i*sin(coordlazos(u,11)*pi/180));
* I(numero,1)=mislazos(numero,8);
*
* numero=numero+1;
*
* % SEGUNDO LAZO SIMPLE
* mislazos(numero,1)=2;
* for(t=2:5)
* mislazos(numero,t)=coordlazos(u,t+2);
* end
*
* % Intensidad
* mislazos(numero,8)=coordlazos(u,12)*(cos(coordlazos(u,13)*pi/180)+i*sin(coordlazos(u,13)*pi/180));
* I(numero,1)=mislazos(numero,8);
*
* numero=numero+1;
*
* end
*
* if(coordlazos(u,1)==4) % Si es de 4 conductores salen 2 lazos simples.
*
* % PRIMER LAZO
* mislazos(numero,1)=2;
* for(t=2:5)
* mislazos(numero,t)=coordlazos(u,t);
* end
*
* % Intensidad
* mislazos(numero,8)=coordlazos(u,10)*(cos(coordlazos(u,11)*pi/180)+i*sin(coordlazos(u,11))*pi/180);
* I(numero,1)=mislazos(numero,8);
*
* numero=numero+1;
*
* % SEGUNDO LAZO
* mislazos(numero,1)=2;
* for(t=2:5)
* mislazos(numero,t)=coordlazos(u,t+4);
* end
*
* % Intensidad
* mislazos(numero,8)=coordlazos(u,12)*(cos(coordlazos(u,13)*pi/180)+i*sin(coordlazos(u,13)*pi/180));
* I(numero,1)=mislazos(numero,8);
*
* numero=numero+1;
*
* end
* end
*/
for (; ; ) {
if (mclEqBool(
mclIntArrayRef2(mclVa(coordlazos, "coordlazos"), v_, 1),
_mxarray3_)) {
mclArrayAssign2(
mislazos,
_mxarray3_,
mclVv(numero, "numero"),
_mxarray1_);
{
int v_0 = mclForIntStart(2);
int e_0 = 5;
if (v_0 > e_0) {
mlfAssign(&t, _mxarray2_);
} else {
for (; ; ) {
mclArrayAssign2(
mislazos,
mclIntArrayRef2(
mclVa(coordlazos, "coordlazos"), v_, v_0),
mclVv(numero, "numero"),
mlfScalar(v_0));
if (v_0 == e_0) {
break;
}
++v_0;
}
mlfAssign(&t, mlfScalar(v_0));
}
}
mclArrayAssign2(
mislazos,
mclMtimes(
mclIntArrayRef2(
mclVa(coordlazos, "coordlazos"), v_, 10),
mclPlus(
mlfCos(
mclMrdivide(
mclMtimes(
mclIntArrayRef2(
mclVa(coordlazos, "coordlazos"), v_, 11),
_mxarray4_),
_mxarray5_)),
mclMtimes(
_mxarray6_,
mlfSin(
mclMrdivide(
mclMtimes(
mclIntArrayRef2(
mclVa(coordlazos, "coordlazos"), v_, 11),
_mxarray4_),
_mxarray5_))))),
mclVv(numero, "numero"),
_mxarray7_);
mclArrayAssign2(
&I,
mclArrayRef2(
mclVv(*mislazos, "mislazos"),
mclVv(numero, "numero"),
_mxarray7_),
mclVv(numero, "numero"),
_mxarray1_);
mlfAssign(
&numero, mclPlus(mclVv(numero, "numero"), _mxarray1_));
}
if (mclEqBool(
mclIntArrayRef2(mclVa(coordlazos, "coordlazos"), v_, 1),
_mxarray8_)) {
mclArrayAssign2(
mislazos,
_mxarray3_,
mclVv(numero, "numero"),
_mxarray1_);
{
int v_1 = mclForIntStart(2);
int e_1 = 5;
if (v_1 > e_1) {
mlfAssign(&t, _mxarray2_);
} else {
for (; ; ) {
mclArrayAssign2(
mislazos,
mclIntArrayRef2(
mclVa(coordlazos, "coordlazos"), v_, v_1),
mclVv(numero, "numero"),
mlfScalar(v_1));
if (v_1 == e_1) {
break;
}
++v_1;
}
mlfAssign(&t, mlfScalar(v_1));
}
}
mclArrayAssign2(
mislazos,
mclMtimes(
mclIntArrayRef2(
mclVa(coordlazos, "coordlazos"), v_, 10),
mclPlus(
mlfCos(
mclMrdivide(
mclMtimes(
mclIntArrayRef2(
mclVa(coordlazos, "coordlazos"), v_, 11),
_mxarray4_),
_mxarray5_)),
mclMtimes(
_mxarray6_,
mlfSin(
mclMrdivide(
mclMtimes(
mclIntArrayRef2(
mclVa(coordlazos, "coordlazos"), v_, 11),
_mxarray4_),
_mxarray5_))))),
mclVv(numero, "numero"),
_mxarray9_);
mclArrayAssign2(
&I,
mclArrayRef2(
mclVv(*mislazos, "mislazos"),
mclVv(numero, "numero"),
_mxarray9_),
mclVv(numero, "numero"),
_mxarray1_);
mlfAssign(
&numero, mclPlus(mclVv(numero, "numero"), _mxarray1_));
mclArrayAssign2(
mislazos,
_mxarray3_,
mclVv(numero, "numero"),
_mxarray1_);
{
int v_2 = mclForIntStart(2);
int e_2 = 5;
if (v_2 > e_2) {
mlfAssign(&t, _mxarray2_);
} else {
for (; ; ) {
mclArrayAssign2(
mislazos,
mclIntArrayRef2(
mclVa(coordlazos, "coordlazos"),
v_,
v_2 + 2),
mclVv(numero, "numero"),
mlfScalar(v_2));
if (v_2 == e_2) {
break;
}
++v_2;
}
mlfAssign(&t, mlfScalar(v_2));
}
}
mclArrayAssign2(
mislazos,
mclMtimes(
mclIntArrayRef2(
mclVa(coordlazos, "coordlazos"), v_, 12),
mclPlus(
mlfCos(
mclMrdivide(
mclMtimes(
mclIntArrayRef2(
mclVa(coordlazos, "coordlazos"), v_, 13),
_mxarray4_),
_mxarray5_)),
mclMtimes(
_mxarray6_,
mlfSin(
mclMrdivide(
mclMtimes(
mclIntArrayRef2(
mclVa(coordlazos, "coordlazos"), v_, 13),
_mxarray4_),
_mxarray5_))))),
mclVv(numero, "numero"),
_mxarray9_);
mclArrayAssign2(
&I,
mclArrayRef2(
mclVv(*mislazos, "mislazos"),
mclVv(numero, "numero"),
_mxarray9_),
mclVv(numero, "numero"),
_mxarray1_);
mlfAssign(
&numero, mclPlus(mclVv(numero, "numero"), _mxarray1_));
}
if (mclEqBool(
mclIntArrayRef2(mclVa(coordlazos, "coordlazos"), v_, 1),
_mxarray10_)) {
mclArrayAssign2(
mislazos,
_mxarray3_,
mclVv(numero, "numero"),
_mxarray1_);
{
int v_3 = mclForIntStart(2);
int e_3 = 5;
if (v_3 > e_3) {
mlfAssign(&t, _mxarray2_);
} else {
for (; ; ) {
mclArrayAssign2(
mislazos,
mclIntArrayRef2(
mclVa(coordlazos, "coordlazos"), v_, v_3),
mclVv(numero, "numero"),
mlfScalar(v_3));
if (v_3 == e_3) {
break;
}
++v_3;
}
mlfAssign(&t, mlfScalar(v_3));
}
}
mclArrayAssign2(
mislazos,
mclMtimes(
mclIntArrayRef2(
mclVa(coordlazos, "coordlazos"), v_, 10),
mclPlus(
mlfCos(
mclMrdivide(
mclMtimes(
mclIntArrayRef2(
mclVa(coordlazos, "coordlazos"), v_, 11),
_mxarray4_),
_mxarray5_)),
mclMrdivide(
mclMtimes(
mclMtimes(
_mxarray6_,
mlfSin(
mclIntArrayRef2(
mclVa(coordlazos, "coordlazos"), v_, 11))),
_mxarray4_),
_mxarray5_))),
mclVv(numero, "numero"),
_mxarray9_);
mclArrayAssign2(
&I,
mclArrayRef2(
mclVv(*mislazos, "mislazos"),
mclVv(numero, "numero"),
_mxarray9_),
mclVv(numero, "numero"),
_mxarray1_);
mlfAssign(
&numero, mclPlus(mclVv(numero, "numero"), _mxarray1_));
mclArrayAssign2(
mislazos,
_mxarray3_,
mclVv(numero, "numero"),
_mxarray1_);
{
int v_4 = mclForIntStart(2);
int e_4 = 5;
if (v_4 > e_4) {
mlfAssign(&t, _mxarray2_);
} else {
for (; ; ) {
mclArrayAssign2(
mislazos,
mclIntArrayRef2(
mclVa(coordlazos, "coordlazos"),
v_,
v_4 + 4),
mclVv(numero, "numero"),
mlfScalar(v_4));
if (v_4 == e_4) {
break;
}
++v_4;
}
mlfAssign(&t, mlfScalar(v_4));
}
}
mclArrayAssign2(
mislazos,
mclMtimes(
mclIntArrayRef2(
mclVa(coordlazos, "coordlazos"), v_, 12),
mclPlus(
mlfCos(
mclMrdivide(
mclMtimes(
mclIntArrayRef2(
mclVa(coordlazos, "coordlazos"), v_, 13),
_mxarray4_),
_mxarray5_)),
mclMtimes(
_mxarray6_,
mlfSin(
mclMrdivide(
mclMtimes(
mclIntArrayRef2(
mclVa(coordlazos, "coordlazos"), v_, 13),
_mxarray4_),
_mxarray5_))))),
mclVv(numero, "numero"),
_mxarray9_);
mclArrayAssign2(
&I,
mclArrayRef2(
mclVv(*mislazos, "mislazos"),
mclVv(numero, "numero"),
_mxarray9_),
mclVv(numero, "numero"),
_mxarray1_);
mlfAssign(
&numero, mclPlus(mclVv(numero, "numero"), _mxarray1_));
}
if (v_ == e_) {
break;
}
++v_;
}
mlfAssign(&u, mlfScalar(v_));
}
}
mclValidateOutput(I, 1, nargout_, "I", "matactivo");
mclValidateOutput(*mislazos, 2, nargout_, "mislazos", "matactivo");
mxDestroyArray(numero);
mxDestroyArray(iiii);
mxDestroyArray(u);
mxDestroyArray(t);
mxDestroyArray(mislazoss);
mxDestroyArray(coordlazos);
mclSetCurrentLocalFunctionTable(save_local_function_table_);
return I;
}
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
