mirror of
https://github.com/gwm17/catima.git
synced 2024-11-22 18:28:51 -05:00
60 lines
1.5 KiB
C++
60 lines
1.5 KiB
C++
#include "lest.hpp"
|
|
#include <math.h>
|
|
#include "catima/catima.h"
|
|
#include "catima/reactions.h"
|
|
#include "testutils.h"
|
|
using namespace std;
|
|
using catima::approx;
|
|
using catima::reaction_rate;
|
|
using catima::nonreaction_rate;
|
|
const lest::test specification[] =
|
|
{
|
|
CASE("reaction"){
|
|
catima::Projectile proj{12,6,6,870};
|
|
catima::Projectile proj2{238,92,92,500};
|
|
auto c = catima::get_material(6);
|
|
auto h = catima::get_material(1);
|
|
catima::Material water({{0,8,2},{0,1,1}},1.0);
|
|
c.thickness(2.0);
|
|
double r,r2;
|
|
|
|
r= catima::nonreaction_rate(proj, c);
|
|
EXPECT(r == approx(0.92,0.02));
|
|
|
|
catima::Layers l;
|
|
l.add(c);
|
|
l.add(c);
|
|
l.add(c);
|
|
|
|
auto res = catima::calculate(proj,l);
|
|
EXPECT(res.total_result.sp == approx(r*r*r,0.02));
|
|
|
|
c.thickness(6.0);
|
|
r= catima::nonreaction_rate(proj, c);
|
|
EXPECT(res.total_result.sp == approx(r,0.001));
|
|
|
|
c.thickness(0.0);
|
|
r= catima::nonreaction_rate(proj, c);
|
|
EXPECT(r == 1.0);
|
|
proj.T = 0.0;
|
|
c.thickness(6);
|
|
r= catima::nonreaction_rate(proj, c);
|
|
EXPECT(r == -1.0);
|
|
|
|
proj.T=870;
|
|
water.thickness(1);
|
|
r= catima::nonreaction_rate(proj2, water);
|
|
r2= catima::nonreaction_rate(proj, water);
|
|
EXPECT( (r > 0 && r<1.0) );
|
|
EXPECT( (r2 > 0 && r2<1.0) );
|
|
EXPECT( r2>r );
|
|
|
|
}
|
|
};
|
|
|
|
int main( int argc, char * argv[] )
|
|
{
|
|
return lest::run( specification, argc, argv );
|
|
}
|
|
|