c81data.h File Reference

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Functions
void	c81_data_destroy (c81_data *data)
c81_data *	c81_data_get (long int jpro)
int	c81_data_set (long int jpro, c81_data *data)
int	c81_data_merge (unsigned ndata, const c81_data **data, const doublereal *upper_bounds, doublereal dCsi, doublereal dcltol, c81_data *i_data)
int	c81_data_do_stall (c81_data *data, const doublereal dcltol)
doublereal	c81_data_get_coef (int nm, doublereal *m, int na, doublereal *a, doublereal alpha, doublereal mach)
int	c81_data_flip (c81_data *data)

Function Documentation

void c81_data_destroy

(

c81_data *

data

)

Definition at line 276 of file c81data.cc.

References c81_data::ad, c81_data::al, c81_data::am, c81_data::md, c81_data::ml, c81_data::mm, c81_data::mstall, and c81_data::stall.

Referenced by main(), C81InterpolatedAeroData::~C81InterpolatedAeroData(), and MBDynParser::~MBDynParser().

{
        delete[] data->ml;
        delete[] data->al;

        delete[] data->md;
        delete[] data->ad;

        delete[] data->mm;
        delete[] data->am;

        delete[] data->stall;
        delete[] data->mstall;
}

int c81_data_do_stall	(	c81_data *	data,
		const doublereal	dcltol
	)

Definition at line 1244 of file c81data.cc.

References c81_data::al, c81_data::am, do_stall(), c81_data::mstall, c81_data::NAL, c81_data::NAM, c81_data::NML, c81_data::NMM, and c81_data::stall.

Referenced by c81_data_fc511_read(), c81_data_merge(), c81_data_nrel_read(), c81_data_read(), and c81_data_read_free_format().

{
        if (data == NULL || data->NML <= 0 || data->NMM <= 0) {
                return -1;
        }

        data->stall = new doublereal[3*data->NML];
        do_stall(data->NML, data->NAL, data->al, data->stall, dcltol);

        data->mstall = new doublereal[3*data->NMM];
        do_stall(data->NMM, data->NAM, data->am, data->mstall, dcltol);

        return 0;
}

Here is the call graph for this function:

int c81_data_flip

(

c81_data *

data

)

Definition at line 1286 of file c81data.cc.

References c81_data::ad, c81_data::al, c81_data::am, flip_one(), c81_data::NAD, c81_data::NAL, c81_data::NAM, c81_data::NMD, c81_data::NML, and c81_data::NMM.

Referenced by MBDynParser::C81Data_int(), and main().

{
        int rc;

        rc = flip_one(data->NML, data->NAL, data->al, -1);
        if (rc) {
                return rc;
        }

        rc = flip_one(data->NMD, data->NAD, data->ad, 1);
        if (rc) {
                return rc;
        }

        rc = flip_one(data->NMM, data->NAM, data->am, -1);
        return rc;
}

Here is the call graph for this function:

c81_data* c81_data_get

(

long int

jpro

)

doublereal c81_data_get_coef	(	int	nm,
		doublereal *	m,
		int	na,
		doublereal *	a,
		doublereal	alpha,
		doublereal	mach
	)

Definition at line 90 of file aerodc81.c.

References c, and get_coef().

Referenced by main().

{
        doublereal c;
        
        get_coef(nm, m, na, a, alpha, mach, &c, NULL);

        return c;
}

Here is the call graph for this function:

int c81_data_merge	(	unsigned	ndata,
		const c81_data **	data,
		const doublereal *	upper_bounds,
		doublereal	dCsi,
		doublereal	dcltol,
		c81_data *	i_data
	)

Definition at line 396 of file c81data.cc.

References c81_data::ad, c81_data::al, c81_data::am, ASSERT, c81_data_do_stall(), for(), c81_data::header, c81_data::md, c81_data::ml, c81_data::mm, c81_data::NAD, c81_data::NAL, c81_data::NAM, c81_data::NMD, c81_data::NML, and c81_data::NMM.

Referenced by C81InterpolatedAeroData::C81InterpolatedAeroData().

