numpy matrix operations | ones() function