bruce/src/BruceFactory.cpp

/* $Id$ */
/*
 ******************************************************************************
 *                                                                            *
 *       COPYRIGHT  ACcESS 2005 -  All Rights Reserved                        *
 *                                                                            *
 * This software is the property of ACcESS. No part of this code              *
 * may be copied in any form or by any means without the expressed written    *
 * consent of ACcESS.  Copying, use or modification of this software          *
 * by any unauthorised person is illegal unless that person has a software    *
 * license agreement with ACcESS.                                             *
 *                                                                            *
 ******************************************************************************
*/

#include "bruce/Bruce/BruceFactory.h"
#include "bruce/Bruce/BruceException.h"

#include <iostream>
#include <sstream>
#include <boost/python/extract.hpp>

using namespace std;
using namespace escript;

namespace bruce {

AbstractContinuousDomain* brick(int n0, int n1, int n2,
                                double l0, double l1, double l2)
{
  double v0_len, v1_len, v2_len;

  Bruce::DimVec v0;
  Bruce::DimVec v1;
  Bruce::DimVec v2;
  Bruce::DimVec origin;

  origin.push_back(0);
  origin.push_back(0);
  origin.push_back(0);

  if (n0<2) {
    v0_len = 0;
  } else {
    v0_len = l0/(n0-1);
  }

  if (n1<2) {
    v1_len = 0;
  } else {
    v1_len = l1/(n1-1);
  }

  if (n2<2) {
    v2_len = 0;
  } else {
    v2_len = l2/(n2-1);
  }

  v0.push_back(v0_len);
  v0.push_back(0);
  v0.push_back(0);

  v1.push_back(0);
  v1.push_back(v1_len);
  v1.push_back(0);

  v2.push_back(0);
  v2.push_back(0);
  v2.push_back(v2_len);

  AbstractContinuousDomain* temp=new Bruce(v0, v1, v2, n0, n1, n2, origin);
  return temp;
}

AbstractContinuousDomain* rectangle(int n0, int n1,
                                    double l0, double l1)
{
  double v0_len, v1_len;

  Bruce::DimVec v0;
  Bruce::DimVec v1;
  Bruce::DimVec v2;
  Bruce::DimVec origin;

  origin.push_back(0);
  origin.push_back(0);

  if (n0<2) {
    v0_len = 0;
  } else {
    v0_len = l0/(n0-1);
  }

  if (n1<2) {
    v1_len = 0;
  } else {
    v1_len = l1/(n1-1);
  }

  v0.push_back(v0_len);
  v0.push_back(0);

  v1.push_back(0);
  v1.push_back(v1_len);

  AbstractContinuousDomain* temp=new Bruce(v0, v1, v2, n0, n1, 0, origin);
  return temp;
}

} // end of namespace
1	jgs	150	/* $Id$ */
2			/*
3			******************************************************************************
4			* *
5			* COPYRIGHT ACcESS 2005 - All Rights Reserved *
6			* *
7			* This software is the property of ACcESS. No part of this code *
8			* may be copied in any form or by any means without the expressed written *
9			* consent of ACcESS. Copying, use or modification of this software *
10			* by any unauthorised person is illegal unless that person has a software *
11			* license agreement with ACcESS. *
12			* *
13			******************************************************************************
14			*/
15
16			#include "bruce/Bruce/BruceFactory.h"
17			#include "bruce/Bruce/BruceException.h"
18
19			#include <iostream>
20			#include <sstream>
21			#include <boost/python/extract.hpp>
22
23			using namespace std;
24			using namespace escript;
25
26			namespace bruce {
27
28	jgs	151	AbstractContinuousDomain* brick(int n0, int n1, int n2,
29			double l0, double l1, double l2)
30	jgs	150	{
31	jgs	153	double v0_len, v1_len, v2_len;
32	jgs	150
33	jgs	151	Bruce::DimVec v0;
34			Bruce::DimVec v1;
35			Bruce::DimVec v2;
36			Bruce::DimVec origin;
37
38			origin.push_back(0);
39			origin.push_back(0);
40			origin.push_back(0);
41
42	jgs	153	if (n0<2) {
43			v0_len = 0;
44			} else {
45			v0_len = l0/(n0-1);
46			}
47
48			if (n1<2) {
49			v1_len = 0;
50			} else {
51			v1_len = l1/(n1-1);
52			}
53
54			if (n2<2) {
55			v2_len = 0;
56			} else {
57			v2_len = l2/(n2-1);
58			}
59
60			v0.push_back(v0_len);
61	jgs	151	v0.push_back(0);
62			v0.push_back(0);
63
64			v1.push_back(0);
65	jgs	153	v1.push_back(v1_len);
66	jgs	151	v1.push_back(0);
67
68			v2.push_back(0);
69			v2.push_back(0);
70	jgs	153	v2.push_back(v2_len);
71	jgs	151
72			AbstractContinuousDomain* temp=new Bruce(v0, v1, v2, n0, n1, n2, origin);
73			return temp;
74	jgs	150	}
75
76	jgs	151	AbstractContinuousDomain* rectangle(int n0, int n1,
77	jgs	150	double l0, double l1)
78			{
79	jgs	153	double v0_len, v1_len;
80	jgs	150
81	jgs	151	Bruce::DimVec v0;
82			Bruce::DimVec v1;
83			Bruce::DimVec v2;
84			Bruce::DimVec origin;
85
86			origin.push_back(0);
87			origin.push_back(0);
88
89	jgs	153	if (n0<2) {
90			v0_len = 0;
91			} else {
92			v0_len = l0/(n0-1);
93			}
94
95			if (n1<2) {
96			v1_len = 0;
97			} else {
98			v1_len = l1/(n1-1);
99			}
100
101			v0.push_back(v0_len);
102	jgs	151	v0.push_back(0);
103
104			v1.push_back(0);
105	jgs	153	v1.push_back(v1_len);
106	jgs	151
107			AbstractContinuousDomain* temp=new Bruce(v0, v1, v2, n0, n1, 0, origin);
108			return temp;
109	jgs	150	}
110
111	jgs	153	} // end of namespace
Name	Value
svn:eol-style	native
svn:keywords	Author Date Id Revision