C / C++ / MFC

<br />
// Q2_new.cpp : Defines the entry point for the console application.<br />
//<br />
<br />
#include <math><br />
#include "stdafx.h"<br />
#include <iostream><br />
using namespace std;<br />
<br />
double theta, phi, t, p;	//azimuth and dip angles respectively<br />
double lambda1 = 0.5;<br />
double lambda2 = 2;<br />
double a11,a12,a21,a22,d11,d13,d31,d33;<br />
const double pi = 3.142857;<br />
const int rows = 3;<br />
const int cols = 3;<br />
<br />
void mult(double *fr_array2);<br />
<br />
<br />
//---------------------------------------------------------<br />
//Initializing the arrays to zero<br />
//---------------------------------------------------------<br />
	//array that holds the dip angle, phi<br />
	double dip_array [3][3] = {0};	<br />
<br />
	//array that holds the azimuth angle, theta<br />
	double azim_array [3][3] = {0};	<br />
<br />
	//array that holds the anisotropy ratios<br />
	double lambda_array [3][3] = {0};	<br />
<br />
	//array that holds the original range values<br />
	double r_array [3] = {0};<br />
<br />
	//array that holds final range values<br />
	double fr_array [3][3] = {0};	<br />
<br />
	//array that holds final range values<br />
	double fr_array2 [3] = {0};	<br />
<br />
<br />
<br />
//---------------------------------------------------------<br />
//Code to convert the angles to radians, acceptable in C++<br />
//---------------------------------------------------------<br />
//Dip Angle,p<br />
double a(double p)<br />
{<br />
	d11 = cos((p*pi)/180);<br />
	return d11;<br />
}<br />
<br />
double b(double p)<br />
{<br />
	d13 = -1*sin((p*pi)/180);<br />
	return d13;<br />
}<br />
<br />
double c(double p)<br />
{<br />
	d31 = sin((p*pi)/180);<br />
	return d31;<br />
}<br />
<br />
double d(double p)<br />
{<br />
	d33 = cos((p*pi)/180);<br />
	return d33;<br />
}<br />
<br />
<br />
//Azimuth Angles, t<br />
double e(double t)<br />
{<br />
	a11 = sin((t*pi)/180);<br />
	return a11;<br />
}<br />
<br />
double f(double t)<br />
{<br />
	a12 = cos((t*pi)/180);<br />
	return a12;<br />
}<br />
<br />
double g(double t)<br />
{<br />
	a21 = -1*cos((t*pi)/180);<br />
	return a21;<br />
}<br />
<br />
double h(double t)<br />
{<br />
	a22 = sin((t*pi)/180);<br />
	return a22;<br />
}<br />
<br />
<br />
void printarray (double arg[][3]) <br />
{<br />
  for (int i=0; i<rows; i++)//for each row<br />
  {<br />
	  for (int j=0; j<cols; j++)	//for each column value<br />
		  cout << arg[i][j] << " ";<br />
	  cout << "\n";<br />
  }<br />
}<br />
<br />
void printarray2(double arg[]) <br />
{<br />
  for (int i=0; i<rows; i++)//for each row<br />
  {<br />
	  //for (int j=0; j<3; j++)	//for each column value<br />
		  cout << arg[i] << " ";<br />
	  cout << "\n";<br />
  }<br />
}<br />
<br />
	//code for multiplication of the dip and azim matrices<br />
/*	void mult(double fr_array2 [3])<br />
	{<br />
		for(int i=0;i<rows;i++)	//for each row<br />
		{<br />
			for(int j=0;j<cols;j++)	//for each column<br />
			{<br />
				for(int k=0;k<cols;k++)	//for each column<br />
				{<br />
					for(int l=0;l<cols;l++)	//for each column<br />
					{<br />
					//	fr_array [i][j] += dip_array [i][k]*azim_array [k][j];<br />
						fr_array [i][j] += lambda_array [i][l]*dip_array [l][k]*azim_array[k][j];<br />
						fr_array2 [i] += fr_array[i][j]*r_array[i];<br />
					}<br />
				}<br />
			}<br />
		}<br />
	}*/<br />
<br />
<br />
<br />
int main(int argc, char* argv[])<br />
{<br />
	<br />
	cout<<"ANISOTROPY MODELLING.\n";<br />
	cout<<"\nPlease enter a real number between 0 and 360 for the azimuth angle, theta, then press the 'Enter' key.\n\n";<br />
	cin>>t;<br />
<br />
	cout<<"\nPlease enter a real number between 0 and 360 for the dip angle, phi, then press the 'Enter' key.\n\n";<br />
	cin>>p;<br />
<br />
	a(p);<br />
	b(p);<br />
	c(p);<br />
	d(p);<br />
	e(t);<br />
	f(t);<br />
	g(t);<br />
	h(t);<br />
<br />
	//-----------------------------------------------------------<br />
	//Initializing the matrices<br />
	//-----------------------------------------------------------<br />
	//array that holds the dip angle, phi<br />
	double dip_array [3][3] = {{d11,0,d13},{0,1,0},{d31,0,d33}};	<br />
<br />
	//array that holds the azimuth angle, theta<br />
	double azim_array [3][3] = {{a11,a12,0},{a21,a22,0},{0,0,1}};	<br />
<br />
	//array that holds the anisotropy ratios<br />
	double lambda_array [3][3] = {{lambda1,0,0},{0,1,0},{0,0,lambda2}};	<br />
<br />
	//array that holds the original range values<br />
	double r_array [3] = {4,8,3.94};<br />
<br />
<br />
	//----------------------------------------------------------------<br />
<br />
	//code to check that the arrays are properly initialized in the matrices for any theta,phi<br />
	cout<<"\nThe dip array is:\n";<br />
	printarray(dip_array);<br />
<br />
	cout<<"\nThe azimuth array is:\n";<br />
	printarray(azim_array);<br />
<br />
	cout<<"\nThe anisotropy array is:\n";<br />
	printarray(lambda_array);<br />
<br />
	cout<<"\nThe original range array is:\n";<br />
	printarray2(r_array);<br />
<br />
	//----------------------------------------------------------------<br />
<br />
	mult(fr_array2);<br />
	cout<<"\nThe final range array is:\n";<br />
	printarray2(fr_array2);<br />
<br />
<br />
<br />
	return 0;<br />
}<br />
<br />
	//code for multiplication of the matrices<br />
	void mult(double *fr_array2)<br />
	{<br />
		for(int j=0;j<cols;j++)	//for each column<br />
		{<br />
			for(int i=0;i<rows;i++)	//for each row<br />
			{<br />
				fr_array2 [j] += r_array[i]*azim_array[i][j];<br />
			}<br />
		}<br />
	}