{
        ASSERT(ndata > 0);
        ASSERT(data != 0);
        ASSERT(upper_bounds != 0);
        ASSERT(dCsi >= -1.);
        ASSERT(dCsi <= 1.);
        ASSERT(dcltol > 0.);
        ASSERT(i_data != 0);

        if (dCsi < upper_bounds[0]) {
                silent_cerr("cannot find C81 data lower bound for point xi=" << dCsi << std::endl);
                return -1;
        }

        int to = 0;
        for (unsigned i = 1; i < ndata; i++) {
                if (upper_bounds[i] > dCsi) {
                        to = i;
                        break;
                }
        }
        if (to == 0) {
                silent_cerr("cannot find C81 data upper bound for point xi=" << dCsi << std::endl);
                return -1;
        }

        int from = to - 1;
        if (unsigned(from) == ndata) {
                silent_cerr("cannot find C81 data upper bound for point xi=" << dCsi << std::endl);
                return -1;
        }

        /* we need to interpolate between data[from]
         * and data[to] */

        /* we only accept homogeneous data sources,
         * i.e. same Mach and alpha patterns */
        if (data[from]->NML != data[to]->NML) {
                silent_cerr("number of Mach points for Cl between airfoils "
                        << from << " (" << data[from]->NML << ") and "
                        << to << " (" << data[to]->NML << ") do not match"
                        << std::endl);
                return -1;
        }

        if (data[from]->NAL != data[to]->NAL) {
                silent_cerr("number of AoA points for Cl between airfoils "
                        << from << " (" << data[from]->NAL << ") and "
                        << to << " (" << data[to]->NAL << ") do not match"
                        << std::endl);
                return -1;
        }

        if (data[from]->NMD != data[to]->NMD) {
                silent_cerr("number of Mach points for Cd between airfoils "
                        << from << " (" << data[from]->NMD << ") and "
                        << to << " (" << data[to]->NMD << ") do not match"
                        << std::endl);
                return -1;
        }

        if (data[from]->NAD != data[to]->NAD) {
                silent_cerr("number of AoA points for Cd between airfoils "
                        << from << " (" << data[from]->NAD << ") and "
                        << to << " (" << data[to]->NAD << ") do not match"
                        << std::endl);
                return -1;
        }

        if (data[from]->NMM != data[to]->NMM) {
                silent_cerr("number of Mach points for Cm between airfoils "
                        << from << " (" << data[from]->NMM << ") and "
                        << to << " (" << data[to]->NMM << ") do not match"
                        << std::endl);
                return -1;
        }

        if (data[from]->NAM != data[to]->NAM) {
                silent_cerr("number of AoA points for Cm between airfoils "
                        << from << " (" << data[from]->NAM << ") and "
                        << to << " (" << data[to]->NAM << ") do not match"
                        << std::endl);
                return -1;
        }

        for (int i = 0; i < data[from]->NML; i++) {
                if (data[from]->ml[i] != data[to]->ml[i]) {
                        silent_cerr("Mach point " << i << "for Cl of airfoils "
                                << from << " (" << data[from]->ml[i] << ") and "
                                << to << " (" << data[to]->ml[i] << ") differs"
                                << std::endl);
                        return -1;
                }
        }

        for (int i = 0; i < data[from]->NAL; i++) {
                if (data[from]->al[i] != data[to]->al[i]) {
                        silent_cerr("AoA point " << i << "for Cl of airfoils "
                                << from << " (" << data[from]->al[i] << ") and "
                                << to << " (" << data[to]->al[i] << ") differs"
                                << std::endl);
                        return -1;
                }
        }

        for (int i = 0; i < data[from]->NMD; i++) {
                if (data[from]->md[i] != data[to]->md[i]) {
                        silent_cerr("Mach point " << i << "for Cd of airfoils "
                                << from << " (" << data[from]->md[i] << ") and "
                                << to << " (" << data[to]->md[i] << ") differs"
                                << std::endl);
                        return -1;
                }
        }

        for (int i = 0; i < data[from]->NAD; i++) {
                if (data[from]->ad[i] != data[to]->ad[i]) {
                        silent_cerr("AoA point " << i << "for Cd of airfoils "
                                << from << " (" << data[from]->ad[i] << ") and "
                                << to << " (" << data[to]->ad[i] << ") differs"
                                << std::endl);
                        return -1;
                }
        }

        for (int i = 0; i < data[from]->NMM; i++) {
                if (data[from]->mm[i] != data[to]->mm[i]) {
                        silent_cerr("Mach point " << i << "for Cm of airfoils "
                                << from << " (" << data[from]->mm[i] << ") and "
                                << to << " (" << data[to]->mm[i] << ") differs"
                                << std::endl);
                        return -1;
                }
        }

        for (int i = 0; i < data[from]->NAM; i++) {
                if (data[from]->am[i] != data[to]->am[i]) {
                        silent_cerr("AoA point " << i << "for Cm of airfoils "
                                << from << " (" << data[from]->am[i] << ") and "
                                << to << " (" << data[to]->am[i] << ") differs"
                                << std::endl);
                        return -1;
                }
        }

        snprintf(i_data->header, sizeof(i_data->header),
                "interpolated (\"%s\"[%e]->\"%s\"[%e]: [%e])",
                data[from]->header, upper_bounds[from],
                data[to]->header, upper_bounds[to], dCsi);

        i_data->NML = data[from]->NML;
        i_data->NAL = data[from]->NAL;
        i_data->NMD = data[from]->NMD;
        i_data->NAD = data[from]->NAD;
        i_data->NMM = data[from]->NMM;
        i_data->NAM = data[from]->NAM;

        doublereal dw = upper_bounds[to] - upper_bounds[from];
        doublereal dw_from = (upper_bounds[to] - dCsi)/dw;
        doublereal dw_to = (dCsi - upper_bounds[from])/dw;

        /* lift */
        i_data->ml = new doublereal[i_data->NML];
        for (int i = 0; i < i_data->NML; i++) {
                i_data->ml[i] = data[from]->ml[i];
        }
        i_data->al = new doublereal[(i_data->NML + 1)*i_data->NAL];
        for (int i = 0; i < i_data->NAL; i++) {
                i_data->al[i] = data[from]->al[i];
        }
        for (int i = i_data->NAL; i < (i_data->NML + 1)*i_data->NAL; i++) {
                i_data->al[i] = dw_from*data[from]->al[i] + dw_to*data[to]->al[i];
        }

        /* drag */
        i_data->md = new doublereal[i_data->NMD];
        for (int i = 0; i < i_data->NMD; i++) {
                i_data->md[i] = data[from]->md[i];
        }
        i_data->ad = new doublereal[(i_data->NMD + 1)*i_data->NAD];
        for (int i = 0; i < i_data->NAD; i++) {
                i_data->ad[i] = data[from]->ad[i];
        }
        for (int i = i_data->NAD; i < (i_data->NMD + 1)*i_data->NAD; i++) {
                i_data->ad[i] = dw_from*data[from]->ad[i] + dw_to*data[to]->ad[i];
        }

        /* moment */
        i_data->mm = new doublereal[i_data->NMM];
        for (int i = 0; i < i_data->NMM; i++) {
                i_data->mm[i] = data[from]->mm[i];
        }
        i_data->am = new doublereal[(i_data->NMM + 1)*i_data->NAM];
        for (int i = 0; i < i_data->NAM; i++) {
                i_data->am[i] = data[from]->am[i];
        }
        for (int i = i_data->NAM; i < (i_data->NMM + 1)*i_data->NAM; i++) {
                i_data->am[i] = dw_from*data[from]->am[i] + dw_to*data[to]->am[i];
        }

        // FIXME: maybe this is not the best place
        c81_data_do_stall(i_data, dcltol);

        return 0;
}

Here is the call graph for this function:

int c81_data_set	(	long int	jpro,
		c81_data *	data
	)

Definition at line 1137 of file c81data.cc.

References __c81_pdata, and c81_ndata.

{
        if (__c81_pdata == NULL || jpro >= c81_ndata) {
                c81_data** pp = NULL;
                int ndata = c81_ndata;

                if (jpro >= c81_ndata) {
                        ndata = jpro+1;
                }

                pp = new c81_data*[ndata];
                if (__c81_pdata != NULL) {
                        for (int i = 0; i < ndata; i++) {
                                pp[i] = __c81_pdata[i];
                        }
                        delete[] __c81_pdata;
                }
                __c81_pdata = pp;
                c81_ndata = ndata;
        } else if (__c81_pdata[jpro] != NULL) {
                silent_cerr("profile " << jpro << "already defined" << std::endl);
                return -1;
        }

        __c81_pdata[jpro] = data;

        return 0;
}