Tutorial013_NumpyExercises.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## EXERCISES\n",
    "\n",
    "### Vector Multiplication\n",
    "1. Produce 2 vectors, one with integers in the range [5,10] and the other [15,20] using the np.arange function\n",
    "1. Without using any functions from the numpy module\n",
    "    1. Compute the outer product of those vectors \n",
    "    1. Compute the \"trace\" of the resulting matrix\n",
    "    \n",
    "    ![trace](https://wikimedia.org/api/rest_v1/media/math/render/svg/3e5b6e82272fc5eeca6d510388e0a2bd0a6c6463)\n",
    "    \n",
    "    \n",
    "Complete the same items using numpy operations"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 45,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "x: [ 5  6  7  8  9 10]\n",
      "y: [15 16 17 18 19 20]\n",
      "outer product [[75, 80, 85, 90, 95, 100], [90, 96, 102, 108, 114, 120], [105, 112, 119, 126, 133, 140], [120, 128, 136, 144, 152, 160], [135, 144, 153, 162, 171, 180], [150, 160, 170, 180, 190, 200]]\n",
      "75\n",
      "diagonal [75, 96, 119, 144, 171, 200]\n",
      "trace: 805\n",
      "with numpy\n",
      "outer product: [[ 75  80  85  90  95 100]\n",
      " [ 90  96 102 108 114 120]\n",
      " [105 112 119 126 133 140]\n",
      " [120 128 136 144 152 160]\n",
      " [135 144 153 162 171 180]\n",
      " [150 160 170 180 190 200]]\n",
      "diagnoal: [75, 96, 119, 144, 171, 200]\n",
      "trace: 805\n"
     ]
    }
   ],
   "source": [
    "import numpy as np\n",
    "\n",
    "x=np.arange(5,11)\n",
    "y=np.arange(15,21)\n",
    "print(\"x:\",x)\n",
    "print(\"y:\",y)\n",
    "\n",
    "c_out = [[x[i]*y[j] for j in range(len(x))]for i in range(len(y))]\n",
    "print('outer product',c_out)\n",
    "\n",
    "print(c_out[0][0])\n",
    "\n",
    "trace=[]\n",
    "for i in range(len(x)):\n",
    "    trace.append(c_out[i][i])\n",
    "    \n",
    "print('diagonal',trace)\n",
    "print('trace:',sum(trace))\n",
    "\n",
    "print(\"with numpy\")\n",
    "\n",
    "c_out_np=np.outer(x,y)\n",
    "trace_np=np.trace(c_out_np)\n",
    "print(\"outer product:\",c_out_np)\n",
    "print(\"diagnoal:\",trace)\n",
    "print(\"trace:\",np.sum(trace))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Matrix Multiplication\n",
    "* Two matricies can be multiplied if their inner dimensions match (eg. 2x3 * 3x5 -> 2x5). The best rule for working with and mutliplying matricies is to remember (rows x columns). This applies to both the dimensions of a matrix (a 2x3 matrix has 2 rows and 3 columns) as well as multiplication (you multiply the rows of the first matrix by the columns of the second). \n",
    "\n",
    "* When multiplying two matricies of dimension (M x N)*(N x P) the resulting matrix is (M x P). The upper element of the reslting matrix is the inner (or dot) product of the first row of the first matrix and the first column of the second matrix\n",
    "\n",
    "![two matrices](https://wikimedia.org/api/rest_v1/media/math/render/svg/16b1644351bc2041175b19cbc65da03ef78130c7)\n",
    "\n",
    "![store product in matrix C](https://wikimedia.org/api/rest_v1/media/math/render/svg/00ac0c831c365b7424cc43239aae8cebea27c56c)\n",
    "\n",
    "![matrix multiply](https://wikimedia.org/api/rest_v1/media/math/render/svg/3cfeccef1c8c7e6da0ddf08daed8dbf3c6f50c5e)\n",
    "\n",
    "for i = 1, ..., n and j = 1, ..., p.\n",
    "\n",
    "1. Make two matrices of random numbers (A and B). A should be a 4x3 matrix and B should be a 3x4 matrix. Multiply A by B using (to a resulting matrix C) using:\n",
    "    1. a conventional for-loop\n",
    "    1. list comprehension\n",
    "    1. numpy operator\n",
    "    \n",
    "2. After you have C, pull out the upper quadrant using fancy indexing, and then replace the main diagonal (upper left to lower right) with 0s. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 73,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "A: [[0.71786291 0.12409882 0.866518  ]\n",
      " [0.3801973  0.37084894 0.68146901]\n",
      " [0.08837354 0.7672099  0.6970603 ]\n",
      " [0.81200587 0.70626258 0.61371151]]\n",
      "B: [[0.96290585 0.86195963 0.24913037 0.01962675]\n",
      " [0.67537133 0.20328856 0.73349844 0.03729123]\n",
      " [0.60301835 0.04334068 0.15324645 0.94436693]]\n",
      "[[0.69123439302205, 0.08381278796600462, 0.522526257312649], [0.3277147210415026, 0.0753893475713185, 0.029535330106429204], [0.022016532637883404, 0.5627472678802744, 0.10682201803875399], [0.015937033378788996, 0.026337399602063924, 0.5795688568591897]]\n",
      "[[1.29757344 0.6815522  0.40265855 0.83702806]\n",
      " [1.02749327 0.4326394  0.47116852 0.66484825]\n",
      " [1.02358713 0.26235049 0.69158582 0.68862538]\n",
      " [1.62895401 0.87009006 0.81438693 0.62184329]]\n"
     ]
    }
   ],
   "source": [
    "A = np.random.rand(12).reshape(4,3)\n",
    "B = np.random.rand(12).reshape(3,4)\n",
    "\n",
    "print(\"A:\",A)\n",
    "print(\"B:\",B)\n",
    "\n",
    "#print(len(A[0,:]))\n",
    "\n",
    "C=[]\n",
    "for i in range(len(A[:,0])):\n",
    "    C_row=[]\n",
    "    for j in range(len(A[0,:])):\n",
    "        C_row.append(A[i,j]*B[j,i])\n",
    "    \n",
    "    C.append(C_row)\n",
    "    \n",
    "print(C)\n",
    "\n",
    "C=np.matmul(A,B)\n",
    "print(C)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.0"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}