Can' t use #pragma in DSL

[tmtoku]

I tried to use #pragma in DSL, but an error has occured.

For example, I inserted #pragma unroll 2 in the second line of
sample/c++/Nbody/kernel.pikg.

F64 eps2
#pragma unroll 2
rij = EPI.pos - EPJ.pos
r2 = rij * rij + eps2
r_inv  = rsqrt(r2)
r2_inv = r_inv * r_inv
mr_inv  = EPJ.mass * r_inv
mr3_inv = r2_inv * mr_inv
FORCE.acc -= mr3_inv * rij
FORCE.pot -= mr_inv

Then, I ran make and got the following error message.

conversion type: reference
epi name: Particle
epj name: Particle
force name: Particle
class file: particle.hpp
output file name: kernel.hpp
input file: kernel.pikg
/home/tom/pikg/src/gen_hash.rb:20:in `block (2 levels) in fusion_iotag': undefined method `fusion_iotag' for #<Pragma:0x0000558cf91a1dd8 @name="unroll", @option=["2"]> (NoMethodError)
	from /home/tom/pikg/src/gen_hash.rb:19:in `each'
	from /home/tom/pikg/src/gen_hash.rb:19:in `block in fusion_iotag'
	from /home/tom/pikg/src/gen_hash.rb:18:in `each'
	from /home/tom/pikg/src/gen_hash.rb:18:in `fusion_iotag'
	from /home/tom/pikg/src/gen_hash.rb:313:in `generate_alias'
	from /home/tom/pikg/src/parserdriver.rb:1681:in `<top (required)>'
	from ../../..//bin/pikg:2:in `require_relative'
	from ../../..//bin/pikg:2:in `<main>'
make: *** [Makefile:50: kernel.hpp] Error 1

If we add the following fusion_iotag method to the Pragma class, the error will no longer occur, but the generated kernel will not be changed by #pragma.

def fusion_iotag(iotag)
  []
end

PIKG/src/intermediate_exp_class.rb

Lines 754 to 780 in 751ef58

	class Pragma
	attr_accessor :name, :option
	def initialize(x)
	@name,@option = x
	end
	def get_related_variable
	[]
	end
	def declare_temporal_var
	[]
	end
	def convert_to_code(conversion_type="reference")
	ret = String.new
	if conversion_type != "reference"
	ret = "#pragma #{@name}"
	if @option != nil then
	@option.each{ |x|
	ret += " #{x}"
	}
	end
	end
	ret
	end
	end

[subarutaro]

Sorry for late reply. Now I can’t try your code because I am in vacation, but please try —conversion-type A64FX option. Pragma only works for A64FX mode.

[tmtoku]

Thank you for your reply.

I tried to run it in A64FX mode, but I still get the same error.

command

../../..//bin/pikg --conversion-type A64FX --epi-name Particle --epj-name Particle --force-name Particle --class-file particle.hpp --output kernel.hpp -i kernel.pikg

result

conversion type: A64FX
epi name: Particle
epj name: Particle
force name: Particle
class file: particle.hpp
output file name: kernel.hpp
input file: kernel.pikg
/home/tom/pikg/src/gen_hash.rb:20:in `block (2 levels) in fusion_iotag': undefined method `fusion_iotag' for #<Pragma:0x0000559c79605e08 @name="unroll", @option=["2"]> (NoMethodError)
	from /home/tom/pikg/src/gen_hash.rb:19:in `each'
	from /home/tom/pikg/src/gen_hash.rb:19:in `block in fusion_iotag'
	from /home/tom/pikg/src/gen_hash.rb:18:in `each'
	from /home/tom/pikg/src/gen_hash.rb:18:in `fusion_iotag'
	from /home/tom/pikg/src/gen_hash.rb:313:in `generate_alias'
	from /home/tom/pikg/src/parserdriver.rb:1681:in `<top (required)>'
	from ../../..//bin/pikg:2:in `require_relative'
	from ../../..//bin/pikg:2:in `<main>'

[subarutaro]

Hi, I am sorry for very late reply.
I finally came back from vacation.

To utilize unroll and loop_fission pragma, see below:

1. adding fusion_iotag for Pragma class (as you did, this modification will be included to next release).
2. use both --conversion-type A64FX and --strip-mining N (N should be size of strip mining, e.g. 24) options.

Please use pragma for unroll and loop_fission below (N should be numbers you want to set):
#pragma unroll N
#pragma statement loop_fission_point

please note that N for loop unrolling must be divisor of strip mining size.

[tmtoku]

Thank you! By adding '--strip-mining' option, it worked.

However, there is a syntax error in the innermost j-loop of the generated code.

"kernel.hpp", line 150: error: no suitable conversion function from "const PIKG::F64vec" to "float64_t" exists
  EPJ_pos_swpl0.v0 = svdup_n_f64(epj[(j+jj_swpl0)].pos).v0;

I think we need to call svdup_n_f64x3 instead of svdup_n_f64.

When type is F64vec, get_type_suffix_a64fx(type) returns f64, but the suffix of svdup_n_ should be f64x3.

PIKG/src/A64FX.rb

Lines 89 to 115 in 751ef58

	def get_type_suffix_a64fx(type)
	case type
	when "F64" then
	suffix = "f64"
	when "F32" then
	suffix = "f32"
	when "S64" then
	suffix = "s64"
	when "S32" then
	suffix = "s32"
	when "U64" then
	suffix = "u64"
	when "U32" then
	suffix = "u32"
	when /F64vec/ then
	suffix = "f64"
	when /F32vec/ then
	suffix = "f32"
	when "B32" then
	suffix = "b32"
	when "B64" then
	suffix = "b64"
	else
	abort "error: unsupported scalar type of #{@type} for A64FX (get_type_suffix)"
	end
	suffix
	end

[subarutaro]

Thanks for the feedback.
In our implementation, the vector operation is divided into scalar operations before actual kernel generation. Thus, we should have "EPJ_pos_swpl0.v0 = svdup_n_f64(eps[(j+jj_swpl0)].pos.x)" (and same expression for v1 and v2 as well) here. Let me check and fix this.

[tmtoku]

Okay, I understand.
For now, I was able to correct the generated code by running the following (vi) commands.

:%s/).v0/.x)
:%s/).v1/.y)
:%s/).v2/.z)