MeshModel* EnvWrap::evalMesh(const QString& nm)
{
int ii = evalInt(nm);
QScriptValue mdsv = env->globalObject().property(ScriptAdapterGenerator::meshDocVarName());
MeshDocumentSI* mdsi = dynamic_cast<MeshDocumentSI*>(mdsv.toQObject());
if (mdsi != NULL)
return mdsi->md->getMesh(ii);
return NULL;
}
int
GusdOBJ_usdcamera::applyInputIndependentTransform(OP_Context& ctx, UT_DMatrix4& mx)
{
mx.identity();
fpreal t = ctx.getTime();
if(UsdGeomCamera cam = _LoadCamera(t, ctx.getThread()))
{
float frame = evalFloat(_frameIdx, 0, t);
GfMatrix4d ctm(1.);
bool stat = true;
bool resetsXformStack = false;
switch(evalInt("xformmode", 0, t))
{
case _POSTMULTCTM_TRANSFORM:
stat = true;
ctm = cam.ComputeLocalToWorldTransform(frame);
break;
case _CTM_TRANSFORM:
stat = true;
ctm = cam.ComputeParentToWorldTransform(frame);
break;
case _OBJ_TRANSFORM:
// XXX: how do we reset xformStack here?
// Is that (or should that
// be) handled by the Compute calls above?
stat = cam.GetLocalTransformation(&ctm, &resetsXformStack, frame);
break;
default: // _IGNORE_TRANSFORM:
stat = true;
ctm.SetIdentity();
break;
}
if(!stat)
{
stealErrors(_errors, /*borrow*/ true);
return 0;
}
mx = GusdUT_Gf::Cast(ctm);
}
return OBJ_Camera::applyInputIndependentTransform(ctx, mx);
}
int main(void){
system("clear");
printf("Programa que calcula el área bajo la curva (1/(σ sqrt(2π)))e^((-1/2)((x-µ)/σ)))² por el método rectangular.\n\n");
int sh=0;
char opc='\0';
double sigma=-1; // σ
double mu=0.0; //µ
double c1;// constante que acompaña a la integral
double auxa=0,auxb=0,a=0,b=0; //límites de integración
double areaBajoCurva=0.0;
long int iteraciones=300;//Númerod e iteraciones por defecto
double x,iteracionesaux=0;
do{
fflush( stdin );fflush( stdin );fflush( stdin );
puts ("Constantes: \nIngrese el valor de σ:");
while(sigma<=0){
getConst(&sigma,&sh);
if(sigma==0)printf("No puedo hacer eso, ya que haré el cálculo:\n (1/(σ sqrt(2π))) y si σ=0 eso se indetermina. Intenta con otro número: \n ");
if(sigma<0.0)printf("Sigma debe ser estrictamente mayor que 0 \n ");
}
puts("Ingrese el valor de µ:");
getConst(&mu,&sh);
while(auxa==auxb){
puts("Límites de integración:\nCalcular área desde (ingrese el valor de a):");
getConst(&auxa,&sh);
printf("Hasta (ingrese el valor de b):");
getConst(&auxb,&sh);
if(auxb>auxa)
{
a=auxa;
b=auxb;
}else if(auxa>auxb){
printf("Parece que tus límites están en el orden incorrecto, a debe ser menor a b por definición, pero puedo manejarlo, invirtiendo..\n");
b=auxa;
a=auxb;
}
else{
printf("Eh!, ahí no hay área!, escribe límites de integración que no sean iguales...Va de nuevo:\n");
}
}
x=a;
printf("Calculando área con σ=%2.2f,µ=%2.2f,a=%2.2f,b=%2.2f\n ",sigma,mu,a,b);
double base=(b-a)/iteraciones;
c1=1/(sigma*(sqrt(2*pi)));
double ev;
for(int i=1;i<=iteraciones;i++){
ev=evalInt(x,mu,sigma);
areaBajoCurva=areaBajoCurva+((c1)*(ev)*(base));
//printf("[%f] x=%f",ev,x);
x=x+base;
}
printf("Area bajo (1/(σ sqrt(2π)))e^((-1/2)((x-µ)/σ)))²:[%f]\n\n\n ",areaBajoCurva);
printf("Detalles:\n Calculé el área desde: %2.2f hasta %2.2f\n ",a,b,iteraciones);
printf("Otro Cálculo? S para si.\n");
//getchar();
scanf("%c",&opc);
a=0;
b=0;
auxa=0;
auxb=0;
sigma=-1;
getchar();
system("clear");
}while(opc=='S'||opc=='s');
}
int EnvWrap::evalEnum( const QString& nm )
{
return evalInt(nm);
}
node_t(char c, node_ptr u = node_ptr(), node_ptr v = node_ptr())
: c_(c), left_(u), right_(v)
{
evalInt();
}
bool PoroElasticity::evalIntMx (LocalIntegral& elmInt,
const MxFiniteElement& fe,
const TimeDomain& time, const Vec3& X) const
{
return evalInt(elmInt, fe, time, X);
}
