ECE 4122/6122 Lab 1: Electric Field Produced by Array of Point Charges (100 pts)

Objective:

To understand and apply the principles of class creation and inheritance in C++ while leveraging std::thread for parallel calculations in a computationally-intensive problem related to the Electric Field calculations.

Description:

Given an 2D grid of N x M point charges, all with the same charge q, in the x-y plane, centered around the origin, your task is to compute the magnitude and direction of the electric field at a specific point (x, y, z) in space due to all the charges.

You will ∙use the formula below to calculate the electric field vector:
 ⃗ =    2        ̂
where, ⃗
is the electric field vector,
is magnitude of the charge of the point, 922

    ̂is the radial distance between the point charge and the (x, y, z) location in space, is the unit vector from the point charge to the (x, y, z) location in space.
Your code should include the following:

Base Class: ECE_PointCharge in files ECE_PointCharge.h & ECE_PointCharge.cpp

Protected members variables:

    • double x; // x-coordinate.

    • double y; // y-coordinate.

    • double z; // z-coordinate.

    • double q; // charge of the point.
Public member functions:

    • A constructor that initializes the above member variables.

    • A function void setLocation (double x, double y, double z);

    • A function void setCharge (double q);

Derived Class: ECE_ElectricField in files ECE_ElectricField.h & ECE_ElectricField.cpp

Inherits from: ECE_PointCharge
Protected additional member variables:

    • double Ex; // Electric field in the x-direction.

    • double Ey; // Electric field in the y-direction.

    • double Ez; // Electric field in the z-direction.

Public member functions:

    • A constructor that initializes base class member variables.

    • void computeFieldAt(double x, double y, double y);
        o Calculates the electric field at the point (x, y, z) due to the charge using the above formula. Updates the Ex, Ey, Ez member variables.

    • void getElectricField(double &Ex, double &Ey, double Ez);

OverView:

Create an application that uses the above classes and multithreading to calculate the electric field at a point in space.

    • Your program should notify the user how many threads can run concurrently when doing the calculation.

    • Your program needs to query the user for the size of the N x M array and the charge q.

    • Your program should continuously prompt the user if a new calculation is wanted.

    • Your program should use multithreading to reduce the processing time taken.

    • After calculating the resultant electric field your program should output the Ex, Ey, Ez values and the resultant electric field strength.

    • Your program should also output how long it takes in microseconds from the time the user inputs the xyz location to just before the results are output to the screen.

    • Spaces separate the values when entering multiple inputs.

    • Make sure you check for valid inputs
    o N and M should be natural numbers.

        o Separation distances should be > 0.0 and be valid numerical values o Charge should be a valid numerical value.

o  Point location should be made up of valid numerical values.
    • Use multithreading to breakup the array so that each thread does the calculations for part of the array.

    • Make sure that your code checks for a location that may be the same as a point charge location.
Sample Program Flow:

    • Your computer supports ## concurrent threads.
    • Please enter the number of rows and columns in the N x M array: 100 100
    • Please enter the x and y separation distances in meters: 0.01 0.03
    • Please enter the common charge on the points in micro C: 0.02


    • Please enter the location in space to determine the electric field (x y z) in meters: 1.0 2.0 3.0

    • The electric field at (1.0, 2.0, 3.0) in V/m is
    • Ex = x.xxxx * 10^y
    • Ey = x.xxxx * 10^y
    • Ez = x.xxxx * 10^y
    • |E| = x.xxxx * 10^y
    • The calculation took x.xxxx microsec!
    • Do you want to enter a new location (Y/N)? Y


    • Please enter the location in space to determine the electric field (x y z) in meters: 2.0 2.0 2.0

    • The electric field at (2.0, 2.0, 2.0) in V/m is
    • Ex = x.xxxx * 10^y
    • Ey = x.xxxx * 10^y
    • Ez = x.xxxx * 10^y
    • |E| = x.xxxx * 10^y
    • The calculation took x.xxxx microsec!
    • Do you want to enter a new location (Y/N)? N
    • Bye!


Turn-In Instructions

Place your files in a zip file called Lab1.zip and upload this zip file on the assignment section of Canvas.

Grading Rubric:
If a student’s program runs correctly and produces the desired output, the student has the potential to get a 100 on his or her homework; however, TA’s will look through your code for other elements needed to meet the lab requirements. The table below shows typical deductions that could occur.


AUTOMATIC GRADING POINT DEDUCTIONS PER PROBLEM:


Element
Percentage

Details


Deduction




Does Not Compile
40%

Code does not compile on PACE-ICE!







Does Not Match Output
Up to 90%

The code compiles but does not produce correct outputs.







Runtime and efficiency of
Up to 20%

The code generates the correct output but runs slower than expected.

code setup











Clear Self-Documenting
Up to 25%

This can include incorrect indentation, using unclear variable names,

Coding Styles


unclear/missing comments, or compiling with warnings. (See




Appendix A)






LATE POLICY





Element
Percentage Deduction
Details







Late Deduction Function
score – 0.5 * H

H = number of hours (ceiling function) passed





deadline







Appendix A: Coding Standards

Indentation:

When using if/for/while statements, make sure you indent 4 spaces for the content inside those. Also make sure that you use spaces to make the code more readable.

For example:

for (int i; i < 10; i++)

{
j = j + i;
}

If you have nested statements, you should use multiple indentions. Each { should be on its own line (like the for loop) If you have else or else if statements after your if statement, they should be on their own line.

for (int i; i < 10; i++)
{

if (i < 5)
{
counter++;

k -= i;
}
else

{
k +=1;
}

j += i;
}

Camel Case:

This naming convention has the first letter of the variable be lower case, and the first letter in each new word be capitalized (e.g. firstSecondThird).

This applies for functions and member functions as well!

The main exception to this is class names, where the first letter should also be capitalized.

Variable and Function Names:

Your variable and function names should be clear about what that variable or function represents. Do not use one letter variables, but use abbreviations when it is appropriate (for example: “imag" instead of “imaginary”). The more descriptive your variable and function names are, the more readable your code will be. This is the idea behind self-documenting code.


File Headers:

Every file should have the following header at the top
/*
Author: your name
Class: ECE4122 or ECE6122 (section)

Last Date Modified: date

Description:

What is the purpose of this file?

*/

Code Comments:

    1. Every function must have a comment section describing the purpose of the function, the input and output parameters, the return value (if any).
    2. Every class must have a comment section to describe the purpose of the class.

    3. Comments need to be placed inside of functions/loops to assist in the understanding of the flow of the code.