mapovani-pole-funkci.c

Go to the documentation of this file.

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

double na2(double x){ return x*x; }
double na3(double x){ return x*x*x; }

double **map(double(*fce[])(double),double *vstup,int pocet_fce, int pocet_vstup) {
  int radku;//kolik radku ma vystup: 1+pocet_fce
  double **vysledky;
  double  *vysledek;
  double (*f)(double);
  int i, j;

  assert(pocet_fce>=0);
  assert(pocet_vstup>0);
  assert(vstup);

  radku = pocet_fce+1;//+1 : prvni (nulty) radek opakuje/obsahuje vstup
  vysledky = (double **)malloc(radku * sizeof(double *));
//naalokujeme radky pro vysledky
  for (i=0; i<radku; i++) {
    vysledky[i] = (double *)malloc(pocet_vstup * sizeof(double));
  }
//vysledky[0] obsahuje kopii vstup:
  for(i=0; i<pocet_vstup; i++) {
     vysledky[0][i] = vstup[i];
  }
  for(i=0; i<pocet_fce; i++) {
     vysledek = vysledky[i+1];//+1, vysledky[0] obsahuje kopii vstupu!
     f = fce[i];
     for(j=0; j<pocet_vstup; j++) {
       vysledek[j] = f(vstup[j]);
     }
  }
  return vysledky;
}

static  vypis (double *vstup, int delka) {
  int i;
  for (i=0; i<delka; i++) {
      printf("%c%8g", i?',':'[', vstup[i]);
  };printf("]\n");
}

#ifndef NELEM
#define NELEM(pole) (sizeof pole/sizeof *pole)
#endif

double vstup[] = {1,2,3,4,5, 10};

//ted vyuzijeme math.h
double ludolf[] = {0.0, M_PI_4, M_PI_2, M_PI};

double (*fce[])(double) = {na2, na3, sqrt};

int main(void) {
   double **vysledky_na2;
   double **vysledky_na3;
   double **vysledky_sin;
   int i;

   vysledky_na2 = map(fce, vstup, NELEM(fce), NELEM(vstup));
   vysledky_na3 = map(fce, vstup, NELEM(fce), NELEM(vstup));

//vyuzijeme vypis() z predchoziho prikladu mapovani-funkce.c
   for(i=0; i<= NELEM(fce); i++) {
   vypis(vysledky_na2[i], NELEM(vstup));
   }
}