Ir ao conteúdo
  • Comunicados

    • Gabriel Torres

      Seja um moderador do Clube do Hardware!   12-02-2016

      Prezados membros do Clube do Hardware, Está aberto o processo de seleção de novos moderadores para diversos setores ou áreas do Clube do Hardware. Os requisitos são:   Pelo menos 500 posts e um ano de cadastro; Boa frequência de participação; Ser respeitoso, cordial e educado com os demais membros; Ter bom nível de português; Ter razoável conhecimento da área em que pretende atuar; Saber trabalhar em equipe (com os moderadores, coordenadores e administradores).   Os interessados deverão enviar uma mensagem privada para o usuário @Equipe Clube do Hardware com o título "Candidato a moderador". A mensagem deverá conter respostas às perguntas abaixo:   Qual o seu nome completo? Qual sua data de nascimento? Qual sua formação/profissão? Já atuou como moderador em algo outro fórum, se sim, qual? De forma sucinta, explique o porquê de querer ser moderador do fórum e conte-nos um pouco sobre você.   OBS: Não se trata de função remunerada. Todos que fazem parte do staff são voluntários.
    • DiF

      Poste seus códigos corretamente!   21-05-2016

      Prezados membros do Fórum do Clube do Hardware, O Fórum oferece um recurso chamado CODE, onde o ícone no painel do editor é  <>     O uso deste recurso é  imprescindível para uma melhor leitura, manter a organização, diferenciar de texto comum e principalmente evitar que os compiladores e IDEs acusem erro ao colar um código copiado daqui. Portanto convido-lhes para ler as instruções de como usar este recurso CODE neste tópico:  
Entre para seguir isso  
osvaldo182

[Resolvido] Ajuda - Melhorar eficiência

Recommended Posts

Opa,

Estou com um problema nesse programa. Na verdade ele funciona bem, mas preciso aumentar sua eficiencia. Minha ideia seria passar as variaveis em forma de ponteiro, porém quanto implementei o código o compilador achou muitos erros.

#include <stdio.h>
#include <math.h>


// Recursive auxiliary function for adaptiveSimpsons() function below
//
double adaptiveSimpsonsAux(double (*f)(double), double a, double b, double epsilon,
double S, double fa, double fb, double fc, int bottom) {
double c = (a + b)/2, h = b - a;
double d = (a + c)/2, e = (c + b)/2;
double fd = f(d), fe = f(e);
double Sleft = (h/12)*(fa + 4*fd + fc);
double Sright = (h/12)*(fc + 4*fe + fb);
double S2 = Sleft + Sright;
if (bottom <= 0 || fabs(S2 - S) <= 15*epsilon)
return S2 + (S2 - S)/15;
else return adaptiveSimpsonsAux(f, a, c, epsilon/2, Sleft, fa, fc, fd, bottom-1) + adaptiveSimpsonsAux(f, c, b, epsilon/2, Sright, fc, fb, fe, bottom-1);
}

//
// Adaptive Simpson's Rule
//
double adaptiveSimpsons(double (*f)(double), // ptr to function
double a, double b, // interval [a,b]
double epsilon, // error tolerance
int maxRecursionDepth) { // recursion cap
double c = (a + b)/2, h = b - a;
double fa = f(a), fb = f(, fc = f(c);
double S = (h/6)*(fa + 4*fc + fb);
return adaptiveSimpsonsAux(f, a, b, epsilon, S, fa, fb, fc, maxRecursionDepth);
}

double f (double argument){
return sin(pow(argument,2));
}


int main(){
double a = 0 , b = 113, epsilon = 0.000000001;
int recDepth = 100;
double I = adaptiveSimpsons(&f, a, b, epsilon, recDepth);

printf("\n\nO valor da integral e: %.15lf",I);
printf("\nO erro absoluto e: %.15lf", fabs(I-0.6266421143805967));
}

Acho que o problema é na função auxiliar que é recursiva.

Alguém tem alguma sugestão para isso?

Desde já, Obrigado a todos.

Compartilhar este post


Link para o post
Compartilhar em outros sites
  • Autor do tópico
  • Tentei de uma forma porca arrumar. ( simplesmente coloquei ponteiros nos argumentos da função f, ainda não pensei nas consequencias :wacko:)

    o codigo ficou assim:


    #include <math.h>
    #include <stdlib.h>

    // Recursive auxiliary function for adaptiveSimpsons() function below
    //
    double adaptiveSimpsonsAux(double (*f)(double *), double *a, double *b, double epsilon,
    double S, double fa, double fb, double fc, int bottom) {

    double *d, *e, h, *c;
    c = (double *) malloc (sizeof(double));
    d = (double *) malloc (sizeof(double));
    e = (double *) malloc (sizeof(double));

    *c = (*a + *b)/2, h = *b - *a;
    *d = (*a + *c)/2, *e = (*c + *b)/2;
    double fd = f(d), fe = f(e);
    double Sleft = (h/12)*(fa + 4*fd + fc);
    double Sright = (h/12)*(fc + 4*fe + fb);
    double S2 = Sleft + Sright;
    if (bottom <= 0 || fabs(S2 - S) <= 15*epsilon)
    return S2 + (S2 - S)/15;
    else return adaptiveSimpsonsAux(f, a, c, epsilon/2, Sleft, fa, fc, fd, bottom-1) + adaptiveSimpsonsAux(f, c, b, epsilon/2, Sright, fc, fb, fe, bottom-1);
    }

    //
    // Adaptive Simpson's Rule
    //
    double adaptiveSimpsons(double (*f)(double *), // ptr to function
    double *a, double *b, // interval [a,b]
    double epsilon, // error tolerance
    int maxRecursionDepth) { // recursion cap

    double *c;
    c = (double *) malloc (sizeof(double));
    *c = (*a + *b)/2;
    double h = *b - *a;
    double fa = f(a), fb = f(B), fc = f(c);
    double S = (h/6)*(fa + 4*fc + fb);
    return adaptiveSimpsonsAux(f, a, b, epsilon, S, fa, fb, fc, maxRecursionDepth);
    }

    double f (double *argument){
    return sin(pow(*argument,2));
    }


    int main(){
    double *a , *b, epsilon = 0.000000001;
    a = (double *) malloc (sizeof(double));
    b = (double *) malloc (sizeof(double));
    *a=0;
    *b=113;
    int recDepth = 100;
    double I = adaptiveSimpsons(&f, a, b, epsilon, recDepth); /* <<<< AQUI É O ERRO*/

    printf("\n\nO valor da integral e: %.15lf",I);
    printf("\nO erro absoluto e: %.15lf", fabs(I-0.6266421143805967));
    }
    #include <stdio.h>

    Agora o codigo esta correto porém a eficiencia nao mudou muito.

    Editado por osvaldo182
    Resolvido

    Compartilhar este post


    Link para o post
    Compartilhar em outros sites

    Caso o autor do tópico necessite, o mesmo será reaberto, para isso deverá entrar em contato com a moderação solicitando o desbloqueio.

    Compartilhar este post


    Link para o post
    Compartilhar em outros sites
    Visitante
    Este tópico está impedido de receber novos posts.
    Entre para seguir isso  





    Sobre o Clube do Hardware

    No ar desde 1996, o Clube do Hardware é uma das maiores, mais antigas e mais respeitadas publicações sobre tecnologia do Brasil. Leia mais

    Direitos autorais

    Não permitimos a cópia ou reprodução do conteúdo do nosso site, fórum, newsletters e redes sociais, mesmo citando-se a fonte. Leia mais

    ×